Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,015,835 C→T 14.3% E78E (GAG→GAA fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,015,8350CT14.3% 121.5 / 11.4 49E78E (GAG→GAAfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base C (14/28);  new base T (4/3);  total (18/31)
Fisher's exact test for biased strand distribution p-value = 3.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

GAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAGAATGTCGGCCACCTGTTGTTGGTGCGCCGGACGCAGCATGAACGACGCATCATGA  >  minE/1015765‑1015904
                                                                      |                                                                     
gAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGcc                                                                       >  1:100230/1‑71 (MQ=255)
    gCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTctg                                                                   >  1:3846492/1‑69 (MQ=255)
    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:1570174/71‑1 (MQ=255)
               cgATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAGa                                                        <  1:2946017/71‑1 (MQ=255)
                 aTGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAGAAt                                                      <  1:428186/71‑1 (MQ=255)
                  tGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCggg                                                                 >  1:340474/1‑59 (MQ=255)
                   gTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTc                                                                     <  1:2011663/54‑1 (MQ=255)
                       gAGTTACGCATGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAGAATGTCGGc                                                <  1:579382/71‑1 (MQ=255)
                       gAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAGAATGTCg                                                  >  1:112155/1‑69 (MQ=255)
                       gAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAGAATGTCGGc                                                >  1:2103028/1‑71 (MQ=255)
                       gAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAGAATGTCGGc                                                <  1:2280586/71‑1 (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:578541/1‑71 (MQ=255)
                                     ttGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTGGt                                  <  1:1844678/71‑1 (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:796910/1‑70 (MQ=255)
                                                          aTTTTCGCTGGCCTCCGGGTCACGCAGAATGTCGGCCACCTGTTGTTGGTGCGCCGGACGCAGCATGAac              <  1:80402/70‑1 (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)
                                                                      |                                                                     
GAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAGAATGTCGGCCACCTGTTGTTGGTGCGCCGGACGCAGCATGAACGACGCATCATGA  >  minE/1015765‑1015904

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: