Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,414,9060TG40.9% 13.0 / 23.7 22V211G (GTG→GGG) aidBputative acyl‑CoA dehydrogenase AidB
Reads supporting (aligned to +/- strand):  ref base T (0/13);  new base G (9/0);  total (9/13)
Fisher's exact test for biased strand distribution p-value = 2.01e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGA  >  NC_000913/4414845‑4414965
                                                             |                                                             
tCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTgggggg                                                            >  1:3757982/1‑65 (MQ=255)
      gTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATaa                                                      <  1:1056703/65‑1 (MQ=255)
          tGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATgg                                                 <  1:703232/65‑1 (MQ=255)
           gAGCAGCGCCACCCGTGCAGAGCGTGTGGAAGATGGCTGTTATCGGCTGGTGGGGCATAAATGGt                                                <  1:648349/65‑1 (MQ=255)
           gAGCAACGCCACCCGTGCAGATCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGt                                                <  1:3692204/65‑1 (MQ=255)
             gCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGGTAAAAGGgtt                                              >  1:2333291/1‑65 (MQ=255)
                  accaccCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGGGAAATTGGttttttt                                          >  1:2913801/1‑63 (MQ=255)
                     accCGTNCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAAAGGGTTTTTTTGGt                                     >  1:2211062/1‑65 (MQ=255)
                     accCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCGTAAATGGTTTTTCTCGGt                                     <  1:2796312/65‑1 (MQ=255)
                       ccGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGGGAAATGGTTTTTTTCGGTTc                                   >  1:2663617/1‑65 (MQ=255)
                               agCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGGTAAAAAGGTTTTTTTGGTTTCGCCAAGc                           >  1:2183551/1‑65 (MQ=255)
                               agCGTCTGGAAGATGGCTCTGATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGc                           <  1:1584331/65‑1 (MQ=255)
                                gCGTCTGGAAGACGGCTCTTATCGGCTGGGGGGGGGTAAATGGGTTTTTTCTGGTCCCCCaaacg                           >  1:2292683/1‑63 (MQ=255)
                                 cGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGGTAAAAGGGTTTTTTCGGGTCCGCCAAACGa                         >  1:3667294/1‑65 (MQ=255)
                                  gTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGAt                        <  1:1624090/65‑1 (MQ=255)
                                     tGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATgcg                     <  1:793562/65‑1 (MQ=255)
                                     tGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATgcg                     <  1:1041603/65‑1 (MQ=255)
                                     tGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATgcg                     <  1:1592419/65‑1 (MQ=255)
                                       gAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCa                   <  1:334503/65‑1 (MQ=255)
                                         aGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATc                 <  1:2684352/65‑1 (MQ=255)
                                         aGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATc                 <  1:221479/65‑1 (MQ=255)
                                           aTGGCTCTTATCGGCTGGGGGGGGGTAAATGGGTTTTTTTGGTTCCGCAAAACGATGCGCATCTg               >  1:1090132/1‑65 (MQ=255)
                                           aTGGCTCTTATCGGCTGGGGGGGGATAAATGGGTTTTTTCCGTTCCGCCAAGCGATGCGCCTCTg               >  1:2400168/1‑65 (MQ=255)
                                                tctTATCGGCTGGGGGGGGGTAAAAGGTTTTTCTCCGTTCCCCAAAACGAAGCGCATCTGGTGcc          >  1:3573477/1‑64 (MQ=255)
                                                  ttATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTgg        <  1:3280218/65‑1 (MQ=255)
                                                        gCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGa  <  1:2980497/65‑1 (MQ=255)
                                                        gCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGa  <  1:3556237/65‑1 (MQ=255)
                                                        gCTAGTGGGGCATAAATGGTTTTTCTGGGTTCCGCAAAGCGAGGCGCATCTGGTGCTGGCGCAGa  <  1:2524999/65‑1 (MQ=255)
                                                             |                                                             
TCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGA  >  NC_000913/4414845‑4414965

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