Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 566,936 T→C 78.2% L271L (CTA→CTG poxB ← pyruvate dehydrogenase (pyruvate oxidase), thiamin‑dependent, FAD‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE566,9360TC78.2% 53.1 / 9.0 41L271L (CTA→CTGpoxBpyruvate dehydrogenase (pyruvate oxidase), thiamin‑dependent, FAD‑binding
Reads supporting (aligned to +/- strand):  ref base T (7/2);  new base C (8/24);  total (15/26)
Fisher's exact test for biased strand distribution p-value = 6.46e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

AATCTGAATGATTTTGGCATCGGTCGGGTAGAAGGCGCGGTAGGGAAATTGCGTGCCGAGTAGCACTAACGTGTCGGCGTTCATCATGGTATGGAAACCTGACGAGAAGCCGATTAACCCGGTCATTCC  >  minE/566876‑567004
                                                            |                                                                    
aaTCTGAATGATTTTGGCATCGGTCGGGTAGAAGGCGCGGTAGGGAAATTGCGTGCCGAGTAGCACTaa                                                              >  1:1111108/1‑69 (MQ=255)
aaTCTGAATGATTTTGGCATCGGTCGGGTAGAAGGCGCGGTAGGGAAATTGCGTGCCGAGTAGCACTaa                                                              >  1:838329/1‑69 (MQ=255)
aaTCTGAATGATTTTGGCATCGGTCGGGTAGAAGGCGCGGTAGGGAAATTGCGTGCCGAGTAGCACTaa                                                              >  1:1787030/1‑69 (MQ=255)
aaTCTGAATGATTTTGGCATCGGTCGGGTAGAAGGCGCGGTAGGGAAATTGCGTGCCGAGTAGCACTaa                                                              >  1:2297283/1‑69 (MQ=255)
                    cggtcggGTAGAACGCGCGGTAGGGAAATTGCGTGCCGAGCAGCACTAACGTATCGGCATTCATCATgg                                          >  1:986470/1‑69 (MQ=255)
                    cggtcggGTAGAACGCGCGGTAGGGAAATTGCGTGCCGAGCAGCACTAACGTATCGGCATTCATCATgg                                          >  1:720410/1‑69 (MQ=255)
                    cggtcggGTAGAACGCGCGGTAGGGAAATTGCGTGCCGAGCAGCACTAACGTATCGGCATTCATCATgg                                          >  1:446939/1‑69 (MQ=255)
                    cggtcggGTAGAACGCGCGGTAGGGAAATTGCGTGCCGAGCAGCACTAACGTATCGGCATTCATCATgg                                          >  1:3219690/1‑69 (MQ=255)
                    cggtcggGTAGAACGCGCGGTAGGGAAATTGCGTGCCGAGCAGCACTAACGTATCGGCATTCATCATgg                                          >  1:3035766/1‑69 (MQ=255)
                    cggtcggGTAGAACGCGCGGTAGGGAAATTGCGTGCCGAGCAGCACTAACGTATCGGCATTCATCATgg                                          >  1:3021591/1‑69 (MQ=255)
                    cggtcggGTAGAACGCGCGGTAGGGAAATTGCGTGCCGAGCAGCACTAACGTATCGGCATTCATCATgg                                          >  1:2827293/1‑69 (MQ=255)
                    cggtcggGTAGAACGCGCGGTAGGGAAATTGCGTGCCGAGCAGCACTAACGTATCGGCATTCATCATgg                                          >  1:1932861/1‑69 (MQ=255)
                            tAGAAGGCGCGGTAGGGAAATTGCGTGCCGAGTAGCACTAACGTGTCGGCGTTCATCATGGTATGGaaa                                  <  1:1987816/69‑1 (MQ=255)
                                         aGGGAAATTGCGTGCCGAGTAGCACTAACGTGTCGGCGTTCATCATGCTATGGAAACCTGACGAGAAg                      >  1:1823821/1‑68 (MQ=255)
                                          gggAAATTGCGTGCCGAGTAGCACTAACGTGTCGGCGTTCATCATGGTATGGAAACCTGACGAGAAGcc                    >  1:3054387/1‑69 (MQ=255)
                                          gggAAATTGCGTGCCGAGTAGCACTAACGTGTCGGCGTTCATCATGGTATGGAAACCTGACGAGAAGcc                    >  1:470475/1‑69 (MQ=255)
                                            gAAATTGCGTGCCTAGCAGCACTAACGTATCGGCATTCATCATGGTATGGAAACCTGACgag                         <  1:3170452/62‑1 (MQ=255)
                                            gAAATTGCGTGCCGAGCTGCACTAACGTATCGGCATTCATCATGGTATGGAAACCTGACgag                         <  1:811185/62‑1 (MQ=255)
                                            gAAATTGCGTGCCGAGCAGCACTAACGTATCGGCATTCATCATTGTATGGAAACCTGACgag                         <  1:2388992/62‑1 (MQ=255)
                                            gAAATTGCGTGCCGAGCAGCACTAACGTATCGGCATTCATCATGGTATGTAAACCTGACgag                         <  1:777726/62‑1 (MQ=255)
                                            gAAATTGCGTGCCGAGCAGCACTAACGTATCGGCATTCATCATGGTATGGAAACCTGACgag                         <  1:1322649/62‑1 (MQ=255)
                                            gAAATTGCGTGCCGAGCAGCACTAACGTATCGGCATTCATCATGGTATGGAAACCTGACgag                         <  1:158949/62‑1 (MQ=255)
                                            gAAATTGCGTGCCGAGCAGCACTAACGTATCGGCATTCATCATGGTATGGAAACCTGACgag                         <  1:3248004/62‑1 (MQ=255)
                                            gAAATTGCGTGCCGAGCAGCACTAACGTATCGGCATTCATCATGGTATGGAAACCTGACgag                         <  1:3303434/62‑1 (MQ=255)
                                            gAAATTGCGTGCCGAGCAGCACTAACGTATCGGCATTCATCATGGTATGGAAACCTGACgag                         <  1:618905/62‑1 (MQ=255)
                                            gAAATTGCGTGCCGAGCAGCACTAACGTATCGGCATTCATCATGGTATGGAAACCTGACgag                         <  1:3163611/62‑1 (MQ=255)
                                            gAAATTGCGTGCCGAGCAGCACTAACGTATCGGCATTCATCATGGTATGGAAACCTGACgag                         <  1:1123677/62‑1 (MQ=255)
                                            gAAATTGCGTGCCGAGCAGCACTAACGTATCGGCATTCATCATGGTATGGAAACCTGACgag                         <  1:857454/62‑1 (MQ=255)
                                            gAAATTGCGTGCCGAGCAGCACTAACGTATCGGCATTCATCATGGTATGGAAACCTGACgag                         <  1:861180/62‑1 (MQ=255)
                                            gAAATTGCGTGCCGAGCAGCACTAACGTATCGGCATTCATCATGGTATGGAAACCTGACgag                         <  1:1062042/62‑1 (MQ=255)
                                            gAAATTGCGTGCCGAGCAGCACTAACGTATCGGCATTCATCATGGTATGGAAACCTGACgag                         <  1:3163586/62‑1 (MQ=255)
                                            gAAATTGCGTGCCGAGCAGCACTAACGTATCGGCATTCATCATGGTATGGAAACCTGACgag                         <  1:3146839/62‑1 (MQ=255)
                                            gAAATTGCGTGCCGAGCAGCACTAACGTATCGGCATTCATCATGGTATGGAAACCTGACgag                         <  1:1691910/62‑1 (MQ=255)
                                            gAAATTGCGTGCCGAGCAGCACTAACGTATCGGCATTCATCATGGTATGGAAACCTGACgag                         <  1:2855362/62‑1 (MQ=255)
                                            gAAATTGCGTGCCGAGCAGCACTAACGTATCGGCATTCATCATGGTATGGAAACCTGACgag                         <  1:2767341/62‑1 (MQ=255)
                                            gAAATTGCGTGCCGAGCAGCACTAACGTATCGGCATTCATCATGGTATGGAAACCTGACgag                         <  1:2612634/62‑1 (MQ=255)
                                            gAAATTGCGTGCCGAGCAGCACTAACGTATCGGCATTCATCATGGTATGGAAACCTGACgag                         <  1:2556271/62‑1 (MQ=255)
                                            gAAATTGCGTGCCGAGCAGCACTAACGTATCGGCATTCATCATGGTATGGAAACCTGACgag                         <  1:2538240/62‑1 (MQ=255)
                                            gAAATTGCGTGCCGAGCAGCACTAACGTATCGGCATTCATCATGGTATGGAAACCTGACgag                         <  1:2204804/62‑1 (MQ=255)
                                            gAAATTGCGTGCCGAGCAGCACTAACGTATCGGCATTCATCATGGTATGGAAACATGACgag                         <  1:2189757/62‑1 (MQ=255)
                                                ttGCGTGCCGAGTAGCACTAACGTGTCGGCGTTCATCATGGTATGGAAACCTGACGAGAAGCCGATTaa              <  1:2716202/69‑1 (MQ=255)
                                                            tAGCACTAACGTGTCGGCGTTCATCATGGTATGGAAACCTGACGAGAAGCCGATTAACCCGGTCATTcc  <  1:2128776/69‑1 (MQ=255)
                                                            tAGCACTAACGTGTCGGCGTTCATCATGGTATGGAAACCTGACGAGAAGCCGATTAACCCGGTCATTcc  <  1:1410407/69‑1 (MQ=255)
                                                            tAGCACTAACGTGTCGGCGTTCATCATGGTATGGAAACCTGACGAGAAGCCGATTAACCCGGTCATTcc  <  1:654316/69‑1 (MQ=255)
                                                            tAGCACTAACGTGTCGGCGTTCATCATGGTATGGAAACCTGACGAGAAGCCGATTAACCCGGTCATTcc  <  1:688510/69‑1 (MQ=255)
                                                            |                                                                    
AATCTGAATGATTTTGGCATCGGTCGGGTAGAAGGCGCGGTAGGGAAATTGCGTGCCGAGTAGCACTAACGTGTCGGCGTTCATCATGGTATGGAAACCTGACGAGAAGCCGATTAACCCGGTCATTCC  >  minE/566876‑567004

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: