Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,984,1880CG37.5% 25.5 / 15.6 24G294G (GGG→GGCaraGL‑arabinose ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base C (5/10);  new base G (9/0);  total (14/10)
Fisher's exact test for biased strand distribution p-value = 2.02e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGTCGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGG  >  NC_000913/1984121‑1984250
                                                                   |                                                              
ggTCGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTAcgc                                                                <  1:185659/68‑1 (MQ=255)
 gTCGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCGAACATGCCTTTCATTAATTCGCTAcgcc                                                               >  1:261112/1‑68 (MQ=255)
 gTCGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTAcgcc                                                               <  1:414237/68‑1 (MQ=255)
 gTCGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTAcgcc                                                               >  2:816285/1‑68 (MQ=255)
  tCGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTAcggcc                                                              >  1:823130/1‑68 (MQ=255)
    gATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTAcgccccg                                                            <  2:244693/68‑1 (MQ=255)
     atatAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTAcgccccgc                                                           <  1:952358/68‑1 (MQ=255)
     atatAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTAcgccccgc                                                           <  2:187626/68‑1 (MQ=255)
         aaaCCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCt                                                       <  2:232821/68‑1 (MQ=255)
            ccTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGGCCCGCCCCTAcc                                                    >  2:854150/1‑68 (MQ=255)
            ccTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGGCCCGCCCCTAcc                                                    >  2:832709/1‑68 (MQ=255)
            ccTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGGCCCGCCCCTAcc                                                    >  2:421379/1‑68 (MQ=255)
            ccTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGGCCCGCCCCTAcc                                                    >  1:968869/1‑68 (MQ=255)
            ccTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGGCCCGCCCCTAcc                                                    >  2:219293/1‑68 (MQ=255)
            ccTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGGACCGGCCCTAcc                                                    >  2:975177/1‑68 (MQ=255)
            ccTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTAcc                                                    >  1:46719/1‑68 (MQ=255)
            ccTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCACTAcc                                                    >  2:790949/1‑68 (MQ=255)
            ccTGAACGGCGGTGAATTGCGTCCCGCCAAACATGCCTGTCAATAATTCGCTGCGGCCAGCCCCTAcc                                                    >  1:607431/1‑68 (MQ=255)
              tGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAg                                                  <  2:155806/68‑1 (MQ=255)
                   ggcggTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGGCCCGCCCCTACCAGACCaa                                             >  2:597018/1‑68 (MQ=255)
                        tGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGc                                        <  2:46719/68‑1 (MQ=255)
                        tGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGc                                        <  2:496448/68‑1 (MQ=255)
                          aTTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGccc                                      <  1:807313/68‑1 (MQ=255)
                           tttGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCACTACCAGACCAACCAGCCca                                     >  2:176038/1‑68 (MQ=255)
                              gCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCcaaca                                  <  1:601502/68‑1 (MQ=255)
                               cGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCcaacaa                                 <  1:994665/68‑1 (MQ=255)
                                gTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGc                                        >  2:575293/1‑60 (MQ=255)
                                gTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGc                                        <  1:575293/60‑1 (MQ=255)
                                                           tCGCTACGGCCCGCCCCAACCAGACCAAACAGCCCAACAATTTCACCACAGCGAACCGCCAGACatat     >  1:390291/1‑68 (MQ=255)
                                                           tCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTAt     >  2:500538/1‑68 (MQ=255)
                                                              cTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTgg  <  1:11869/68‑1 (MQ=255)
                                                              cTACGCCCCGCCCCAACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTgg  <  2:1014626/68‑1 (MQ=255)
                                                                   |                                                              
GGTCGATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACTGCGAACCGCCAGACTTATTGG  >  NC_000913/1984121‑1984250

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: