Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,374,0400CT25.0% 34.4 / 12.3 24F273F (TTC→TTT) ‡ycjOputative ABC transporter membrane subunit YcjO
Reads supporting (aligned to +/- strand):  ref base C (7/11);  new base T (0/6);  total (7/17)
Fisher's exact test for biased strand distribution p-value = 1.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAATAAA  >  NC_000913/1373984‑1374104
                                                        |                                                                
aCCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCattat                                                          >  1:697337/1‑65 (MQ=255)
          cgtttAGTGATCTAGCAAAATCCGGTGCGTCCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTc                                                <  1:3067899/63‑1 (MQ=255)
          cgtttAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTc                                                <  1:2414410/63‑1 (MQ=255)
          cctttAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTc                                                <  1:3139075/65‑1 (MQ=255)
          cctttAATGATGTAGGCAAGTCGGCGGCGATCTTGGTGGTGCTCTTCATCATTATTTTCGCTGTc                                                <  1:3067298/65‑1 (MQ=255)
           ctttAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATTATTTTTTTCGCTGTCa                                               >  1:1857254/1‑65 (MQ=255)
                   gtttAGGGGAAGCCGCGTCGATCTCGGTGGTCGCCGTATCATTATTTTCGCTGTCATTCTGCTGa                                      <  1:1465332/64‑1 (MQ=255)
                   ggttAGGAAGAGCCGCGTAGATCTCGGTGGTCGCCGTATCATTATTTTCGCTGTCATTCTGCTGa                                      <  1:2231062/63‑1 (MQ=255)
                   aTTTAGGAAAAGCGGCGGCGGTCTCCGCGGTGCTCTTTATCATTATTTTCTCTGTCATTcttctg                                       >  1:1683221/1‑65 (MQ=255)
                        ggAAAAGGGGCGTAGATCTCGGTGGTCGCCGTATCATTATTTTCGCTGTCATTCTGCTGACCAgg                                 <  1:2789594/65‑1 (MQ=255)
                         gAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAgg                                 <  1:1359307/65‑1 (MQ=255)
                         gAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAgg                                 <  1:1667022/65‑1 (MQ=255)
                         gAAAAGCGGCGGCGATCTCGGTGGTACTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAgg                                 <  1:1474628/65‑1 (MQ=255)
                            aaGGGGCGGCGATCTCGGTGGTCGCCGTATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAAc                             <  1:1789855/65‑1 (MQ=255)
                                    tcgATCTCGGTGGTCGCCGTATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAAc                     <  1:1838305/64‑1 (MQ=255)
                                    gcgATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGaa                      <  1:2359243/65‑1 (MQ=255)
                                     cgATCTCGGTGGTCGCCGTATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAAcc                    <  1:2132290/65‑1 (MQ=255)
                                     cgATCTCGGTGGTCGCCGTATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAAcc                    <  1:1916273/65‑1 (MQ=255)
                                     agATCTCGGTGGTCGCCGTATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAAcc                    <  1:1887546/64‑1 (MQ=255)
                                      gATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAAcc                    <  1:2302806/65‑1 (MQ=255)
                                         ctcGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCa                 <  1:3137535/65‑1 (MQ=255)
                                          tcGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCaa                >  1:339400/1‑65 (MQ=255)
                                             gtggtgCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATgg             <  1:559870/65‑1 (MQ=255)
                                               ggtgCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGca           >  1:1983477/1‑65 (MQ=255)
                                               ggtgCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGca           <  1:1822405/65‑1 (MQ=255)
                                                 tgCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGcaac         >  1:309583/1‑65 (MQ=255)
                                                  gCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGcaaca        >  1:673180/1‑65 (MQ=255)
                                                    tctTCATCATTATTTTCGCTGTCATTCTGCTGACTAGGAAACGGGTGAACCTCAATGGCAACaaa      >  1:2142399/1‑65 (MQ=255)
                                                    tctTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACaaa      >  1:687469/1‑65 (MQ=255)
                                                    tctTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGAGAACCTCAATGGCAACaaa      >  1:1573045/1‑65 (MQ=255)
                                                        catcatTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAATaaa  <  1:2318679/65‑1 (MQ=255)
                                                        catcatTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAATaaa  <  1:2549287/65‑1 (MQ=255)
                                                        |                                                                
ACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAATAAA  >  NC_000913/1373984‑1374104

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: