Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,848,301	0	C	T	79.3%	116.8 / 28.8	59	G126R (GGG→AGG) ‡	ydjA	putative oxidoreductase YdjA
Reads supporting (aligned to +/- strand):  ref base C (4/8);  new base T (25/21);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 3.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.95e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GATAAGTGAGAGTGTCTGAATTCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGC‑TGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGAAACGCTCGCGCCCTTCCCCTTCAATCACAAAAAAATGCCACGGTTGC  >  NC_000913/1848061‑1848551                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                            gagaaGTGAGAGTGTCTGAATTCCTGCGCCTTTGCGCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAAct                                                                                                                                                                                                                                                   <  2:171670/248‑1 (MQ=255)            gtgtCTGAATTCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGC‑TGCCATTTg                                                                                                                                                                                                                                       >  1:8707/1‑251 (MQ=255)            gtgtCTGAATTCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGC‑TGCCATTTg                                                                                                                                                                                                                                       <  1:296643/251‑1 (MQ=255)                  gAATTCCTGCGCCTTTGCGCACAATCCATACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTTC‑TGCCATTTGCATCGc                                                                                                                                                                                                                                 <  2:54401/251‑1 (MQ=255)                   aaTTCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGC‑TGCCATTTGCATCGcc                                                                                                                                                                                                                                >  2:226743/1‑251 (MQ=255)                    aTTCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGC‑TGCCATTTGCATCGcc                                                                                                                                                                                                                                >  1:59704/1‑250 (MQ=255)                      tCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGC‑TGCCATTTGCATCGCCAt                                                                                                                                                                                                                              <  2:8707/250‑1 (MQ=255)                            gcCTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGGGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGC‑TGCCATTTGCATCGCCATGACcgcg                                                                                                                                                                                                                       >  1:102210/1‑251 (MQ=255)                             ccTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGC‑TGCCATTTGCATCGCCATGACCgcgc                                                                                                                                                                                                                      >  1:59406/1‑251 (MQ=255)                                     cacaATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGC‑TGCCATTTGCATCGCCATGACCGCGCATCCGGcc                                                                                                                                                                                                              >  2:152536/1‑250 (MQ=255)                                     cacaATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGC‑TGCCATTTGCATCGCCATGACCGCGCATCCGGCa                                                                                                                                                                                                              <  2:20388/251‑1 (MQ=255)                                       caagccAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGC‑TGCCATTTGCATCGCCATGACCGCGCATCCGGCAGa                                                                                                                                                                                                            <  2:102210/247‑1 (MQ=255)                                           ccAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGC‑TGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACAtt                                                                                                                                                                                                        <  1:181530/251‑1 (MQ=255)                                             aGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGC‑TGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTc                                                                                                                                                                                                      <  1:405356/251‑1 (MQ=255)                                                   tttCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGC‑TGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTc                                                                                                                                                                                                >  1:27625/1‑251 (MQ=255)                                                      cgcgACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGC‑TGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCa                                                                                                                                                                                             <  1:56113/251‑1 (MQ=255)                                                             aTTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGC‑TGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCggg                                                                                                                                                                                      >  1:66178/1‑251 (MQ=255)                                                                aTCAGAAATAAGTTACAAACGGCGTCGGGTCCTGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAAGTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTCCCGGACTTTCAGTTAATGCTCCACTGCACCCGATGCCGCCAAACTGCTGGGCATCTGC‑TGCCATTTGCATCGCCATGACCCCGCATCCGGCAGGCATTTCCTGTACCCCGCGCGGGAc                                                                                                                                                                                    <  2:213698/250‑1 (MQ=255)                                                                aTCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGC‑TGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGAc                                                                                                                                                                                    <  2:113828/250‑1 (MQ=255)                                                                  cAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGC‑TGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACttt                                                                                                                                                                                 <  2:27625/251‑1 (MQ=255)                                                                  cAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGC‑TGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACttt                                                                                                                                                                                 >  2:45429/1‑251 (MQ=255)                                                                       aTAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGC‑TGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGAt                                                                                                                                                                            >  1:339343/1‑251 (MQ=255)                                                                         aaGTTACACACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGC‑TGCCATTTGCATCGCCATGACGGCGCATCCGGCAGACATTTACTGTTCCCAGCGCGGGACTTTATGATtc                                                                                                                                                                          >  2:333359/1‑251 (MQ=255)                                                                              aCAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGC‑TGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTc                                                                                                                                                                     <  2:362989/251‑1 (MQ=255)                                                                                 aaCGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGC‑TGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGc                                                                                                                                                                   >  1:107014/1‑250 (MQ=255)                                                                                 aaCGGCGTCGGGTCCGGGACGGTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAAATTCTGGGCAACTGC‑TGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGc                                                                                                                                                                   >  2:407249/1‑250 (MQ=255)                                                                                         cGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCGCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGC‑TGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGcc                                                                                                                                                           <  2:343556/250‑1 (MQ=255)                                                                                          gggTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGC‑TGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGccac                                                                                                                                                         >  1:303887/1‑251 (MQ=255)                                                                                               cGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGC‑TGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCg                                                                                                                                                    >  1:42035/1‑251 (MQ=255)                                                                                                  gACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCTCGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGC‑TGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGtga                                                                                                                                                 <  1:466927/251‑1 (MQ=255)                                                                                                            gACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACTGCAACCGAATGCGTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAAC‑GCGTGCCATGTGCATCGCCATGACCGCTCATCCTGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCgg                                                                                                                                       <  1:458826/251‑1 (MQ=255)                                                                                                                        ttCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGACTGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGC‑TGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCg                                                                                                                                                    >  1:69894/1‑226 (MQ=255)                                                                                                                        ttCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGC‑TGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCg                                                                                                                                                    <  2:69894/226‑1 (MQ=255)                                                                                                                            gCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGC‑TGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGC‑ACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCAt                                                                                                                      <  2:208707/251‑1 (MQ=255)                                                                                                                                    gTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGTCAACTGC‑TGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGc                                                                                                               >  2:81145/1‑251 (MQ=255)                                                                                                                                          aGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGC‑TGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTc                                                                                                         <  1:383957/251‑1 (MQ=255)                                                                                                                                           ggTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGC‑TGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCa                                                                                                                       <  2:19280/236‑1 (MQ=255)                                                                                                                                           ggTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGC‑TGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCa                                                                                                                       >  1:19280/1‑236 (MQ=255)                                                                                                                                               gagaAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGC‑TGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAg                                                                                                    <  1:265630/251‑1 (MQ=255)                                                                                                                                                   aaaCCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGC‑TGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCttt                                                                                                <  1:272232/251‑1 (MQ=255)                                                                                                                                                      ccGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCTCCAGATGCCGCCAAACTTCTGGGCAACTGC‑TGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGtc                                                                                             <  1:97707/251‑1 (MQ=255)                                                                                                                                                       cGACAATTGTATACTGATCAGGGCAACCCAATGCTTCACGCACTAGCGGAATTTCAGTTACTGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGC‑TGCCAGTTGCATCGCCATGACCGCGCATCCGCCAGACATATCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTTATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCATTGTCGATAGCTTTGtc                                                                                             >  2:470726/1‑250 (MQ=255)                                                                                                                                                        gaccaTTTTATCCTGCTCACGGCAACCGAATGCTTCCCGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGC‑TGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGtcat                                                                                           <  2:330260/247‑1 (MQ=255)                                                                                                                                                               ttATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGC‑TGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTAcc                                                                                    >  2:285620/1‑251 (MQ=255)                                                                                                                                                                         cACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGC‑TGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATc                                                                          >  1:414447/1‑251 (MQ=255)                                                                                                                                                                         cACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGC‑TGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATc                                                                          >  2:461061/1‑251 (MQ=255)                                                                                                                                                                              cAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGC‑TGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGcccc                                                                     >  2:151045/1‑251 (MQ=255)                                                                                                                                                                                   gAATGCTTCACGCAATACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGC‑TGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTAGATCGCTTTGTCATCACTCCCGGCCGCAATCGCCCCCTGtt                                                                >  2:69937/1‑251 (MQ=255)                                                                                                                                                                                                   aCCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGC‑TGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAGCGGCGCATTACGGGCGTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGCTCCCTGTTCCAGTACGTCGCTGaa                                                <  2:59406/251‑1 (MQ=255)                                                                                                                                                                                                        aCTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGC‑TGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGAAACGCt                                           <  2:42035/251‑1 (MQ=255)                                                                                                                                                                                                           ttCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGC‑TGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGAAACGCTc                                          <  1:251304/249‑1 (MQ=255)                                                                                                                                                                                                            tCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGC‑TGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCAGCACTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGAAACGCTcgc                                        <  2:59704/250‑1 (MQ=255)                                                                                                                                                                                                               gTTAATACGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTCC‑TGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTACGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGTTGAAACGATCGCGccc                                    >  2:308858/1‑251 (MQ=255)                                                                                                                                                                                                                  aaTGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGC‑TGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGAAACGCTCGCGcccttc                                 <  2:66178/251‑1 (MQ=255)                                                                                                                                                                                                                               gcCAGATGCCGCCAAACTTCTGGGCAACTGC‑TGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTTTACCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCATCAATCGCGCCCTGTTCCAGTACGGCGCTGAAACGCTCGCGCCCTTCCCCTTCAATcaca                    <  1:313349/251‑1 (MQ=255)                                                                                                                                                                                                                               gcCAGATGCCGCCAAACTTCTGGGCAACTGC‑TGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGAAACGCTCGCGCCCTTCCCCTTCAATcac                     >  1:274058/1‑250 (MQ=255)                                                                                                                                                                                                                                        cgccAAACTTCTGGGCAACTGC‑TGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGAAACGCTCGCGCCCTTCCCCTTCAATCACAAAAAAATGc           >  2:214213/1‑251 (MQ=255)                                                                                                                                                                                                                                        cgccAAACTTCTGGGCAACTGC‑TGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGAAACGCTCGCGCCCTTCCCCTTCAATCACAAAAAAATGc           >  1:285376/1‑251 (MQ=255)                                                                                                                                                                                                                                        cgccAAACCCCTGGGCAACTGC‑TGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGAAACGCTCGCGCCCTTCCCCTTCAATCACAAAAAAATGc           >  1:378312/1‑251 (MQ=255)                                                                                                                                                                                                                                                 ccTGGGCAACTGC‑TGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGAAACGCTCGCGCCCTTCCCCTTCAATCACAAAAAAATGCCACGGTTGc  <  1:461061/251‑1 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                            GATAAGTGAGAGTGTCTGAATTCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGC‑TGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGAAACGCTCGCGCCCTTCCCCTTCAATCACAAAAAAATGCCACGGTTGC  >  NC_000913/1848061‑1848551

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