Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,348,968	0	T	A	72.7%	36.0 / 10.1	22	intergenic (‑224/+347)	dcuB/dcuR	anaerobic C4‑dicarboxylate transporter DcuB/DNA‑binding transcriptional activator DcuR
Reads supporting (aligned to +/- strand):  ref base T (3/3);  new base A (8/8);  total (11/11)
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.41e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AAATAACATGTGTGAACCCTCGCGATAATCCTATTTAAATTTTTGCTGAATAGATCACAGTCACGTTCTGTTTTGTATGAACTGTTTTCAGAGCCTGACTAAATATCCGTTTGGTCGTTACCGGCTTTAGCAAATACCTCACAGTGAATATTGGCTGATTAAAGCAGCGGCAAATTTATGCACTGTGATAAGCGGCTTTTTCAGAGAGAAAGACTCTCTGGCGATGGGTTTATAGTTTGTACTGCGCAGATGTTAATGTTTTGTTTGCTGTTATCTTATTGATATTTATCGGTCTATTTTTTGTTCTCAGGGCAGATGTCATTAGGTTTATAGATTAATCTGATCTACCCATTTGTGGGTAAAAATACACATAATGCGGGTGACATAATAGTTAATTAACTTTTGTTAGCGTTTTGAAATTAAAAACACCGTTCACCTGAAGAGATATTAATTTTTAGCGATGATGG  >  NC_000913/4348736‑4349202                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                           aaaTAACATGTGTGAACCCTCGCGATAATCCTATTTAAATTTTTGCTGAATAGATCACAGTCACGTTCTGTTTTGTATGAACTGTTTTCAGAGCCTGACTAAATATCCGTTTGGTCGTTACCGGCTTTAGCAAATACCTCACAGTGAATATTGGCTGATTAAAGCAGCGGCAAATTTATGCACTGTGATAAGCGGCTTTTTCAGAGAGAAAGACTCTCTGGCGATGGGTTTAAAGTTTGTACTgcgc                                                                                                                                                                                                                              >  2:352132/1‑247 (MQ=255)     aaCATGTGTGAACCCTCGCGATAATCCTATTTAAATTTTTGCTGAATAGATCACAGTCACGTTCTGTTTTGTATGAACTGTTTTCAGAGCCTGACTAAATATCCGTTTGGTCGTTACCGGCTTTAGCAAATACCTCACAGTGAATATTGGCTGATTAAAGCAGCGGCAAATTTATGCACTGTGATAAGCGGCTTTTTCAGAGAGAAAGACTCTCTGGCGATGGGTTTAAAGTTTGTACTGCGCAGATGtt                                                                                                                                                                                                                       <  1:352132/250‑1 (MQ=255)          gtgtgAACCCTCGCGATAATCCTATTTAAATTTTTGCTGAATAGATCACAGTCACGTTCTGTTTTGTATGAACTGTTTTCAGAGCCTGACTAAATATCCGTTTGGTCGTTACCGGCTTTAGCAAATACCTCACAGTGAATATTGGCTGATTAAAGCAGCGGCAAATTTATGCACTGTGATAAGCGGCTTTTTCAGAGAGAAAGACTCTCTGGCGATGGGTTTAAAGTTTGTACTGCGCAGATGTTAAtgtt                                                                                                                                                                                                                 <  1:342017/251‑1 (MQ=255)                     cgcgATAATCCTATTTAAATTTTTGCTGAATAGATCACAGTCACGTTCTGTTTTGTATGAACTGTTTTCAGAGCCTGACTAAATATCCGTTTGGTCGTTACCGGCTTTAGCAAATACCTCACAGTGAATATTGGCTGATTAAAGCAGCGGCAAATTTATGCACTGTGATAAGCGGCTGTTTCAGAGAGAAAGACTCTCTGGCGATGGGTTTATAGTTTGTACTGCTCAGATGTTAATGTTTTGGTTGCt                                                                                                                                                                                                        <  2:415803/249‑1 (MQ=255)                              ccTATTTAAATTTTTGCTGAATAGATCACAGTCACGTTCTGTTTTGTATGAACTGTTTTCAGAGCCTGACTAAATATCCGTTTGGTCGTTACCGGCTTTAGCAAATACCTCACAGTGAATATTGGCTGATTAAAGCAGCGGCAAATTTATGCACTGTGATAAGCGGCTTTTTCAGAGAGAAAGACTCTCTGGCGATGGGTTTATAGTTTGTACTGCGCAGATGTTAATGTTTTGTTTGCTGTTATCTTAtt                                                                                                                                                                                             >  1:158264/1‑251 (MQ=255)                                 aTTTAAATTTTTGCTGAATAGATCACAGTCACGTTCTGTTTTGTATGAACTGTTTTCAGAGCCTGACTAAATATCCGTTTGGTCGTTACCGGCTTTAGCAAATTCCTCACAGTGAATATTGGCTGATTAAAGCAGCGGCAAATTTATGCACTGTGATAAGCGGCTTTTTCAGAGAGAAAGACTCTCTGGCGATGGTTTTATAGTTTGTACTGCGCAGATGTTAATGTTTTGTTTGCTGTTATCTTATTGat                                                                                                                                                                                          <  1:178386/251‑1 (MQ=255)                                     aaaTTTTTGCTGAATAGATCACAGTCACTTTCTGTTTTGTATGAACTGTTTTCAGAGCCTGACTAAATATCCGTTTGGTCGTTACCGGCTTTAGCAAATACCTCACAGTGAATATTGGCTGATTAAAGCAGCGGCAAATTTATGCACTGTGATAAGCGGCTTTTTCAGAGAGAAAGACTCTCTGGCGATGGGTTTAAAGTTTGTACTGCGCAGATGTTAATGTTTTGTTTGCTGTTATCTTATTGATAttt                                                                                                                                                                                      <  2:131440/251‑1 (MQ=255)                                                    aGATCACAGTCACGTTCTGTTTTGTATGAACTGTTTTCAGAGCCTGACTAAATATCCGTTTGGTCGTTACCGGCTTTAGCAAATACCTCACAGTGAATATTGGCTGATTAAAGCAGCGGCAAATTTATGCACTGTGATAAGCGGCTTTTTCAGAGAGAAAGACTCTCTGGCGATGGGTTTAAAGTTTGTACTGCGCAGATGTTAATGTTTTGTTTGCTGTTATCTTATTGATATTTATCGGTCTAtttttt                                                                                                                                                                       >  1:129163/1‑251 (MQ=255)                                                                cGTTCTTTATTGTATGAGCTGTTTGCAGAGCCTGACTAAATATCCGTGTGGTCGTTTCCGGCTTTAGCAAATACATCACATTGAATATTGGCTGATTAAAGCAGCGGCAAATTTATGCACTGTGAGAAGCGGCTTTTTCAGAGAGAAAGTCTTTCTGGCGATGTGTTTATAGTTTGTACTGCGTCGATGTTAATGTTTTGTTTGCTGTTATCTTATTGATATTTATCTTTCTATTTTTTGTTCTCAGGGCa                                                                                                                                                           <  2:120161/251‑1 (MQ=255)                                                                         ttGTATGAACTGTTTTCAGAGCCTTACTAAATATCCGTTTTGTCGTTACCGGCTTTAGCAAATACCTCACAGTGAATATTGGCTGATTAAAGCAGCGGCAAATTTATGCACTGTGATAAGCGGCTTTTTCAGAGAGAAAGACTCTCTGGCGATGGGTTTAAAGTTTGTACTGCGCAGATGTTAATGTTTTGTTTGCTGTTATCTTAGTGATATTTATCGGTCTATTTTTTGTTCTCAGGGCAGATGTCAt                                                                                                                                                   <  2:298884/250‑1 (MQ=255)                                                                              tGAACTGTTTTCAGCCCCAGACTCAATATCCGTTTGGTCGTTACCGGCTTTAGCACATACCTCACAGTGAATCTTGGCTGATTAAAGCAGCGGCAAATTTATGCACTGTGATAAGCGGCTTTTTCACACAGAAAGCCTCTCTGGCGATGGGTTTATAGTTTGTACTGCGCAGATGATACTGTTTTGTTTGCTGTTATCTTATTGATATTTATCGTTCTATTTTTTGTTCTCAGGGCAGATGTCATTAGGtt                                                                                                                                             >  2:160618/1‑251 (MQ=255)                                                                                           gagCCTAACTAAATCTCGGTTTGGTCGTTACCCGATTTAGGCCATACCTCACACTGAATATTGGCTGATTCAAGCAGCGGCCAATTTATGCACTGTGATAAGCGGCTTTTTCAGAGAGAAAGACTCTCTGGCGATGGGTTTCAAGTTTGTACTGCGCAGATGTTAATGTTTTGTGTGCTGTTATCTTATTGATATTTATCGGTCTATTTTTTGTTCTCAGGGCAGATGTCATTAGGTTTATAGATTAATCt                                                                                                                                >  2:186068/1‑251 (MQ=255)                                                                                                                                 aGCAAATACCTCACAGTGAATATTGGCTGATTAAAGCAGCGGCAAATTTATGCACTGTGATAAGCGGCTTTTTCAGAGAGAAAGACTCTCTGGCGATGGGTTTAAAGTTTGTACTGCGCAGATGTTAATGTTTTGTTTGCTGTTATCTTATTGATATTTATCGGTCTATTTTTTGTTCTCAGGGCAGATGTCATTAGGTTTATAGATTAATCTGATCTACCCATTTGTGGGTAAAAATACACATAATGCgg                                                                                          >  1:102060/1‑251 (MQ=255)                                                                                                                                   cAAATACCTCACAGTGAATATTGGCTGATTAAAGCAGCGGCAAATTTATGCACTGTGATAAGCGGCTTTTTCAGAGAGAAAGACTCTCTGGCGATGGGTTTAACGTTTGTACTGCGCAGATGTTAATGTTTTGTTTGCTGTTATCTTATTGATATTTATCGGTCTATTTTTTGTTCTCAGGGCAGATGTCATTAGGTTTATAGATTAATCTGATCTACCCATTTGTGGGTAAAAATACACATAATGCGGGt                                                                                        >  1:421385/1‑251 (MQ=255)                                                                                                                                              caGTGAATATTGGCTGATTAAAGCAGCGGCAAATTTATGCACTGTGATAAGCGGCTTTTTCAGAGAGAAAGACTCTCTGGCGATGGGTTTAAAGTTTGTTCTGCGCAGATGTTAATGTTTTGTTTGCTGTTATCTTATTGATATTTATCGGTCTATTTTTTGTTCTCAGGGCAGATGTCATTAGGTTTATAGATTAATCTGATCTACCCATTTGTGGGTAAAAATACACATAATGCGGGTGACATAATAGt                                                                             <  2:102060/251‑1 (MQ=255)                                                                                                                                                                                ttATGCACTGTGATAAGCGGCTTTTTCAGAGAGAAAGACTCTCTGGCGATGGGTTTAAAGTTTGTACTGCGCAGATGTTAATGTTTTGTTTGCTGTTATCTTATTGATATTTATCGGTCTATTTTTTGTTCTCAGGGCAGAGGTCATTAGGTTTATAGATTAATCTGATCTACCCATTTGTGGGTAAAAATACACATAATGCGGGTGACATAATAGTTAATTAACTTTTGTTAGCGTTTTGAAATTaaaa                                            <  1:408375/250‑1 (MQ=255)                                                                                                                                                                                                 cGGCTTTTTCAGAGAGAAAGACTCTCTGGCGATGGGTTTAAAGTTTGTACTGCGCAGTTGTTAATGTTTTGTTTGCTGTTATCTTAGTGATATTTATCGGTCTATTTTTTGTTCTCAGGGCAGATGTCATTAGGTTTATAGATTAATCTGATCTTCCCATTTGTGGGTAAAAATACACATAATGCGGGTGACATAATAGTTAATTAACTTTTGTTAGCGTTTTGAAAATAAAAACTCCGTTCTCCTGAaga                          <  2:419995/251‑1 (MQ=255)                                                                                                                                                                                                        ttCAGAGAGAAAGACTCTCTGACGATGGGTTTATAGTTTGTACTGAGCAGATGTTAATGTTTTGTTTGCTGTTATCTTATTGATATTTATCGGACTATTTTTTGTTCTCAGGGCAGATGTCATTAGGTTTATAGATTAATCTGATCTACCCATTTGTGGGTAAAAATACACATAATGCGGGTGACATAATAGTTAATTAACTTTTGTTAGCGTTTTGAAATTAAAAACACCGTTCACCTGAAGAGATAtt                    >  2:233551/1‑250 (MQ=255)                                                                                                                                                                                                                gaAAGCCTCTCTGGCGATGGGTTTAAAGTTTGTACTGCGCAGATGTTAATGTTTTGTTTGCTGTTATCTTATTGATATTTATCGGTCTATTTTTTGTTCTCAGGGCAGATGTCATTAGGTTTATAGATTAATCTGATCTACCCATTTGTGGGTAAAAATACACATAATGCGGGTGACATAATAGTTAATTAACTTTTGTTAGCGTTTTGAAATTAAAAACACCGTTCACCTGAAGAGATATTAAttttt             >  1:130599/1‑249 (MQ=255)                                                                                                                                                                                                                     aCTCTCTGGCGATGGGTTTAAAGTTTGTACTCCGCAGAAGTTACTGATTTGCTTACTGTTATCTTATTGATATTTATCGGTCTATTTTTTGTTCTCAGGGCAGATGTCATTAGGTTTATAGATTACTCTGATCTACCCCTTTGTGGGTAAAAATACACATAATGCGGGTGACATAATAGTTAATTAACTTTTGTTAGCGTTTTGAAATTAAAAACACCGTTCACCTGAAGAGATATTAATTTTTAGCgatg      >  1:381092/1‑251 (MQ=255)                                                                                                                                                                                                                          ctGGCGATGGGTTTAAAGTTTGTACTGCGCAGATGTTAATGTTTTGTTTGCTGTTATCTTATTGATATTTATCGGTCTATTTTTTGTTCTCAGGGCAGATGTCATTAGGTTTATAGATTAATCTGATCTACCCATTTGTGGGTAAAAATACACATAATGCGGGTGACATAATAGTTAATTAACTTTTGTTAGCGTTTTGAAATTAAAAACACCGTTCACCTGAAGAGATATTAATTTTTAGCGATGATgg  >  1:334256/1‑250 (MQ=255)                                                                                                                                                                                                                            ggCGATGGGTTTAAAGTTTGTACTGCGCAGATGTTAATGTTTTGTTTGCTGTTATCTTATTGATATTTATCGGTCTATTTTTTGTTCTCAGGGCAGATGTCATTAGGTTTATAGATTAATCTGATCTACCCATTTGTGGGTAAAAATACACATAATGCGGGTGACATAATAGTTAATTAACTTTTGTTAGCGTTTTGAAATTAAAAACACCGTTCACCTGAAGAGATATTAATTTTTAGCGATGATgg  <  1:78092/248‑1 (MQ=255)                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                           AAATAACATGTGTGAACCCTCGCGATAATCCTATTTAAATTTTTGCTGAATAGATCACAGTCACGTTCTGTTTTGTATGAACTGTTTTCAGAGCCTGACTAAATATCCGTTTGGTCGTTACCGGCTTTAGCAAATACCTCACAGTGAATATTGGCTGATTAAAGCAGCGGCAAATTTATGCACTGTGATAAGCGGCTTTTTCAGAGAGAAAGACTCTCTGGCGATGGGTTTATAGTTTGTACTGCGCAGATGTTAATGTTTTGTTTGCTGTTATCTTATTGATATTTATCGGTCTATTTTTTGTTCTCAGGGCAGATGTCATTAGGTTTATAGATTAATCTGATCTACCCATTTGTGGGTAAAAATACACATAATGCGGGTGACATAATAGTTAATTAACTTTTGTTAGCGTTTTGAAATTAAAAACACCGTTCACCTGAAGAGATATTAATTTTTAGCGATGATGG  >  NC_000913/4348736‑4349202

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