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.7%	140.0 / ‑4.5	43	coding (13/534 nt)	rrrQ	putative lysozyme
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (16/26);  total (16/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑T‑TTTTGTGTTCATAGATATTTATGATGAGGACGCTCGTGCTTATTGGCAGGATTTTCAATCTTAAAGGAGTACTGATGCTGCAGATAAGACTCAACTTTTTCTGACAATTTTTCTGCTACTTCCAGGAAGACTTGCCGGACGCTCCTTCTGGCTGCTGCCTCATAAAACTCCAGCGCAGCTCCTTCAACACGGTCCATGGCGACATCCAGGTCAAAAATTTCACCGTCAAAGCGTTC  >  NC_000913/1639821‑1640281                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                             ggCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgt                                                                                                                                                                                                                     <  1:101341‑M1/251‑15 (MQ=255)      tGACCGGGTCGCAGTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGGCAGGTATCGGAGGAGCAGAGGCGCCAGCGCCAATCGGCGCCAGAACGGCGGCCGACATTCCGTATCTGA‑TCGTTTGTGTTCAcccgcccgcttcgtctta                                                                                                                                                                                                                 <  2:70435‑M1/250‑19 (MQ=255)                 cATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatc                                                                                                                                                                                                      >  2:110563‑M1/1‑221 (MQ=255)                    ttttCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcata                                                                                                                                                                                                   >  2:391196‑M1/1‑218 (MQ=255)                        ccTTCGACAGGTTCATATTGGGAAAAACGGTTTTTCTATCCACCACTGTGGCACCCCGACATATTGTCCATATGCAAGAACCAGCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCGGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgctcgcttcgtcttacgatattcatcataaact                                                                                                                                                                                               <  2:106258‑M1/250‑37 (MQ=255)                             gACAGTTTCATATTGGGAAAAACGGATTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgtg                                                                                                                                                                                          >  2:266537‑M1/1‑209 (MQ=255)                                aGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatc                                                                                                                                                                                       >  1:340172‑M1/1‑206 (MQ=255)                                    tCATATTGGGAAAAACGGTTTTTCCATCCCCCACTGTGGCACCCCGAAAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTCCCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatac                                                                                                                                                                                   <  1:407229‑M1/250‑49 (MQ=255)                                          ggggcaaaacggtaTTTCCATCCACCCGTGTGGCACCCCGACAGATGGTCCATCTGCCATACCCATCGCGGTATGCCATTGTGTGGGTACCTTCTTTTTCGTACAGAAACTGGTCAAGTACTTAAGGAGCAAATGCGCCAGCACCAATCAGCGCCAGAACGGCTGCCGACAGTCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataatctttgggatcatactgaagct                                                                                                                                                                            <  2:1633‑M1/246‑56 (MQ=255)                                            ggAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGAATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctcc                                                                                                                                                                          >  2:201735‑M1/1‑194 (MQ=255)                                             gAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccc                                                                                                                                                                         <  1:312569‑M1/251‑59 (MQ=255)                                                   cGGTTTTTCCTTCCACCACTGTGCAACCCCGACAGATGGTCCATATGCAAGATCCATCGCGGTATGCCTTTGTTTGGTTGCCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTAAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGAAAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcggcttacgatattcatcataaactttgggatcatactgaagctccccgccag                                                                                                                                                                   <  1:317886‑M1/251‑65 (MQ=255)                                                      tttttCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTCTGTGTTCAcccgcccccttcgacttacgatcttcatcataaactttgggatcatactgaacctccccctcagat                                                                                                                                                                 >  2:227734‑M1/1‑184 (MQ=255)                                                                       gtgGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctcccagccagatgcctcctgtagacgcatc                                                                                                                                               >  1:152953‑M1/1‑167 (MQ=255)                                                                            aCCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcg                                                                                                                                          >  1:416482‑M1/1‑162 (MQ=255)                                                                                cGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAACCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgacct                                                                                                                                      <  2:194603‑M1/251‑94 (MQ=255)                                                                                     gATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcggg                                                                                                                                 <  1:167419‑M1/251‑99 (MQ=255)                                                                                           ccaattgcccgAACCATCGCGTGATGCCATTGTGTGGTTACCGTCTTGTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCATCACCAATCAGCGCCTGAACGGCACCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCCccctcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagctgcctcctgtagacgcatcgcgcgaccttcgggaactaa                                                                                                                           <  2:272186‑M1/246‑105 (MQ=255)                                                                                                   cAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctt                                                                                                                   <  2:340172‑M1/251‑113 (MQ=255)                                                                                                    aGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatccccttt                                                                                                                  >  1:119478‑M1/1‑138 (MQ=255)                                                                                                      aaCCATCGCGGTATGCCATTGTGTGGTTTCCTTCTTTTTCGTCCAGAAACTGGTCAAGTATGTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGTCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttcc                                                                                                                <  2:221404‑M1/251‑116 (MQ=255)                                                                                                        ccATAGCGGTATGCCATTGTGTGGTTCCCTGCTTTTGCGTCCAGAAACGGGTCAAGTATCGGAGGAGCAGATGCGCCAGCACCAATCAGCTCCAGAACGGCGGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccag                                                                                                              <  2:9722‑M1/251‑118 (MQ=255)                                                                                                          aTCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcatacactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagct                                                                                                            >  2:206053‑M1/1‑132 (MQ=255)                                                                                                          aTCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatactcatcataaacgttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagct                                                                                                            <  1:161666‑M1/251‑120 (MQ=255)                                                                                                              cgGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataa                                                                                                        <  1:110563‑M1/251‑124 (MQ=255)                                                                                                                            gtgGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacg                                                                                          >  1:190578‑M1/1‑114 (MQ=255)                                                                                                                                    cttcttTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctg                                                                                  >  1:154670‑M1/1‑106 (MQ=255)                                                                                                                                         ttttCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgct                                                                                <  1:266537‑M1/248‑148 (MQ=255)                                                                                                                                                 cAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaag                                                                     <  1:391196‑M1/251‑159 (MQ=255)                                                                                                                                                 cAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaacttggggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgcttcttgtgcaag                                                                     <  1:253187‑M1/251‑159 (MQ=255)                                                                                                                                                   aaacctgttcagggatctgcggagcagatGCTCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGTCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaag                                                                     <  2:301460‑M1/235‑159 (MQ=255)                                                                                                                                                    aaaCTGGTCAAGTATCTTAGGAGCAGATGCGCCAGCCCCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcacgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttt                                                                  <  2:439697‑M1/251‑162 (MQ=255)                                                                                                                                                    aaaCTGGACAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgctaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgttgcttgtgcaagttt                                                                  >  2:382094‑M1/1‑90 (MQ=255)                                                                                                                                                          gTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttt                                                            >  1:266543‑M1/1‑84 (MQ=255)                                                                                                                                                             aaggatcggaggagcagatgcgCCAGCACCATTCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttga                                                         <  2:416482‑M1/243‑171 (MQ=255)                                                                                                                                                              aGTATCTGAGGAGCGGGTGCGCCAGCACCATCCAGCGCCAGACCGGCAGCCGACAGGCCTTATCTGA‑TCTTTTGTTTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaataccggctgcttgtgcaagttttgtttttgaa                                                        <  2:339689‑M1/251‑172 (MQ=255)                                                                                                                                                                 aTCTGAGGAGCATATGCTCCAGCACCAATCACCGCCAGAACTGCAGCCTACCGGCTGTATCTGACC‑TTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgtgatcatactgaagctccccgccagatgcctcctgtagacgcatcgctcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatccctgctgcttgcgcaagttttgtttttgaacc                                                      <  2:329188‑M1/250‑174 (MQ=255)                                                                                                                                                                             gATGCGACAGCACCAATCAGCGCCAGAACGGCATCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaaga                                         <  2:26772‑M1/251‑187 (MQ=255)                                                                                                                                                                                gcggcagcaCAAATCATCGCCAGACCGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacggtattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagc                                       <  2:443133‑M1/246‑189 (MQ=255)                                                                                                                                                                                   ccggcaacaatccgcGCCGTAACGGCGGCCAACAGGTCGTATCTGA‑TCTTTTTTGTTCCcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctgtccagctataaagcgaagccaaacgaatacctgatgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtca                                   <  2:213951‑M1/244‑193 (MQ=255)                                                                                                                                                                                       cACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagttt                               >  1:323539‑M1/1‑55 (MQ=255)                                                                                                                                                                                                            cAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagttttaagcatttaaaacaccttta          <  1:14404‑M1/251‑218 (MQ=255)                                                                                                                                                                                                                    aGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagttttaagcatttaaaacacctttattgttagt  >  1:166783‑M1/1‑26 (MQ=255)                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                             GGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑T‑TTTTGTGTTCATAGATATTTATGATGAGGACGCTCGTGCTTATTGGCAGGATTTTCAATCTTAAAGGAGTACTGATGCTGCAGATAAGACTCAACTTTTTCTGACAATTTTTCTGCTACTTCCAGGAAGACTTGCCGGACGCTCCTTCTGGCTGCTGCCTCATAAAACTCCAGCGCAGCTCCTTCAACACGGTCCATGGCGACATCCAGGTCAAAAATTTCACCGTCAAAGCGTTC  >  NC_000913/1639821‑1640281

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