Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,841,988 A→G 70.0% G64G (GGT→GGC frdB ← fumarate reductase (anaerobic), Fe‑S subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,841,9880AG70.0% 82.3 / 35.7 84G64G (GGT→GGCfrdBfumarate reductase (anaerobic), Fe‑S subunit
Reads supporting (aligned to +/- strand):  ref base A (19/6);  new base G (22/36);  total (41/43)
Fisher's exact test for biased strand distribution p-value = 1.84e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACGCAGGAAGGTTTTACATGCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCT  >  minE/2841930‑2842045
                                                          |                                                         
cACGCAGGAAGGTTTTACATGCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATc                                                  >  1:23984/1‑68 (MQ=255)
cACGCAGGAAGGTTTTACATGCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATc                                                  >  1:28537/1‑68 (MQ=255)
 aCGCAGGAAGGTTTTACATGCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGc                                                <  1:875413/69‑1 (MQ=255)
         aGGTTTTACATGCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGc                                        >  1:563039/1‑69 (MQ=255)
                   tGCCAGTTTTGGCACGTTTTAACCATCATGCCGCAGGAGCCACAAATCGCCATACGGCAGGACCAGCgg                             <  1:2892156/69‑1 (MQ=255)
                    gCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAGCCACAAATCGCCATACGGCAGGACCAGCgg                             <  1:2020885/69‑1 (MQ=255)
                    gCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAGCCACAAATCGCCATACGGCAGGACCAGCgg                             <  1:2577514/69‑1 (MQ=255)
                    gCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAGCCACAAATCGCCATACGGCAGGACCAGCgg                             <  1:2545048/69‑1 (MQ=255)
                    gCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAGCCACAAATCGCCATACGGCAGGACCAGCgg                             <  1:274951/69‑1 (MQ=255)
                    gCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAGCCACAAATCGCCATACGGCAGGACCAGCgg                             <  1:281115/69‑1 (MQ=255)
                    gCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAGCCACAAATCGCCATACGGCAGGACCAGCgg                             <  1:3152857/69‑1 (MQ=255)
                    gCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAGCCACAAATCGCCATACGGCAGGACCAGCgg                             <  1:3220300/69‑1 (MQ=255)
                    gCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAGCCACAAATCGCCATACGGCAGGACCAGCgg                             <  1:3263271/69‑1 (MQ=255)
                    gCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAGCCACAAATCGCCATACGGCAGGACCAGCgg                             <  1:3345050/69‑1 (MQ=255)
                    gCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAGCCACAAATCGCCATACGGCAGGACCAGCgg                             <  1:3364026/69‑1 (MQ=255)
                    gCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAGCCACAAATCGCCATACGGCAGGACCAGCgg                             <  1:3414257/69‑1 (MQ=255)
                    gCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAGCCACAAATCGCCATACGGCAGGACCAGCgg                             <  1:503442/69‑1 (MQ=255)
                    gCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAGCCACAAATCGCCATACGGCAGGACCAGCgg                             <  1:827379/69‑1 (MQ=255)
                    gCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAGCCACAAATCGCCATACGGCAGGACCAGCgg                             <  1:163400/69‑1 (MQ=255)
                    gCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAGCCACAAATCGCCATACGGCAGGACCAGCgg                             <  1:1438489/69‑1 (MQ=255)
                    gCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAGCCACAAATCGCCATACGGCAGGACCAGCgg                             <  1:1344187/69‑1 (MQ=255)
                     ccAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAGCCACAAATCGCCATACGGCAGGACCAGCgg                             <  1:3046731/68‑1 (MQ=255)
                           ttGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAgg                                     >  1:3240921/1‑54 (MQ=255)
                                    tgttAACCATCATGCCGCAGGAGCCACAAATCGCCATACGGCGGGACCAGCGGTAGCTCAGGTCCGGt              >  1:732098/1‑68 (MQ=255)
                                    tgttAACCATCATGCCGCAGGAGCCACAAATCGCCATACGGCAGGAc                                   >  1:8559/1‑47 (MQ=255)
                                    tgttAACCATCATGCCGCAGGAGCCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGt              >  1:3084120/1‑68 (MQ=255)
                                    tgttAACCATCATGCCGCAGGAGCCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGt              >  1:2809497/1‑68 (MQ=255)
                                    tgttAACCATCATGCCGCAGGAGCCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGt              >  1:2348648/1‑68 (MQ=255)
                                    tgttAACCATCATGCCGCAGGAGCCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGt              >  1:1530288/1‑68 (MQ=255)
                                    tgttAACCATCATGCCGCAGGAGCCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGt              >  1:1481245/1‑68 (MQ=255)
                                    tgttAACCATCATGCCGCAGGAGCCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGt              >  1:67687/1‑68 (MQ=255)
                                    tgttAACCATCATGCCGCAGGAGCCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGt              >  1:2860916/1‑68 (MQ=255)
                                    tgttAACCATCATGCCGCAGGAGCCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGt              >  1:2668304/1‑68 (MQ=255)
                                    tgttAACCATCATGCCGCAGGAGCCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGt              >  1:2736609/1‑68 (MQ=255)
                                    tgttAACCATCATGCCGCAGGAGCCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGt              >  1:144809/1‑68 (MQ=255)
                                    tgttAACCATCATGCCGCAGGAGCCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGt              >  1:1433237/1‑68 (MQ=255)
                                    tgttAACCATCATGCCGCAGGAGCCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGt              >  1:2116325/1‑68 (MQ=255)
                                    tgttAACCATCATGCCGCAGGAGCCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGt              >  1:467696/1‑68 (MQ=255)
                                    tgttAACCATCATGCCGCAGGAGCCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGt              >  1:2003235/1‑68 (MQ=255)
                                    tgttAACCATCATGCCGCAGGAGCCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGt              >  1:853828/1‑68 (MQ=255)
                                    tgttAACCATCATGCCGCAGGAGCCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGt              >  1:1955073/1‑68 (MQ=255)
                                    tgttAACCATCATGCCGCAGGAGCCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGt              >  1:1946420/1‑68 (MQ=255)
                                    tgttAACCATCATGCCGCAGGAGCCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGt              >  1:1670806/1‑68 (MQ=255)
                                    tgttAACCATCATGCCGCAGGAGCCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGt              >  1:1637802/1‑68 (MQ=255)
                                    tgttAACCATCATGCCGCAGGAGCCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTg             <  1:3012820/69‑1 (MQ=255)
                                    tgttAACCATCATGCCGCAGGAGCCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTg             <  1:335037/69‑1 (MQ=255)
                                    tgttAACCATCATGCCGCAGGAGCCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTg             <  1:44418/69‑1 (MQ=255)
                                    tgttAACCATCATGCCGCAGGAGCCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTg             <  1:2864890/69‑1 (MQ=255)
                                    tgttAACCATCATGCCGCAGGAGCCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTg             <  1:1217143/69‑1 (MQ=255)
                                    tgttAACCATCATGCCGCAGGAGCCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTg             <  1:76836/69‑1 (MQ=255)
                                    tgttAACCATCATGCCGCAGGAGCCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTg             <  1:2663467/69‑1 (MQ=255)
                                    tgttAACCATCATGCCGCAGGAGCCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTg             <  1:2603998/69‑1 (MQ=255)
                                    tgttAACCATCATGCCGCAGGAGCCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTg             <  1:2571342/69‑1 (MQ=255)
                                    tgttAACCATCATGCCGCAGGAGCCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTg             <  1:2468115/69‑1 (MQ=255)
                                    tgttAACCATCATGCCGCAGGAGCCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTg             <  1:2465178/69‑1 (MQ=255)
                                    tgttAACCATCATGCCGCAGGAGCCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTg             <  1:2266517/69‑1 (MQ=255)
                                    tgttAACCATCATGCCGCAGGAGCCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTg             <  1:2258659/69‑1 (MQ=255)
                                    tgttAACCATCATGCCGCAGGAGCCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTg             <  1:208700/69‑1 (MQ=255)
                                    tgttAACCATCATGCCGCAGGAGCCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTg             <  1:2235651/69‑1 (MQ=255)
                                    tgttAACCATCATGCCGCAGGACCCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTg             <  1:444276/69‑1 (MQ=255)
                                    tgctAACCATCATGCCGCAGGAGCCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGt              >  1:1174270/4‑68 (MQ=255)
                                        aaCCATCATGCCGCAGGAGCCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTg             <  1:50741/65‑1 (MQ=255)
                                        aaCCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAg         <  1:2761213/69‑1 (MQ=255)
                                        aaCCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAg         <  1:1104797/69‑1 (MQ=255)
                                        aaCCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAg         <  1:63291/69‑1 (MQ=255)
                                        aaCCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAg         <  1:115400/69‑1 (MQ=255)
                                        aaCCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAg         <  1:2800739/69‑1 (MQ=255)
                                            atcatGCCGCAGGAGCCACAAATCGCCATACGGCAGGACCAGCgg                             <  1:1968220/45‑1 (MQ=255)
                                               atGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGTTAGCTCAGGTCCGGTGCCAGGTTGTCt  >  1:1411295/1‑69 (MQ=255)
                                               atGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCt  >  1:1096931/1‑69 (MQ=255)
                                               atGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCt  >  1:3393967/1‑69 (MQ=255)
                                               atGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCt  >  1:1262078/1‑69 (MQ=255)
                                               atGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCt  >  1:3318702/1‑69 (MQ=255)
                                               atGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCt  >  1:1373864/1‑69 (MQ=255)
                                               atGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCt  >  1:1453434/1‑69 (MQ=255)
                                               atGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCt  >  1:1479457/1‑69 (MQ=255)
                                               atGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCt  >  1:1043021/1‑69 (MQ=255)
                                               atGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCt  >  1:2142472/1‑69 (MQ=255)
                                               atGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCt  >  1:1536165/1‑69 (MQ=255)
                                               atGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCt  >  1:1980709/1‑69 (MQ=255)
                                               atGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCt  >  1:188460/1‑69 (MQ=255)
                                               atGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCt  >  1:1853898/1‑69 (MQ=255)
                                               atGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGAGGTAGCTCAGGACCGGTGCCAGGTTGTCt  >  1:1978637/1‑69 (MQ=255)
                                                tGCCGCAGGAGCCACAAATCGCCATACGGCAGGACCAGCgg                             <  1:3047083/41‑1 (MQ=255)
                                                          |                                                         
CACGCAGGAAGGTTTTACATGCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCT  >  minE/2841930‑2842045

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: