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	96.7%	222.2 / ‑6.2	62	I5V (ATC→GTC)	rrrQ	putative lysozyme
Reads supporting (aligned to +/- strand):  ref base . (1/0);  major base C (29/30);  minor base A (2/0);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

ATGCCAGCGCCTTATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGAT‑TTTTG‑TGTTCATAGATATTTATGATGAGGACGCTCGTGCTTATTGGCAGGATTTTCAATCTTAAAGGAGTACTGATGCTGCAGATAAGACTCAACTTTTTCTGACAATTTTTCTGCTACTTCCAGGAAGACTTGCCGGACGCTCCTTCTGGCTGCTGCCTCATAAAACTCCAGCGCAGCTCCTTCAACACGGTCCATGGCGACATCCAGGTCAAAAATTTCACCGTCAAAGCGTTCTTTGTC  >  NC_000913/1639796‑1640287                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                    aTGCCAGCGCCTTATCACGCTCAATGGCGTTGACCTTGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGCCAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGAtc                                                                                                                                                                                                                                                     <  2:395934/251‑2 (MQ=255)  tGCCAGCGCCTTATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCCCTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGTTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCCGCCTACAGTCCGTATCTGAtct                                                                                                                                                                                                                                                    >  1:83788/1‑249 (MQ=255)    ccagctacTTATCACTCTCAAGGGCGTTGACCTGGGCGCATTTTTCCTTCGACAGTTTCATTTTGGGAAAATCGGTCTTTCCATCCACCCCTGGGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTTTGTGGTTACCTTCTTTTTCGTCCATAAACTGGTCAATTATCTGAGGAGCAGAGTCGCCAGCACCAATCAGCTCCAGAACGGCAGCCGACAGGCCGTATCTGAt‑cttt                                                                                                                                                                                                                                                 <  2:655386/244‑1 (MQ=255)              aTCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTG‑TGTTCAc                                                                                                                                                                                                                                        <  1:345215‑M1/250‑2 (MQ=255)              aTCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTG‑TGTTCAc                                                                                                                                                                                                                                        <  1:175315‑M1/250‑2 (MQ=255)               tCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTG‑TGTTCAccc                                                                                                                                                                                                                                      >  2:607778‑M1/1‑248 (MQ=255)                cACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTG‑TGTTCAcccg                                                                                                                                                                                                                                     <  1:306035‑M1/251‑5 (MQ=255)                cACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTG‑TGTTCAcc                                                                                                                                                                                                                                       >  1:536189‑M1/1‑247 (MQ=255)                cACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTG‑TGTTCAcc                                                                                                                                                                                                                                       <  2:536189‑M1/249‑3 (MQ=255)                cACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGGTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCAATTGAGTGGTTACCTACTTTTTCGTCCAGCACCTGGACAAGTATCGGAGGAGCAGATGCGCCAGAACCAATCAGCGCCAGAACGGCATCCGACAGGCCGGATCTGAT‑ATTTGATGTTCAcc                                                                                                                                                                                                                                       >  2:257866‑M1/1‑247 (MQ=255)                 aCGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTG‑TGTTCAcccgc                                                                                                                                                                                                                                    >  1:120414‑M1/1‑246 (MQ=255)                   gCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTG‑TGTTCAcccgccc                                                                                                                                                                                                                                  <  2:623013‑M1/251‑8 (MQ=255)                    cTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTG‑TGTTCAcccgcccg                                                                                                                                                                                                                                 >  1:668026‑M1/1‑243 (MQ=255)                        aTGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTG‑TGTTCAcccgcccgcttc                                                                                                                                                                                                                             <  1:339521‑M1/251‑13 (MQ=255)                         tGGCGTTGACCTGGACGCATTTTTCCTTCGCCAGCTACATATTGGGCAAAACGGATTTTCCATCCAACACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGAGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGCGGATCAGATGCGCCAGCACCAACAAGCGCCAGCAAGGCACCCGACAGGCCGAACCTGAt‑ctttcgtcttcaccagcccccttcg                                                                                                                                                                                                                            >  2:211581/1‑230 (MQ=255)                             gTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAACACGGGTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATTCGTCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGCTATTTTG‑TGTTCAcccgcccgcctagtctt                                                                                                                                                                                                                        >  2:515103‑M1/1‑234 (MQ=255)                                  ccTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGTTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTG‑TGTTCAcccgcccgcttcgtcttacgat                                                                                                                                                                                                                   <  1:152711‑M1/251‑23 (MQ=255)                                                    tCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTG‑TGTTCAcccgcccgcttcgtcttacgatattcatcataaacttt                                                                                                                                                                                                   <  1:396461‑M1/249‑39 (MQ=255)                                                     cGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTG‑TGTTCAcccgcccgcttcgtcttacgatattcatcataaactttggg                                                                                                                                                                                                >  1:9044‑M1/1‑210 (MQ=255)                                                          ggttcaTAATGGGACAAACGGATTTTCCCTCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTAACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAAAGGCAGCCGACAGGCGGTATCTGATATTTTG‑TGTTCAcccgcccgcttcgtctaacgatcttcatcttcacctttgggctca                                                                                                                                                                                            >  2:698100‑M1/3‑205 (MQ=255)                                                          gTTTCATATTGGGAAAAACGGTTTTTCCATCCCCCACTGTGGCTCCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCGGAGAAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCATCCGACAGGCCGTATCTGATCTTTTG‑TGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcat                                                                                                                                                                                           <  2:386336‑M1/251‑47 (MQ=255)                                                            ttCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTG‑TGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcata                                                                                                                                                                                          <  1:397173‑M1/250‑48 (MQ=255)                                                                 atTGGGAAAAACGGTTTTTCCATCCACCCCTGTGGAACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCAATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGGATCTGAGGAGCAGATGCGCAAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGAATCTGATCTTTTG‑TGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactga                                                                                                                                                                                      >  2:176370‑M1/1‑198 (MQ=255)                                                                 atTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTG‑TGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactga                                                                                                                                                                                      >  1:199584‑M1/1‑198 (MQ=255)                                                                 atTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTG‑TGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactga                                                                                                                                                                                      >  1:626281‑M1/1‑198 (MQ=255)                                                                   tGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTG‑TGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagct                                                                                                                                                                                  >  1:150811‑M1/1‑196 (MQ=255)                                                                    gggAAAAACTGTTTTTCCATCCCCCACTGTGGCACCCCGACATGTGTTCCATATGCCAGATCCATCGCGTTATGCTATTGTGTGGTCACCTGCTTTTTCGTCCAAAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACTGCAGCCGTCAGGCTGTATCTGATCTTTTG‑TGTTCAcccgcccgctgcgtctgacgatattcatcataaactttgtgatcatactgaagct                                                                                                                                                                                  <  2:549087‑M1/250‑56 (MQ=255)                                                                        aaaaCGGTTTTTCCATCCCCCACTGTGGCCCCCCGCCAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTG‑TGTTCAcccgcccgcttcgtcttacgatattcatcataacctttgggatcatactgaagctccccg                                                                                                                                                                             <  1:131985‑M1/251‑61 (MQ=255)                                                                             ggTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTTCCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCAGATCTTTTG‑TGTTCAccc                                                                                                                                                                                                                                      <  2:595687‑M1/189‑4 (MQ=255)                                                                             ggTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTG‑TGTTCAccc                                                                                                                                                                                                                                      >  1:595687‑M1/1‑186 (MQ=255)                                                                                      atccaCCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTG‑TGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaatctccccgccagatgcctcctt                                                                                                                                                               >  2:133488‑M1/1‑177 (MQ=255)                                                                                      atccaCCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTG‑TGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctg                                                                                                                                                               >  1:166730‑M1/1‑177 (MQ=255)                                                                                        ccaccaCTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTG‑TGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgt                                                                                                                                                              >  1:269613‑M1/1‑175 (MQ=255)                                                                                          ccccctgtgGCACCCCGACAGATGGTCCATATGCCAGACCCATCGCGGTATGCCATGGTGTGGTTACCTTCTTTTTCGTCCAGAAACGGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTG‑TGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtaga                                                                                                                                                           <  2:539713‑M1/247‑79 (MQ=255)                                                                                          accacTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTG‑TGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgt                                                                                                                                                              >  2:271526‑M1/1‑173 (MQ=255)                                                                                             acTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTG‑TGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactnaagctccccgccagatgcctcctgtagacgc                                                                                                                                                        >  2:357172‑M1/1‑170 (MQ=255)                                                                                               ggtggCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCAATTGTGGGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAAGCAGCGCAAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTG‑TGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcat                                                                                                                                                      <  2:405413‑M1/250‑84 (MQ=255)                                                                                                gtgGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTG‑TGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatc                                                                                                                                                     >  2:204629‑M1/1‑167 (MQ=255)                                                                                                        ccGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTG‑TGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgacc                                                                                                                                             >  1:67130‑M1/1‑159 (MQ=255)                                                                                                          gACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTG‑TGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgacctt                                                                                                                                           <  1:696923‑M1/251‑95 (MQ=255)                                                                                                              gATGGTCCATATGCCATAACCATCGCGTTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAATTATCTGAGGATCAGATGCGCCAGCCCCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTG‑TGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcggg                                                                                                                                       <  2:567446‑M1/251‑99 (MQ=255)                                                                                                                  gTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTCCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCTCCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTG‑CGTTCAcccgcccgcttcgtcttacgatattcatcataaactgtgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaact                                                                                                                                   <  1:406426‑M1/251‑103 (MQ=255)                                                                                                                        atGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTG‑TGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatc                                                                                                                              <  1:102663‑M1/250‑108 (MQ=255)                                                                                                                                   aTCGGGGTTTGCCATTGTGTGGTTACCTTCTTTTTCGTCCATAACCTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCCCCAATCAGCGCCAGACCGGCAGCCGACAGGCCGTATCTGATCTTTTG‑TGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagct                                                                                                                  <  1:571543‑M1/251‑120 (MQ=255)                                                                                                                                      gcgGTATGCCATTGTGGTTTTACCTTCTTTTGCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGTCAGCACCAACCATCGCCAGAACGGCAGCCTTCAGGCCGTATCTGATCTTTTG‑TTTTCGcccgcccgcttcgtcttacgatattcatcataaactttgggatcattctgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccatctata                                                                                                               <  1:206690‑M1/251‑123 (MQ=255)                                                                                                                                         gTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTG‑TGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcctcgcgcgaccttcgggaactacatagtaagccaatcgtggcgc                                                                                                            >  2:658363‑M1/1‑126 (MQ=255)                                                                                                                                           aTGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTG‑TGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatcctgacgctccccg                                                                                                                                                                             >  2:288711‑M1/1‑124 (MQ=255)                                                                                                                                           aTGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTG‑TGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcg                                                                                                          >  1:46627‑M1/1‑124 (MQ=255)                                                                                                                                           aTGCCATTGTGTGGTTACCTTCGTTTTCGGCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTG‑TGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctacccg                                                                                                                                                                             <  1:288711‑M1/184‑61 (MQ=255)                                                                                                                                                          tACCTTCTTTTTCGTCCCGAAACTGGTCCAGTATCTGAGGAGCAGATGCGCCAGCACCACTCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTG‑TGTTCCcccgcccgcttcgtcttacgatattcatcatacactttgggatcatactgaagctccccgccagatgcctcctgtagccgcatcgctcgaccttcggacgctaaatcccctttccagctctaaagcgaagccacacgaatac                                                                                           >  2:29659‑M1/1‑109 (MQ=255)                                                                                                                                                               tcttTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTG‑CGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgactcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttcccgctataaagcgaagccaaacgaatacctgcc                                                                                      >  2:2520‑M1/1‑104 (MQ=255)                                                                                                                                                                         ccAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTG‑TGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgca                                                                             >  1:282532‑M1/1‑94 (MQ=255)                                                                                                                                                                             aaaCTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTG‑TGTTCAcccgcccgcttcgtcttacgctattcatcataaactttgggataatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcacgttt                                                                        >  2:452533‑M1/1‑90 (MQ=255)                                                                                                                                                                             aaaCGGGTCAAGTATCTTAGGAGCAGTTCCGCCAGCCCCAACCAGCTCCGGAACGGCAGCCGACAGGCCGTATTTGATCTTTTG‑TGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttt                                                                        <  2:97233‑M1/251‑162 (MQ=255)                                                                                                                                                                                  gtgcaagTACCTGAGGAGCAGAGGCGCCACCACCAATCAGCGCCAGAACGGCAGCCGACAGTCCGTATCTGATCTTTTG‑TGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgttt                                                                   <  2:548497‑M1/248‑167 (MQ=255)                                                                                                                                                                                     cAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTG‑TGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgttttt                                                                 <  1:607778‑M1/250‑169 (MQ=255)                                                                                                                                                                                               aggagCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTG‑TGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaat                                                       >  1:439836‑M1/1‑72 (MQ=255)                                                                                                                                                                                                    cAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTG‑TGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaa                                                 <  1:163958‑M1/251‑185 (MQ=255)                                                                                                                                                                                                     aGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTG‑TGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaa                                                  >  1:166207‑M1/1‑66 (MQ=255)                                                                                                                                                                                                                           gcCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTG‑TGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagttttaagcatttaa                          <  1:366328‑M1/251‑208 (MQ=255)                                                                                                                                                                                                                                accggCAGCCGACAGGCCGTATCTGATCTTTTG‑TGTTCAcccgcccgctgcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagttttaagcatttaaaacac                     <  1:701368‑M1/249‑213 (MQ=255)                                                                                                                                                                                                                                   ggcggccaacaggccgtATCTGATCTGTTG‑TTTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagttttaagcatttaaaacacctt                  <  2:626281‑M1/243‑216 (MQ=255)                                                                                                                                                                                                                                   ggCAGCCGACAGGCCGTATCTGATCTTTTG‑TGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagttttaagcatttaaaacacctt                  >  1:443752‑M1/1‑36 (MQ=255)                                                                                                                                                                                                                                         cgcaagtccgtatctgATCTTTTG‑TGTTCAcccgcccgcttcttcttacgatattcatcataaactttgggatcatactaaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaagcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagttttaagcatttaaaacacctttattgt            <  2:199584‑M1/244‑222 (MQ=255)                                                                                                                                                                                                                                                ccGTATCTGATCTTTTG‑TGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagttttaagcatttaaaacacctttattgttagtcat     >  1:558375‑M1/1‑23 (MQ=255)                                                                                                                                                                                                                                                   tATCTGATCTTTTG‑TGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagttttaagcatttaaaacacctttattgttagtcataac  >  1:687586‑M1/1‑20 (MQ=255)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                    ATGCCAGCGCCTTATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGAT‑TTTTG‑TGTTCATAGATATTTATGATGAGGACGCTCGTGCTTATTGGCAGGATTTTCAATCTTAAAGGAGTACTGATGCTGCAGATAAGACTCAACTTTTTCTGACAATTTTTCTGCTACTTCCAGGAAGACTTGCCGGACGCTCCTTCTGGCTGCTGCCTCATAAAACTCCAGCGCAGCTCCTTCAACACGGTCCATGGCGACATCCAGGTCAAAAATTTCACCGTCAAAGCGTTCTTTGTC  >  NC_000913/1639796‑1640287

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