Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,283,233	0	C	T	25.0%	44.3 / 14.2	28	intergenic (‑158/+20)	yjcF/actP	pentapeptide repeats protein/acetate transporter
Reads supporting (aligned to +/- strand):  ref base C (9/12);  new base T (3/4);  total (12/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ATAATTAACAAGAAAGAAAAGGCATATTCTCTGCATTATCATTTTCTACGACTGTCAAAAATCGCTCATTTTTTAATGAGTTTATTTGTTTAATATTATGGGAAAAGGTGATGCATTTGGGAGAGGAAGAGTATTCCCCGGTCAGACGACCGGGGAAGGGGAGAGATTAATGCGCGCGGCCTTGCTCAACGCCAAAGCCGGTCTGGGAGCGGATAAACTGCGCGCGGAACAGTTCACGCTCACGCGCGCCTTCCGCTGAGTTATCGGTTGCCGAGAAGAACC  >  NC_000913/4283087‑4283368                                                                                                                                                   |                                                                                                                                        ataATTAACAAGAAAGAAAAGGCATATTCTCTGCATTATCATTTTCTACGACTGTCAAAAATCGCTCATTTTTTAATGAGTTTATTTGTTTAATATTATGGGAAAAGGTGATGCATTTGGGAGAGGAAGAGTATTCCCCGGTCAgacga                                                                                                                                       >  2:105817/1‑149 (MQ=255)  taATTAACAAGAAAGAAAAGGCATATTCTCTGCATTATCATTTTCTACGACTGTCAAAAATCGCTCATTTTTTAATGAGTTTATTTGTTTAATATTATGGGAAAAGGTGATGCATTTGGGAGAGGAAGAGTATTCCCCGGTCAgacgac                                                                                                                                      <  6:175135/149‑1 (MQ=255)     ttAACAAGAAAGAAAAGGCATATTCTCTGCATTATCATTTTCTACGACTGTCAAAAATCGCTCATTTTTTAATGAGTTTATTTGTTTAATATTATGGGAAAAGGTGATGCATTTGGGAGAGGAAGAGTATTCCCCGGTCAGACGACCgg                                                                                                                                   <  1:105817/149‑1 (MQ=255)          aagaaagaAAAGGCATATTCTCTGCATTATCATTTTCTACGACTGTCAAAAATCGCTCATTTTTTAATGAGTTTATTTGTTTATTATTATGGGAAAAGGTGATGCATTTGGGAGAGGAAGAGTATTCCCCGGTCGTCTGACCGGGGaa                                                                                                                               >  3:4581/1‑148 (MQ=255)          aagaaagaAAAGGCATATTCTCTGCATTATCATTTTCTACGACTGTCAAAAATCGCTCATTTTTTAATGAGTTTATTTGTTTAATATTATGGGAAAAGGTGATGCATTTGGGAGAGGAAGAGTATTCCCCGGTCAGACGACCGGGGAAg                                                                                                                              >  4:143343/1‑149 (MQ=255)           agaaagaaAAGGCATATTCTCTGCATTATCATTTTCTACGACTGTCAAAAATCGCTCATTTTTTAATGAGTTTATTTGTTTATTATTATGGGAAAAGGTGATGCATTTGGGAGAGGAAGAGTATTCCCCGGTCGTCTGACCGGGGaata                                                                                                                             <  4:4581/149‑3 (MQ=255)            gaaagaaaAGGCATATTCTCTGCATTATCATTTTCTACGACTGTCAAAAATCGCTCATTTTTTAATGAGTTTATTTGTTTAATATTATGGGAAAAGGTGATGCATTTGGGAGAGGAAGAGTATTCCCCGGTCGTCTGACCGGGGa                                                                                                                                >  3:10271/1‑145 (MQ=255)            gaaagaaaAGGCATATTCTCTGCATTATCATTTTCTACGACTGTCAAAAATCGCTCATTTTTTAATGAGTTTATTTGTTTAATATTATGGGAAAAGGTGATGCATTTGGGAGAGGAAGAGTATTCCCCGGTCGTCTGACCGGGGa                                                                                                                                <  4:10271/145‑1 (MQ=255)            gaaagaaaAGGCATATTCTCTGCATTATCATTTTCTACGACTGTCAAAAATCGCTCATTTTTTAATGAGTTTATTTGTTTAATATTATGGGAAAAGGTGATGCATTTGGGAGAGGAAGAGTATTCCCCGGTCAGACGACCGGGGAAggg                                                                                                                            <  6:169655/149‑1 (MQ=255)                gaaaAGGCATATTCTCTGCATTATCATTTTCTACGACTGTCAAAAATCGCTCATTTTTTAATGAGTTTATTTGTTTAATATTATGGGAAAAGGTGATGCATTTGGGAGAGGAAGAGTATTCCCCGGTCGTCTGACCGGGGaatactctt                                                                                                                        <  7:17692/149‑8 (MQ=255)                      gCATATTCTCTGCATTATCATTTTCTACGACTGTCAAAAATCGCTCATTTTTTAATGAGTTTATTTGTTTAATATTATGGGAAAAGGTGATGCATTTGGGAGAGGAAGAGTATTCCCCGGTCAGACGACCGGGGAAGGGGAGAGATTaa                                                                                                                  >  4:132572/1‑149 (MQ=255)                       cATATTCTCTGCATTATCATTTTCTACGACTGTCAAAAATCGCTCATTTTTTAATGAGTTTATTTGTTTAATATTATGGGAAAAGGTGATGCATTTGGGAGAGGAAGAGTATTCCCCGGTCAGACGACCGGGGAAGGGGAGAGATTAAt                                                                                                                 <  6:6941/149‑1 (MQ=255)                       cATATTCTCTGCATTATCATTTTCTACGACTGTCAAAAATCGCTCATTTTTTAATGAGTTTATTTGTTTAATATTATGGGAAAAGGTGATGCATTTGGGAGAGGAAGAGTATTCCCCGGTCAGACGACCGGGGAAGGGGAGAGATTAAt                                                                                                                 <  3:143343/149‑1 (MQ=255)                        atatTCTCTGCATTATCATTTTCTACGACTGTCAAAAATCGCTCATTTTTTAATGAGTTTATTTGTTTAATATTATGGGAAAAGGTGATGCATTTGGGAGAGGAAGAGTATTCCCCGGTCAGACGACCGGGGAAggg                                                                                                                            <  8:43853/137‑1 (MQ=255)                        atatTCTCTGCATTATCATTTTCTACGACTGTCAAAAATCGCTCATTTTTTAATGAGTTTATTTGTTTAATATTATGGGAAAAGGTGATGCATTTGGGAGAGGAAGAGTATTCCCCGGTCAGACGACCGGGGAAggg                                                                                                                            >  7:43853/1‑137 (MQ=255)                                                    cTGTCAAAAATCGCTCATTTTTTAATGAGTTTATTTGTTTAATATTATGGGAAAAGGTGATGCATTTGGGAGAGGAAGAGTATTCCCCGGTCAGACGACCGGGGAAGGGGAGAGATTAATGCGCGCGGCCTTGCTCAACGCCAAAGCCg                                                                                    >  3:60866/1‑149 (MQ=255)                                                                                        gTTTAATATTATGGGAAAAGGTGCTGCATTTGGGAGAGGAAGAGTATGCCCCGGTCAGACGACCGGGGAAGGGGAGAGATTAATGCGCGCGGCCTTGCTCAACGCCAAAGCCGGTCTGGGAGCGGATAAACTGCGCGCGGAACAGTtc                                                 <  7:86739/148‑1 (MQ=255)                                                                                            aaTATTATGGGAAAAGGTGATGCATTTGGGAGAGGAAGAGTATTCCCCGGGCAGACGACCGGGGAAGGGGAGAGATTAATGCGCGCGGCCTTGCTCAACGCCAAAGCCGGTCTGGGAGCGGATAAACTGCGCGCGGAACAGTtcacgc                                             <  3:132572/148‑1 (MQ=255)                                                                                              tattatGGGAAAAGGTGATGCATTTGGGAGAGGAAGAGTATTCCCCGGTCAGACGACCGGGGAAGGGGAGAGATTAATGCGCGCGGCCTTGCTCAACGCCAAAGCCGGTCTGGGAGCGGATAAACTGCGCGCGGAACAGTtcacgctca                                          >  6:65804/1‑149 (MQ=255)                                                                                                     ggAAAAGGTGATGCATTTGGGAGAGGAAGAGTATTCCCCGGTCAGACGACCGGGGAAGGGGAGAGATTAATGCGCGCGGCCTTGCTCAACGCCAAAGCCGGTCTGGGAGCGGATAAACTGCGCGCGGAACAGTTCACGCTCACgcgcgc                                   <  5:65804/149‑1 (MQ=255)                                                                                                                            agGAAGAGTATTCCCCGGTCAGACGACCGGGGAAGGGGAGAGATTAATGCGCGCGGCCTTGCTCAACGCCAAAGCCGGTCTGGGAGCGGATAAACTGCGCGCGGAACAGTTCACGCTCACGCGCGCCTTCCGCTGAGTTATCGGTTGcc            <  4:60866/149‑1 (MQ=255)                                                                                                                              gAAGAGTATTCCCCGGTCAGACGACCGGGGAAGGGGAGAGATTAATGCGCGCGGCCTTGCTCAACGCCAAAGCCGGTCTGGGAGc                                                                          >  2:148201/1‑85 (MQ=255)                                                                                                                              gAAGAGTATTCCCCGGTCAGACGACCGGGGAAGGGGAGAGATTAATGCGCGCGGCCTTGCTCAACGCCAAAGCCGGTCTGGGAGc                                                                          <  1:148201/85‑1 (MQ=255)                                                                                                                              gAAGAGTATTCCCCGGTCAGACGACCGGGGAAGGGGAGAGATTAATGCGCGCGGCCTTGCTCAACGCCAAAGCCGGTCTGGGAGCGGATAAACTGCGCGCGGAACAGTTCACGCTCACGCGCGCCTTCCGCTGAGTTATCGGTTGcc            >  4:113224/1‑147 (MQ=255)                                                                                                                               aaGAGTATTCCCCGGTCAGACGACCGGGGAAGGGGAGAGATTAATGCGCGCGGCCTTGCTCAACGCCAAAGCCGGTCTGGGAGCGGATAAACTGCGCGCGGAACAGTTCACGCTCACGCGCGCCTTCCGCTGAGTTATCGGTTGcc            >  6:15582/1‑146 (MQ=255)                                                                                                                                      ttCCCCGGTCGTCTGACCGGGGAAGGGGAGAGATTAATGCGCGCGGCCTTGCTCAACGCCAAAGCCGGTCTGGGAGCGGATAAACTGCGCGCGGAACAGTTCACGCTCACGCGCGcc                                  >  2:26927/1‑117 (MQ=255)                                                                                                                                      ttCCCCGGTCGTCTGACCGGGGAAGGGGAGAGATTAATGCGCGCGGCCTTGCTCAACGCCAAAGCCGGTCTGGGAGCGGATAAACTGCGCGCGGAACAGTTCACGCTCACGCGCGcc                                  <  1:26927/117‑1 (MQ=255)                                                                                                                                      ttCCCCGGTCAGACGACCGGGGAAGGGGAGAGATTAATGCGCGCGGCCTTGCTCAACGCCAAAGCCGGTCTGGGAGCGGATAAACTGCGCGCGGAACAGTTCACGCTCACGCGCGCCTTCCGCTGAGTTATCGGTTGCCGAGAAGAAcc  <  5:15582/149‑1 (MQ=255)                                                                                                                                       tCCCCGGTCAGACGACCGGGGAAGGGGAGAGATTAATGCGCGCGGCCTTGCTCAACGCCAAAGCCGGTCTGGGAGCGGATAAACTGCGCGCGGAACAGTTCACGCTCACGCGCGCCTTCCGCTGAGTTATCGGTTGcc            <  2:150970/138‑1 (MQ=255)                                                                                                                                       tCCCCGGTCAGACGACCGGGGAAGGGGAGAGATTAATGCGCGCGGCCTTGCTCAACGCCAAAGCCGGTCTGGGAGCGGATAAACTGCGCGCGGAACAGTTCACGCTCACGCGCGCCTTCCGCTGAGTTATCGGTTGcc            >  1:150970/1‑138 (MQ=255)                                                                                                                                                   |                                                                                                                                        ATAATTAACAAGAAAGAAAAGGCATATTCTCTGCATTATCATTTTCTACGACTGTCAAAAATCGCTCATTTTTTAATGAGTTTATTTGTTTAATATTATGGGAAAAGGTGATGCATTTGGGAGAGGAAGAGTATTCCCCGGTCAGACGACCGGGGAAGGGGAGAGATTAATGCGCGCGGCCTTGCTCAACGCCAAAGCCGGTCTGGGAGCGGATAAACTGCGCGCGGAACAGTTCACGCTCACGCGCGCCTTCCGCTGAGTTATCGGTTGCCGAGAAGAACC  >  NC_000913/4283087‑4283368

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