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	100.0%	243.9 / NA	66	I5V (ATC→GTC)	rrrQ	putative lysozyme
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (36/30);  total (36/30)

ATGCCAGCGCCTTATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGAT‑TTTTGTGTTCATAGATATTTATGATGAGGACGCTCGTGCTTATTGGCAGGATTTTCAATCTTAAAGGAGTACTGATGCTGCAGATAAGACTCAACTTTTTCTGACAATTTTTCTGCTACTTCCAGGAAGACTTGCCGGACGCTCCTTCTGGCTGCTGCCTCATAAAACTCCAGCGCAGCTCCTTCAACACGGTCCATGGCGACATCCAGGTCAAAAATTTCACCGTCAAAGCGTTCTTTGTCC  >  NC_000913/1639796‑1640288                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                    aTGCCAGCGCCTTATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGGATCTGAtc                                                                                                                                                                                                                                                     >  2:370861/1‑250 (MQ=255)   gCCAGCGCCTTATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGCAAGGCCGTATCTGAtctt                                                                                                                                                                                                                                                   >  2:232364/1‑248 (MQ=255)      aGCGCCTTATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCAAGAACGGCAGCCGACAGGCCGTATCTGATCTTTtg                                                                                                                                                                                                                                                >  2:393516/1‑251 (MQ=255)             tATCACGCTCCATGGCGTTCACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAc                                                                                                                                                                                                                                         >  2:563460‑M1/1‑250 (MQ=255)                           gCGTTGACCTGGTCGCATTTTTCCTTAGACAGTTTCATATTGGGATCCACGGTTTTTCCATCCACCACTGTGGCACCACGACAGATGGTACATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtc                                                                                                                                                                                                                           >  1:159022‑M1/1‑236 (MQ=255)                             gTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTGCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtctt                                                                                                                                                                                                                         <  2:29267‑M1/251‑18 (MQ=255)                                          cATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatc                                                                                                                                                                                                             >  2:387234‑M1/1‑221 (MQ=255)                                                    tCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgg                                                                                                                                                                                                  >  1:512863‑M1/1‑211 (MQ=255)                                                        cAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatc                                                                                                                                                                                              >  2:319570‑M1/1‑207 (MQ=255)                                                         aGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggctc                                                                                                                                                                                              >  2:434424‑M1/1‑206 (MQ=255)                                                              cATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCTCACCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactg                                                                                                                                                                                        <  2:407789‑M1/251‑51 (MQ=255)                                                              cATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactg                                                                                                                                                                                        >  1:178419‑M1/1‑201 (MQ=255)                                                                 atTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactga                                                                                                                                                                                       >  2:57616‑M1/1‑198 (MQ=255)                                                                      gAAAAACGGTTTTTCCATCCACCCCTGTGGCACCCCGACAGATGGTCCAAATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCGcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccc                                                                                                                                                                                <  2:574161‑M1/251‑59 (MQ=255)                                                                       aaaaaCGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGAAGATGCGCCAGCACCAAGCAGCGCCAGAACGGCAGCCGACAGGCCGTAGCTGATCTTTGGTGTTCAcccgcccgctgcgtcttacgatattcatcataaactttgggatcatactgaagctcccc                                                                                                                                                                               <  1:181956‑M1/251‑60 (MQ=255)                                                                             ggTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgcc                                                                                                                                                                            >  2:205544‑M1/1‑186 (MQ=255)                                                                                  ttCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcc                                                                                                                                                                     <  1:95135‑M1/250‑70 (MQ=255)                                                                                       tccaCCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgcccgctgcctcctgt                                                                                                                                                               >  2:103023‑M1/1‑176 (MQ=255)                                                                                            cacTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTGCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGACCGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactgtgggatcatactgaagctccccgccagatgcctcctgtagacg                                                                                                                                                          <  2:51321‑M1/251‑81 (MQ=255)                                                                                              cTGTGGCACCCCGACAGATGGTCCATAAACAAGAACCATCGCGGTAAGCCATTGTGTGGCTACATTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGACGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccaacacagatgcctcctgtagacgc                                                                                                                                                         >  1:465119‑M1/1‑169 (MQ=255)                                                                                               tgtgGCACTCCGACAGATGGGCCATATGCCAGAACCATCGCGGAATGCCATTGTGTGGTTACATTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcatacactttgggatcataccgaagctcccctccagatgcctcctgtatacgcat                                                                                                                                                       >  2:527084‑M1/1‑168 (MQ=255)                                                                                                gtgGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatc                                                                                                                                                      <  1:84721‑M1/251‑85 (MQ=255)                                                                                                gtgGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatc                                                                                                                                                      <  1:255487‑M1/251‑85 (MQ=255)                                                                                                gtgGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcat                                                                                                                                                       <  1:57616‑M1/250‑84 (MQ=255)                                                                                                  ggCACCCCGACAGATGGTCCTTAGGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCTGCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcg                                                                                                                                                     <  2:430063‑M1/250‑86 (MQ=255)                                                                                                           aCAGATGGTCCATAGGCCAGAACCATCGCGTTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTTTCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttc                                                                                                                                           <  2:541272‑M1/251‑96 (MQ=255)                                                                                                                               aaCCCTCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttcc                                                                                                                       <  2:573585‑M1/251‑116 (MQ=255)                                                                                                                                       cgGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataa                                                                                                               <  2:388771‑M1/251‑124 (MQ=255)                                                                                                                                        ggTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataa                                                                                                               >  1:425595‑M1/1‑127 (MQ=255)                                                                                                                                           aTGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcg                                                                                                           >  1:145841‑M1/1‑124 (MQ=255)                                                                                                                                              ccATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaag                                                                                                        <  1:370861‑M1/251‑131 (MQ=255)                                                                                                                                              ccATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcga                                                                                                          >  1:420307‑M1/1‑121 (MQ=255)                                                                                                                                                        gTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaat                                                                                              >  1:88334‑M1/1‑111 (MQ=255)                                                                                                                                                           aCCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACAAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctctaaagcgaagccaaacgaatacc                                                                                           >  1:227661‑M1/1‑108 (MQ=255)                                                                                                                                                           aCCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGAGGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatac                                                                                            <  1:205544‑M1/250‑143 (MQ=255)                                                                                                                                                              ttcttTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgc                                                                                        >  2:228508‑M1/1‑105 (MQ=255)                                                                                                                                                                 tttttCGTCCAGAAACTGGTCAAGTCTCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttggtatcatactgaagctccctgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgc                                                                                     >  2:565460‑M1/1‑102 (MQ=255)                                                                                                                                                                   tttCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgct                                                                                    <  1:319570‑M1/250‑151 (MQ=255)                                                                                                                                                                    ttCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttggggtcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgcatacctgcggcttg                                                                                  >  2:319271‑M1/1‑99 (MQ=255)                                                                                                                                                                       gTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgc                                                                               >  1:182325‑M1/1‑96 (MQ=255)                                                                                                                                                                            gAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagtt                                                                          >  1:83872‑M1/1‑91 (MQ=255)                                                                                                                                                                            gAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttggaatcatactgaagctcccctccagatgcctcctgtagacgcatctcgcgaccttcgggacctaaatcccctttccagctataacgcgaagccaaacgaatacctgctacttgtgcacgtt                                                                          >  2:375756‑M1/1‑91 (MQ=255)                                                                                                                                                                              aaCTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagtttt                                                                        <  1:62743‑M1/251‑163 (MQ=255)                                                                                                                                                                                 tGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtt                                                                     >  1:604345‑M1/1‑86 (MQ=255)                                                                                                                                                                                  ggtaaagtatCTGAGGAGCGGATGCGCCACCACCACTAAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcagcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtt                                                                     <  2:512863‑M1/246‑166 (MQ=255)                                                                                                                                                                                  ggTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgttt                                                                    >  1:470426‑M1/1‑85 (MQ=255)                                                                                                                                                                                  ggTCAAGGATCGGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgttt                                                                    <  2:234415‑M1/251‑167 (MQ=255)                                                                                                                                                                                      aaGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctctaaagcgaagccaaacgcatacctgctgcttgtgcaagttttgtttttg                                                                 >  2:543262‑M1/1‑81 (MQ=255)                                                                                                                                                                                       aGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACCGGCAGTATCTGATCTTTTGTGTTCAcccgcccgattcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaaataaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttga                                                                >  1:499591‑M1/1‑80 (MQ=255)                                                                                                                                                                                          aTCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccg                                                            >  2:638031‑M1/1‑77 (MQ=255)                                                                                                                                                                                               gggagcAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaata                                                       <  1:380829‑M1/250‑180 (MQ=255)                                                                                                                                                                                                 gagCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatac                                                      >  1:168764‑M1/1‑70 (MQ=255)                                                                                                                                                                                                       aTGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagag                                               <  1:398796‑M1/251‑188 (MQ=255)                                                                                                                                                                                                         gcgcCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagc                                              <  1:264631‑M1/250‑189 (MQ=255)                                                                                                                                                                                                                 aCCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagtttt                                     >  1:91280‑M1/1‑54 (MQ=255)                                                                                                                                                                                                                   cAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagttttaa                                   <  1:387234‑M1/251‑200 (MQ=255)                                                                                                                                                                                                                       aagCGCCAGAACGGCAGCCGACAGGCCGTAGCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcatgtagacgcatcgcgcgaccttcgggaacaaaatcccctttccagctataaagcgaagccaaacgaatacatgctgcttgtgcaagttttgtttttgacccgaaatacaaaagagcggcagttttaagcat                               >  2:133202‑M1/2‑48 (MQ=255)                                                                                                                                                                                                                        aGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagttttaagcatt                              <  1:232364‑M1/251‑205 (MQ=255)                                                                                                                                                                                                                              aGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgtttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagttttaagcatttaaaac                        <  2:133993‑M1/251‑211 (MQ=255)                                                                                                                                                                                                                                 ccgGCAGCCGACAGGCCGTATCTGATCTTTTGGTTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagttttaagcatttaaaacacc                     <  2:170542‑M1/250‑214 (MQ=255)                                                                                                                                                                                                                                 aCGGCCGCCGACAGGCCGTATCTGATCTTTTGGGTCCAcccgcccgcttcgtcttacgagattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcaccgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatagctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagttttaagcatttaaaccagc                     <  2:159022‑M1/251‑214 (MQ=255)                                                                                                                                                                                                                                   ggCAGCCGCCAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaa                                                    <  1:165678‑M1/218‑183 (MQ=255)                                                                                                                                                                                                                                   ggCAGCCGACAGGCCGTATCTGCTCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgattacctgctgcttgtgcaagttttgtttttgaaccgaaatacaa                                                    >  2:165678‑M1/1‑36 (MQ=255)                                                                                                                                                                                                                                           aCAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcc                                                                                                                              <  1:349077‑M1/136‑109 (MQ=255)                                                                                                                                                                                                                                           aCAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcc                                                                                                                              >  2:349077‑M1/1‑28 (MQ=255)                                                                                                                                                                                                                                             aGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagttttaagcatttaaaacacctttattgttagt         >  2:186999‑M1/1‑26 (MQ=255)                                                                                                                                                                                                                                              ggCCGTATCTGATCTTGTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagttttaagcatttaaaacacctttattgttagtc        <  2:26467‑M1/251‑227 (MQ=255)                                                                                                                                                                                                                                                    aTCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagttttaagcatttaaaacacctttattgttagtcataact  >  1:380258‑M1/1‑19 (MQ=255)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                    ATGCCAGCGCCTTATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGAT‑TTTTGTGTTCATAGATATTTATGATGAGGACGCTCGTGCTTATTGGCAGGATTTTCAATCTTAAAGGAGTACTGATGCTGCAGATAAGACTCAACTTTTTCTGACAATTTTTCTGCTACTTCCAGGAAGACTTGCCGGACGCTCCTTCTGGCTGCTGCCTCATAAAACTCCAGCGCAGCTCCTTCAACACGGTCCATGGCGACATCCAGGTCAAAAATTTCACCGTCAAAGCGTTCTTTGTCC  >  NC_000913/1639796‑1640288

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