Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,374,0400CT32.1% 24.9 / 22.7 28F273F (TTC→TTT) ‡ycjOputative ABC transporter membrane subunit YcjO
Reads supporting (aligned to +/- strand):  ref base C (9/10);  new base T (0/9);  total (9/19)
Fisher's exact test for biased strand distribution p-value = 2.61e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.29e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAA  >  NC_000913/1373984‑1374100
                                                        |                                                            
aCCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCCGGGGTGCTCTTCATCatttt                                                      >  1:1865104/1‑63 (MQ=255)
aCCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCCAGCTTGGTGGGGGTCTTCATCattat                                                      >  1:649825/1‑65 (MQ=255)
   gcctttgcctttAATGATTTAGGAAAAGCGGCGGCGATCTCGGCGGTGCTCTTCATCATTAtttt                                                   >  1:1400406/1‑65 (MQ=255)
        tgcttttAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTCGCCGTATCATTATTTTCGCTGt                                             <  1:1755169/65‑1 (MQ=255)
          cctttAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTc                                            <  1:2955138/65‑1 (MQ=255)
            tttAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCAt                                          <  1:2468439/65‑1 (MQ=255)
                   aTTTAGGAAAAGCGGCGGAGATCTCGGTGGTCGCCGTATCATTATTTTCGCTGTCATTCTGCTGa                                  <  1:496659/65‑1 (MQ=255)
                    tttAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGa                                  <  1:2848552/65‑1 (MQ=255)
                        ggAAAAGGCGTGTCGATCTCGGTGGTCGCCGTATCATTATTTTCGCTGTCATTCTGCTGACCAgg                             <  1:1473200/65‑1 (MQ=255)
                        ggAAAAGCGGCGGCGATCTCGGTGGTCGCCGTATCATTATTTTCGCTGTCATTCTGCTGACCAgg                             <  1:2804024/65‑1 (MQ=255)
                         gAAAAGCGGCGGCGATCTCGGTGGTTCTCCGCATCATTATTTTCGCTGTCATTCTGCTGACCAgg                             <  1:347420/65‑1 (MQ=255)
                         gAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAgg                             <  1:769098/65‑1 (MQ=255)
                         gAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAgg                             <  1:1652054/65‑1 (MQ=255)
                         gAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAgg                             <  1:2955208/65‑1 (MQ=255)
                                ggcggcgATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACggg                      <  1:2028739/65‑1 (MQ=255)
                                    tagATCTCGGTGGTCGCCGTATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAAc                 <  1:380128/63‑1 (MQ=255)
                                    gcgATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGaa                  <  1:546826/65‑1 (MQ=255)
                                     cgATCTCGGTGGTCGCCGTATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAAcc                <  1:7319/65‑1 (MQ=255)
                                     cgATCTCGGTGGTCGCCGTATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAAcc                <  1:2586655/65‑1 (MQ=255)
                                     agATCTCGGTGGTCGCCGTATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAAcc                <  1:2977990/64‑1 (MQ=255)
                                     agATCTCGGTGGTCGCCGTATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAAcc                <  1:2853750/64‑1 (MQ=255)
                                      gATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAAcc                <  1:1201442/65‑1 (MQ=255)
                                       aTCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCt               <  1:3687247/65‑1 (MQ=255)
                                               ggtgCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGca       >  1:3368311/1‑65 (MQ=255)
                                               ggtgCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGca       >  1:200063/1‑65 (MQ=255)
                                               ggtgCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGca       >  1:2230211/1‑65 (MQ=255)
                                                    tctTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACaaa  >  1:1902265/1‑65 (MQ=255)
                                                    tctTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACaaa  >  1:2469201/1‑65 (MQ=255)
                                                    tctTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACaaa  >  1:2021057/1‑65 (MQ=255)
                                                    tctTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACaaa  >  1:618299/1‑65 (MQ=255)
                                                    tctTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACaaa  >  1:1477023/1‑65 (MQ=255)
                                                    tctTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACaaa  >  1:2197242/1‑65 (MQ=255)
                                                        |                                                            
ACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAA  >  NC_000913/1373984‑1374100

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