Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	78,761	G→A	V381I (GTC→ATC)	setA →	sugar exporter SetA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	78,761	0	G	A	98.6%	212.4 / NA	72	V381I (GTC→ATC)	setA	sugar exporter SetA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (34/37);  minor base T (0/1);  total (34/38)
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.45e-01

TTCTTTAATAGCCGTATGGCGTTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAA‑TTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTC  >  NC_000913/78518‑79000                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                ttatttcaTAGCCGTCTGGTGTTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGGTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGGTATTCAGTGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGGCTACTGGGTAA‑TTGCGGTTATTTCTGTTATCGCat                                                                                                                                                                                                                                           <  1:913607/243‑1 (MQ=255)  tatttcatagttgtcttgCGTTGATGACGCTGTAACTTGTTAACGCTGTATTTAGCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTTAGGATTTAATGCTTGGAAGAGCGGGGGCATCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCGGGCGTTACTCAGGGAGCAATTGCATAAAGTTGTGGGCACTTTGCATTCTAATGGGTAAGTTG‑GGTTATTTCTGTTTTCGCatt                                                                                                                                                                                                                                          <  2:111820/233‑1 (MQ=255)  tatttcaTAGCCGTCTGGCGTTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAA‑TTGCGGTTATTTCTGTTATCGCatt                                                                                                                                                                                                                                          <  2:194276/244‑1 (MQ=255)       caTAGCCGTCTGGCGTTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAA‑TTGCGGTTATTTCTGTTATCGCATTAttttt                                                                                                                                                                                                                                    >  1:898630/2‑251 (MQ=255)          aGCCGTCTGGCGTTGATGACGCTGCACCTTTTTAACGCTGTATTTATCGGCATTGATGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGCAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAA‑TTGCGGTTATTTCTGTTATCGCCTTATTTTTAAc                                                                                                                                                                                                                                 >  1:857693/1‑251 (MQ=255)          aGCCGTCTGGCGTTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGATATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGTTGTCTACTGGGTAA‑TTGCGGTTATTTCTGTTATCGAATTATTTTTAAc                                                                                                                                                                                                                                 >  2:802764/1‑251 (MQ=255)                  tgCGTTTATGACGCTTCAACGTTTTAACGTTGTATTTTTCGGCATTTTTGCGGGGATTTGGATGCTTTGTATTCCGGATGTAATGTCTGTAAGTGCTGGTGCAGCTTCCCCCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTTGTTGGCGTTATTCATGGTGCATGTTCACAATGTTGTGGGCACTTTGCTGTCTACTGGGTAA‑TTGCGGTTATCTCTGTTATCGCATTATGTTTAACCGCTAAg                                                                                                                                                                                                                          <  2:544481/249‑1 (MQ=255)                    ggTTGATGACGCTGCAACTTTTTAACGTTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAA‑TTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGtt                                                                                                                                                                                                                       <  2:43112/250‑1 (MQ=255)                        gatgaCACTGCAACTTGTTAACTCTGTATTGAACGGCAGTTTTGCGGGTATTGGGAAGCAATGGTTGCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATGTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGTGCACTTTGCTGTCTACTGGGTAA‑TTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTaaa                                                                                                                                                                                                                    <  1:890898/250‑1 (MQ=255)                         atgaCGCTTCCACTTGTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTAACACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAA‑TTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGa                                                                                                                                                                                                                  <  1:232087/251‑1 (MQ=255)                                 gCAACTTTTTTACGCGGTATGTATCTGCTATGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGTCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTGTGGCTGGCGTTATTCAGGGAGCAATTTCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAA‑TTTCGGTTATTTCTGTTATCGCATTATTTTTAACTGCAAAGGTTAAAGACTTTtgat                                                                                                                                                                                                          <  2:949142/251‑1 (MQ=255)                                    aCTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAA‑TTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGAc                                                                                                                                                                                                       >  1:405516/1‑251 (MQ=255)                                      tttttAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGCGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAA‑TTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGTTTAAGGACGTTTGATGACGt                                                                                                                                                                                                     >  2:181758/1‑251 (MQ=255)                                           aaCGCTGTATTTATCGGCATTGTTGCGGGGATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAA‑TTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACg                                                                                                                                                                                                <  2:187131/251‑1 (MQ=255)                                              gCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTTGAAGAGCGGGGGCAGCTTCCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAA‑TTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATa                                                                                                                                                                                             <  1:433770/251‑1 (MQ=255)                                                tGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAA‑TTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTACAGACGTTTGATGACGTGGCCGATAgg                                                                                                                                                                                           >  2:962487/1‑250 (MQ=255)                                                      tATCGGCAATGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAA‑TTGCGGTTATTTCTGTTATCGCATTAGTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGaaa                                                                                                                                                                                      <  2:906667/250‑1 (MQ=255)                                                       ataggCATTGGTGCGTGTATTGAGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGATACCACCTTATTTACTAACAGATTTTCTACCGGGGGAATTCTGGCCGGCGTTATTCAGGAAGAAATTACACAAAGTTGGGGTACCTTTGCTGTATATTGGGTAA‑TTGCGGTTCTTTCTATTATCGAATTCTTTTTAACCGCAAAGGTTAAAGAAGTTTGTTGATGAGGACGCTATCGGAAAGc                                                                                                                                                                                    >  2:527773/4‑251 (MQ=255)                                                          ggCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAA‑TTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCg                                                                                                                                                                                 >  1:84541/1‑251 (MQ=255)                                                                       gtagtGGGATGCTATGGTTTCATGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCTCCTTATGTTCTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAA‑TTGCGGTTATTTCTGTTATCTCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCgg                                                                                                                                                                    <  2:265185/247‑1 (MQ=255)                                                                       gTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAA‑TTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCgg                                                                                                                                                                    >  1:669078/1‑251 (MQ=255)                                                                        tATTGGGATGCTATGGGTTCAGGATTTAATGCCTGGAAGAGCTGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAACGTTGGGGGCACTTTGCTGTCTACTGGGTAA‑TTGCGGTTATTTCTGTTATCGCATTATTTTTCACCGCAACGGTTAAAGACGTTTGATGACGTGGACGAGTGCGGAAAGCCCGGTCATTTGACAggg                                                                                                                                                                   >  2:655617/1‑251 (MQ=255)                                                                           tGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAA‑TTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCa                                                                                                                                                                 <  1:696081/250‑1 (MQ=255)                                                                             ggATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAA‑TTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAgg                                                                                                                                                              <  1:474392/251‑1 (MQ=255)                                                                             ggATGCTATGGTTTAAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAA‑TTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAgg                                                                                                                                                              <  2:351126/251‑1 (MQ=255)                                                                              gATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGAAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAA‑TTGAGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAggg                                                                                                                                                             >  1:57269/1‑251 (MQ=255)                                                                                 gCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAA‑TTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGAt                                                                                                                                                          <  1:659957/251‑1 (MQ=255)                                                                                 gCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAA‑TTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGAt                                                                                                                                                          >  1:334953/1‑251 (MQ=255)                                                                                 gCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGCCTACTGGGTAA‑TTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGAt                                                                                                                                                          >  2:295691/1‑251 (MQ=255)                                                                                                   aattccTGGTAGAGCGGGGGTAGTTACCATCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCTTTATTCAGTGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGAGTAA‑TTGCGGTTATTTCTGTTATCGCATTATATTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAgg                                                                                                                                         <  1:667593/246‑1 (MQ=255)                                                                                                    aTGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGAAATTCTGGCTGGCGTTGTTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAA‑TTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGCTGACGTGTACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAgtgt                                                                                                                                       >  2:954441/1‑248 (MQ=255)                                                                                                    aTGCCTGGAAGAGCGAGGGCAGCTCCCACCTTATTTACTAACAGTATTTCTACCGGGTGAATTCTGGCTGGCGTTATTCAAGGAGCAATTGCACAGAGTTGGGGGCACTTTGCTGTCTACTGGGTAA‑ATGCGTTTATTTCTGTTATCGCATTCTTTTTAACCGCAAAGGTTAAATACGTTTGATGGCGTGGACGATAGCGGCAAGCCCAGTCTTTATACCGGGCCGGGGGATTCATTAATAAAAGCCGgtt                                                                                                                                       >  2:30916/1‑251 (MQ=255)                                                                                                      gCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAA‑TTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGgttgt                                                                                                                                     <  1:490176/251‑1 (MQ=255)                                                                                                      gCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGCATTTCTACCGGGGTAATTCTGGCCGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGAACTTTGCTGTCTACTGGGTAA‑TTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAACGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCACGGTCATTTGACCGGGCAAGGGGATTAATTCATACACGCAGgtgtt                                                                                                                                     >  2:376425/1‑248 (MQ=255)                                                                                                       ccTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAA‑TTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGgttgtt                                                                                                                                    <  1:320573/251‑1 (MQ=255)                                                                                                          agAAGAGAGGGGGCAGCAACCACATTATTTACTAACAGTATTTCTACCGGGGGAATTCTGGCTGGCGGTATTCAGGGAGAAATTGCACACAGTTGGGGGCGCTTTGCTGTCTACTGGGTAA‑TTGCGGATATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGTTGACGGGGACGATAGCGGTAAGCCCGGGCATTTGACCGGGCAAGGGGATTAATTCATAAAAGCAGgttgttctg                                                                                                                                 >  2:762068/2‑248 (MQ=255)                                                                                                           gAAGAGCATGGGCTGCTACCAACTGATTTACTAACAGAATTTAAAAAGGGGTAATTCTGGCTGGCGTTAATAAGAGAGCATTTGCACCAAGTTGGGGGCACTTTGCTGTCTACTGGGTAA‑TTGCGGTTATTTCTGTTATCGCATTATATTAAAACGAACACGTTACAGCAGATTGATGACGTGGACGCTAGCGGAAAACCCGGTCCTTTGCCCAGACAAGGGGTATAGTGCATAATCGCAGGGTGtttttc                                                                                                                                >  2:959083/1‑249 (MQ=255)                                                                                                            aaGATCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAA‑TTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCt                                                                                                                               <  1:523781/251‑1 (MQ=255)                                                                                                                gCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAA‑TTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCAt                                                                                                                           <  1:655617/251‑1 (MQ=255)                                                                                                                    gggCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAA‑TTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGc                                                                                                                       <  1:296131/251‑1 (MQ=255)                                                                                                                         gCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAA‑TTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCa                                                                                                                   >  2:192978/1‑250 (MQ=255)                                                                                                                                      ttACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAA‑TTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGTcgtcgtcgt                                                                                                      <  1:328308/250‑1 (MQ=255)                                                                                                                                               gTATGTCTACCTGGGGAATTCTGGCTGGCGTTATTCAGGGAGTAAGTGCTGAAATTTTGGAGCACTTTGCTGTCTATTGGGTAA‑ATTAGGTTATTTCTGTTAGAGCATTATTTTTAATCTCAAAGGATAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGATTATTTGACCGGGCAAGGGGAGTAATGCATAAACACAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGt                                                                                            <  2:487652/251‑1 (MQ=255)                                                                                                                                                tATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAA‑TTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGt                                                                                            >  1:366184/1‑250 (MQ=255)                                                                                                                                                     cTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAA‑TTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGACGTGCTGCAAGGTAAAccc                                                                                      >  2:543920/1‑251 (MQ=255)                                                                                                                                                            ggTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAA‑TTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCCTAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTg                                                                               >  2:716355/1‑251 (MQ=255)                                                                                                                                                               aagtCTGGCTGGAGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCTATTTGCTGTCTACTGGGTAA‑TTGCGTTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGGGCTGCAAGGTAAGCCCAATACTGTc                                                                             <  1:30916/247‑1 (MQ=255)                                                                                                                                                               aaTTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAA‑TTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGGAAGCCCAATACTGTcc                                                                            >  2:384151/1‑251 (MQ=255)                                                                                                                                                                    gggctggcGTTATTCAGGGAGCTATTTCCCAAAGTTTGGGGTACTTTGCTGTCTACTGGGTAA‑TTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTGTGATGACGTTGACGATAGCGGAAAGCCCGGTCAGTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTTTCCAGAcc                                                                       <  2:698149/250‑1 (MQ=255)                                                                                                                                                                     ggctggcGTTATTCAGGAAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAA‑TTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAACCCGGTAATTTGACCGGGCACGGGGGTAAATTCATAAACACAGGTTGTTTTGCTTCATAAGCAGCAATGGTGTAGTAGTGGTGCAAGGTACGCCCAATAAGGTCCAGACCg                                                                      >  2:770948/1‑251 (MQ=255)                                                                                                                                                                           cGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAA‑TTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGTAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTGAGCCCAATACTGTCCAGACCGTtcatc                                                                >  2:679301/1‑251 (MQ=255)                                                                                                                                                                               aTTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAA‑TTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATg                                                             <  1:181758/250‑1 (MQ=255)                                                                                                                                                                                    gggagacatagcaacaaGTGGGGGACACGTTTCTGTCTACAGGGAAA‑TTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGATAAAGACGTTTGATGACGTGGACGATAGCGGAATGCCCGGTCCTTTGACCGGTCACGGGGATTAATTACTAAACGCCGGTTGTTTTGATTCATAAGCGGCAAAGGCGTAGTCGTTCTGCAAGGTACACCCCAGACTGTCCAGACCGTccatcctgactagg                                                       >  2:944744/16‑242 (MQ=255)                                                                                                                                                                                      gAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAA‑TTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTggcg                                                     <  1:578226/251‑1 (MQ=255)                                                                                                                                                                                              gcaacaaGTTGGGGGCCATTTGCTGTCTACTGGGTAA‑TTGTGGTTATTGCGGTTATCGCAGTATTTGTAACCGCAAAGGTTAAAGACGGTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGTGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATATGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGTCCGTTCATCATGCAGTGGCGTCGGAAg                                              <  2:41737/245‑1 (MQ=255)                                                                                                                                                                                              gCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAA‑TTGCTGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAgg                                             <  1:619117/251‑1 (MQ=255)                                                                                                                                                                                              gCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAA‑TTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCATGCCGTTCATCATGCAGTGGCGGCGGAAg                                              <  1:376425/250‑1 (MQ=255)                                                                                                                                                                                                acaAAGTTGGGGGCCCTTTGCTGTCTACTGGGTAA‑TTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGc                                            >  1:761099/1‑250 (MQ=255)                                                                                                                                                                                                  aaaGTTGGGGGCACTTTGCTGTCTACTGGGTAA‑TTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATc                                         <  1:192978/251‑1 (MQ=255)                                                                                                                                                                                                              aCTTTGCTGTCTACTGGGTAA‑TTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGcgtcg                                                                                                          <  1:287950/174‑1 (MQ=255)                                                                                                                                                                                                              aCTTTGCTGTCTACTGGGTAA‑TTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGcgtcg                                                                                                          >  2:287950/1‑174 (MQ=255)                                                                                                                                                                                                                tttGCTGTCTACTGGGTAA‑TTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGATAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATACACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCACGGTACGCCCAATACTGTCCAGACCGCTCATCATACAGTGCCGGCGGAAGGCATCGATGGGCAAGCGa                            >  2:586853/1‑250 (MQ=255)                                                                                                                                                                                                                   gCTGTCTACTGGGTAA‑TTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCACGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCCTCGATGGTAAAGCGATAgg                        >  1:53284/1‑251 (MQ=255)                                                                                                                                                                                                                          aCTGGGTAA‑TTTCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTtctc                 <  2:787733/251‑1 (MQ=255)                                                                                                                                                                                                                          aCTGGGTAA‑TGGCGGTTATTTCTGATATCGCATGATTTTTAAACGCAAAGGTTAAAGACATTTGATGACGTGGACGATGGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTGGTTTTGCTTCATAAGCGTCAATGGCTTCGTCGTGCTGTAAGGTAAGCGCAATACTGTCCAGACCGTTCATCTTGCAGTGGCGGCGGAAGTCATCGATTGTAAAGCGATAGGTTTtctc                 <  2:365023/251‑1 (MQ=255)                                                                                                                                                                                                                             gggTAA‑TTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGCACGTTCATCATGCAGAGGCGGCGGAAGGAATAGATGATAAAGCGATAGGTTTTTTCTcc              >  2:792423/1‑251 (MQ=255)                                                                                                                                                                                                                              ggTAA‑TTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTccc             <  2:898630/251‑1 (MQ=255)                                                                                                                                                                                                                                 aa‑TTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGGCTCGATGGGAAAACGAATGGTTTTTTCTCCCGcc          >  2:148537/1‑250 (MQ=255)                                                                                                                                                                                                                                        ggTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGACCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCAAATACTGTCCAGACCGTTCATCATGCAGTTGCGGCGGAAGGCACCGATGGTAAAGAGAGTGGTTTTTTCTCCCCCTTTCAcc    >  2:863553/1‑251 (MQ=255)                                                                                                                                                                                                                                           tATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACctc  <  1:641228/250‑1 (MQ=255)                                                                                                                                                                                                                                              ttCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGAt                                      <  2:424000/211‑1 (MQ=255)                                                                                                                                                                                                                                              ttCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGAt                                      >  1:424000/1‑211 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                TTCTTTAATAGCCGTATGGCGTTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAA‑TTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTC  >  NC_000913/78518‑79000

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