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

GCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGC  >  NC_000913/4414872‑4414940
                                  |                                  
gCAGAGCGTCTGGAAGATGGCTCTTATCGGCtggt                                    <  1:1615581/35‑1 (MQ=255)
  agagCGTCTGGAAGATGGCTCTTATCGGCtggtgg                                  <  1:1637282/35‑1 (MQ=255)
   gagCGTCTGGAAGATGGCTCTTATCGGCTgggggg                                 >  1:3729400/1‑35 (MQ=255)
    agCGTCTGGAAGATGGCTCTTATCGGCTGGTgggg                                <  1:4092850/35‑1 (MQ=255)
      cGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCa                              <  1:1913603/35‑1 (MQ=255)
       gTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCAt                             <  1:4380481/35‑1 (MQ=255)
       gTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCAt                             <  1:3800796/35‑1 (MQ=255)
       gTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCAt                             <  1:2515560/35‑1 (MQ=255)
        tCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATa                            >  1:5773740/1‑35 (MQ=37)
        tCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATa                            >  1:5331306/1‑35 (MQ=37)
         cTGGAAGATGGCTCTTATCGGCTGGTGGGGCATaa                           <  1:7147317/35‑1 (MQ=255)
         cTGGAAGATGGCTCTTATCGGCTGGTGGGGCATaa                           <  1:5480772/35‑1 (MQ=255)
         cTGGAAGATGGCTCTTATCGGCTGGTGGGGCATaa                           <  1:469858/35‑1 (MQ=255)
          tGGAAGATGGCTCTTATCGGCTGGTGGGGCATaaa                          <  1:6988010/35‑1 (MQ=255)
          tGGAAGATGGCTCTTATCGGCTGGTGGGGCATaaa                          <  1:448059/35‑1 (MQ=255)
           ggAAGATGGCTCTTATCGGCTGGGGGGGGATaaaa                         >  1:255957/1‑34 (MQ=25)
              aGATGGCTCTTATCGGCTGGGGGGGGATAAATggg                      >  1:7230377/1‑34 (MQ=25)
              aGATGGCTCTTATCGGCTGGGGGGGCATAAATggg                      >  1:4052176/1‑34 (MQ=38)
                aTGGCTCTTATCGGCTGGTGGGGCATAAATGGttt                    <  1:5499113/35‑1 (MQ=255)
                aTGGCTCTTATCGGCTGGTGGGGCATAAATGGttt                    <  1:6113851/35‑1 (MQ=255)
                aTGGCTCTTATCGGCTGGTGGGGCATAAATGGttt                    <  1:4462553/35‑1 (MQ=255)
                aTGGCTCTTATCGGCTGGGGGGGGGTAAATGGgtt                    >  1:694069/1‑35 (MQ=16)
                 tGGCTCTTATCGGCTGGTGGGGCATAAATGGtttt                   <  1:7176079/35‑1 (MQ=255)
                 tGGCTCTTATCGGCTGGTGGGGCATAAATGGtttt                   <  1:6025452/35‑1 (MQ=255)
                 tGGCTCTTATCGGCTGGGGGGGGATAAATGGgttt                   >  1:4190148/1‑35 (MQ=25)
                  ggCTCTTATCGGCTGGTGGGGCATAAATGGttttt                  <  1:1191327/35‑1 (MQ=255)
                   gCTCTTATCGGCTGGTGGGGCATAAATGGTTTTtc                 <  1:6253102/35‑1 (MQ=255)
                   gCTCTTATCGGCTGGGGGGGCATAAATGGGTTTtc                 >  1:2113918/1‑35 (MQ=37)
                    ctctTATCGGCTGGTGGGGCATAAATGGTTTTtct                <  1:7801235/35‑1 (MQ=255)
                     tctTATCGGCTGGTGGGGCATAAATGGTTTTtctc               <  1:5624899/35‑1 (MQ=255)
                       ttATCGGCTGGTGGGGCATAAATGGTTTTTCTCgg             >  1:935638/1‑35 (MQ=255)
                        tATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGt            <  1:5719931/35‑1 (MQ=255)
                          tCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTc          <  1:1403475/35‑1 (MQ=255)
                          tCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTc          <  1:3239805/35‑1 (MQ=255)
                            ggCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCg        <  1:7479183/35‑1 (MQ=255)
                             gCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGc       <  1:7051373/35‑1 (MQ=255)
                             gCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGc       <  1:3908436/35‑1 (MQ=255)
                             gCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGc       >  1:2000711/1‑35 (MQ=255)
                              cTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCa      >  1:6082152/1‑35 (MQ=255)
                              cTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCa      >  1:1786840/1‑35 (MQ=255)
                                ggtggGGCATAAATGGTTTTTCTCGGTTCCGCaaa    >  1:5130864/1‑35 (MQ=255)
                                ggtggGGCATAAATGGTTTTTCTCGGTTCCGCaaa    >  1:5067/1‑35 (MQ=255)
                                ggtggGGCATAAATGGTTTTTCTCGGTTCCGCaaa    <  1:4970331/35‑1 (MQ=255)
                                ggtggGGCATAAATGGTTTTTCTCGGTTCCGCaaa    >  1:4741020/1‑35 (MQ=255)
                                 gtggGGCATAAATGGTTTTTCTCGGTTCCGCAAAg   >  1:4440602/1‑35 (MQ=255)
                                  tggGGCATAAATGGTTTTTCTCGGTTCCGCAAAGc  <  1:718996/35‑1 (MQ=255)
                                  tggGGCATAAATGGTTTTTCTCGGTTCCGCAAAGc  >  1:3213224/1‑35 (MQ=255)
                                  |                                  
GCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGC  >  NC_000913/4414872‑4414940

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: