Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,423,2870GC45.5% 14.9 / 14.4 22intergenic (+93/+136)yhdZ/rrfFputative ABC transporter ATP‑binding subunit YhdZ/5S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (6/6);  new base C (10/0);  total (16/6)
Fisher's exact test for biased strand distribution p-value = 1.52e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTT  >  NC_000913/3423240‑3423304
                                               |                 
cTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTTCCTGcc                 >  1:628207/1‑50 (MQ=18)
cTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCGcccccc                 >  1:1833902/1‑50 (MQ=18)
 tttCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGCCCTTCCCGCCa                >  1:2506808/1‑50 (MQ=17)
 tttCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCa                >  1:3207604/1‑50 (MQ=17)
  ttCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGcccccccc               >  1:1892820/1‑48 (MQ=17)
  ttCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCAc               >  1:1531466/1‑50 (MQ=255)
     gCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCt            <  1:1382818/50‑1 (MQ=255)
      cTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCTc           >  1:297691/1‑50 (MQ=14)
      cTTGATTTGATGTCTGGCAGTTTATGGCGGGCGCCCCGCCCCCCAccccc           >  1:2406789/1‑48 (MQ=14)
       ttgatttgatGTCTGGCAGTTTATGGCGGGCGGCCTGCCCGCcccccccc          >  1:1676639/1‑50 (MQ=14)
       ttgatttgatGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCTcc          >  1:88896/1‑50 (MQ=14)
        tgatttgatGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCAccccccc         >  1:1216517/1‑49 (MQ=14)
         gatttgatGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCcccccccccg        >  1:1575254/1‑48 (MQ=14)
         gatttgatGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCCCCCTcc          <  1:2373061/48‑1 (MQ=255)
         gatttgatGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCg         <  1:2837395/49‑1 (MQ=40)
          atttgatGTCTGGCAGTTTATGGCGGGCGTCCTTCCCGCCCCCCCCCcgg       >  1:2428532/1‑50 (MQ=11)
           tttgatGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGc      <  1:515545/50‑1 (MQ=38)
           tttgatGTCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCACCCCCCCGGc      >  1:2123516/1‑50 (MQ=11)
            ttgatGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGcc     >  1:1132330/1‑50 (MQ=11)
            ttgatGTCTGGCAGTTTATGGCGGGCGGTCTGCCCCCCACCCCCCCGGcc     >  1:3886866/1‑50 (MQ=11)
             tgatGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGcc     <  1:516894/49‑1 (MQ=38)
              gatGTCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCCCCCCCCCGGccct   >  1:3838866/1‑48 (MQ=11)
              gatGTCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCACCCCCCCGGCCgc   >  1:277479/1‑49 (MQ=11)
               atGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGtt  <  1:674001/50‑1 (MQ=38)
                                               |                 
CTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTT  >  NC_000913/3423240‑3423304

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