Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,015,838 C→T 12.3% P77P (CCG→CCA fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,015,8380CT12.3% 125.9 / 8.4 49P77P (CCG→CCAfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base C (13/29);  new base T (3/3);  total (17/32)
Fisher's exact test for biased strand distribution p-value = 3.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.47e-01

GCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAGAATGTCGGCCACCTGTTGTTGGTGCGCCGGACGCAGCATGAACGACGCATCATGAAAC  >  minE/1015769‑1015907
                                                                     |                                                                     
gCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCg                                                                      <  1:1622153/71‑1 (MQ=255)
 ccTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCgg                                                                     <  1:2916760/71‑1 (MQ=255)
 ccTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCgg                                                                     <  1:3641654/71‑1 (MQ=255)
   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)
                                 ttGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTgga                                                                    <  1:2691527/40‑2 (MQ=12)
                                 ttGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTGGt                                     <  1:1844678/71‑1 (MQ=255)
                                 ttGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTGGt                                     >  1:578541/1‑71 (MQ=255)
                                 ttGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAGAATGTCGGCCACCtgttgt                                         >  1:1546073/1‑67 (MQ=255)
                                   gCAGCGCCACATATTTATCATTTTCGCTGGCCTCCgg                                                                     >  1:1952346/1‑37 (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:3807331/70‑1 (MQ=255)
                                                            gCTGGCCTCCGGGTCACGCAGAATGTCGGCCACCTGTTGTTGGTGCGCCGGACGCAGCATGAACGACGca           <  1:158835/70‑1 (MQ=255)
                                                             cTGGCCTCCGGGTCACGCAGAATGTCGGCCACCTGTTGTTGGTGCGCCGGACGCAGCATGAACGACGcatc         <  1:2308537/71‑1 (MQ=255)
                                                             cTGGCCTCCGGGTCACGCAGAATGTCGGCCACCTGTTGTTGGTGCGCCGGACGCAGCATGAACGACGcatc         <  1:1524610/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:2086710/1‑60 (MQ=255)
                                                                     cGGGTCACGCAGAATGTCGGCCACCTGTTGTTGGTGCGCCGGACGCAGCATGAACGACGc            >  1:3310553/1‑60 (MQ=255)
                                                                     cGGGTCACGCAGAATGTCGGCCACCTGTTGTTGGTGCGCCGGACGCAGCATGAACGACGCATCATGAAAc  <  1:108142/70‑1 (MQ=255)
                                                                     |                                                                     
GCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAGAATGTCGGCCACCTGTTGTTGGTGCGCCGGACGCAGCATGAACGACGCATCATGAAAC  >  minE/1015769‑1015907

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: