Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,948,8090GA43.2% 19.0 / 56.2 44N3N (AAC→AATgcvHglycine cleavage complex lipoylprotein
Reads supporting (aligned to +/- strand):  ref base G (4/21);  new base A (19/0);  total (23/21)
Fisher's exact test for biased strand distribution p-value = 6.43e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TGTAAGTGCCGTCGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGCACAAAAACAGGTTT  >  minE/1948742‑1948879
                                                                   |                                                                      
tgtAAGTGCCGTCGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTACGTTg                                                                     <  1:1453246/71‑1 (MQ=255)
tgtAAGTGCCGTCGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTACGTTg                                                                     <  1:55459/71‑1 (MQ=255)
tgtAAGTGCCGTCGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTACGTTg                                                                     <  1:441377/71‑1 (MQ=255)
 gtAAGTGCCGTCGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTACGTTg                                                                     <  1:4313588/70‑1 (MQ=255)
  tAAGTGCCGTCGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTACGTTGCt                                                                   <  1:1114403/71‑1 (MQ=255)
  tAAGTGCCGTCGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTACGTTGCt                                                                   <  1:4304432/71‑1 (MQ=255)
  tAAGTGCCGTCGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTACGTTGCt                                                                   <  1:3763639/71‑1 (MQ=255)
  tAAGTGCCGTCGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTACGTTGCt                                                                   <  1:3595063/71‑1 (MQ=255)
  tAAGTGCCGTCGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTACGTTGCt                                                                   <  1:3297079/71‑1 (MQ=255)
                     tACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTACGTTGCTCATCAATCAATCTCCaaa                                                 >  1:779797/1‑70 (MQ=255)
                     tACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTACGTTGCTCATCAATCAATCTCCaaa                                                 >  1:442052/1‑70 (MQ=255)
                     tACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTACGTTGCTCATCAATCAATCTCCaaa                                                 >  1:3348998/1‑70 (MQ=255)
                     tACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTACGTTGCTCATCAATCAATCTCCaaa                                                 >  1:4197278/1‑70 (MQ=255)
                          aGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTACGTTGCTCATCAATCAATCTCCAAAAAAGTa                                           <  1:3313395/71‑1 (MQ=255)
                          aGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTACGTTGCTCATCAATCAATCTCCAAAAAAGTa                                           <  1:3078255/71‑1 (MQ=255)
                          aGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTACGTTGCTCATCAATCAATCTCCAAAAAAGTa                                           <  1:3060595/71‑1 (MQ=255)
                          aGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTACGTTGCTCATCAATCAATCTCCAAAAAAGTa                                           <  1:4251787/71‑1 (MQ=255)
                          aGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTACGTTGCTCATCAATCAATCTCCAAAAAAGTa                                           <  1:3042648/71‑1 (MQ=255)
                          aGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTACGTTGCTCATCAATCAATCTCCAAAAAAGTa                                           <  1:1289872/71‑1 (MQ=255)
                          aGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTACGTTGCTCATCAATCAATCTCCAAAAAAGTa                                           <  1:1557011/71‑1 (MQ=255)
                                       cTTTGCTGTATTTCAGTTCTGCTGGTACATTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGAc                                >  1:3271006/1‑69 (MQ=255)
                                       cTTTGCTGTATTTCAGTTCTGCTGGTACATTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACgg                              >  1:1826601/1‑71 (MQ=255)
                                       cTTTGCTGTATTTCAGTTCTGCTGGTACATTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACgg                              >  1:4369304/1‑71 (MQ=255)
                                       cTTTGCTGTATTTCAGTTCTGCTGGTACATTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACgg                              >  1:4295614/1‑71 (MQ=255)
                                       cTTTGCTGTATTTCAGTTCTGCTGGTACATTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACgg                              >  1:4289399/1‑71 (MQ=255)
                                       cTTTGCTGTATTTCAGTTCTGCTGGTACATTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACgg                              >  1:4260294/1‑71 (MQ=255)
                                       cTTTGCTGTATTTCAGTTCTGCTGGTACATTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACgg                              >  1:1636946/1‑71 (MQ=255)
                                       cTTTGCTGTATTTCAGTTCTGCTGGTACATTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACgg                              >  1:3930298/1‑71 (MQ=255)
                                       cTTTGCTGTATTTCAGTTCTGCTGGTACATTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACgg                              >  1:875667/1‑71 (MQ=255)
                                       cTTTGCTGTATTTCAGTTCTGCTGGTACATTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACgg                              >  1:2993764/1‑71 (MQ=255)
                                       cTTTGCTGTATTTCAGTTCTGCTGGTACATTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACgg                              >  1:708899/1‑71 (MQ=255)
                                       cTTTGCTGTATTTCAGTTCTGCTGGTACATTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACgg                              >  1:2581953/1‑71 (MQ=255)
                                       cTTTGCTGTATTTCAGTTCTGCTGGTACATTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACgg                              >  1:641518/1‑71 (MQ=255)
                                       cTTTGCTGTATTTCAGTTCTGCTGGTACATTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACgg                              >  1:2721017/1‑71 (MQ=255)
                                       cTTTGCTGTATTTCAGTTCTGCTGGTACATTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACgg                              >  1:3112788/1‑71 (MQ=255)
                                       cTTTGCTGTATTTCAGTTCTGCTGGTACATTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACgg                              >  1:2925516/1‑71 (MQ=255)
                                       cTTTGCTGTATTTCAGTTCTGCTGGTACATTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACg                               >  1:3796552/1‑70 (MQ=255)
                                       cTTTGCTGTATTTCAGTTCTGCTGGTACATTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACg                               >  1:3855369/1‑70 (MQ=255)
                                       cTTTGCTGTATTTCAGTTCTGCTGGTACATTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACg                               >  1:892002/1‑70 (MQ=255)
                                                   tCAGTTCTGCTGGTACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTa                   <  1:3343071/70‑1 (MQ=255)
                                                       ttCTGCTGGTACGTTGCTCATCAATCAATCTCCAAAAAAGTa                                           <  1:275564/42‑1 (MQ=255)
                                                            ctgGTACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGCACaaa           <  1:76109/69‑1 (MQ=255)
                                                            ctgGTACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGCACaaa           <  1:940203/69‑1 (MQ=255)
                                                            ctgGTACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGCACaaa           <  1:155030/69‑1 (MQ=255)
                                                                   gTTGCTCATCTATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGCACAAAAACAGGttt  <  1:1096727/71‑1 (MQ=255)
                                                                   gTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGCACAAAAACAGGttt  <  1:3440718/71‑1 (MQ=255)
                                                                   gTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGCACAAAAACAGGttt  <  1:3205504/71‑1 (MQ=255)
                                                                   gTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGCACAAAAACAGGttt  <  1:444947/71‑1 (MQ=255)
                                                                   gTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGCACAAAAACAGGttt  <  1:2891999/71‑1 (MQ=255)
                                                                   gTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGCACAAAAACAGGttt  <  1:2595646/71‑1 (MQ=255)
                                                                   gTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGCACAAAAACAGGttt  <  1:2295259/71‑1 (MQ=255)
                                                                   gTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGCACAAAAACAGGttt  <  1:228578/71‑1 (MQ=255)
                                                                   gTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGCACAAAAACAGGttt  <  1:2135531/71‑1 (MQ=255)
                                                                   gTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGCACAAAAACAGGttt  <  1:848993/71‑1 (MQ=255)
                                                                   gTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGCACAAAAACAGGttt  <  1:851699/71‑1 (MQ=255)
                                                                   gTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGCACAAAAACAGGttt  <  1:1751513/71‑1 (MQ=255)
                                                                   gTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGCACAAAAACAGGttt  <  1:885567/71‑1 (MQ=255)
                                                                   gTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGCACAAAAACAGGttt  <  1:890639/71‑1 (MQ=255)
                                                                   gTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGCACAAAAACAGGttt  <  1:914950/71‑1 (MQ=255)
                                                                   |                                                                      
TGTAAGTGCCGTCGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGCACAAAAACAGGTTT  >  minE/1948742‑1948879

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