Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,015,049 A→G 6.0% H340H (CAT→CAC fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,015,0490AG 6.0% 192.8 / 3.6 67H340H (CAT→CACfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base A (32/31);  new base G (2/2);  total (34/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.95e-01

TAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGC  >  minE/1014979‑1015116
                                                                      |                                                                   
tAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTa                                                                     >  1:899932/1‑71 (MQ=255)
tAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTa                                                                     >  1:2330732/1‑71 (MQ=255)
 aaGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTa                                                                     <  1:141893/70‑1 (MQ=255)
      cAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTc                                                                >  1:2154091/1‑70 (MQ=255)
         cgcgCACCGCTTCGCCTTCTCCAGCTTTTAGCCGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGt                                                            >  1:68083/1‑71 (MQ=255)
           ttcACCGCTTCGCCTTCACCGGCCTGGCGCAGTTCTTGTGGAATGTACTGGCCTGGGTTGTGTTCCAGttt                                                          >  1:1674801/3‑71 (MQ=255)
            gcACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGtttt                                                         <  1:2405387/71‑1 (MQ=255)
            gcACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGtttt                                                         <  1:1440185/71‑1 (MQ=255)
            gcACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGtttt                                                         <  1:2237799/71‑1 (MQ=255)
            gcACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGtttt                                                         <  1:1504462/71‑1 (MQ=255)
            gcACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGtttt                                                         >  1:165342/1‑71 (MQ=255)
              aCCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTc                                                       >  1:99078/1‑71 (MQ=255)
              aCCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTc                                                       <  1:869609/71‑1 (MQ=255)
                cGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCg                                                      >  1:423275/1‑70 (MQ=255)
                cGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGa                                                     >  1:2705091/1‑71 (MQ=255)
                cGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGa                                                     >  1:3842882/1‑71 (MQ=255)
                  cTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATc                                                   <  1:3447399/71‑1 (MQ=255)
                  cTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATc                                                   <  1:3260101/71‑1 (MQ=255)
                     cGCCTTCACCGGCCTGGCGCAGTTCTTGTGGAATGTACTGGCCTGGGTTGTGTTCCAGTTTTTCGATCCa                                                 <  1:326395/70‑1 (MQ=255)
                      gCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTa                                                                     <  1:1672012/49‑1 (MQ=255)
                         ttCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATc                                                   <  1:3598271/64‑1 (MQ=255)
                         ttCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATcc                                            >  1:274429/1‑71 (MQ=255)
                          tctcCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATccc                                           >  1:417230/1‑71 (MQ=255)
                           ctcCCGCTTTGCGGAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGttttt                                                        >  1:3477920/1‑57 (MQ=255)
                           ctcCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGttttt                                                        >  1:727207/1‑57 (MQ=255)
                           ctcCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGttttt                                                        >  1:3289995/1‑57 (MQ=255)
                           ctcCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATcccc                                          >  1:2102682/1‑71 (MQ=255)
                              ccGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGAt                                                    >  1:1201972/1‑58 (MQ=255)
                              ccGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATc                                                   <  1:1437471/59‑1 (MQ=255)
                              ccGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATc                                                   <  1:1035314/59‑1 (MQ=255)
                              ccGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTg                                        <  1:1254966/70‑1 (MQ=255)
                               cGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGAc                                      >  1:1971624/1‑71 (MQ=255)
                               cGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGAc                                      >  1:1914875/1‑71 (MQ=255)
                                 cTTTGCGCAGCTCTTCCGGGATATATTTGCCTTGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACgg                                    <  1:4101072/71‑1 (MQ=255)
                                 cTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGtttt                                                         <  1:3314699/50‑1 (MQ=255)
                                 cTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACgg                                    <  1:3294437/71‑1 (MQ=255)
                                 cTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACgg                                    <  1:2770773/71‑1 (MQ=255)
                                 cTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACgg                                    <  1:1286753/71‑1 (MQ=255)
                                  tttGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACgg                                    <  1:2962426/70‑1 (MQ=255)
                                   ttGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGttt                                                          >  1:2446283/1‑47 (MQ=255)
                                   ttGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGtt                                  >  1:763344/1‑71 (MQ=255)
                                   tggcgcAGTTCTTGTGGAATGTACTTGCCTGGGTTGTGTTCCAGtttt                                                         <  1:2002911/46‑1 (MQ=255)
                                         aGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCtt                            >  1:3143774/1‑71 (MQ=255)
                                          gCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCg                                                      <  1:570136/44‑1 (MQ=255)
                                             ctTCCGGGATATATTTGCCTGGATTATGTTCCAGttttt                                                        >  1:3933054/1‑39 (MQ=255)
                                             ctTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCt                        >  1:3407522/1‑71 (MQ=255)
                                                ccGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCttt                      <  1:550704/70‑1 (MQ=255)
                                                ccGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCttt                      >  1:2997828/1‑70 (MQ=255)
                                                    gATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGatat                 <  1:1667321/71‑1 (MQ=255)
                                                     atatatTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCttt                      >  1:1598816/1‑65 (MQ=255)
                                                      tatatTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTa               <  1:1176771/71‑1 (MQ=255)
                                                      tatatTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTa               <  1:2056677/71‑1 (MQ=255)
                                                      tatatTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTa               <  1:1990288/71‑1 (MQ=255)
                                                       atatTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTa               <  1:452900/70‑1 (MQ=255)
                                                        tatTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGAt                                              >  1:3636079/1‑38 (MQ=255)
                                                        tatTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACgg                                    <  1:3788921/48‑1 (MQ=255)
                                                        tatTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGAt                               >  1:21074/1‑53 (MQ=255)
                                                        tatTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGAt                               >  1:2017608/1‑53 (MQ=255)
                                                        tatTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGc                         <  1:1472120/59‑1 (MQ=255)
                                                         atTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACgg            >  1:1734537/1‑71 (MQ=255)
                                                          tttGCCTGGATTATGTTCCAGTTTTTCGATCCAGAt                                              >  1:1158612/1‑36 (MQ=255)
                                                          tttGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGtt                                  >  1:418821/1‑48 (MQ=255)
                                                          tttGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCttt                      <  1:1879492/60‑1 (MQ=255)
                                                          tttGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGata                  >  1:2838542/1‑64 (MQ=255)
                                                          tttGCCTGGATTATGTTCCAGTTTTTCGATCCAGAACCCCTGACGGTTGATCTTCGCTTTGATATTACgg            >  1:310023/1‑70 (MQ=255)
                                                            ggCCTGGGTTGTGTTCCAGTTTTTCGATCCAGATACCTTCGCGGTTGa                                >  1:3579025/2‑48 (MQ=19)
                                                              ccTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTcag        <  1:1662314/70‑1 (MQ=255)
                                                                tGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTg                                 <  1:739414/43‑1 (MQ=255)
                                                                  gATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCagcagag   <  1:1117211/71‑1 (MQ=255)
                                                                   aTTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCagcagagc  <  1:1152466/71‑1 (MQ=255)
                                                                     tATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCagcagag   >  1:1901724/1‑68 (MQ=255)
                                                                      |                                                                   
TAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGC  >  minE/1014979‑1015116

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: