Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,640,045	+C	coding (13/534 nt)	rrrQ ←	putative lysozyme

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,640,045	1	.	C	97.8%	164.5 / ‑3.3	46	coding (13/534 nt)	rrrQ	putative lysozyme
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (24/21);  total (24/22)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

CATGCCAGCGCCTTATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGT‑TACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑T‑TTTTGTGTTCATAGATATTTATGATGAGGACGCTCGTGCTTATTGGCAGGATTTTCAATCTTAAAGGAGTACTGATGCTGCAGATAAGACTCAACTTTTTCTGACAATTTTTCTGCTACTTCCAGGAAGACTTGCCGGACGCTCCTTCTGGCTGCTGCCTCATAAAACTCCAGCGCAGCTCCTTCAACACGGTCCATGGCGACATCCAGGTCAAAAATTTCACCGTCAAAGCGTTCTTTGTCC  >  NC_000913/1639795‑1640288                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                    cATGCCAGCGCCTTATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGT‑TACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑t                                                                                                                                                                                                                                                      <  1:11577/251‑1 (MQ=255)              tATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGT‑TACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAc                                                                                                                                                                                                                                         >  1:293738‑M1/1‑250 (MQ=255)               atcgcgctcAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGT‑TACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcc                                                                                                                                                                                                                                        <  1:42933‑M1/247‑3 (MQ=255)                   cGCTCAATGGCGTGGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTCCCATCCACCCCTGTGGCCCCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGT‑TACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcc                                                                                                                                                                                                                                    <  2:141207‑M1/251‑7 (MQ=255)                         aTGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATCGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTTTGTGGT‑TACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCATCACCAATCAGCGCCATAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttc                                                                                                                                                                                                                              <  1:488612‑M1/251‑13 (MQ=255)                           ggCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGT‑TACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgt                                                                                                                                                                                                                            >  1:286456‑M1/1‑237 (MQ=255)                                 gACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGT‑TACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacg                                                                                                                                                                                                                      >  1:301564‑M1/1‑231 (MQ=255)                                    cTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGT‑TACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGCACGGCAGCCGACCGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcatctgtcgct                                                                                                                                                                                                                    >  2:138376‑M1/1‑228 (MQ=255)                                      ggTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGT‑TACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatatt                                                                                                                                                                                                                 <  1:298227‑M1/251‑26 (MQ=255)                                         cGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGT‑TACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGNTCAcccgcccgcttcgtcttacgatattcat                                                                                                                                                                                                              >  2:237587‑M1/1‑223 (MQ=255)                                                 tCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGT‑TACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcctaaact                                                                                                                                                                                                      >  2:287398‑M1/1‑215 (MQ=255)                                                               cATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGT‑TACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactg                                                                                                                                                                                        <  1:299232‑M1/251‑51 (MQ=255)                                                                  atTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGT‑TACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactga                                                                                                                                                                                       >  1:532694‑M1/1‑198 (MQ=255)                                                                  atTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGAATGCCATTGTGTGGT‑TACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGATCAccctcccgcttcgtcgtacgctattcatcctaaactttcggaacatactgaat                                                                                                                                                                                     >  1:325792‑M1/1‑198 (MQ=255)                                                                     gggAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGT‑TACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctc                                                                                                                                                                                  >  1:269081‑M1/1‑195 (MQ=255)                                                                           aaCGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGT‑TACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGTTGCGCAAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgcc                                                                                                                                                                            <  1:24502‑M1/251‑63 (MQ=255)                                                                              gaattttcaaTCCACCACGGTGGCACCCCGAAAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGAGTGGT‑TACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcatcaactttgggatcatactgaagctccccgcc                                                                                                                                                                            >  2:447829‑M1/4‑186 (MQ=255)                                                                                   ttCCATCCACCCCTGCGGCCCCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGT‑TACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcc                                                                                                                                                                     <  1:250200‑M1/250‑70 (MQ=255)                                                                                      catccaCCACTGTGGCCCCCCGACAGATGGTCCATATGCCAGACCCATCGCGGTATGCCATTGTGTGGT‑TACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcct                                                                                                                                                                 <  2:301564‑M1/251‑74 (MQ=255)                                                                                            ccacTGTGGCCCCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCAATGGTGTGGT‑TCCCTTCTTTTGCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCGTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggagcatactgaagctccccgccagatgcctcctgtagac                                                                                                                                                           <  2:164772‑M1/251‑80 (MQ=255)                                                                                                 gtgGCCCCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGT‑TACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatc                                                                                                                                                      <  1:500977‑M1/251‑85 (MQ=255)                                                                                                  tgGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGT‑TACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGACAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatacagaagctccccgccagatgcctcctgtagacgcatcc                                                                                                                                                     >  2:92865‑M1/1‑166 (MQ=255)                                                                                                          cGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGT‑TACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgacct                                                                                                                                             <  1:287398‑M1/251‑94 (MQ=255)                                                                                                           gACAGTTGGTCCATATGCCAGAACTATCGCGGTAGGCCCTTGT‑TGGTGTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGTTGCGCCAGCACCAATCGGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTTTTCCcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgacctt                                                                                                                                            <  2:356602‑M1/251‑95 (MQ=255)                                                                                                                 tGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGG‑TACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcggga                                                                                                                                       >  2:107212‑M1/1‑151 (MQ=255)                                                                                                                           gCCAGAACCATCGCGGTATGCCATTGTGTGGT‑TACCTTCTTTTTCGTCCAGAACCTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGACGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgcaagttgactcctgtaganacatcgcgagcccttcgggaactaaaacccc                                                                                                                            >  2:263364‑M1/1‑141 (MQ=255)                                                                                                                                   cATCGGGGTATGCCATTGTGTGGT‑TACCTTCTTTTTCGTCCAAAAACGGGTCAAGTAGCTGAGGAGCAGATGCTCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGGCTTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagc                                                                                                                    <  2:89854‑M1/251‑119 (MQ=255)                                                                                                                                    agcgcggttggcctttgtgtggt‑tacgttcTTTTCCTTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagct                                                                                                                   <  2:332614‑M1/237‑120 (MQ=255)                                                                                                                                                      gtgGT‑TACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGAGGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaac                                                                                                  <  1:55341‑M1/250‑137 (MQ=255)                                                                                                                                                        ggT‑TACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacga                                                                                                >  1:56303‑M1/1‑112 (MQ=255)                                                                                                                                                          t‑tACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCCCCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaa                                                                                               <  1:519225‑M1/249‑140 (MQ=255)                                                                                                                                                          t‑tACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaa                                                                                               >  2:519225‑M1/1‑110 (MQ=255)                                                                                                                                                                    ttttCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgct                                                                                       <  1:92865‑M1/248‑148 (MQ=255)                                                                                                                                                                         gTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgta                                                                                                                                                              >  2:399578‑M1/1‑96 (MQ=255)                                                                                                                                                                         gTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgta                                                                                                                                                              <  1:399578‑M1/172‑77 (MQ=255)                                                                                                                                                                                     gTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgcagccaaacgaatacctgctgattgtgccagttttgtttt                                                                   >  2:220283‑M1/1‑84 (MQ=255)                                                                                                                                                                                      tCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGACGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgcccgctgcctcctgtagacgcatagcgcgaccttcgggaa                                                                                                                                      >  2:296002‑M1/1‑83 (MQ=255)                                                                                                                                                                                       cAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgttttt                                                                  <  1:58426‑M1/250‑169 (MQ=255)                                                                                                                                                                                            aTCTGAGGGGCATATGCGCCAGCACCAACCAGGGCAATAACGTCAGCCGACAGGACGTCTCTGT‑TCTTTTTTGTTCAcccgcccgcgtcgtcttacgatattcatcataaacgttgcgatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcggtaactaaatcccctttccagctataaatcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaacc                                                             <  2:190095‑M1/250‑174 (MQ=255)                                                                                                                                                                                               tGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTACAcccgcccgcttcgtcttacgctattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacacatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgatgcttgtgcaagttttgtttttgaaccgaa                                                          >  1:499360‑M1/1‑74 (MQ=255)                                                                                                                                                                                                             gcCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtc                                           >  1:429294‑M1/1‑60 (MQ=255)                                                                                                                                                                                                                                        aGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagttttaagcatttaaaacacctttat                >  1:54107‑M1/1‑33 (MQ=255)                                                                                                                                                                                                                                           cGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagttttaagcatttaaaacacctttattgt             >  1:241919‑M1/1‑30 (MQ=255)                                                                                                                                                                                                                                                ggCCTTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagttttaagcatttaaaacacctttattgttagtc        <  2:493305‑M1/251‑227 (MQ=255)                                                                                                                                                                                                                                                   cGTATCTGAGC‑TTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagcca                                                                                                     <  2:150516‑M1/155‑134 (MQ=255)                                                                                                                                                                                                                                                   cGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagcca                                                                                                     >  1:150516‑M1/1‑22 (MQ=255)                                                                                                                                                                                                                                                      aTCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagttttaagcatttaaaacacctttattgttagtcataact  >  2:186534‑M1/1‑19 (MQ=255)                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                    CATGCCAGCGCCTTATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGT‑TACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑T‑TTTTGTGTTCATAGATATTTATGATGAGGACGCTCGTGCTTATTGGCAGGATTTTCAATCTTAAAGGAGTACTGATGCTGCAGATAAGACTCAACTTTTTCTGACAATTTTTCTGCTACTTCCAGGAAGACTTGCCGGACGCTCCTTCTGGCTGCTGCCTCATAAAACTCCAGCGCAGCTCCTTCAACACGGTCCATGGCGACATCCAGGTCAAAAATTTCACCGTCAAAGCGTTCTTTGTCC  >  NC_000913/1639795‑1640288

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑