Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,374,0360TG27.8% 28.1 / 13.1 18L272R (CTC→CGC) ycjOputative ABC transporter membrane subunit YcjO
Reads supporting (aligned to +/- strand):  ref base T (4/9);  new base G (0/5);  total (4/14)
Fisher's exact test for biased strand distribution p-value = 2.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.15e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TATCTCTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAA  >  NC_000913/1373972‑1374100
                                                                |                                                                
tatCTCTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGct                                                                  <  1:1248927/65‑1 (MQ=255)
  tctcTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGctct                                                                >  1:1601188/1‑65 (MQ=255)
  tctcTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGGGgtct                                                                >  1:2810777/1‑65 (MQ=255)
  tctcTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGGGctct                                                                >  1:1380827/1‑65 (MQ=255)
  tctcTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGAGCTCGGTGGTGctct                                                                >  1:2537332/1‑65 (MQ=255)
    tctACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGGTCTCGGGGGGGCTTTtc                                                              >  1:1703201/1‑65 (MQ=255)
    tctACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTtc                                                              <  1:2800701/65‑1 (MQ=255)
     ctACCAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTtca                                                             <  1:1317849/65‑1 (MQ=255)
            aCCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGGTCTCGGTGGTGCTCTTCATCattat                                                      >  1:1172573/1‑65 (MQ=255)
            aCCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGGGCTCGGTGGGGCTCTTCATCagtat                                                      >  1:3028620/1‑65 (MQ=255)
            aCCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGGAATCGGTGGCTTTTTTCATCatttt                                                      >  1:658114/1‑63 (MQ=255)
            aCCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGGGCTCTTCATCattat                                                      >  1:756964/1‑65 (MQ=255)
             ccGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCattatt                                                     <  1:551574/65‑1 (MQ=255)
                 ctttgcctttAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCg                                                 >  1:1875703/1‑65 (MQ=255)
                     gcctttAATGATTTGGGAGAAGGGGTGTAGATCTCGGTGGTCGCCGTATCATTATTTTCGCTGTc                                            <  1:783282/65‑1 (MQ=255)
                     gcctttAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGt                                             <  1:1572020/65‑1 (MQ=255)
                                  tAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGAcc                                <  1:329671/65‑1 (MQ=255)
                                          ggcgcgtcgATCTCGGTGGTCGCCGTATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACgg                       <  1:2376886/62‑1 (MQ=255)
                                          gcggcggcgATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACg                        <  1:963523/65‑1 (MQ=255)
                                            ggcggcgATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACggg                      <  1:1567794/65‑1 (MQ=255)
                                                tagATCTCGGTGGTCGCCGTATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAAc                 <  1:2428665/63‑1 (MQ=255)
                                                 cgATCTCGGTGGTCGCCGTATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAAcc                <  1:2703502/65‑1 (MQ=255)
                                                 agATCTCGGTGGTCGCCGTATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAAcc                <  1:1765329/64‑1 (MQ=255)
                                                  gATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAAcc                <  1:2741376/65‑1 (MQ=255)
                                                         gtggtgCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATgg         >  1:22714/1‑65 (MQ=255)
                                                         gtggtgCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATgg         >  1:360425/1‑65 (MQ=255)
                                                         gtggtgCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATgg         >  1:1811831/1‑65 (MQ=255)
                                                           ggtgCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGca       <  1:2540856/65‑1 (MQ=255)
                                                              gCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGcaaca    >  1:2049585/1‑65 (MQ=255)
                                                                tctTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACaaa  >  1:1110655/1‑65 (MQ=255)
                                                                tctTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACaaa  >  1:2107156/1‑65 (MQ=255)
                                                                tctTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACaaa  >  1:2043374/1‑65 (MQ=255)
                                                                tctTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACaaa  >  1:926678/1‑65 (MQ=255)
                                                                |                                                                
TATCTCTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAA  >  NC_000913/1373972‑1374100

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