Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,374,0340G.30.0% 42.9 / 20.8 20coding (813/882 nt)ycjOputative ABC transporter membrane subunit YcjO
Reads supporting (aligned to +/- strand):  ref base G (1/13);  new base . (0/6);  total (1/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.80e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TATCTCTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCA  >  NC_000913/1373972‑1374095
                                                              |                                                             
tatCTCTACAAACCCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGct                                                             <  1:2973080/65‑1 (MQ=255)
 atCTCTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGAACTCCGGTGTGctc                                                            >  1:375244/1‑65 (MQ=255)
  tctcTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGGTCTCGGTGGGGgtct                                                           >  1:251688/1‑65 (MQ=255)
  tctcTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGGTCTCGGGGGGGctct                                                           >  1:2540643/1‑65 (MQ=255)
  tctcTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGGTCTCGGGGGGGctct                                                           >  1:1509983/1‑65 (MQ=255)
  tctcTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGGTCTCGGGGGGGctcc                                                           >  1:362143/1‑64 (MQ=255)
  tctcTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGCTCTCGGTGGTGctct                                                           >  1:449037/1‑65 (MQ=255)
  tctcTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATTTCGGCGGGGcttt                                                           >  1:1225577/1‑63 (MQ=255)
  tctcTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGctct                                                           >  1:2560289/1‑65 (MQ=255)
  tctcTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGctct                                                           >  1:2547946/1‑65 (MQ=255)
  tctcTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGGGgtct                                                           >  1:2872590/1‑65 (MQ=255)
  tctcTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGGGgtct                                                           >  1:1288710/1‑65 (MQ=255)
  tctcTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGATGttcc                                                           >  1:320693/1‑64 (MQ=255)
  tctcTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGGGGGGgtct                                                           >  1:3468581/1‑65 (MQ=255)
  tctcTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGGGGGGctct                                                           >  1:2172818/1‑65 (MQ=255)
  tctcTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGGGGGGctct                                                           >  1:609204/1‑65 (MQ=255)
    tctACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTtc                                                         >  1:1783070/1‑65 (MQ=255)
      tACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTtcat                                                       >  1:517852/1‑65 (MQ=255)
         aaaaCCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGGGCTCTTcatcat                                                    >  1:257622/1‑65 (MQ=255)
            aCCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGGGCTCTTCATCattat                                                 >  1:3078270/1‑65 (MQ=255)
                cctttgcctttAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTc                                             <  1:2885741/65‑1 (MQ=255)
                    tccctttAATGATGTAGAGGAAGGCGTGTAGATCTCGGTGGTCGCCGTATCATTATTTTCGCTGt                                        <  1:775754/63‑1 (MQ=255)
                      cctttAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTc                                       <  1:3000231/65‑1 (MQ=255)
                      cctttAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTc                                       <  1:913171/65‑1 (MQ=255)
                      cctttAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTc                                       <  1:1647721/65‑1 (MQ=255)
                            aTGGTTTGGGAACAGCGGGGTCTATCTCGGTGGTCGCCGTATCATTATTTTCGCTGTCATTctgc                                <  1:3636034/65‑1 (MQ=255)
                            aTGATTTAGGAAAAGCGGCGGCGGTCTCGGTGGGGCTCTTCATCATTATTTTCGCTGTCATTctg                                 >  1:3390791/1‑65 (MQ=255)
                            aTGATTTAGGAAAAGCGGCGGCGATCTCGGGGGTGCTCTTTATCATTATTTTTTCTGTCCTTctg                                 >  1:1372727/1‑65 (MQ=255)
                                  tAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGAcc                           <  1:2110860/65‑1 (MQ=255)
                                    gggaaaGCCGTGTCGATCTCGGTGGTCGCCGTATCATTATTTTCGCTGTCATTCTGCTGACCAgg                        <  1:3305366/62‑1 (MQ=255)
                                     gAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAgg                        <  1:2546199/65‑1 (MQ=255)
                                     gAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAgg                        <  1:2032715/65‑1 (MQ=255)
                                     gAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAgg                        <  1:2724001/65‑1 (MQ=255)
                                          ggggcggcgATCTCGGTGGTCGCCGTATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACgg                  <  1:1897720/63‑1 (MQ=255)
                                               gtagATCTCGGTGGTCGCCGTATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGaa             <  1:2964375/62‑1 (MQ=255)
                                                 cgATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGaaa            >  1:3567092/1‑64 (MQ=255)
                                                 cgATCTCGGTGGTCGCCGTATCATTATTTTCGCTGTGATTCTGCTGACCAGGAAACGGGTGAAcc           <  1:2695256/65‑1 (MQ=255)
                                                 agATCTCGGTGGTCGCCGTATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAAcc           <  1:1968438/64‑1 (MQ=255)
                                                 agATCTCGGTGGTCGCCGTATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAAcc           <  1:3381143/64‑1 (MQ=255)
                                                  gATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAAcc           <  1:3129776/65‑1 (MQ=255)
                                                  gATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAAcc           <  1:52208/65‑1 (MQ=255)
                                                  gATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAAcc           <  1:534383/65‑1 (MQ=255)
                                                   aTCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCt          <  1:1232952/65‑1 (MQ=255)
                                                         gtggtgCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATgg    >  1:1245228/1‑65 (MQ=255)
                                                         gtggtgCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAACgg    >  1:1256593/1‑65 (MQ=255)
                                                           ggtgCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGca  >  1:11931/1‑65 (MQ=255)
                                                           gctgCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGca  <  1:2352283/63‑1 (MQ=255)
                                                              |                                                             
TATCTCTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCA  >  NC_000913/1373972‑1374095

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