Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,984,1880CG70.6% 11.5 / 13.7 17G294G (GGG→GGCaraGL‑arabinose ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base C (0/5);  new base G (12/0);  total (12/5)
Fisher's exact test for biased strand distribution p-value = 1.62e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.34e-02
Rejected as consensus: Frequency below/above cutoff threshold.

CGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTAC  >  NC_000913/1984124‑1984255
                                                                |                                                                   
cGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTAcggccc                                                                  >  2:715446/1‑68 (MQ=255)
  atatAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTAcgccccgc                                                                <  2:761451/68‑1 (MQ=255)
  atatAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTAcgccccgc                                                                <  1:656192/68‑1 (MQ=255)
  atatAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTAcgccccgc                                                                <  1:715445/68‑1 (MQ=255)
  atatAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTAcgccccgc                                                                <  2:330949/68‑1 (MQ=255)
  atatAAACCTGACCGGCGGTGAGTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTAcgccccgc                                                                <  1:709602/68‑1 (MQ=255)
         ccTGGCCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAAATCGCTACGGCCCGCCCCTagc                                                         >  1:571161/1‑66 (MQ=255)
         ccTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGGCCCGCCCCTAcc                                                         >  2:695808/1‑68 (MQ=255)
         ccTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGGCCCGCCCCTAcc                                                         >  1:651346/1‑68 (MQ=255)
         ccTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGGCCCGCCCCTAcc                                                         >  1:708203/1‑68 (MQ=255)
         ccTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGGCCCGCCCCTAcc                                                         >  2:59949/1‑68 (MQ=255)
         ccTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGGCCCGCCCCTAcc                                                         >  2:360009/1‑68 (MQ=255)
         ccTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGGACCGCCACTAcc                                                         >  2:525433/1‑68 (MQ=255)
         ccTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTAAGGCCCGCCCCTAcc                                                         >  2:690389/1‑68 (MQ=255)
                     tgtttttCGTCCCGACAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGc                                             <  1:230004/65‑1 (MQ=255)
                     tGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGc                                             <  2:655188/68‑1 (MQ=255)
                                            ccTTTCATTAATTCGCTACGGCCCGCCCCTACCAGACCAAACAGCCcaacaa                                      >  2:582852/1‑52 (MQ=255)
                                            ccTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCcaacaa                                      <  1:582852/52‑1 (MQ=255)
                                                   ttAATTCGCTAGGGCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGa               >  2:187021/1‑68 (MQ=255)
                                                     aaTTCGCTACGGCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGGCt             >  1:366020/1‑68 (MQ=255)
                                                       ttCGCTACGGCCCGCCCCTACCAGACCAAACAGCCCAACAATTTcaccac                             >  2:460776/1‑50 (MQ=255)
                                                       ttCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTcaccac                             <  1:460776/50‑1 (MQ=255)
                                                           cTACTCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCTAACCGCCAGACTTATTgg       <  2:228756/68‑1 (MQ=255)
                                                                ccccgccccTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTac  <  1:187021/68‑1 (MQ=255)
                                                                |                                                                   
CGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGGCGTAC  >  NC_000913/1984124‑1984255

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: