Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 940,715 G→A 39.2% L20L (CTG→TTG)  sfcA ← malate dehydrogenase, NAD‑requiring

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE940,7150GA39.2% 24.3 / 30.4 52L20L (CTG→TTG) sfcAmalate dehydrogenase, NAD‑requiring
Reads supporting (aligned to +/- strand):  ref base G (14/17);  new base A (6/14);  total (20/32)
Fisher's exact test for biased strand distribution p-value = 3.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAAGTTACGGCGTTCTTCCATGCTGAAGGCACTGCCTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTC  >  minE/940652‑940769
                                                               |                                                      
tGAAGTTACGGCGTTCTTCCATGCTGAAGGCACTGCCTTTATTCAACAACGGAAATTCCAGCAGTACAg                                                   >  1:1250851/1‑69 (MQ=255)
       aCGGCGTTCTTCCATGCTGAAGGCACTGCCTTTATTCAACAACGGAAATTCCAGCAGTAc                                                     >  1:1415447/1‑60 (MQ=255)
       aCGGCGTTCTTCCATGCTGAAGGCACTGCCTTTACTCAACAACGGAAATTCCAGCAGTAc                                                     >  1:1219284/1‑60 (MQ=255)
          gCGTTCTTCCATGCTGAAGGCACTGCCTTTATTCAACAACGGAAATTCCAGTAGTACAGGGCCAGCGTa                                         <  1:2003010/69‑1 (MQ=255)
          gCGTTCTTCCATGCTGAAGGCACTGCCTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTa                                         <  1:1226528/69‑1 (MQ=255)
          gCGTTCTTCCATGCTGAAGGCACTGCCTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTa                                         <  1:1227966/69‑1 (MQ=255)
              tcttcCATGCTGAAGGCACTGCCTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAggg                                     <  1:1851835/69‑1 (MQ=255)
              tcttcCATGCTGAAGGCACTGCCTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAggg                                     <  1:1182075/69‑1 (MQ=255)
              tcttcCATGCTGAAGGCACTGCCTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAggg                                     <  1:3178454/69‑1 (MQ=255)
               cttcCATGCTGAAGGCACTGCCTTTATTCAACAACTGAAATTCCAGCAGTACAGGGCCAGCGTAAggg                                     >  1:470538/1‑68 (MQ=255)
               cttcCATGCTGAAGGCACTGCCTTTATTCAACAACTGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGa                                    >  1:1457852/1‑69 (MQ=255)
               cttcCATGCTGAAGGCACTGCCTTTATTCAACAACGGAAATTCcagcag                                                        >  1:1971650/1‑49 (MQ=255)
               cttcCATGCTGAAGGCACTGCCTTTATTCAACAACGGAAATTCcagcag                                                        >  1:284493/1‑49 (MQ=255)
               cttcCATGCTGAAGGCACTGCCTTTATTCAACAACGGAAATTCcagcag                                                        >  1:2934446/1‑49 (MQ=255)
               cttcCATGCTGAAGGCACTGCCTTTATTCAACAACGGAAATTCcagcag                                                        >  1:21577/1‑49 (MQ=255)
               cttcCATGCTGAAGGCACTGCCTTTATTCAACAACGGAAATTCcagcag                                                        >  1:265285/1‑49 (MQ=255)
               cttcCATGCTGAAGGCACTGCCTTTATTCAACAACGGAAATTCcagcag                                                        >  1:2434362/1‑49 (MQ=255)
               cttcCATGCTGAAGGCACTGCCTTTATTCAACAACGGAAATTCcagcag                                                        >  1:548668/1‑49 (MQ=255)
               cttcCATGCTGAAGGCACTGCCTTTATTCAACAACGGAAATTCcagcag                                                        >  1:691446/1‑49 (MQ=255)
               cttcCATGCTGAAGGCACTGCCTTTATTCAACAACGGAAATTCcagcag                                                        >  1:1173506/1‑49 (MQ=255)
               cttcCATGCTGAAGGCACTGCCTTTATTCAACAACGGAAATTCcagcag                                                        >  1:767847/1‑49 (MQ=255)
               cttcCATGCTGAAGGCACTGCCTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGa                                    <  1:2697870/69‑1 (MQ=255)
               cttcCATGCTGAAGGCACTGCCTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGa                                    >  1:1146007/1‑69 (MQ=255)
               cttcCATGCAGAAGGCACTGCCTTTATTCAACAACGGAAATTCcagcag                                                        >  1:44997/1‑49 (MQ=255)
                         gAAGGCACTGCCTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCg                           <  1:3366235/68‑1 (MQ=255)
                         gAAGGCACTGCCTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCg                           <  1:3226474/68‑1 (MQ=255)
                              cACTGCCTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCt                     >  1:2255138/1‑69 (MQ=255)
                                cTGCCTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTg                    >  1:405617/1‑68 (MQ=255)
                                cTGCCTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGt                   <  1:2994491/69‑1 (MQ=255)
                                cTGCCTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGt                   <  1:3062568/69‑1 (MQ=255)
                                cTGCCTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGt                   >  1:866034/1‑69 (MQ=255)
                                cTGCCTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGt                   >  1:2616406/1‑69 (MQ=255)
                                cTGCCTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGt                   >  1:2597628/1‑69 (MQ=255)
                                cTGCCTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGt                   <  1:2231301/69‑1 (MQ=255)
                                cTGCCTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGt                   >  1:1430252/1‑69 (MQ=255)
                                cTGCCTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGt                   >  1:2070275/1‑69 (MQ=255)
                                cTGCCTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGt                   <  1:1156004/69‑1 (MQ=255)
                                cTGCCTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGt                   <  1:1881948/69‑1 (MQ=255)
                                cTGCCTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGt                   <  1:1224468/69‑1 (MQ=255)
                                cTGCCTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGt                   <  1:163464/69‑1 (MQ=255)
                                cTGCCTTTATTCAACAACGGAAATTCAAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGt                   >  1:1550251/1‑69 (MQ=255)
                                   ccTTTATTCAACAACGGAAATTCCAGCAATACAGGGCCa                                              <  1:1336702/39‑1 (MQ=255)
                                   ccTTTATTCAACAACGGAAATTCCAGCAATACAGGGCCa                                              <  1:1363000/39‑1 (MQ=255)
                                   ccTTTATTCAACAACGGAAATTCCAGCAATACAGGGCCa                                              <  1:678133/39‑1 (MQ=255)
                                   ccTTTATTCAACAACGGAAATTCCAGCAATACAGGGCCa                                              <  1:3301453/39‑1 (MQ=255)
                                   ccTTTATTCAACAACGGAAATTCCAGCAATACAGGGCCa                                              <  1:3297450/39‑1 (MQ=255)
                                   ccTTTATTCAACAACGGAAATTCCAGCAATACAGGGCCa                                              <  1:3249817/39‑1 (MQ=255)
                                   ccTTTATTCAACAACGGAAATTCCAGCAATACAGGGCCa                                              <  1:1800713/39‑1 (MQ=255)
                                   ccTTTATTCAACAACGGAAATTCCAGCAATACAGGGCCa                                              <  1:3000495/39‑1 (MQ=255)
                                   ccTTTATTCAACAACGGAAATTCCAGCAATACAGGGCCa                                              <  1:1813153/39‑1 (MQ=255)
                                   ccTTTATTCAACAACGGAAATTCCAGCAATACAGGGCCa                                              <  1:1921271/39‑1 (MQ=255)
                                   ccTTTATTCAACAACGGAAATTCCAGCAATACAGGGCCa                                              <  1:2807266/39‑1 (MQ=255)
                                   ccTTTATTCAACAACGGAAATTCCAGCAATACAGGGCCa                                              <  1:2698897/39‑1 (MQ=255)
                                   ccTTTATTCAACAACGGAAATTCCAGCAATACAGGGCCa                                              <  1:2678357/39‑1 (MQ=255)
                                   ccTTTATTCAACAACGGAAATTCCAGCAATACAGGGCCa                                              <  1:2619930/39‑1 (MQ=255)
                                        aTTCAACAACGGAAATTCCAGCAATACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTtttttgt           >  1:758693/1‑69 (MQ=255)
                                        aTTCAACAACGGAAATTCCAGCAATACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTtttttgt           >  1:1131435/1‑69 (MQ=255)
                                        aTTCAACAACGGAAATTCCAGCAATACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTtttttgt           >  1:539322/1‑69 (MQ=255)
                                        aTTCAACAACGGAAATTCCAGCAATACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTtttttgt           >  1:1553231/1‑69 (MQ=255)
                                        aTTCAACAACGGAAATTCCAGCAATACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTtttttgt           >  1:2897027/1‑69 (MQ=255)
                                        aTTCAACAACGGAAATTCCAGCAATACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTtttttgt           >  1:2573294/1‑69 (MQ=255)
                                         ttCTACAACGGAAATTCCAGCATTACAGGGCCAGCGTa                                         <  1:3248237/38‑1 (MQ=38)
                                                 cGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTc  <  1:1707213/69‑1 (MQ=255)
                                                               |                                                      
TGAAGTTACGGCGTTCTTCCATGCTGAAGGCACTGCCTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTC  >  minE/940652‑940769

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: