Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,014,752 C→G 19.6% T439T (ACG→ACC fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,014,7520CG19.6% 116.8 / 20.7 56T439T (ACG→ACCfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base C (27/18);  new base G (7/4);  total (34/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01

TCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATG  >  minE/1014682‑1014820
                                                                      |                                                                    
tCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTc                                                                      >  1:4617686/1‑71 (MQ=255)
  gCCAGCATGATCATGCTGCCGCCGTGGGATTGCAGCAGATCCACGTAGGAGTCCATACGGCCTGCggtggt                                                                    <  1:208187/71‑1 (MQ=255)
   ccAGCATGATCATGCTGCCGCCGTGGGATTGCAGCAGATCCACGTAGGAGTCCATACGGCCTGCGGTGGtt                                                                   <  1:5185377/71‑1 (MQ=255)
   ccAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGt                                                                    <  1:3198675/70‑1 (MQ=255)
     aGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTgg                                                                 <  1:3163240/71‑1 (MQ=255)
     aGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTgg                                                                 >  1:1936462/1‑71 (MQ=255)
       cATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGt                                                                    <  1:5507567/66‑1 (MQ=255)
       cATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGc                                                               >  1:2295277/1‑71 (MQ=255)
         tgttCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCa                                                             <  1:3217070/68‑1 (MQ=255)
         tGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTg                                                                  >  1:1444583/1‑66 (MQ=255)
         tGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGcc                                                              <  1:4603298/70‑1 (MQ=255)
          gATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGtt                                                                   <  1:2871881/64‑1 (MQ=255)
           aTCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGcc                                                              >  1:520481/1‑68 (MQ=255)
            tCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCg                                                                     >  1:4616898/1‑60 (MQ=255)
              aTACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAagag                                                         >  1:5288592/1‑70 (MQ=255)
                 ctgccgccgTGGGATTGCAGCAGATCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGGCCAAGTGAAc                                                      >  1:1768185/4‑70 (MQ=255)
                   tCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCg                                                    >  1:3303521/1‑70 (MQ=255)
                    ccgccgTGGGATTGCAGCAGATCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGGCCAAGTGAAcctg                                                   <  1:2496049/70‑3 (MQ=255)
                    ccGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCAGTCCGGCGGTCGTTGGGCCAAGAGAACCg                                                    >  1:532879/1‑69 (MQ=255)
                         cTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAaga                                                          >  1:5226654/1‑58 (MQ=255)
                         cTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAaga                                                          >  1:3919429/1‑58 (MQ=255)
                         cTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAaga                                                          >  1:3284340/1‑58 (MQ=255)
                         cTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAaga                                                          >  1:73698/1‑58 (MQ=255)
                         cTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAaga                                                          >  1:2758279/1‑58 (MQ=255)
                               ttGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGt                                                                    <  1:1779207/42‑1 (MQ=255)
                               ttGCAGCAGATCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGGCCAAGTGAAcctg                                                   <  1:2145841/59‑3 (MQ=255)
                                tGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTc                                      >  1:1543561/1‑71 (MQ=255)
                                tGCAGCAGATCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGGCCAAgtga                                                        <  1:5144256/53‑1 (MQ=255)
                                tGCAGCAGATCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAAccgg                                       >  1:3406880/1‑68 (MQ=255)
                                 gCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCg                                                                     <  1:3510167/39‑1 (MQ=255)
                                 gCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCg                                                                     <  1:2092348/39‑1 (MQ=255)
                                 gCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGt                                                                    >  1:1098524/1‑40 (MQ=255)
                                 gCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGt                                                                    >  1:1125083/1‑40 (MQ=255)
                                 gCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGt                                                                    >  1:61690/1‑40 (MQ=255)
                                 gCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTcc                                     >  1:2276558/1‑71 (MQ=255)
                                     ttGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCaa                                                            <  1:2636206/44‑1 (MQ=255)
                                     ttGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCaa                                                            <  1:3044223/44‑1 (MQ=255)
                                      agATCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAAccg                                        >  1:1166034/2‑62 (MQ=255)
                                      agATCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAAcc                                         >  1:4609802/2‑62 (MQ=255)
                                       gATCGACATAAGAATCCATCCGTCCGGCGGTCGTTg                                                                  >  1:1441623/1‑36 (MQ=255)
                                       gATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCaa                                                            >  1:4746208/1‑42 (MQ=255)
                                       gATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAaga                                                          >  1:5664825/1‑44 (MQ=255)
                                       gATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAag                                                           >  1:63370/1‑43 (MQ=255)
                                       gATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGt                                <  1:2935683/70‑1 (MQ=255)
                                         tCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGttt                              >  1:116620/1‑70 (MQ=255)
                                            aCATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGc                          <  1:2207995/71‑1 (MQ=255)
                                              aTAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCg                        >  1:1045288/1‑71 (MQ=255)
                                                aaGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGa                      >  1:3023545/1‑71 (MQ=255)
                                                  gAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCAt                                             >  1:429418/1‑46 (MQ=255)
                                                    atccatccGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGAccc                   >  1:2789702/1‑70 (MQ=255)
                                                         tccgtccgGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGttt                              <  1:3929613/54‑1 (MQ=255)
                                                           cgtcccGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCgtagta           <  1:3608595/71‑1 (MQ=255)
                                                               ctgcggTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCgtag             >  1:2669235/3‑65 (MQ=255)
                                                                tgcggTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCgtag             >  1:3835326/2‑64 (MQ=255)
                                                                ggcggTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCg      <  1:1963437/71‑1 (MQ=255)
                                                                  cggTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCggg    >  1:3979551/1‑70 (MQ=255)
                                                                   ggTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGgat   <  1:2927941/71‑1 (MQ=255)
                                                                    gTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGgatg  <  1:1281739/71‑1 (MQ=255)
                                                                      |                                                                    
TCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATG  >  minE/1014682‑1014820

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: