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	95.6%	137.5 / NA	46	V381I (GTC→ATC)	setA	sugar exporter SetA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (18/25);  minor base C (2/0);  total (20/26)
Fisher's exact test for biased strand distribution p-value = 1.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

ATAGCCGTATGGCGTTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTT  >  NC_000913/78525‑79002                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                  aTAGCCGTCTGGCGTTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCAAAAAGTTGGGGGCACTTTGCTGTCTCCTGGGAAATTGAGGTTATTTCTGTTCTCGCATTCttttt                                                                                                                                                                                                                                      >  1:440470/1‑250 (MQ=255)   aGCCGTCTGGCGTTGATGACGCTGCAACGTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTATCGCATTATTTTTAAc                                                                                                                                                                                                                                   <  1:389704/251‑1 (MQ=255)      cgtctGGCGTTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGc                                                                                                                                                                                                                                >  1:405768/5‑251 (MQ=255)             cGTTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAACTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGtt                                                                                                                                                                                                                         <  1:639004/251‑1 (MQ=255)                                   tAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGa                                                                                                                                                                                                    >  1:373072/1‑250 (MQ=255)                                   tAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAAGATTTAAGGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGAACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGa                                                                                                                                                                                                    <  2:603353/250‑1 (MQ=255)                                    aaCGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTCTTCAGGGAGCAATTGCACAACGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTATCGCATTATTTTTAACTGCAAAGGTTAAAGACGTTTGATGAAGTGGACg                                                                                                                                                                                                  >  2:239085/1‑251 (MQ=255)                                    aaCGCTGTATTTATAGGCATTGTTGCGGGTATTGAGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGAAGTTTGATGACGTGGACg                                                                                                                                                                                                  >  2:103774/1‑251 (MQ=255)                                            aTTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCCATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCgg                                                                                                                                                                                           >  1:461462/1‑250 (MQ=255)                                                           tgCGGGTTTTGGGATGCTATGGGTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTAATTCTGCCGGGGTAATTCTTGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGAAGTGGACGATAGCGGAAAGCCCGGTCAttt                                                                                                                                                                            <  2:4211/250‑1 (MQ=255)                                                                      ggATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAgg                                                                                                                                                                >  1:595785/1‑251 (MQ=255)                                                                       gATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAggg                                                                                                                                                               <  1:254640/251‑1 (MQ=255)                                                                          gCTATTGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCTCTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGa                                                                                                                                                             <  2:133544/250‑1 (MQ=255)                                                                             aTGGTGTCAGTATTTAAGGCCTGGAAGAGCGGTGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACGTGGCTGTCTACTGGGTAAGTGCGTTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATTACGGGGACGATAGCGGAAAGCCCGGGCATTTTACCGGTCAAGGGGAt                                                                                                                                                            <  2:351186/248‑1 (MQ=255)                                                                               ggTTTTAGGAGTTAATGCCAGGAAGAGCTGGGGCAGCTACCTCCTTATTTGCTAACAGTAGTTCTACCTGGGTAATTCTGGCTGGCGTGATTCAGGGAGCAATTGCACAAAGTTGGGGGCGCTTTGCTGTCTACTGGGTAATTGCGGTTATGTCTGTTATCTCATTAGTTTTAATCGCAAAGGGTATAGAAGTTTGATGATGTGTACGATAGAGGAAATCGCGGTCATTGGACCGGGCAAGTGGATTAAt                                                                                                                                                        <  2:689017/250‑1 (MQ=255)                                                                               ggTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCAAATGACCGGGCTAGGGGAtgaagg                                                                                                                                                       <  1:650345/251‑6 (MQ=255)                                                                                    cAGGATTTAATGCCTGGTAGAGCGGGGGCAGCTACCACCTTATCTACTAACAGTATTTCTAACGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTAGTGGGCAATTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTATCGCATTAATTTTAAACGCAACGGTTAAAGACGTTTGATGACGTGGATGATAGCGGAAAGCCAGGTCATATAACCGGGCAAGGTGTTTAATTCATa                                                                                                                                                   >  2:473427/1‑250 (MQ=255)                                                                                                   ggAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTg                                                                                                                                   <  2:512694/251‑1 (MQ=255)                                                                                                       gagCGGGGGCAGCGACCACCTTATTTACGAACAGTATTTCTGCCGGGGTAATTCTTGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGTGTAATTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTGTGATGACTTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGTGGATTAATTCATAAACGCAGGTTGTTTTGCTTc                                                                                                                               <  1:666275/251‑1 (MQ=255)                                                                                                        agCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTc                                                                                                                               >  1:708783/1‑250 (MQ=255)                                                                                                           gggggCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATa                                                                                                                            >  2:328144/1‑250 (MQ=255)                                                                                                             tggCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAAGTCTGGCTGGCGTTATTCAGGGAGCATTTGCACAAAGTTGGGGGCATTTTGCTGTCTACTGGGTAATTGTGGTTATTTCTGTTATCGCATTTTTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGAc                                                                                                                                                                                                   <  2:502805/176‑1 (MQ=255)                                                                                                             gggCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGAc                                                                                                                                                                                                   >  1:502805/1‑177 (MQ=255)                                                                                                                             aTTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGcgtcgtc                                                                                                          <  1:328144/250‑1 (MQ=255)                                                                                                                                 aCTAACAGTATTTCTACCGGGGTAATTCTGGCTGTCGTTATTCAGGGAGCAATTGCACAAATTTGGGTGCTCTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCTCGGTCATTTGACCGGGCAAGGGGATTAATTCAGAAACGCATGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGtgct                                                                                                     <  1:249447/251‑1 (MQ=255)                                                                                                                                   tAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGtgctgc                                                                                                   >  1:329523/1‑251 (MQ=255)                                                                                                                                     aCAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCACAGGCGTCGTCGTGCTGCa                                                                                                  >  2:335972/1‑250 (MQ=255)                                                                                                                                          aTTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTa                                                                                             >  2:145026/1‑250 (MQ=255)                                                                                                                                                aCCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCa                                                                                       <  1:103774/250‑1 (MQ=255)                                                                                                                                                aCCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTAGCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGcc                                                                                         <  2:530143/248‑1 (MQ=255)                                                                                                                                                     ggTAATTATGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCAGGTAATTTGACCGGGCAAGGGGCTTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGACGAGCGGCAAGGTAAGCCCAAtcatg                                                                                 >  2:605645/1‑247 (MQ=255)                                                                                                                                                     ggTAAATCTGGCTGGCGTTATTCAGGGAGCCATTTCCCAAAGTTGTGGGCTCTTTGCTGTCTACGGGGTAATTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGCCGTTTGATTACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGACGTCGTCGTGCTGCAAGGTAAGCCCAATACTg                                                                                 <  2:530888/251‑1 (MQ=255)                                                                                                                                                           tCTGGATGGCGTTATTCAGGGAGCAATTGCTCACAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTATCGCCTTATTTTTAAATGCAAAGGTTACAGACGTTTGATGACGTGGCCGTTAGCGGAAAGCCCGGTAATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTTGTGCTGCAAGGTACGCCCAATACTGtac                                                                              >  2:641211/1‑246 (MQ=255)                                                                                                                                                                    cGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCaa                                                                                      <  1:667255/231‑1 (MQ=255)                                                                                                                                                                    cGTTATGCAGGGAGCAATTGCAGAAAGTTGGGGACACTTTGCTGTCTACTGGGACATTGCGGTTATTTCTGTTCTCGCATTCTTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGAAAGGAAAGccccaagatcggaagggaggcgtgg                                                                  >  2:667255/1‑229 (MQ=255)                                                                                                                                                                        aTTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGATCGTTCATCATGc                                                              <  1:481458/251‑1 (MQ=255)                                                                                                                                                                                         agaaaGTTGGGGGCACTTTGCTGTTTACTGGGTAATTGCGGTTAGTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGAAGTGGCGGCGGAAGGCa                                             <  2:570446/249‑1 (MQ=255)                                                                                                                                                                                                      cACTTTGCTGTCTACTGGGCAATTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTACAGACGTTTGCTGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGCTTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGCCCGTTCATCATGCAGTGGCGGCGGGAGGCATCgtagggaaagc                                >  2:296023/1‑241 (MQ=255)                                                                                                                                                                                                       aCTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCg                               <  1:296023/251‑1 (MQ=255)                                                                                                                                                                                                            gcattCTACAGGGTAATTGCGGTTATTTCGGGTAGCGCATTATTGTTAACCGCAAATGTTAAAGACGTGTGTGGACGTGGACGATAGCGGAAAGCTCGGTCATTTGACCTGTCAAGGGGATTAATTCATAAACGCAGGTGGTTTTGCGTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACGGTCCAAACTGTTTAGCATGCAGTGTTTGCTGAATGCGTCGATGGTAAAGCGATAgg                          <  2:89697/247‑1 (MQ=255)                                                                                                                                                                                                            gCTGTCTACTGGGTAATTGCTGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTTACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAgg                          <  1:40859/251‑1 (MQ=255)                                                                                                                                                                                                            gCTGTCTACTGGGTAATTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGt                                                                                >  1:426939/1‑197 (MQ=255)                                                                                                                                                                                                            gCTGTCTACTGGGTAAGTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGGCGTCGTGCTGCAAGGTAAGCCCAATACTGt                                                                                <  2:426939/197‑1 (MQ=255)                                                                                                                                                                                                               gTCTACTGGGTAATTGCGGTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGt                         <  1:100761/249‑1 (MQ=255)                                                                                                                                                                                                                                ggTTATTTCTGTTATCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCTCTTTcccc      >  1:658333/1‑248 (MQ=255)                                                                                                                                                                                                                                    aTTTCTGTTTGCGTATTATGTTTAACCGCAAAGGTTAAAGACGTCTGATGACGTGGACGATAGCGGAAAGCCCGGTCATGTGACCGGGCAAGGGGATTAATTCATAAACTCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTTGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTTGCGGCGGAAGGCAGCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCtt  <  2:34123/251‑1 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                  ATAGCCGTATGGCGTTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCGCATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCTT  >  NC_000913/78525‑79002

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