Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	632,141	0	C	T	21.2%	63.3 / 13.1	33	intergenic (+142/‑41)	cstA/ybdD	carbon starvation protein involved in peptide utilization; APC peptide transporter family protein/DUF466 family protein
Reads supporting (aligned to +/- strand):  ref base C (14/12);  new base T (3/4);  total (17/16)
Fisher's exact test for biased strand distribution p-value = 6.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTAC  >  NC_000913/631999‑632256                                                                                                                                               |                                                                                                                    aaGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCgg                                                                                                               <  1:24104/149‑1 (MQ=39)    tCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAc                                                                                                            >  3:75801/1‑149 (MQ=255)           gAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCa                                                                                                     <  4:75801/149‑1 (MQ=255)            aGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGcctac                                                                                                          <  3:132523/143‑5 (MQ=38)            aGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGcctac                                                                                                          >  4:132523/1‑139 (MQ=38)                          gCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCg                                                                                      >  7:75888/1‑149 (MQ=255)                                     cTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCttt                                                                                                       <  6:57891/121‑1 (MQ=255)                                     cTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCttt                                                                                                       >  5:57891/1‑121 (MQ=255)                                         gTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGttt                                                                       <  5:111238/149‑1 (MQ=255)                                            cccAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGAt                                                                    >  5:113017/1‑149 (MQ=255)                                                gggTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCAc                                                                <  6:113017/149‑1 (MQ=255)                                                               agtccccTCGCCCCACCGGGGAGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGa                                                  >  8:76311/5‑148 (MQ=38)                                                                tAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGaa                                                 >  5:166414/1‑148 (MQ=255)                                                                                        gggACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCt                        <  7:76311/149‑1 (MQ=255)                                                                                        gggACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCt                        >  4:87029/1‑149 (MQ=255)                                                                                        gggACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATATTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCt                        >  2:75347/1‑149 (MQ=255)                                                                                             tGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCtgatga                   >  5:42450/1‑149 (MQ=255)                                                                                               aCCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGAtt                 <  6:42450/149‑1 (MQ=255)                                                                                                  gAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGt              <  8:75888/149‑1 (MQ=255)                                                                                                    gcgcGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGTCAAGATGTTTGATTCACTg                                                              >  4:21954/1‑99 (MQ=255)                                                                                                    gcgcGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGTCAAGATGTTTGATTCACTg                                                              <  3:21954/99‑1 (MQ=255)                                                                                                      gcgTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGt              <  1:20656/145‑1 (MQ=255)                                                                                                      gcgTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGt              >  2:20656/1‑145 (MQ=255)                                                                                                        gTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATAtt                                           >  3:36059/1‑114 (MQ=255)                                                                                                        gTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATAtt                                           <  4:36059/114‑1 (MQ=255)                                                                                                               gCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACt                                                               <  8:75962/87‑1 (MQ=255)                                                                                                               gCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACt                                                               >  7:75962/1‑87 (MQ=255)                                                                                                                        tgCCGGATGCGGCGTGAACGCCTTATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGAtt    <  7:57095/136‑1 (MQ=255)                                                                                                                        tgCCGGATGCGGCGTGAACGCCTTATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGAtt    >  8:57095/2‑137 (MQ=255)                                                                                                                                       gaaCGCCTTATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTAc  >  4:28972/3‑124 (MQ=255)                                                                                                                                       gaaCGCCTTATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTAc  <  3:28972/122‑1 (MQ=255)                                                                                                                                         aCGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGc                         >  1:46339/1‑99 (MQ=255)                                                                                                                                         aCGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGc                         <  2:46339/99‑1 (MQ=255)                                                                                                                                               |                                                                                                                    AAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTAC  >  NC_000913/631999‑632256

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