Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,015,655 A→G 39.1% T138T (ACT→ACC fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,015,6550AG39.1% 10.6 / 21.4 23T138T (ACT→ACCfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base A (5/9);  new base G (5/4);  total (10/13)
Fisher's exact test for biased strand distribution p-value = 4.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

TTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCT  >  minE/1015593‑1015713
                                                              |                                                          
ttCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTtata                                                    <  1:1956898/71‑1 (MQ=255)
    cttcttTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAg                                                   >  1:306349/1‑68 (MQ=255)
      tcttTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGcc                                              >  1:1308063/1‑71 (MQ=255)
      tcttTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGcc                                              >  1:3695163/1‑71 (MQ=255)
      tcttTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACAcc                                              >  1:1078024/1‑71 (MQ=255)
          tgtaCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGcct                                             <  1:228787/66‑2 (MQ=255)
              cATATCCAGCGGTGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCa                                      <  1:3592687/71‑1 (MQ=255)
              cATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCa                                      <  1:1279535/71‑1 (MQ=255)
               atatCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAg                                                           <  1:762288/49‑1 (MQ=255)
                tatCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAgc                                    <  1:3927639/71‑1 (MQ=255)
                  tCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGcctt                                            <  1:184665/61‑3 (MQ=255)
                         gcgcGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTtc                            <  1:1202475/70‑1 (MQ=255)
                          cgcaTTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAg                                                   <  1:3494254/46‑1 (MQ=255)
                           gcgTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTtata                                                    <  1:1139577/44‑1 (MQ=255)
                               tttGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAg                                                   >  1:912376/1‑41 (MQ=255)
                               tttGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAg                                     >  1:745678/1‑55 (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)
                                               aTTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGcgcc                                 <  1:472427/43‑1 (MQ=255)
                                                  aTCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCt  >  1:1774856/1‑71 (MQ=255)
                                                  aTCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCt  >  1:969383/1‑71 (MQ=255)
                                                              aGTGTTATAGACACCGCGCGCCAGCGCCGCTtcatc                         >  1:3317241/1‑36 (MQ=255)
                                                              |                                                          
TTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCT  >  minE/1015593‑1015713

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: