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	98.3%	221.2 / ‑4.7	60	coding (13/534 nt)	rrrQ	putative lysozyme
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (24/35);  total (24/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGCGCCTTATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑T‑TTTTGTGTTCATAGATATTTATGATGAGGACGCTCGTGCTTATTGGCAGGATTTTCAATCTTAAAGGAGTACTGATGCTGCAGATAAGACTCAACTTTTTCTGACAATTTTTCTGCTACTTCCAGGAAGACTTGCCGGACGCTCCTTCTGGCTGCTGCCTCATAAAACTCCAGCGCAGCTCCTTCAACACGGTCCATGGCGACATCCAGGTCAAAAATTTC  >  NC_000913/1639799‑1640266                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                              ccAGCGCCTTATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑t‑cttt                                                                                                                                                                                                                           >  1:289731/1‑251 (MQ=255)  cAGCGCCTTATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑t‑ctttt                                                                                                                                                                                                                          <  2:633289/251‑2 (MQ=255)         ttATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑t‑ctttt                                                                                                                                                                                                                          <  2:223188/244‑2 (MQ=255)         ttATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑t‑ctttt                                                                                                                                                                                                                          >  1:223188/1‑243 (MQ=255)           aTCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcc                                                                                                                                                                                                                  <  1:33607‑M1/251‑3 (MQ=255)             cACGCTCAATGGCGTTGGCCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGCGTGGTTGCCTTCTTTTTCGTCCAGAAACGGGTCAAGTATCTGAGTAGCAGATTCGCCATCATCAATCATCGCCAGAACGGCAGCCGACAGGCTGTATCTGA‑TCTTTTGTGGTCAcccg                                                                                                                                                                                                                <  2:537059‑M1/251‑5 (MQ=255)              aCGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgc                                                                                                                                                                                                               >  1:589531‑M1/1‑246 (MQ=255)                gCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgccc                                                                                                                                                                                                             <  2:124552‑M1/251‑8 (MQ=255)                   cAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGTTTTTTCCCACCACCACTGTGGCGCCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACGGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAGCGGCAGCCTACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccg                                                                                                                                                                                                            <  2:37389‑M1/249‑9 (MQ=255)                   cAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgct                                                                                                                                                                                                          <  1:326928‑M1/251‑11 (MQ=255)                     aTGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttc                                                                                                                                                                                                        <  2:620431‑M1/251‑13 (MQ=255)                       ggCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGTTTTTTCCATCCACCACTGTGGCACCCCGACGTATGGTCCATATGCCAGATCCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAATCTGGGCAAGTATCTGAGGAGCAGATGTGCCAGCCCCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcg                                                                                                                                                                                                       <  1:249693‑M1/250‑14 (MQ=255)                        gCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtc                                                                                                                                                                                                     >  1:497082‑M1/1‑236 (MQ=255)                             gACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCCCTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTCCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTAGCTGA‑TCTTGTGTGTTCAcccgcccgcttcgtcttacg                                                                                                                                                                                                <  2:3702‑M1/251‑21 (MQ=255)                             gACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacg                                                                                                                                                                                                <  1:534054‑M1/251‑21 (MQ=255)                                     cGCATTTTTCCTTCGACAGTTTCATATTGGGAAAACCGGGTTTTCCATCCACCCCTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCCGAACGGCAGCCGACAGGCCGTATCTGC‑TCTTTTGTGTTCAcccgcacgcttcgtcttacgatattcat                                                                                                                                                                                        >  2:574895‑M1/1‑223 (MQ=255)                                      gCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgttattcatc                                                                                                                                                                                       >  2:395088‑M1/1‑222 (MQ=255)                                      gCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatc                                                                                                                                                                                       <  2:11560‑M1/251‑30 (MQ=255)                                        aTTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcat                                                                                                                                                                                     >  1:500203‑M1/1‑220 (MQ=255)                                             tCCTTCGCCAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTGCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaac                                                                                                                                                                                 <  2:340898‑M1/250‑36 (MQ=255)                                             tCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaac                                                                                                                                                                                 <  1:405717‑M1/250‑36 (MQ=255)                                              ccTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaact                                                                                                                                                                                <  1:171318‑M1/250‑37 (MQ=255)                                               cTTCGACAGTTTCATATTGGGAAAAATGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCCCCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattaatcataaacttt                                                                                                                                                                              >  2:265878‑M1/1‑213 (MQ=255)                                                 tCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGCTCAcccgccctcttcgtcttacgatattcatcataaactttgg                                                                                                                                                                            >  1:168995‑M1/1‑211 (MQ=255)                                                    aCAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggat                                                                                                                                                                         <  1:98183‑M1/251‑44 (MQ=255)                                                     cAGTTTCATATTGGGAAAAACGGTTTTTCCATCCCCCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGTCCGTATCTGA‑TCTTTTGTGTTCAcccgccctcttcgtcttacgatattcatcataaactttgggatc                                                                                                                                                                        <  2:504282‑M1/251‑45 (MQ=255)                                                       gTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcat                                                                                                                                                                      >  1:345114‑M1/1‑205 (MQ=255)                                                          tCATATTGGGAAAAACGGTTTTTCCATCCACCAGTTTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGCGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAACCGCCAGGCCGTATCTGA‑TCTCTTCCGTTCAcccgcccgcttcgtcttacgatattcatcatacactttgggaccagact                                                                                                                                                                   >  2:410109‑M1/1‑202 (MQ=255)                                                           cATATTGGGAAAAACGGTTTTTCCATCCACCACTGTTGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactg                                                                                                                                                                  <  2:309888‑M1/251‑51 (MQ=255)                                                                 gggAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataagctttgggat                                                                                                                                                                         <  1:397335‑M1/238‑44 (MQ=255)                                                                 gggAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgaaattcatcataacctttgggat                                                                                                                                                                         >  2:397335‑M1/1‑195 (MQ=255)                                                                  ggAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAACCTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctcc                                                                                                                                                           >  2:304053‑M1/1‑194 (MQ=255)                                                                   gAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccc                                                                                                                                                          >  1:96593‑M1/1‑193 (MQ=255)                                                                   gAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccc                                                                                                                                                          >  1:378524‑M1/1‑193 (MQ=255)                                                                               ttCCATCCCCCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTGCCTTCTTTTGCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcc                                                                                                                                               <  2:597113‑M1/250‑70 (MQ=255)                                                                                      caccacTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgt                                                                                                                                         >  1:107540‑M1/1‑174 (MQ=255)                                                                                        ccacTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagac                                                                                                                                     >  1:619132‑M1/1‑172 (MQ=255)                                                                                             gtgGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatc                                                                                                                                <  2:96593‑M1/251‑85 (MQ=255)                                                                                             gtgGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatc                                                                                                                                >  2:20564‑M1/1‑167 (MQ=255)                                                                                                       gACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccaaatgcctcctgtagacgcatcgcgcgacctt                                                                                                                      >  2:399730‑M1/1‑157 (MQ=255)                                                                                                             tgctccaTATGCCAGAACCATCGCGGTAGGCCATTGTGTGGTTACCTTCTGTTTCGTCCAGAAACGGGTCAAGTATCTGAGGAGCAGAGGCGCCAGCCCCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgccctcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcggga                                                                                                                 <  2:321855‑M1/247‑100 (MQ=255)                                                                                                                     atGCCAGAACCATCGCGGTATGCCATTGTGTGGTTCCCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatc                                                                                                         <  2:349525‑M1/250‑108 (MQ=255)                                                                                                                            aaCCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactcaaccccctttcc                                                                                                 >  2:598761‑M1/1‑136 (MQ=255)                                                                                                                                  cgcgGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagcta                                                                                            <  1:546042‑M1/250‑121 (MQ=255)                                                                                                                                      gTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaacaaaatcccctttccagctataa                                                                                         >  2:331271‑M1/1‑126 (MQ=255)                                                                                                                                       tATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATTCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAacaacccgcttcgtcttacgctattcatcataaactttgggatgatactgacgctccccgccagatgccccctgtagacgaatcgagcgaccttcggaaacaaacacccctttccagatataaagc                                                                                      >  2:74310‑M1/1‑125 (MQ=255)                                                                                                                                          gCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgtgatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcga                                                                                    >  1:631351‑M1/1‑122 (MQ=255)                                                                                                                                          gCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaa                                                                                   <  1:331271‑M1/251‑130 (MQ=255)                                                                                                                                               tgtgtgtgtacattcgtgttcgtccagaatctggtcaaGTATCTGAGGAGCAGATCCGCCAGCACAAAGCTGCGCCAGACCGCCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTGCAcccgcccgcttcgtcttactatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagcca                                                                               <  2:278981‑M1/232‑134 (MQ=255)                                                                                                                                                gtgtgGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcggcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaa                                                                              >  2:567018‑M1/1‑116 (MQ=255)                                                                                                                                                  gtgGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacg                                                                           <  1:353751‑M1/251‑138 (MQ=255)                                                                                                                                                                     tCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgt                                                           <  2:619132‑M1/248‑154 (MQ=255)                                                                                                                                                                       cAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaag                                                      <  1:34591‑M1/251‑159 (MQ=255)                                                                                                                                                                          aaaCGGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttt                                                   <  2:176781‑M1/251‑162 (MQ=255)                                                                                                                                                                                   aaGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGCCAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttg                                           >  1:270303‑M1/1‑81 (MQ=255)                                                                                                                                                                                    agaatcggaggagcagatgcGCCAGCCCCAATCATCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaa                                         <  2:133583‑M1/244‑172 (MQ=255)                                                                                                                                                                                       aTCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGACCGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaacc                                       <  2:497082‑M1/250‑174 (MQ=255)                                                                                                                                                                                        cctgagggtccgatgcccctacaccaatcagccccagaacggcAGCCGTCAGGCCGTATCTGAGC‑TTTTGTTTTCCcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccatatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccga                                     <  2:591679‑M1/230‑176 (MQ=255)                                                                                                                                                                                             ggagCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatac                                <  1:395088‑M1/251‑181 (MQ=255)                                                                                                                                                                                              gagCAGATGCGCCAGCCCCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaattcctgctgcttgtgcaagttttgtttttgaaccgaaataca                               <  1:104737‑M1/251‑182 (MQ=255)                                                                                                                                                                                                               ccAATCAGCTCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagtttta              <  2:46147‑M1/251‑199 (MQ=255)                                                                                                                                                                                                                     atcgccaGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagttttaagcatt        <  2:252520‑M1/249‑205 (MQ=255)                                                                                                                                                                                                                      gcgcCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagttttaagcattt       <  1:37745‑M1/251‑206 (MQ=255)                                                                                                                                                                                                                           aGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagttttaagcatttaaaac  >  1:83150‑M1/1‑41 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                              CCAGCGCCTTATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑T‑TTTTGTGTTCATAGATATTTATGATGAGGACGCTCGTGCTTATTGGCAGGATTTTCAATCTTAAAGGAGTACTGATGCTGCAGATAAGACTCAACTTTTTCTGACAATTTTTCTGCTACTTCCAGGAAGACTTGCCGGACGCTCCTTCTGGCTGCTGCCTCATAAAACTCCAGCGCAGCTCCTTCAACACGGTCCATGGCGACATCCAGGTCAAAAATTTC  >  NC_000913/1639799‑1640266

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