Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,842,244 G→A 73.1% R579R (CGC→CGT frdA ← fumarate reductase (anaerobic) catalytic and NAD/flavoprotein subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,842,2440GA73.1% 35.8 / 11.6 37R579R (CGC→CGTfrdAfumarate reductase (anaerobic) catalytic and NAD/flavoprotein subunit
Reads supporting (aligned to +/- strand):  ref base G (3/7);  new base A (16/11);  total (19/18)
Fisher's exact test for biased strand distribution p-value = 1.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCGTAGTAATCTTCACGTCGCTGTA  >  minE/2842181‑2842285
                                                               |                                         
cGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAAcgc                                          >  1:3312610/1‑65 (MQ=255)
cGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAAcgc                                          >  1:1548656/1‑65 (MQ=255)
cGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAAcgc                                          >  1:735334/1‑65 (MQ=255)
cGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGttt                                      >  1:2372073/1‑69 (MQ=255)
     tctccttctTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAAcacg                                         <  1:1004136/61‑1 (MQ=255)
     tctccttctTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAAcacg                                         <  1:398381/61‑1 (MQ=255)
     tctccttctTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAAcacg                                         <  1:631909/61‑1 (MQ=255)
     tctccttctTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAAcacg                                         <  1:2703279/61‑1 (MQ=255)
     tctccttctTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAAcacg                                         <  1:2699688/61‑1 (MQ=255)
     tctccttctTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAAcacg                                         <  1:266620/61‑1 (MQ=255)
     tctccttctTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAAcacg                                         <  1:2663740/61‑1 (MQ=255)
     tctccttctTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAAcacg                                         <  1:2396295/61‑1 (MQ=255)
     tctccttctTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAAcacg                                         <  1:22932/61‑1 (MQ=255)
     tctccttctTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAAcacg                                         <  1:1162018/61‑1 (MQ=255)
     tctccttctTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAAcacg                                         <  1:1483376/61‑1 (MQ=255)
          ttcttATTGTCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAAcacg                                         <  1:531507/56‑1 (MQ=255)
                  ggCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCgtagt                    >  1:3191207/1‑69 (MQ=255)
                    ctgctTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCgtagta                   <  1:1964307/68‑1 (MQ=255)
                         tCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCGTAGTAATCtt              <  1:2504462/68‑1 (MQ=255)
                         tCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCGTAGTAATCtt              <  1:701078/68‑1 (MQ=255)
                         tCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCGTAGTAATCtt              <  1:793831/68‑1 (MQ=255)
                         tCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCGTAGTAATCtt              <  1:1608659/68‑1 (MQ=255)
                         tCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCGTAGTAATCtt              <  1:2029431/68‑1 (MQ=255)
                         tCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCGTAGTAATCtt              <  1:1741277/68‑1 (MQ=255)
                               ttATCGGCTGCATCCGCTTCGCCACCGTAAACACGTTTAGCTGGCGGCAGCGTAGTAATCTTCACGTCg       <  1:1665736/69‑1 (MQ=255)
                               ttATCGGCTGCATCCGCTTCGCCACCGTAAACACGTTTAGCTGGCGGCAGCGTAGTAATCTTCACGTCg       <  1:1586192/69‑1 (MQ=255)
                               ttATCGGCTGCATCCGCTTCGCCACCGTAAACACGTTTAGCTGGCGGCAGCGTAGTAATCTTCACGTCg       <  1:3434974/69‑1 (MQ=255)
                               ttATCGGCTGCATCCGCTTCGCCACCGTAAACACGTTTAGCTGGCGGCAGCGTAGTAATCTTCACGTCg       <  1:1532542/69‑1 (MQ=255)
                               ttATCGGCTGCATCCGCTTCGCCACCGTAAACACGTTTAGCTGGCGGCAGCGTAGTAATCTTCACGTCg       <  1:1288914/69‑1 (MQ=255)
                               ttATCGGCTGCATCCGCTTCGCCACCGTAAACACGTTTAGCTGGCGGCAGCGTAGTAATCTTCACGTCg       <  1:1191011/69‑1 (MQ=255)
                               ttATCGGCTGCATCCGCTTCGCCACCGTAAACACGTTTAGCTGGCGGCAGCGTAGTAATCTTCACGTCg       <  1:285514/69‑1 (MQ=255)
                               ttATCGGCTGCATCCGCTTCGCCACCGTAAACACGTTTAGCTGGCGGCAGCGTAGTAATCTTCACGTCg       <  1:444172/69‑1 (MQ=255)
                               ttATCGGCTGCATCCGCTTCGCCACCGTAAACACGTTTAGCTGGCGGCAGCGTAGTAATCTTCACGTCg       <  1:3148767/69‑1 (MQ=255)
                               ttATCGGCTGCATACGCTTCGCCACCGTAAACACGTTTAGCTGGCGGCAGCGTAGTAATCTTCACGTCg       <  1:2019459/69‑1 (MQ=255)
                                tATCGGCTGCAGCCGCTTCGCCACCGTAAACACGTTTAGCTGGCGGCAGCGTAGTAATCTTCACGTCg       <  1:195775/68‑1 (MQ=255)
                                    ggCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCGTAGTAATCTTCACGTCGCTGTa  >  1:1811547/1‑69 (MQ=255)
                                                    ccACCGTAAACACGTTTAGCTGGCGGCAGCGTAGTaa                  >  1:805424/1‑37 (MQ=255)
                                                    ccACCGTAAACACGTTTAGCTGGCGGCAGCGTAGTaa                  >  1:785876/1‑37 (MQ=255)
                                                    ccACCGTAAACACGTTTAGCTGGCGGCAGCGTAGTaa                  >  1:982665/1‑37 (MQ=255)
                                                    ccACCGTAAACACGTTTAGCTGGCGGCAGCGTAGTaa                  >  1:3228804/1‑37 (MQ=255)
                                                    ccACCGTAAACACGTTTAGCTGGCGGCAGCGTAGTaa                  >  1:3228551/1‑37 (MQ=255)
                                                    ccACCGTAAACACGTTTAGCTGGCGGCAGCGTAGTaa                  >  1:3224571/1‑37 (MQ=255)
                                                    ccACCGTAAACACGTTTAGCTGGCGGCAGCGTAGTaa                  >  1:3216640/1‑37 (MQ=255)
                                                    ccACCGTAAACACGTTTAGCTGGCGGCAGCGTAGTaa                  >  1:3001549/1‑37 (MQ=255)
                                                    ccACCGTAAACACGTTTAGCTGGCGGCAGCGTAGTaa                  >  1:2976073/1‑37 (MQ=255)
                                                    ccACCGTAAACACGTTTAGCTGGCGGCAGCGTAGTaa                  >  1:2830892/1‑37 (MQ=255)
                                                    ccACCGTAAACACGTTTAGCTGGCGGCAGCGTAGTaa                  >  1:2523402/1‑37 (MQ=255)
                                                    ccACCGTAAACACGTTTAGCTGGCGGCAGCGTAGTaa                  >  1:2397553/1‑37 (MQ=255)
                                                    ccACCGTAAACACGTTTAGCTGGCGGCAGCGTAGTaa                  >  1:2292323/1‑37 (MQ=255)
                                                    ccACCGTAAACACGTTTAGCTGGCGGCAGCGTAGTaa                  >  1:1589662/1‑37 (MQ=255)
                                                    ccACCGTAAACACGTTTAGCTGGCGGCAGCGTAGTaa                  >  1:1203825/1‑37 (MQ=255)
                                                    ccACCGTAAACACGTTTAGCTGGCGGCAGCGTAGTaa                  >  1:1036276/1‑37 (MQ=255)
                                                               |                                         
CGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCGTAGTAATCTTCACGTCGCTGTA  >  minE/2842181‑2842285

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: