Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,937,909	G→A	A206T (GCG→ACG)	rbsK →	ribokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,937,909	0	G	A	98.0%	159.0 / NA	51	A206T (GCG→ACG)	rbsK	ribokinase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (19/30);  total (20/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

TGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAGCGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGGTGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGT  >  NC_000913/3937664‑3938150                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  tggtggAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGa                                                                                                                                                                                                                                                >  2:428804/1‑249 (MQ=255)       aaGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGCAGATGCAACGAAGgtggcg                                                                                                                                                                                                                                        >  1:302946/1‑251 (MQ=255)           gcAACGTGAGCGTATTGCCAACGCGTCAGCATTATGAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCAGCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAgg                                                                                                                                                                                                                                    <  2:373962/251‑1 (MQ=255)                gTGATCGTATGGACAACGCGACAGCATTATGAATGCAGTTTGAACGACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTGGGGCTTAACCCGGCGCCTGCTCTCGAACTGCCTGACGAACTGGTTTCGCTGGTGGACATTTTTACGCCTAGCGAAACGGAATCAGAAAAGCTCTCCGGTATTCGTGTTGAAAATGATGAAGATGCAGCGGAggtgtcacagttactg                                                                                                                                                                                                                               <  2:414431/251‑15 (MQ=255)                  gAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGc                                                                                                                                                                                                                              >  1:344009/1‑250 (MQ=255)                        aTTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGaaaa                                                                                                                                                                                                                       <  1:272195/251‑1 (MQ=255)                             cAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTAGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTTTTACGCCAAACGAAACGGAAGCGGAAAAGATCACCGGTATTCGAGTTGAAAATGATGAAGATGCAACGCAGGCGGCGCCGTAACTGAATGAAAAAGGtg                                                                                                                                                                                                                  >  1:2848/1‑250 (MQ=255)                              aaCGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTTGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGGGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTAt                                                                                                                                                                                                                 <  1:280667/251‑1 (MQ=255)                                      aGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCgtactg                                                                                                                                                                                                         <  2:152922/251‑1 (MQ=255)                                                             caccacTCGACAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGATTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGATCACCGGTATTAGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGCAGTCgg                                                                                                                                                                                  >  2:405904/1‑250 (MQ=255)                                                                tacTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCCGGTGGACATTATTACGCACAACGAAACGGAAGCAGAAAAGCTCACCGGAATTCGTGTTGAAAATGATGAAGATGCAACGAAGGAGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTAATTTAGGAAGTCggggt                                                                                                                                                                               >  2:343870/2‑251 (MQ=255)                                                                 acTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCgtggtg                                                                                                                                                                              >  2:301818/1‑251 (MQ=255)                                                                  cTCGACAGTGTGATGGCAGCGGCGAACATCGCCCATCAAAATAAGCCTATCGTTGCGCTTAACCCGGCTCCGGCACGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAACTGATGAAGATGCAACGAAGGCGGCGCAGGAACTGAATGAAAAAGGTATCCGTACTGTACTGTTTACTTTAGGTAGTCgtggtgg                                                                                                                                                                             >  2:399258/1‑250 (MQ=255)                                                                           gtgATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTa                                                                                                                                                                     <  2:306060/250‑1 (MQ=255)                                                                                         gAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAg                                                                                                                                                      <  1:301818/251‑1 (MQ=255)                                                                                          aaaaTCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATGATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCTTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGGACTTCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGATGTTGTGGGGTATGGGCGAGCGTGAATGATGAAgg                                                                                                                                                     <  2:399392/251‑1 (MQ=255)                                                                                                  ccATCAAAAGAAGTCTAGCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTTACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGcgcg                                                                                                                                             <  2:379120/251‑1 (MQ=255)                                                                                                           tAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGAATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGAt                                                                                                                                    >  2:398251/1‑251 (MQ=255)                                                                                                               aCTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATGACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGGCAGCGCGTTCCTGGATTCCg                                                                                                                                <  1:45919/251‑1 (MQ=255)                                                                                                                atATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGGGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTAGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGGCGTGGTGGATGGGCTAGCGTGAATGATGAAGGTAAGCGCGTTCCTGGATTCCgg                                                                                                                               >  2:245157/2‑251 (MQ=255)                                                                                                                  aTCGTTGCGCTTAACCCGGCTCCTGCTCGCGAACGTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGGATGGGCTAGCGGGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGt                                                                                                                             <  2:28810/251‑1 (MQ=255)                                                                                                                    cGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGc                                                                                                                           >  2:97415/1‑251 (MQ=255)                                                                                                                     gTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGc                                                                                                                           >  1:138501/1‑250 (MQ=255)                                                                                                                                cccGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATa                                                                                                               <  1:27513/251‑1 (MQ=255)                                                                                                                                 ccggctccggctCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATAc                                                                                                              >  1:150411/1‑251 (MQ=255)                                                                                                                                   ggctccggcacGCGAACTTCACGCCGCACTGCTGGCGCTGGTGGGCATTATTAGGCCAAACGAAACGGACGCAGAAAAGATTACCGGTATTCGAGTTGAACATGATGAAGATGCACAGAAGGCGGCGCAGGTACTGCATGAAAAAGGAATCCGTACTGTACTGATTACATTCGGTAGTCGTGGTGTATGGGCTAGCGTGCATGCTGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCGTGCTGTCGATAcc                                                                                                             >  2:214790/1‑250 (MQ=255)                                                                                                                                       ccggctcGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCAtt                                                                                                          <  1:97415/249‑1 (MQ=255)                                                                                                                                        cggctcGCGAACTTGCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATtgct                                                                                                       <  2:344927/251‑1 (MQ=255)                                                                                                                                                gAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGACGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGGGCAGGCTGTCGATACCATTGCTGCCgg                                                                                                  >  2:153265/1‑248 (MQ=255)                                                                                                                                                  aCTTCCTGTCGAACTGCTTGCGCTGGTGGACATTATTTCGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGCATTCGTGTTGAAAATGATGAAGATGCTACGAAGGCGGTGCAGGTACAGCATGAAAAAGGTATCCGGACTGTACTGATTACGTTAGGAAGTCGTTGTGTATGGGCTAGCGTGAATGATGAAGGTCAGTGCTGTCCTGGATTCCGGGTTCAGGCTGTCGATACCATTGCTGCCGGAGATa                                                                                             <  1:259547/251‑1 (MQ=255)                                                                                                                                                  aCTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATa                                                                                             <  2:122885/251‑1 (MQ=255)                                                                                                                                                      ccTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAGCGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGTCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCtt                                                                                         <  1:235629/251‑1 (MQ=255)                                                                                                                                                          aCGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAAc                                                                                     <  2:344009/251‑1 (MQ=255)                                                                                                                                                                   tgGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGtt                                                                             >  1:55251/1‑250 (MQ=255)                                                                                                                                                                    ggCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTa                                                                            <  2:55251/250‑1 (MQ=255)                                                                                                                                                                       gcTGGTGTACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATGCGTGTTGAAAATGATGAAGATGCAACGAATTCGGCGCAGGTTCTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCa                                                                        <  1:371318/251‑1 (MQ=255)                                                                                                                                                                          ggtggACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACgg                                                                     <  2:368612/251‑1 (MQ=255)                                                                                                                                                                             ggACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCac                                                                  <  1:191608/251‑2 (MQ=255)                                                                                                                                                                             ggACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCAt                                                                  >  1:352398/1‑251 (MQ=255)                                                                                                                                                                                    attaCGCCAAATGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTgg                                                            >  1:362002/1‑250 (MQ=255)                                                                                                                                                                                    attaCGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTgg                                                            <  1:167227/250‑1 (MQ=255)                                                                                                                                                                                      taCGCCAAACCAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGGGGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTAATGCATGAAAAAGGTAGCCGTACTGGACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGGGAATGATGAAGGTCAGTGCGTTCCTGGATTCCGGGTGCAGGCTTTCGATACCATTGCTGCCGGGGATTCCTTTAACGGTGAATTAATCACgtgacggctcgaa                                                          <  2:44814/250‑13 (MQ=255)                                                                                                                                                                                                        aGCAGAAAAGCTCACCGGTATTCGTGTTGACAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTGTCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAgg                                       >  2:117105/1‑251 (MQ=255)                                                                                                                                                                                                                         ggaTTCGGGTTGAAAATGATGAAGATGCAACGAAGGCGTCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGAGTACTTTAGGAAGTCGTGGTGTATGGGCGAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGTGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCAGTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCAt                      <  1:27705/249‑1 (MQ=255)                                                                                                                                                                                                                             tCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCAtgct                   <  2:144751/250‑1 (MQ=255)                                                                                                                                                                                                                               gtgtTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTCCTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGACGGACAGCGCGCTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCGTGCTGcc                >  1:186361/1‑251 (MQ=255)                                                                                                                                                                                                                               gtgtTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGTCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGcc                <  2:149981/251‑1 (MQ=255)                                                                                                                                                                                                                                   tGAAAATGATGAAGATGCAACGAAGGCGGAGCAGGTACTGCATGAAAAAGGTATCCGTTCTGTACTGATTACTTTAGGAAGTCGTTGTGTATGGGCTAGCGTGAATGAGGAAGGTCAGCGCGGTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCGTTAACGTTGCGGTAATCACGTCATTGCTTGAATAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGcctct             <  1:193368/250‑4 (MQ=255)                                                                                                                                                                                                                                          gatgaAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTg      <  1:347644/250‑1 (MQ=255)                                                                                                                                                                                                                                             gaAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGt  <  1:153265/251‑1 (MQ=255)                                                                                                                                                                                                                                             gaAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGt  <  1:90685/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  TGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAGCGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGGTGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGT  >  NC_000913/3937664‑3938150

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