Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,015,683 G→T 23.7% A129E (GCG→GAG) ‡ fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,015,6830GT23.7% 39.6 / 4.7 21A129E (GCG→GAG) ‡fumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base G (9/7);  new base T (3/2);  total (12/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTAT  >  minE/1015618‑1015749
                                                                 |                                                                  
gcgcGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTtc                                                                <  1:1202475/70‑1 (MQ=255)
         gCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATcaccacca                                                      >  1:3236232/1‑71 (MQ=255)
         gCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATcaccac                                                        <  1:3645300/69‑1 (MQ=255)
             atAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCgg                                                  <  1:997455/70‑1 (MQ=255)
                 cgcAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTtcatca                                                            >  1:649220/1‑57 (MQ=255)
                         aTCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCt                                      >  1:969383/1‑71 (MQ=255)
                         aTCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCt                                      >  1:1774856/1‑71 (MQ=255)
                                     aGTGTTATAGACACCGCGCGCCAGCGCCGCTtcatc                                                             >  1:3317241/1‑36 (MQ=255)
                                      gtgtTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTAcc                         <  1:3023536/71‑1 (MQ=255)
                                        gtTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACcaa                       >  1:3869561/1‑71 (MQ=255)
                                         ttATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACcaa                       <  1:3554671/70‑1 (MQ=255)
                                           ataGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTc                                                <  1:3671896/43‑1 (MQ=255)
                                             aGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATaca                                          >  1:2230814/1‑47 (MQ=255)
                                             aGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCtttt                              >  1:2434756/1‑59 (MQ=255)
                                             aGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACaat                   >  1:449900/1‑70 (MQ=255)
                                               acacCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTCTTACCAACaata                  <  1:612243/69‑1 (MQ=255)
                                                 acCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACca                        <  1:670534/61‑1 (MQ=255)
                                                   cgcgcgCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATcgcg            <  1:1651817/71‑1 (MQ=255)
                                                        gcCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACaataat                >  1:3655628/1‑62 (MQ=255)
                                                        gcCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACaataat                >  1:3707824/1‑62 (MQ=255)
                                                        gcCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGcc       >  1:3082898/1‑71 (MQ=255)
                                                             cgccgcTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTAt  >  1:3474463/1‑71 (MQ=255)
                                                             cgccgcTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTAt  >  1:2402544/1‑71 (MQ=255)
                                                              gccgcTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTa   >  1:1417215/1‑69 (MQ=255)
                                                                cgcTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACaat                   <  1:701938/51‑1 (MQ=255)
                                                                 |                                                                  
GCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTAT  >  minE/1015618‑1015749

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: