Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,015,841 G→A 10.6% D76D (GAC→GAT fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,015,8410GA10.6% 130.1 / 5.7 47D76D (GAC→GATfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base G (15/27);  new base A (3/2);  total (18/29)
Fisher's exact test for biased strand distribution p-value = 3.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

TTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAGAATGTCGGCCACCTGTTGTTGGTGCGCCGGACGCAGCATGAACGACGCATCATGAAACGCC  >  minE/1015772‑1015910
                                                                     |                                                                     
tttCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGt                                                                      >  1:821981/1‑71 (MQ=255)
      cgcgATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGAACTTGTCGTTTTCGCTGGCTTCTGGATcgtg                                                                  >  1:1904568/1‑66 (MQ=255)
       gcgATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATcgtg                                                                  >  1:50103/1‑65 (MQ=255)
       gcgATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAg                                                               <  1:1130330/71‑1 (MQ=255)
       gcgATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAg                                                               <  1:284531/71‑1 (MQ=255)
        cgATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAGa                                                              <  1:2946017/71‑1 (MQ=255)
        cgATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAGa                                                              <  1:1570174/71‑1 (MQ=255)
          aTGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAGAAt                                                            <  1:428186/71‑1 (MQ=255)
           tGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCggg                                                                       >  1:340474/1‑59 (MQ=255)
                gAGTTACGCATGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAGAATGTCGGc                                                      <  1:579382/71‑1 (MQ=255)
                gAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAGAATGTCg                                                        >  1:112155/1‑69 (MQ=255)
                gAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAGAATGTCGGc                                                      <  1:2280586/71‑1 (MQ=255)
                gAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAGAATGTCGGc                                                      >  1:2103028/1‑71 (MQ=255)
                  gTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAGAATGTCGGCCa                                                    >  1:50623/1‑71 (MQ=255)
                   ttACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTACGGGTCACGCAGAATGTCGGCCAc                                                   <  1:1390487/71‑1 (MQ=255)
                      cGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAGAATGTCGGCCACCtg                                                <  1:1286828/71‑1 (MQ=255)
                        cAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAGAATGTCGGCCACCtgtt                                              <  1:2228343/71‑1 (MQ=255)
                        cAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCAGGTCACGCAGAATGTCGGCCAcc                                                  <  1:2561581/67‑1 (MQ=255)
                          ggAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAGAATGTCGGcc                                                     <  1:196257/62‑1 (MQ=255)
                          ggAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAGAATGTCGGCCACCtgttgt                                            <  1:3369008/71‑1 (MQ=255)
                           gAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCtgtttct                                           <  1:1757065/71‑4 (MQ=255)
                            aaTTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAGAATGTCGGCCACCtgttg                                             <  1:477566/68‑1 (MQ=255)
                            aaTTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAGAATGTCGGCCACCTGttgttg                                          <  1:277602/71‑1 (MQ=255)
                              ttGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTGGt                                        <  1:1844678/71‑1 (MQ=255)
                              ttGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTGGt                                        >  1:578541/1‑71 (MQ=255)
                              ttGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAGAATGTCGGCCACCtgttgt                                            >  1:1546073/1‑67 (MQ=255)
                                 cAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAGAATGTCGGCCACCTGTTGTTGGTgcg                                     <  1:3243852/71‑1 (MQ=255)
                                 cAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAGAATGTCGGCCACCTGTTGTTGGTgcg                                     <  1:246284/71‑1 (MQ=255)
                                   gcgcCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAGAATGTCGGCCACCTGTTGTTGGTgcgcg                                   >  1:3648842/1‑70 (MQ=255)
                                       cacaTATTTATCATTTTCGCTGGCCTCGGGGTCACGCAGAATGTCGGCCACCTGTTGTTGGTGCGCCGGAc                               >  1:3078885/1‑71 (MQ=255)
                                        acaTATTTATCATTTTCGCTGGCCTCCGGGTCACGCAGAATGTCGGCCACCTGTTGTTGGTGCGCCGGACg                              >  1:3829596/1‑71 (MQ=255)
                                          atatTTATCATTTTCGCTGGCCTCCGGGTCACGCAGAATGTCGGCCACCTGTTGTTGGTGCGCCGGACgca                            <  1:2138846/71‑1 (MQ=255)
                                             tttATCATTTTCGCTGGCCTCCGGGTCACGCAGAATGTc                                                         >  1:1970993/1‑39 (MQ=255)
                                                   aTTTTCGCTGGCCTCCGGGTCACGCAGAATGTCGGCCACCTGTTGTTGGTGCGCCGGACGCAGCATGAac                    <  1:80402/70‑1 (MQ=255)
                                                   aTTTTCGCTGGCCTCCGGGTCACGCAGAATGTCGGCCACCTGTTGTTGGTGCGCCGGACGCAGCATGAac                    >  1:796910/1‑70 (MQ=255)
                                                     tttCGCTGGCCTCCGGGTCACGCAGAATGTCGGCCACCTGTTGTTGGTGCGCCGGACGCAGCATGAacga                  >  1:811236/1‑70 (MQ=255)
                                                         gCTGGCCTCCGGGTCACGCAGAATGTCGGCCACCTGTTGTTGGTGCGCCGGACGCAGCATGAACGACGca              <  1:158835/70‑1 (MQ=255)
                                                         gCTGGCCTCCGGGTCACGCAGAATGTCGGCCACCTGTTGTTGGTGCGCCGGACGCAGCATGAACGACGca              <  1:3807331/70‑1 (MQ=255)
                                                          cTGGCCTCCGGGTCACGCAGAATGTCGGCCACCTGTTGTTGGTGCGCCGGACGCAGCATGAACGACGcatc            <  1:1524610/71‑1 (MQ=255)
                                                          cTGGCCTCCGGGTCACGCAGAATGTCGGCCACCTGTTGTTGGTGCGCCGGACGCAGCATGAACGACGcatc            <  1:2308537/71‑1 (MQ=255)
                                                            ggCCTCCGGGTCACGCAGAATGTCGGCCACCTGTTGTTGGTGCGCCGGACGCAGCATGAACGACGcatcat          >  1:1629260/1‑71 (MQ=255)
                                                             gCCTCCGGGTCACGCAGAATGTCGGCCACCTGTTGTTGGTGCGCCGGACGCAGCATGAACGACGCATCATg         <  1:245173/71‑1 (MQ=255)
                                                               cTCCGGGTCACGCAGAATGTCGGCCACCTGTTGTTGGTGCGCCGGACGCAGCATGAacgacg                <  1:638584/62‑1 (MQ=255)
                                                               cTCCGGGTCACGCAGAATGTCGGCCACCTGTTGTTGGTGCGCCGGACGCAGCATGAACGACGCATCATGa        <  1:1565932/70‑1 (MQ=255)
                                                                 ccGGGTCACGCAGAATGTCGGCCACCTGTTGTTGGTGCGCCGGACGCAGCATGAACGACGcatcat          >  1:1060536/1‑66 (MQ=255)
                                                                  cGGGTCACGCAGAATGTCGGCCACCTGTTGTTGGTGCGCCGGACGCAGCATGAACGACGc               >  1:3310553/1‑60 (MQ=255)
                                                                  cGGGTCACGCAGAATGTCGGCCACCTGTTGTTGGTGCGCCGGACGCAGCATGAACGACGc               >  1:2086710/1‑60 (MQ=255)
                                                                  cGGGTCACGCAGAATGTCGGCCACCTGTTGTTGGTGCGCCGGACGCAGCATGAACGACGCATCATGAAAc     <  1:108142/70‑1 (MQ=255)
                                                                    ggTCACGCAGAATGTCGGCCACCTGTTGTTGGTGCGCCGGACGCAGCATGAACGACGCATCATGAAACGcc  >  1:2369084/1‑71 (MQ=255)
                                                                     |                                                                     
TTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAGAATGTCGGCCACCTGTTGTTGGTGCGCCGGACGCAGCATGAACGACGCATCATGAAACGCC  >  minE/1015772‑1015910

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: