Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,336,397 G→A 60.8% A261A (GCG→GCA fdoH → formate dehydrogenase‑O, Fe‑S subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,336,3970GA60.8% 32.3 / 49.8 74A261A (GCG→GCAfdoHformate dehydrogenase‑O, Fe‑S subunit
Reads supporting (aligned to +/- strand):  ref base G (10/19);  new base A (20/25);  total (30/44)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

GGCTTGCCGGAGAACCCGGAAATCAGCGAAACCGTGAAATTCTGGAAAGGCATCTGGAAACCGCTCGCGGCTGTTGGCTTTGCGGCTACCTTTGCGGCCAGTATCTTCCACTACGTGGGTGTCGGT  >  minE/2336329‑2336454
                                                                    |                                                         
ggCTTGCCGGAGAACCCGGAAATCAGCGAAACCGTGAAATTCTGGAAAGGCATCTGGAAACCGCTcgcg                                                           <  1:417941/69‑1 (MQ=255)
ggCTTGCCGGAGAACCCGGAAATCAGCGAAACCGTGAAATTCTGGAAAGGCATCTGGAAACCGCTcgcg                                                           <  1:1717038/69‑1 (MQ=255)
ggCTTGCCGGAGAACCCGGAAATCAGCGAAACCGTGAAATTCTGGAAAGGCATCTGGAAACCGCTcgcg                                                           <  1:1725139/69‑1 (MQ=255)
ggCTTGCCGGAGAACCCGGAAATCAGCGAAACCGTGAAATTCTGGAAAGGCATCTGGAAACCGCTcgcg                                                           <  1:1529487/69‑1 (MQ=255)
ggCTTGCCGGAGAACCCGGAAATCAGCGAAACCGTGAAATTCTGGAAAGGCATCTGGAAACCGCTcgcg                                                           <  1:663533/69‑1 (MQ=255)
ggCTTGCCGGAGAACCCGGAAATCAGCGAAACCGTGAAATTCTGGAAAGGCATCTGGAAACCGCTcgcg                                                           <  1:845602/69‑1 (MQ=255)
ggCTTGCCGGAGAACCCGGAAATCAGCGAAACCGTGAAATTCTGGAAAGGCATCTGGAAACCGCTcgcg                                                           <  1:488915/69‑1 (MQ=255)
ggCTTGCCGGAGAACCCGGAAATCAGCGAAACCGTGAAATTCTGGAAAGGCATCTGGAAACCGCTcgcg                                                           <  1:3360305/69‑1 (MQ=255)
 gCTTGCCGGAGAACCCGGAAATCAGCGAAACCGTGAAATTCTGGAAAGGCATCTGGAAACCGCTcgcg                                                           <  1:1295050/68‑1 (MQ=255)
 gCTTGCCGGAGAACCCGGAAATCAGCGAAACCGTGAAATTCTGGAAAGGCATCTGGAAACCGCTcgcg                                                           <  1:2842591/68‑1 (MQ=255)
 gCTTGCCGGAGAACCCGGAAATCAGCGAAACCGTGAAATTCTGGAAAGGCATCTGGAAACCGCTCGCgg                                                          >  1:460867/1‑69 (MQ=255)
                  gAAATCAGCGAAACCGTGAAATTCTGGAAAGGCATCTGGAAACCGCTCGCGGCTGTTGGCTTTGCGGCt                                         <  1:1230288/69‑1 (MQ=255)
                  gAAATCAGCGAAACCGTGAAATTCTGGAAAGGCATCTGGAAACCGCTCGCGGCTGTTGGCTTTGCGGCt                                         <  1:837222/69‑1 (MQ=255)
                  gAAATCAGCGAAACCGTGAAATTCTGGAAAGGCATCTGGAAACCGCTCGCGGCTGTTGGCTTTGCGGCt                                         <  1:1114473/69‑1 (MQ=255)
                  gAAATCAGCGAAACCGTGAAATTCTGGAAAGGCATCTGGAAACCGCTCGCGGCTGTTGGCTTTGCGGCt                                         <  1:1881198/69‑1 (MQ=255)
                          cGAAACCGTGAAATTCTGGAAAGGCATCTGGAAACCGCTcgcg                                                           <  1:2795076/43‑1 (MQ=255)
                           gAAACCGTGAAATTCTGGAAAGGCATCTGGAAACCGCTCGCGGCTGTTGGCTTTGCGGCTACCTTTGCg                                <  1:1579527/69‑1 (MQ=255)
                           gAAACCGTGAAATTCTGGAAAGGCATCTGGAAACCGCTCGCGGCTGTTGGCTTTGCGGCTACCTTTGCg                                <  1:1865011/69‑1 (MQ=255)
                               ccGTGAAATTCTGGAAAGGCATCTGGAAACCGCTCGCGGCTGTTGGCTTTGCGGCTACCTTTGCGGCCa                            <  1:275293/69‑1 (MQ=255)
                               ccGTGAAATTCTGGAAAGGCATCTGGAAACCGCTCGCGGCTGTTGGCTTTGCGGCTACCTTTGCGGCCa                            <  1:2720948/69‑1 (MQ=255)
                                 gTGAGATTCTGGAAAGGCATCTGGAAACCGCTCGCGGCTGTTGGCTTTGCGGCTACCTTTGCGGCCAGt                          <  1:2263320/69‑1 (MQ=255)
                                 gTGAAATTCTGGAAAGGCATCTGGAAACCGCTCGCGGCTGTTGGCTTTGCGGCTACCTTTGCGGCCAGt                          <  1:375768/69‑1 (MQ=255)
                                    aaaTTCTGGAAAGGCATCTGGAAACCGCTCGCGGCTGTTGGCTTTGCGGCTACCTTTGCGGCCAGTa                         <  1:3382081/67‑1 (MQ=255)
                                    aaaTTCTGGAAAGGCATCTGGAAACCGCTCGCGGCTGTTGGCTTTGCGGCTACCTTTGCGGCCAGTa                         <  1:582154/67‑1 (MQ=255)
                                    aaaTTCTGGAAAGGCATCTGGAAACCGCTCGCGGCTGTTGGCTTTGCGGCTACCTTTGCGGCCAGTa                         <  1:2691380/67‑1 (MQ=255)
                                    aaaTTCTGGAAAGGCATCTGGAAACCGCTCGCGGCTGTTGGCTTTGCGGCTACCTTTGCGGCCAGTa                         <  1:392941/67‑1 (MQ=255)
                                    aaaTTCTGGAAAGGCATCTGGAAACCGCTCGCGGCTGTTGGCTTTGCGGCTACCTTTGCGGCCAGTa                         <  1:1894145/67‑1 (MQ=255)
                                    aaaTTCTGGAAAGGCATCTGGAAACCGCTCGCGGCTGTTGGCTTTGCGGCTACCTTTGCGGCCAGTa                         <  1:3005629/67‑1 (MQ=255)
                                    aaaTTCTGGAAAGGCATCTGGAAACCGCTCGCGGCTGTTGGCTTTGCGGCTACCTTTGCGGCCAGTa                         <  1:3022816/67‑1 (MQ=255)
                                    aaaTTCTGGAAAGGCATCTGGAAACCGCTCGCGGCTGTTGGCTTTGCGGCTACCTTTGCGGCCAGTa                         <  1:159047/67‑1 (MQ=255)
                                           ggAAAGGCATCTGGAAACCGCTCGCGGCTGTTGGCTTTGCGGCTACCTTTGCGGCCAGTATCTTCCAc                 <  1:1243648/68‑1 (MQ=255)
                                                 gCATCTGGAAACCGCTCGCGGCTGTTGGCTTTGCGGCTACCTTTGCGGCCAGTATCTTCCACTACGTgg          >  1:3017379/1‑69 (MQ=255)
                                                 gCATCTGGAAACCGCTCGCGGCTGTTGGCTTTGCGGCTACCTTTGCGGCCAGTATCTTCCACTACGTgg          >  1:304547/1‑69 (MQ=255)
                                                 gCATCTGGAAACCGCTCGCGGCTGTTGGCTTTGCGGCTACCTTTGCGGCCAGTATCTTCCACTACGTgg          >  1:404976/1‑69 (MQ=255)
                                                 gCATCTGGAAACCGCTCGCGGCTGTTGGCTTTGCGGCTACCTTTGCGGCCAGTATCTTCCACTACGTgg          >  1:308853/1‑69 (MQ=255)
                                                 gCATCTGGAAACCGCTCGCGGCTGTTGGCTTTGCGGCTACCTTTGCGGCCAGTATCTTCCACTACGTgg          >  1:1156832/1‑69 (MQ=255)
                                                   aTCTGGATACCGCTCGCAGCTGTTGGCTTTGCGGCTACCTTCGCAGCCAGTATCTTCCACTACGTGGgt        <  1:1260425/69‑1 (MQ=255)
                                                   aTCTGGAAACCGCTCGCAGCTGTTGGCTTTGCGGCTACCTTCGCAGCCAGTATCTTCCACTACgtgtgt        <  1:2068845/69‑1 (MQ=255)
                                                   aTCTGGAAACCGCTCGCAGCTGTTGGCTTTGCGGCTACCTTCGCAGCCAGTATCTTCCACTACGTGGgt        <  1:47951/69‑1 (MQ=255)
                                                   aTCTGGAAACCGCTCGCAGCTGTTGGCTTTGCGGCTACCTTCGCAGCCAGTATCTTCCACTACGTGGgt        <  1:114890/69‑1 (MQ=255)
                                                   aTCTGGAAACCGCTCGCAGCTGTTGGCTTTGCGGCTACCTTCGCAGCCAGTATCTTCCACTACGTGGgt        <  1:1062612/69‑1 (MQ=255)
                                                   aTCTGGAAACCGCTCGCAGCTGTTGGCTTTGCGGCTACCTTCGCAGCCAGTATCTTCCACTACGTGGgt        <  1:286319/69‑1 (MQ=255)
                                                   aTCTGGAAACCGCTCGCAGCTGTTGGCTTTGCGGCTACCTTCGCAGCCAGTATCTTCCACTACGTGGgt        <  1:2875917/69‑1 (MQ=255)
                                                   aTCTGGAAACCGCTCGCAGCTGTTGGCTTTGCGGCTACCTTCGCAGCCAGTATCTTCCACTACGTGGgt        <  1:163248/69‑1 (MQ=255)
                                                   aTCTGGAAACCGCTCGCAGCTGTTGGCTTTGCGGCTACCTTCGCAGCCAGTATCTTCCACTACGTGGgt        <  1:2980596/69‑1 (MQ=255)
                                                   aTCTGGAAACCGCTCGCAGCTGTTGGCTTTGCGGCTACCTTCGCAGCCAGTATCTTCCACTACGTGGgt        <  1:2483140/69‑1 (MQ=255)
                                                   aTCTGGAAACCGCTCGCAGCTGTTGGCTTTGCGGCTACCTTCGCAGCCAGTATCTTCCACTACGTGGgt        <  1:1698998/69‑1 (MQ=255)
                                                   aTCTGGAAACCGCTCGCAGCTGTTGGCTTTGCGGCTACCTTCGCAGCCAGTATCTTCCACTACGTGGgt        <  1:1801054/69‑1 (MQ=255)
                                                   aTCTGGAAACCGCTCGCAGCTGTTGGCTTTGCGGCTACCTTCGCAGCCAGTATCTTCCACTACGTGGgt        <  1:180159/69‑1 (MQ=255)
                                                   aTCTGGAAACCGCTCGCAGCTGTTGGCTTTGCGGCTACCTTCGCAGCCAGTATCTTCCACTACGTGGgt        <  1:1815689/69‑1 (MQ=255)
                                                   aTCTGGAAACCGCTCGCAGCTGTTGGCTTTGCGGCTACCTTCGCAGCCAGTATCTTCCACTACGTGGgt        <  1:1874963/69‑1 (MQ=255)
                                                   aTCTGGAAACCGCTCGCAGCTGTTGGCTTTGCGGCTACCTTCGCAGCCAGTATCTTCCACTACGTGGgt        <  1:1878838/69‑1 (MQ=255)
                                                   aTCTGGAAACCGCTCGCAGCTGTTGGCTTTGCGGCTACCTTCGCAGCCAGTATCTTCCACTACGTGGgt        <  1:1504288/69‑1 (MQ=255)
                                                   aTCTGGAAACCGCTCGCAGCTGTTGGCTTTGCGGCTACCTTCGCAGCCAGTATCTTCCACTACGTGGgt        <  1:1881921/69‑1 (MQ=255)
                                                   aTCTGGAAACCGCTCGCAGCTGTTGGCTTTGCGGCTACCTTCGCAGCCAGTATCTTCCACTACGTGGgt        <  1:1909653/69‑1 (MQ=255)
                                                   aTCTGGAAACCGCTCGCAGCTGTTGGCTTTGCGGCTACCTTCGCAGCCAGTATCTTCCACTACGTGGgt        <  1:767389/69‑1 (MQ=255)
                                                   aTCTGGAAACCGCTCGCAGCTGTTGGCTTTGCGGCTACCTTCGCAGCCAGTATCTTCCACTACGTGGgt        <  1:613357/69‑1 (MQ=255)
                                                   aTCTGGAAACCGCTCGCAGCTGTTGGCTTTGCGGCTACCTTCGCAGCCAGTATCTTCCACTACGTGGgt        <  1:603679/69‑1 (MQ=255)
                                                   aTCTGGAAACCGCTCGCAGCTGTTGGCTTTGCGGCTACCTTCGCAGCCAGTATCTTCCACTACGTGGgt        <  1:216826/69‑1 (MQ=255)
                                                   aTCTGGAAACCGCTCGCAGCTGTTGGCTTTGCGGCTACCTTCGCAGCCAGTATCTTCCACTACGTGGgt        <  1:567379/69‑1 (MQ=255)
                                                      tGGAAACCGCTCGCGGCTGTTGGCTTTGCGGCTACCTTTGCGGCCAGTTTCTTCCACTACGTGGGTGTc     >  1:2210874/1‑69 (MQ=255)
                                                      tGGAAACCGCTCGCGGCTGTTGGCTTTGCGGCTACCTTTGCGGCCAGTATCTTc                    >  1:94952/1‑54 (MQ=255)
                                                      tGGAAACCGCTCGCGGCTGTTGGCTTTGCGGCTACCTTTGCGGCCAGTATCTTc                    >  1:3382803/1‑54 (MQ=255)
                                                      tGGAAACCGCTCGCGGCTGTTGGCTTTGCGGCTACCTTTGCGGCCAGTATCTTCCACTACGTGGGTGTc     >  1:1183819/1‑69 (MQ=255)
                                                      tGGAAACCGCTCGCGGCTGTTGGCTTTGCGGCTACCTTTGCGGCCAGTATCTTCCACTACGTGGGTGTc     >  1:363255/1‑69 (MQ=255)
                                                      tGGAAACCGCTCGCAGCTGTTGGCTTTGCGGCTACCTTCGCAGCCAGTATCTTCCACTACGTGGgt        <  1:2607448/66‑1 (MQ=255)
                                                          aaCCGCTCGCAGCTGTTGGCTTTGCGGCTACCTTCGCAGCCAGTATCTTCCACTACGTGGGTGTCgg   >  1:1582813/1‑67 (MQ=255)
                                                          aaCCGCTCGCAGCTGTTGGCTTTGCGGCTACCTTCGCAGCCAGTATCTTCCACTACGTGGGTGTCGGt  >  1:70449/1‑68 (MQ=255)
                                                          aaCCGCTCGCAGCTGTTGGCTTTGCGGCTACCTTCGCAGCCAGTATCTTCCACTACGTGGGTGTCGGt  >  1:2469354/1‑68 (MQ=255)
                                                          aaCCGCTCGCAGCTGTTGGCTTTGCGGCTACCTTCGCAGCCAGTATCTTCCACTACGTGGGTGTCGGt  >  1:882988/1‑68 (MQ=255)
                                                          aaCCGCTCGCAGCTGTTGGCTTTGCGGCTACCTTCGCAGCCAGTATCTTCCACTACGTGGGTGTCGGt  >  1:98409/1‑68 (MQ=255)
                                                          aaCCGCTCGCAGCTGTTGGCTTTGCGGCTACCTTCGCAGCCAGTATCTTCCACTACGTGGGTGTCGGt  >  1:3228630/1‑68 (MQ=255)
                                                          aaCCGCTCGCAGCTGTTGGCTTTGCGGCTACCTTCGCAGCCAGTATCTTCCACTACGTGGGTGTCGGt  >  1:2801665/1‑68 (MQ=255)
                                                          aaCCGCTCGCAGCTGTTGGCTTTGCGGCTACCTTCGCAGCCAGTATCTTCCACTACGTGGGTGTCGGt  >  1:265994/1‑68 (MQ=255)
                                                          aaCCGCTCGCAGCTGTTGGCTTTGCGGCTACCTTCGCAGCCAGTATCTTCCACTACGTGGGTGTCGGt  >  1:2473295/1‑68 (MQ=255)
                                                          aaCCGCTCGCAGCTGTTGGCTTTGCGGCTACCTTCGCAGCCAGTATCTTCCACTACGTGGGTGTCGGt  >  1:1019586/1‑68 (MQ=255)
                                                          aaCCGCTCGCAGCTGTTGGCTTTGCGGCTACCTTCGCAGCCAGTATCTTCCACTACGTGGGTGTCGGt  >  1:2119399/1‑68 (MQ=255)
                                                          aaCCGCTCGCAGCTGTTGGCTTTGCGGCTACCTTCGCAGCCAGTATCTTCCACTACGTGGGTGTCGGt  >  1:2101183/1‑68 (MQ=255)
                                                          aaCCGCTCGCAGCTGTTGGCTTTGCGGCTACCTTCGCAGCCAGTATCTTCCACTACGTGGGTGTCGGt  >  1:1996880/1‑68 (MQ=255)
                                                          aaCCGCTCGCAGCTGTTGGCTTTGCGGCTACCTTCGCAGCCAGTATCTTCCACTACGTGGGTGTCGGt  >  1:1923212/1‑68 (MQ=255)
                                                          aaCCGCTCGCAGCTGTTGGCTTTGCGGCTACCTTCGCAGCCAGTATCTTCCACTACGTGGGTGTCGGt  >  1:1484520/1‑68 (MQ=255)
                                                          aaCCGCTCGCAGCTGTTGGCTTTGCGGCTACCTTCGCAGCCAGTATCTTCCACTACGTGGGTGTCGGt  >  1:1403917/1‑68 (MQ=255)
                                                          aaCCGCTCGCAGCTGTTGGCTTTGCGGCTACCTTCGCAGCCAGTATCTTCCACTACGTGGGTGTCGGt  >  1:1403755/1‑68 (MQ=255)
                                                          aaCCGCTCGCAGCTGTTGGCTTTGCGGCTACCTTCGCAGCCAGTATCTTCCACTACGTGGGTGTCGGt  >  1:1270016/1‑68 (MQ=255)
                                                          aaCCGCTCGCAGCTGTTGGCTTTGCGGCTACCTTCGCAGCCAGTATCTTCCACTACGTGGGTGTCGGt  >  1:107059/1‑68 (MQ=255)
                                                          aaCCGCTCGCAGCTGTTGGCTTTGCGGCTACCTTCGCAGCCAGTATCTTCCACTACGTGGGTGTCGGt  >  1:1049945/1‑68 (MQ=255)
                                                                    |                                                         
GGCTTGCCGGAGAACCCGGAAATCAGCGAAACCGTGAAATTCTGGAAAGGCATCTGGAAACCGCTCGCGGCTGTTGGCTTTGCGGCTACCTTTGCGGCCAGTATCTTCCACTACGTGGGTGTCGGT  >  minE/2336329‑2336454

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: