Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,015,046 A→G 5.5% N341N (AAT→AAC fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,015,0460AG 5.5% 213.0 / 3.0 73N341N (AAT→AACfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base A (33/36);  new base G (2/2);  total (35/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

GTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGC  >  minE/1014977‑1015116
                                                                     |                                                                      
gTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGAt                                                                       <  1:1378592/71‑1 (MQ=255)
 ttAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGAtt                                                                      <  1:1909983/71‑1 (MQ=255)
  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)
         aaCGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGAt                                                                       <  1:1811078/62‑1 (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:869609/71‑1 (MQ=255)
                aCCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTc                                                       >  1:99078/1‑71 (MQ=255)
                  cGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCg                                                      >  1:423275/1‑70 (MQ=255)
                  cGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGa                                                     >  1:3842882/1‑71 (MQ=255)
                  cGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGa                                                     >  1:2705091/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)
                           ttCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGAt                                                                       <  1:4140745/44‑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:1286753/71‑1 (MQ=255)
                                   cTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACgg                                    <  1:3294437/71‑1 (MQ=255)
                                   cTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACgg                                    <  1:2770773/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:1990288/71‑1 (MQ=255)
                                                        tatatTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTa               <  1:2056677/71‑1 (MQ=255)
                                                        tatatTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTa               <  1:1176771/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:2017608/1‑53 (MQ=255)
                                                          tatTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGAt                               >  1:21074/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)
                                                                     |                                                                      
GTTAAGGTCAACGCGCACCGCTTCGCCTTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGC  >  minE/1014977‑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: