Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,092,627	0	T	C	24.3%	58.1 / 20.1	37	intergenic (+448/+470)	PP_1865/aroH	transposase/3‑deoxy‑7‑phosphoheptulonate synthase
Reads supporting (aligned to +/- strand):  ref base T (14/14);  new base C (7/2);  total (21/16)
Fisher's exact test for biased strand distribution p-value = 2.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGG  >  NC_002947/2092498‑2092747                                                                                                                                  |                                                                                                                         gAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCg                                                                                                                    <  1:234666/136‑1 (MQ=255) gAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCg                                                                                                                    <  1:69879/136‑1 (MQ=255)   ggTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTCCAGGCGAc                                                                                                                  >  1:46171/1‑136 (MQ=255)   ggTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTCCAGGCGAc                                                                                                                  >  1:43408/1‑136 (MQ=255)                             atatCGCCGGCATGGCTGCCCTCTGCCTGGTGTACCCGGAAGAACTCCGTTGCTCTTGCTCTTGCTATTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTCCAGGCGACGCAGATTGCGACTTCAGGACGCCTAg                                                                                        >  1:329555/1‑136 (MQ=25)                                            cTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTCCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGGAGGCATTGCgga                                                                         <  2:351204/136‑1 (MQ=21)                                                   cGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGt                                                                  >  1:232975/1‑136 (MQ=255)                                                   cGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTCCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGGAGGCATTGCggagggagat                                                                  >  1:104400/1‑134 (MQ=18)                                                    ggCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTCCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGGAGGCATTGCGGAGGGAGATg                                                                 >  2:41256/1‑136 (MQ=18)                                                    ggCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTCCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGGAGGCATTGCGGAGGGAGATg                                                                 >  2:364051/1‑136 (MQ=18)                                                            gtACCTGGAAGAACGCCGTTGCTGTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCAt                                                         <  2:162162/136‑1 (MQ=255)                                                              aCCTGGAAGAACGCCGTTGCTCATGCTCTTGCTCTTGATTTGGCTTTTGCTTCTAAGCGCGCGGTAGATCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTAGCGGAGGGAGGTGACCGGCATgg                                                       >  1:339118/1‑136 (MQ=17)                                                                           ccGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCTGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCATgc                                          <  2:26600/136‑1 (MQ=255)                                                                             gTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAgcgc                                        <  2:160314/136‑1 (MQ=255)                                                                                             ttctcgctTTTGATCTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATggg                        >  1:306773/2‑136 (MQ=255)                                                                                                 tgctTTGTCTTTTGCTTCTCAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGTGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCt                    <  2:86882/136‑1 (MQ=12)                                                                                                 tgctTAGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCt                    <  2:382702/136‑1 (MQ=17)                                                                                                 acagttGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCt                    <  1:180508/132‑1 (MQ=17)                                                                                                     ttGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCt                    <  2:248853/132‑1 (MQ=32)                                                                                                     ttGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCt                    >  1:248853/1‑132 (MQ=32)                                                                                                                 tCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCt                    >  1:273800/1‑120 (MQ=32)                                                                                                                 tCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCt                    <  2:273800/120‑1 (MQ=32)                                                                                                                 tCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCAGCAGCGCGGTCCTccc    >  1:47257/1‑136 (MQ=17)                                                                                                                   tAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGCTGGGGGCTGCCGCGCGGTCCAccccg  >  1:245881/1‑134 (MQ=11)                                                                                                                   tAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCt                    >  1:4539/1‑118 (MQ=32)                                                                                                                   tAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCt                    <  2:68719/118‑1 (MQ=32)                                                                                                                   tAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCt                    <  2:4539/118‑1 (MQ=32)                                                                                                                   tAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCt                    <  1:15922/118‑1 (MQ=32)                                                                                                                   tAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCt                    >  2:348866/1‑118 (MQ=32)                                                                                                                   tAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCt                    >  2:15922/1‑118 (MQ=32)                                                                                                                   tAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCt                    >  1:68719/1‑118 (MQ=32)                                                                                                                   tAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCgg  >  1:301998/1‑136 (MQ=32)                                                                                                                   tAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCgg  >  1:343382/1‑136 (MQ=32)                                                                                                                   tAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCCTCAAGCGCTGGGCCCCCGGATGGGGCCTGCCGCCCGGTCCTCCCgg  >  1:154611/1‑136 (MQ=255)                                                                                                                   tAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCAGGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCt                    <  1:348866/118‑1 (MQ=17)                                                                                                                      gcgcgcGGTAGTCCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGc                                                                            <  2:357240/59‑1 (MQ=16)                                                                                                                      gcgcgcGGTAGTCCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGc                                                                            >  1:357240/1‑59 (MQ=16)                                                                                                                                  |                                                                                                                         GAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGG  >  NC_002947/2092498‑2092747

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