Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,015,675 C→A 23.7% A132S (GCG→TCG)  fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,015,6750CA23.7% 100.3 / 17.4 46A132S (GCG→TCG) fumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base C (25/10);  new base A (3/8);  total (28/18)
Fisher's exact test for biased strand distribution p-value = 1.38e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGG  >  minE/1015605‑1015740
                                                                      |                                                                 
taCATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCgcgcgc                                                                   >  1:4243451/1‑71 (MQ=255)
 aCATATCCAGCGGCGCGTTTTGCGACTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCgcgcgc                                                                   >  1:3526384/1‑70 (MQ=255)
  cATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCa                                                                 <  1:4582949/71‑1 (MQ=255)
     atCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAgcg                                                              >  1:3766262/1‑71 (MQ=255)
       ccAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGcgcc                                                            >  1:1071204/1‑71 (MQ=255)
       ccAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGcgcc                                                            >  1:5688093/1‑71 (MQ=255)
       ccAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGcgcc                                                            >  1:1919440/1‑71 (MQ=255)
        cAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGcgccg                                                           >  1:967218/1‑71 (MQ=255)
        cAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGcgccg                                                           >  1:4531729/1‑71 (MQ=255)
           cggcgCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCtt                                                        <  1:2243750/71‑1 (MQ=255)
                   tttGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATcacca                                                >  1:163289/1‑71 (MQ=255)
                   tttGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATcacca                                                >  1:327459/1‑71 (MQ=255)
                    ttGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATcaccac                                               >  1:4547975/1‑71 (MQ=255)
                         aGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAgcgc                                                             <  1:3941957/52‑1 (MQ=255)
                          gTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCgg                                         >  1:1569488/1‑71 (MQ=255)
                              cgcAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCa                                     <  1:2678046/71‑1 (MQ=255)
                                cAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCaca                                   >  1:5239332/1‑69 (MQ=255)
                                cAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATa                                   >  1:707321/1‑71 (MQ=255)
                                  ggttATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCacac                                  >  1:2203505/3‑70 (MQ=255)
                                   aTTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTcc                                      <  1:141177/65‑1 (MQ=255)
                                    ttATCCTCGATATAAGTGTTATAGACACCGCGCGcc                                                                  >  1:2181820/1‑36 (MQ=255)
                                    ttATCCTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATcacc                                                 >  1:2000193/1‑53 (MQ=255)
                                     tATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACgcgc                              >  1:2326450/1‑71 (MQ=255)
                                      aTCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTtca                                                      >  1:4562658/1‑46 (MQ=255)
                                       tCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTg                            >  1:3654795/1‑71 (MQ=255)
                                        cTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATcaccgc                                               <  1:963907/51‑3 (MQ=255)
                                        cTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATcaccaccac                                            <  1:2301525/54‑1 (MQ=255)
                                            gATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCtt                       <  1:1704682/71‑1 (MQ=255)
                                            gATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTtcat                                                     >  1:4994944/1‑41 (MQ=255)
                                            gATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTtcat                                                     >  1:5581721/1‑41 (MQ=255)
                                            gATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCAt                                    >  1:1070353/1‑58 (MQ=255)
                                             atatAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCttt                      <  1:2118962/71‑1 (MQ=255)
                                             atatAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCttt                      <  1:1647556/71‑1 (MQ=255)
                                             atatAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCttt                      <  1:1357801/71‑1 (MQ=255)
                                              tataGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCttt                      <  1:1981250/70‑1 (MQ=255)
                                               ataAGTGTTATAGACACCGCGCGCCAGCGCCGCTCCATCACCACCACCGGTCCATACACGCTGCCCttttt                    >  1:4085079/1‑71 (MQ=255)
                                                taAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCa                                     <  1:4237096/53‑1 (MQ=255)
                                                  aGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTAc                 >  1:2649153/1‑71 (MQ=255)
                                                   gtgtTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTAcc                <  1:5075048/71‑1 (MQ=255)
                                                        ataGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACcaaca            >  1:5127351/1‑70 (MQ=255)
                                                         taGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACcaacaa           >  1:2528372/1‑70 (MQ=255)
                                                         taGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACcaacaa           >  1:3402511/1‑70 (MQ=255)
                                                         taGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACaat          <  1:4727397/71‑1 (MQ=255)
                                                           gACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACaataa        >  1:3811583/1‑71 (MQ=255)
                                                            acacCGCGCGCCAGCGCCGCTTCATCACCACCCCCGGTCCATACCCGCTGCCCTTTTTTACCAACaataat       >  1:520643/1‑71 (MQ=255)
                                                             cacCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATc      <  1:3447993/71‑1 (MQ=255)
                                                              acCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATcg     <  1:2695996/71‑1 (MQ=255)
                                                               ccGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATcgc    >  1:3982361/1‑71 (MQ=255)
                                                                cgcgcgCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATcg     <  1:3658047/69‑1 (MQ=255)
                                                                 gcgcgcCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCgg  >  1:232508/1‑71 (MQ=255)
                                                                  cgcgcCAGCGCCGCTTCATCACCACCACCGGTCCATa                                   >  1:435346/1‑37 (MQ=255)
                                                                  cgcgcCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCttttt                    >  1:5052715/1‑52 (MQ=255)
                                                                     gcCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGccc                         <  1:3002665/44‑1 (MQ=255)
                                                                      ccAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGccc                         <  1:4759638/43‑1 (MQ=255)
                                                                      |                                                                 
TACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGG  >  minE/1015605‑1015740

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: