Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,014,482 A→G 20.0% D529D (GAT→GAC fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,014,4820AG20.0% 119.3 / 24.3 60D529D (GAT→GACfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base A (29/19);  new base G (7/5);  total (36/24)
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.67e-01

CTGTTTATTTCACACAGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCC  >  minE/1014418‑1014544
                                                                |                                                              
ctGTTTATTTCACACAGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTatcatc                                                           >  1:5478054/1‑70 (MQ=255)
ctGTTTATTTCACACAGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTatc                                                              <  1:5480882/67‑1 (MQ=255)
 tGTTTATTTCACACAGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTat                                                               <  1:200366/65‑1 (MQ=255)
 tGTTTATTTCACACAGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTat                                                               <  1:641805/65‑1 (MQ=255)
  gTTTATTTCACACAGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTatcat                                                            >  1:3974054/1‑67 (MQ=255)
  gTTTATTTCACACAGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTatcat                                                            >  1:1669720/1‑67 (MQ=255)
   tttATTTCACACAGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAATAAGTCATTTCCTTTATCATCcac                                                        >  1:4075677/1‑70 (MQ=255)
    ttATTTCACACAGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCcac                                                        <  1:4162374/69‑1 (MQ=255)
    ttATTTCACACAGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCcac                                                        >  1:1964758/1‑69 (MQ=255)
     tATTTCACACAGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCAGAAg                                                     >  1:651327/1‑71 (MQ=255)
     tATTTCACACAGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAg                                                     >  1:3445596/1‑71 (MQ=255)
     tATTTCACACAGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAg                                                     >  1:5554675/1‑71 (MQ=255)
      aTTTCACACAGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAg                                                     >  1:1624172/1‑70 (MQ=255)
         tCACACAGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTat                                                               <  1:4830794/57‑1 (MQ=255)
         tCACACAGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTat                                                               <  1:2043783/57‑1 (MQ=255)
             acaGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATaa                                                >  1:1598061/1‑68 (MQ=255)
             acaGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATaa                                                >  1:4695406/1‑68 (MQ=255)
             acaGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATaa                                                >  1:5647785/1‑68 (MQ=255)
                 cGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCgg                                         <  1:1554401/71‑1 (MQ=255)
                     tGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCcac                                                        <  1:4179032/52‑1 (MQ=255)
                      gCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTatcat                                                            >  1:236799/1‑47 (MQ=255)
                        aTTGTGTGAGTTGTATCTGGTGGAAGATGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCt                                  >  1:3113571/1‑71 (MQ=255)
                         ttGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAgg                                                    >  1:1743950/1‑52 (MQ=255)
                         ttGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATa                                                 <  1:1742119/55‑1 (MQ=255)
                         ttGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATc                                   <  1:1179722/69‑1 (MQ=255)
                         ttGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATc                                   <  1:4372053/69‑1 (MQ=255)
                          tgtgtgAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAAt                                    >  1:1728731/1‑67 (MQ=255)
                          tgtgtgAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAAt                                    >  1:1273370/1‑67 (MQ=255)
                          tgtgtgAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCtt                                 >  1:5460216/1‑70 (MQ=255)
                          tgtgtgAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTc                                >  1:3988852/1‑71 (MQ=255)
                          tgtgtgAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTc                                >  1:4682293/1‑71 (MQ=255)
                           gtgtgAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTc                                >  1:3169026/1‑70 (MQ=255)
                            tgtgAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATa                                                 <  1:3197063/52‑1 (MQ=255)
                            tgtgAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATa                                                 <  1:5010572/52‑1 (MQ=255)
                            tgtgAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATa                                                 <  1:4757560/52‑1 (MQ=255)
                            tgtgAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATa                                                 <  1:3759953/52‑1 (MQ=255)
                            tgtgAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATa                                                 <  1:1142094/52‑1 (MQ=255)
                            tgtgAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTcc                               >  1:5068504/1‑70 (MQ=255)
                             gtgAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTatcat                                                            >  1:1972261/1‑40 (MQ=255)
                             gtgAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTatcat                                                            >  1:3816463/1‑40 (MQ=255)
                             gtgAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTatcat                                                            >  1:3516100/1‑40 (MQ=255)
                             gtgAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCAc                             >  1:4427468/1‑71 (MQ=255)
                               gAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATa                                                 <  1:1465795/49‑1 (MQ=255)
                               gAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATa                                                 <  1:1792050/49‑1 (MQ=255)
                                aGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATa                                                 <  1:813566/48‑1 (MQ=255)
                                   tGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAAt                                    >  1:4483828/1‑58 (MQ=255)
                                   tGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAAt                       >  1:2839048/1‑71 (MQ=255)
                                    gTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAAt                       >  1:565999/1‑70 (MQ=255)
                                      aTCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAAtttt                    >  1:535536/1‑71 (MQ=255)
                                       tttgctgGAAGAAGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATc                                   >  1:5242002/3‑55 (MQ=25)
                                        ttgctgGAAGAAGTCGTTACCTTTGTCATCGACCAGGATAAACGc                                            >  1:3080222/2‑45 (MQ=19)
                                        ttgctgGAAGAAGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTc                                >  1:443686/2‑57 (MQ=25)
                                        ttgctgGAAGAAGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTc                                >  1:3661546/2‑57 (MQ=25)
                                           ctgGAAGAAGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAg                <  1:3389495/70‑1 (MQ=255)
                                              gaagaaGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGa               >  1:4708865/1‑68 (MQ=255)
                                              gaagaaGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGAtag            >  1:4840571/1‑69 (MQ=255)
                                              gaagaaGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATg             >  1:1281765/1‑70 (MQ=255)
                                              gaagaaGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATg             <  1:3396548/70‑1 (MQ=255)
                                              caagaagTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGAtag            >  1:4840576/2‑69 (MQ=255)
                                                  aagTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTc        >  1:3585077/1‑71 (MQ=255)
                                                   agTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTc        <  1:1851522/70‑1 (MQ=255)
                                                     tCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCAt     >  1:1554285/1‑71 (MQ=255)
                                                      cATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCAtt    >  1:3199116/1‑71 (MQ=255)
                                                       aTTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAAtt                      <  1:3480619/52‑1 (MQ=255)
                                                       aTTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGAt              >  1:5655099/1‑60 (MQ=255)
                                                        tttCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTc                                <  1:1187718/41‑1 (MQ=255)
                                                        tttCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTc                                <  1:3929580/41‑1 (MQ=255)
                                                        tttCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTcc  >  1:489388/1‑71 (MQ=255)
                                                         ttCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTc                                <  1:4017569/40‑1 (MQ=255)
                                                         taccTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTc                   <  1:384221/51‑1 (MQ=25)
                                                            cTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTcga                        <  1:4102925/45‑3 (MQ=21)
                                                             tttGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCAt     <  1:5383765/63‑1 (MQ=255)
                                                                |                                                              
CTGTTTATTTCACACAGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCC  >  minE/1014418‑1014544

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: