Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,374,0380TC42.9% 7.3 / 21.2 21F273L (TTC→CTC) ‡ycjOputative ABC transporter membrane subunit YcjO
Reads supporting (aligned to +/- strand):  ref base T (3/9);  new base C (0/9);  total (3/18)
Fisher's exact test for biased strand distribution p-value = 2.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TCTCTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAA  >  NC_000913/1373974‑1374100
                                                                |                                                              
tctcTACAAATCCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGctct                                                                <  1:2366826/65‑1 (MQ=255)
tctcTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGGGATCTCGGGGGTGctct                                                                >  1:793655/1‑65 (MQ=255)
tctcTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGCTCTCGGTGGGGctct                                                                >  1:2301846/1‑65 (MQ=255)
tctcTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGctct                                                                >  1:2533041/1‑65 (MQ=255)
tctcTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGGGGTGctct                                                                >  1:1288575/1‑65 (MQ=255)
tctcTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCCCGGGGGGGctct                                                                >  1:3585110/1‑65 (MQ=255)
tctcTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGAGCCCGGGGGTGctct                                                                >  1:3603001/1‑65 (MQ=255)
tctcTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGAGATCTCGGTGGGGatct                                                                >  1:522567/1‑65 (MQ=255)
          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:1652054/65‑1 (MQ=255)
                                   gAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAgg                             <  1:2955208/65‑1 (MQ=255)
                                   gAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAgg                             <  1:769098/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:2586655/65‑1 (MQ=255)
                                               cgATCTCGGTGGTCGCCGTATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAAcc                <  1:7319/65‑1 (MQ=255)
                                               agATCTCGGTGGTCGCCGTATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAAcc                <  1:2853750/64‑1 (MQ=255)
                                               agATCTCGGTGGTCGCCGTATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAAcc                <  1:2977990/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:2230211/1‑65 (MQ=255)
                                                         ggtgCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGca       >  1:200063/1‑65 (MQ=255)
                                                              tctTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACaaa  >  1:2469201/1‑65 (MQ=255)
                                                              tctTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACaaa  >  1:1477023/1‑65 (MQ=255)
                                                              tctTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACaaa  >  1:2197242/1‑65 (MQ=255)
                                                              tctTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACaaa  >  1:618299/1‑65 (MQ=255)
                                                              tctTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACaaa  >  1:2021057/1‑65 (MQ=255)
                                                              tctTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACaaa  >  1:1902265/1‑65 (MQ=255)
                                                                |                                                              
TCTCTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAA  >  NC_000913/1373974‑1374100

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