Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,015,667 A→G 36.0% G134G (GGT→GGC fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,015,6670AG36.0% 17.3 / 20.9 25G134G (GGT→GGCfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base A (8/8);  new base G (5/4);  total (13/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01

TCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATA  >  minE/1015599‑1015733
                                                                    |                                                                  
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)
          tatCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAgc                                                        <  1:3927639/71‑1 (MQ=255)
            tCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGcctt                                                                <  1:184665/61‑3 (MQ=255)
                   gcgcGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTtc                                                <  1:1202475/70‑1 (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)
                                                         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)
                                                                    |                                                                  
TCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATA  >  minE/1015599‑1015733

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: