Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,374,0390TG25.0% 38.8 / 12.8 20F273C (TTC→TGC) ycjOputative ABC transporter membrane subunit YcjO
Reads supporting (aligned to +/- strand):  ref base T (8/7);  new base G (0/5);  total (8/12)
Fisher's exact test for biased strand distribution p-value = 5.47e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.01e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TCTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAA  >  NC_000913/1373976‑1374100
                                                               |                                                             
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:2107156/1‑65 (MQ=255)
                                                            tctTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACaaa  >  1:2043374/1‑65 (MQ=255)
                                                            tctTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACaaa  >  1:926678/1‑65 (MQ=255)
                                                            tctTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACaaa  >  1:1110655/1‑65 (MQ=255)
                                                               |                                                             
TCTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAA  >  NC_000913/1373976‑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: