Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,014,677 G→A 14.5% G464G (GGC→GGT fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,014,6770GA14.5% 161.4 / 17.6 62G464G (GGC→GGTfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base G (28/25);  new base A (6/3);  total (34/28)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

ACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTC  >  minE/1014611‑1014744
                                                                  |                                                                   
gccGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCt                                                                  >  1:320690/2‑70 (MQ=255)
accGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCtt                                                                 <  1:3607258/71‑1 (MQ=255)
 ccGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCttt                                                                >  1:1037621/1‑71 (MQ=255)
 ccGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCttt                                                                >  1:3542282/1‑71 (MQ=255)
  cGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTc                                                               <  1:2295753/71‑1 (MQ=255)
    aTACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCt                                                                  <  1:1736244/66‑1 (MQ=255)
    aTACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGc                                                             <  1:1788725/71‑1 (MQ=255)
     tACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGcc                                                            <  1:1873513/71‑1 (MQ=255)
       cTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAg                                                          <  1:448507/71‑1 (MQ=255)
           cAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGcc                                                            <  1:1702727/65‑1 (MQ=255)
            aaGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTg                                                                     <  1:1357149/55‑1 (MQ=255)
             aGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGAt                                                    >  1:1908248/1‑71 (MQ=255)
              ggTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATc                                                   <  1:1931665/71‑1 (MQ=255)
                  gAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCTTTCGCCAGCATGATCAtg                                                <  1:2703189/70‑2 (MQ=255)
                          cgTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCTTTCGCCAGCATGATCATGCTgccgccg                                       <  1:3049470/71‑2 (MQ=255)
                           gtgtTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCTTTCGCCAGCATGATCATGCTgccgccg                                       <  1:1691351/70‑2 (MQ=255)
                              ttttttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTgcg                                   >  1:3358151/1‑71 (MQ=255)
                               tttttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCt                                      >  1:1178129/1‑67 (MQ=255)
                               tttttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCt                                      >  1:2549758/1‑67 (MQ=255)
                                ttttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCt                                 >  1:2443948/1‑71 (MQ=255)
                                ttttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCt                                 >  1:1725149/1‑71 (MQ=255)
                                ttttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCt                                 >  1:117427/1‑71 (MQ=255)
                                 tttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTgc                                    >  1:66361/1‑67 (MQ=255)
                                 tttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCtt                                >  1:3504993/1‑71 (MQ=255)
                                  ttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCtt                                <  1:63543/70‑1 (MQ=255)
                                    aCAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTg                               >  1:1779442/1‑69 (MQ=255)
                                    aCAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTg                               >  1:3463334/1‑69 (MQ=255)
                                        gCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTg                         <  1:1115746/71‑1 (MQ=255)
                                          aTCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGAt                       >  1:256349/1‑71 (MQ=255)
                                          aTCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGAt                       >  1:688963/1‑71 (MQ=255)
                                          aTCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCACCATGATCATACTTCCGCCCTGCGCTTGCAGTTGa                        >  1:1528744/1‑70 (MQ=255)
                                          aTCCGTCACCTGCTGGCTGCGGTTGCCTTTAGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGAt                       >  1:1630231/1‑71 (MQ=255)
                                           tCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACtt                                             <  1:842147/48‑1 (MQ=255)
                                           tCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATc                      >  1:2316786/1‑71 (MQ=255)
                                           tCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATc                      >  1:3298645/1‑71 (MQ=255)
                                             cGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCg                     <  1:3010248/70‑1 (MQ=255)
                                             cGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGa                    <  1:96148/71‑1 (MQ=255)
                                              gTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTc                                            >  1:3073524/1‑46 (MQ=255)
                                               tCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGtt                          <  1:1729560/63‑1 (MQ=255)
                                                cACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACa                  >  1:1543197/1‑70 (MQ=255)
                                                 aCCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACAt                 >  1:3553830/1‑70 (MQ=255)
                                                 aCCTGCTGACTGCGGTTACCTTTCGCCAGCATGATCATGCTGCCGCCGTGGGATTGCAGCAGATCCAcg                  >  1:1696151/1‑68 (MQ=255)
                                                    tgctgGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGa                    >  1:3612077/1‑64 (MQ=255)
                                                    tgctgGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAg              >  1:1654041/1‑70 (MQ=255)
                                                    tgctgGCTGCGGTTGCCTTTCGCCAGCATGATCAAACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAg              >  1:2290286/1‑70 (MQ=255)
                                                      ctgactgCGGTTACCTTTCGCCAGCATGATCATGCTGCCGCCGTGGGATTGCAGCAGATCCACGTAGGagt           >  1:3391991/1‑69 (MQ=255)
                                                      ctgactgCGGTTACCTTTCGCCAGCATGATCATGCTGCCGCCGTGGGATTGCAGCAGATCCACGTAGGagt           >  1:1372482/1‑69 (MQ=255)
                                                        ggctgCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGc                              >  1:1696675/1‑50 (MQ=255)
                                                        ggctgCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGc                              >  1:860589/1‑50 (MQ=255)
                                                         gctgCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGa                        <  1:3372384/55‑1 (MQ=255)
                                                         gctgCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGa                        <  1:3376674/55‑1 (MQ=255)
                                                         gctgCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAatcc         <  1:1438703/70‑1 (MQ=255)
                                                         cctgCGGTTGCCTTTCGCCAGCATGATCATACTTCCTCCCTGCGCTTGCAGTTGATCGACATAAGAatcca        <  1:2763591/70‑1 (MQ=255)
                                                          ctgCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGa                        <  1:1423252/54‑1 (MQ=255)
                                                          ctgCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAat           <  1:354527/67‑1 (MQ=255)
                                                          ctgCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAat           <  1:2274862/67‑1 (MQ=255)
                                                           tgCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATc                      >  1:516534/1‑55 (MQ=255)
                                                           tgCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATc                      >  1:2548659/1‑55 (MQ=255)
                                                           tgCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAatccatc      >  1:2617295/1‑71 (MQ=255)
                                                               gTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCAtccgtc  <  1:2894135/71‑1 (MQ=255)
                                                               gTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCAtccgtc  <  1:3149736/71‑1 (MQ=255)
                                                                ttGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGc                              <  1:1810521/42‑1 (MQ=255)
                                                                ttGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGAc                   <  1:666455/53‑1 (MQ=255)
                                                                  |                                                                   
ACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTC  >  minE/1014611‑1014744

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: