Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,794,7870CT27.8% 39.3 / 19.3 36P179S (CCG→TCG) gabP4‑aminobutyrate:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base C (13/13);  new base T (10/0);  total (23/13)
Fisher's exact test for biased strand distribution p-value = 5.85e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GTGCAAAGTCATCGCTATCCTGGCCTTTATTTTCCTTGGTGCAGTCGCAATTAGCGGTTTTTACCCGTATGCCGAAGTGAGCGGGATCTCAAGATTGTGGGATAGCGGCGGCTTTATGCCCAACGGT  >  NC_000913/2794723‑2794849
                                                                |                                                              
gtgCAAAGTCATCGCTATCCTGGCCTTTATTTTCCTTGGTGCAGTCGCAATTAGCGGTTTTTAcc                                                                <  1:3688993/65‑1 (MQ=255)
gtgCAAAGTCATCGCTATCCTGGCCTTTATTTTCCTTGGTGCAGTCGCAATTAGCGGTTTTTAcc                                                                <  1:1827662/65‑1 (MQ=255)
  gCAAAGTCATCGCTATCCTGGCCTTTATTTTCCTTGGTGCAGTCGCAATTAGCGGTTTTTACTCg                                                              >  1:122736/1‑65 (MQ=255)
  gCAAAGTCATCGCTATCCTGGCCTTTATTTTCCTTGGTGCAGTCGCAATTAGCGGTTTTTACTCg                                                              >  1:427678/1‑65 (MQ=255)
   cAAAGTCATCGCTATCCTGGCCTTTATTTTCCTTGGTGCAGTCGCAATTAGCGGTTTTTACCCGt                                                             <  1:2435340/65‑1 (MQ=255)
       gTCATCTCTATCCTGGCCTTTATTTTCCTTGGTGCAGTCGCAATTAGCGGTTTTTACTCGTATga                                                         >  1:2225595/1‑64 (MQ=255)
       gTCATCGCTATCCTGGCCTTTATTTTCCTTGGTGCAGTCGCAATTAGCGGTTTTTACTCGTATGc                                                         >  1:1176554/1‑65 (MQ=255)
       gTCATCGCTATCCTGGCCTTTATTTTCCTTGGTGCAGTCGCAATTAGCGGTTTTTACTCGTATGc                                                         >  1:1457100/1‑65 (MQ=255)
       gTCATCGCTATCCTGGCCTTTATTTTCCTTGGTGCAGTCGCAATTAGCGGTTTTTACTCGTATGc                                                         >  1:918594/1‑65 (MQ=255)
       gTCATCGCTATCCTGGCCTTTATTTTCCTTGGTGCAGTCGCAATTAGCGGTTTTTACTCGTATGc                                                         >  1:127192/1‑65 (MQ=255)
       gTCATCGCTATCCTGGCCTTTATTTTCCTTGGTGCAGTCGCAATTAGCGGTTTTTACTCGTATGc                                                         >  1:11841/1‑65 (MQ=255)
       gTCATCGCTATCCTGGCCTTTATTTTCCTTGGTGCAGTCGCAATTAGCGGTTTTTACCCGTATGc                                                         >  1:822976/1‑65 (MQ=255)
       gTCATCGCTATCCTGGCCTTTATTTTCCTTGGTGCAGTCGCAATTAGCGGTTTTTACCCGTATGc                                                         >  1:749831/1‑65 (MQ=255)
       gTCATCGCTATCCTGGCCTTTATTTTCCTTGGTGCAGTCGCAATTAGCGGTTTTTACCCGTATGc                                                         >  1:1591706/1‑65 (MQ=255)
       gTCATCGCTATCCTGGCCTTTATTTTCCTTGGTGCAGTCGCAATTAGCGGTTTTTACCCGTATGc                                                         >  1:3745549/1‑65 (MQ=255)
       gTCATCGCTATCCTGGCCTTTATTTTCCTTGGTGCAGTCGCAATTAGCGGTTTTTACCCGTATGc                                                         >  1:1726426/1‑65 (MQ=255)
       gTCATCGCTATCCTGGCCTTTATTTTCCTTGGTGCAGTCGCAATTAGCGGTTTTTACCCGTATGc                                                         >  1:2290700/1‑65 (MQ=255)
           tCGCTATCCTGACCTTTATTTTCCTTGGTGCAGTCGCAATTAGCGCTTTTTACCCGTATGCCGaa                                                     <  1:3741150/65‑1 (MQ=255)
                     ggCCTTTATTTTCCTTGGTGCAGTCGCAATTAGCGGTTTTTACCCGTATGCCGAAGTGAGCGGGa                                           <  1:963230/65‑1 (MQ=255)
                       ccTTTATTTTCCTTGGTGCAGTCGCAATTAGCGGTTTTTACTCGTATGCCGTATTCAGCTGGAtc                                         >  1:3045107/1‑65 (MQ=255)
                       ccTTTATTTTCCTTGGTGCAGTCGCAATTAGCGGTTTTTACTCGTATGCCGACGTGAGCGGGAtc                                         >  1:168148/1‑65 (MQ=255)
                            aTTTTCCTTGGTGCAGTCGCAATTAGCGGTTTTTACCCGTATGCCGAAGTGAGCGGGATCTCAAg                                    <  1:2879883/65‑1 (MQ=255)
                             ttttCCTTGGTGCAGTCGCAATTAGCGGTTTTTACCCGTATGCCGAAGTGAGCGGGATCTCAAGa                                   <  1:1586656/65‑1 (MQ=255)
                                    tGGTGCAGTCGCAATTAGCGGTTTTTACCCGTATGCCGAAGTGAGCGGGATCTCAAGATTGTggg                            <  1:3443186/65‑1 (MQ=255)
                                        gaaGTCGCAATTAGCGGTTTTTACCCGTATGCCGAAGTGAGCGGGATCTCAAGATTGTGGGATAg                        >  1:1006110/3‑65 (MQ=255)
                                        gCAGTCGCAATTAGCGGTTTTTACCCGTATGCCGAATTGTGCTTGAACTCAAGATTGTGGGATAg                        >  1:1715417/1‑65 (MQ=255)
                                         cAGTCGCAATTAGCGGTTTTTACCCGTATGCCGAAGTGAGCGGGATCTCAAGATTGTGGGATAgc                       <  1:935820/65‑1 (MQ=255)
                                             cGCAATTAGCGGTTTTTACCCGTATGCCGAAGTGAGCGGGATCTCAAGATTGTGGGATAgcggcg                   >  1:448294/1‑65 (MQ=255)
                                               cAATTAGCGGGTTTTACCCGTATGCCGAAGTGAGCGGGATCTCAAGATTGTGGGATAGCggcggc                 <  1:1324732/65‑1 (MQ=255)
                                                aaTTAGCGGTTTTTACCCGTATGCCGAAGTGAGCGGGATCTCAAGATTGTGGGATAGCGGCGGCt                <  1:2566850/65‑1 (MQ=255)
                                                 aTTAGCGGTTTTTACCCGTATGCCGAAGTGAGCGGGATCTCAAGATTGTGGGATAGCGGCGGCtt               <  1:1548970/65‑1 (MQ=255)
                                                  ttAGCGGTTTTTACTCGTATGCCGAAGTGTGCTTGATATCAAGATTGTGGGATAGCggcggcgct              >  1:1808067/1‑62 (MQ=255)
                                                  ttAGCGGTTTTTACCCGTATGCCGAAGTGTGCGTGATCTCAAGATTGTGGGATAGCGGCGGCttt              >  1:2016088/1‑65 (MQ=255)
                                                  ttAGCGGTTTTTACCCGTATGCCGAAGTGAGCGGGATCTCAAGATTGTGGGATAGCGGCGGGttt              >  1:2485271/1‑65 (MQ=255)
                                                  ttAGCGGTTTTTACCCGTATGCCGAAGTGAGCGGGATCTCAAGATTGTGGGATAGCGGCGGCttt              >  1:645856/1‑65 (MQ=255)
                                                      cGGTTTTTACCCGTATGCCGAAGTGAGCGGGATCTCAAGATTGTGGGATAGCGGCGGCTTTATGc          <  1:3370644/65‑1 (MQ=255)
                                                      cGGTTTTTACCCGTATGCCGAAGTGAGCGGGATCTCAAGATTGTGGGATAGCGGCGGCTCTATGc          >  1:3555001/1‑65 (MQ=255)
                                                       ggTTTTTGCCCGTATGCCGAAGTGAGCGGGATCTCAAGATTGTGGGATAGCGGCGGCTTTATGcc         <  1:946150/65‑1 (MQ=255)
                                                              aCCCGTATGTCGAAGTGAGCGGGATCTCAAGATTGTGGGATAGCGGCGGCTTTATGCCCAACGGt  <  1:3210782/65‑1 (MQ=255)
                                                              aCCCGTATGCCGAAGTGAGCGGGATCTCAAGATTGTGGGATAGCGGCGGCTTTATGCCCAGCGGt  <  1:2421050/65‑1 (MQ=255)
                                                              aCCCGTATGCCGAAGTGAGCGGGATCTCAAGATTGTGGGATAGCGGCGGCTTTATGCCCAACGGt  <  1:3540463/65‑1 (MQ=255)
                                                              aCCCGTATGCCGAAGTGAGCGGGATCTCAAGATTGTGGGATAGCGGCGGCTTTATGCCCAACGGt  <  1:3396836/65‑1 (MQ=255)
                                                                |                                                              
GTGCAAAGTCATCGCTATCCTGGCCTTTATTTTCCTTGGTGCAGTCGCAATTAGCGGTTTTTACCCGTATGCCGAAGTGAGCGGGATCTCAAGATTGTGGGATAGCGGCGGCTTTATGCCCAACGGT  >  NC_000913/2794723‑2794849

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