Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,015,079 A→G 9.8% R330R (CGT→CGC fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,015,0790AG 9.8% 194.4 / 12.0 71R330R (CGT→CGCfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base A (28/36);  new base G (2/5);  total (30/41)
Fisher's exact test for biased strand distribution p-value = 6.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.35e-01

CGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGC  >  minE/1015010‑1015149
                                                                     |                                                                      
cGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGAc                                                                       >  1:1914875/1‑71 (MQ=255)
cGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGAc                                                                       >  1:1971624/1‑71 (MQ=255)
  cTTTGCGCAGCTCTTCCGGGATATATTTGCCTTGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACgg                                                                     <  1:4101072/71‑1 (MQ=255)
  cTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACgg                                                                     <  1:3294437/71‑1 (MQ=255)
  cTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACgg                                                                     <  1:1286753/71‑1 (MQ=255)
  cTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACgg                                                                     <  1:2770773/71‑1 (MQ=255)
   tttGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACgg                                                                     <  1:2962426/70‑1 (MQ=255)
    ttGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGtt                                                                   >  1:763344/1‑71 (MQ=255)
          aGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCtt                                                             >  1:3143774/1‑71 (MQ=255)
              ctTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCt                                                         >  1:3407522/1‑71 (MQ=255)
                 ccGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCttt                                                       >  1:2997828/1‑70 (MQ=255)
                 ccGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCttt                                                       <  1:550704/70‑1 (MQ=255)
                     gATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGatat                                                  <  1:1667321/71‑1 (MQ=255)
                      atatatTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCttt                                                       >  1:1598816/1‑65 (MQ=255)
                       tatatTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTa                                                <  1:2056677/71‑1 (MQ=255)
                       tatatTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTa                                                <  1:1990288/71‑1 (MQ=255)
                       tatatTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTa                                                <  1:1176771/71‑1 (MQ=255)
                        atatTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTa                                                <  1:452900/70‑1 (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)
                           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)
                                       gtGTTCCAGTTTTTCGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTCagcggagcag                                 >  1:3940202/2‑70 (MQ=255)
                                        tGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGga                               <  1:4049330/71‑1 (MQ=255)
                                        tGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGga                               <  1:398089/71‑1 (MQ=255)
                                        tGTTCCAGTTTTTCGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGga                               <  1:2483611/71‑1 (MQ=255)
                                         gTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGga                               <  1:1760336/70‑1 (MQ=255)
                                          ttCCAGTTTTTCGATCCAGATACCTTCGCGGTTGATTTTCGCttt                                                       <  1:1576478/45‑1 (MQ=25)
                                            ccAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAgg                                >  1:2332628/1‑66 (MQ=255)
                                             cAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGa                                                     >  1:3175627/1‑44 (MQ=255)
                                              aGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGcc                         <  1:3148770/71‑1 (MQ=255)
                                               gTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCttt                                                       <  1:2168970/40‑1 (MQ=255)
                                               gTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCttt                                                       <  1:3742768/40‑1 (MQ=255)
                                               gTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCttt                                                       <  1:2242600/40‑1 (MQ=255)
                                               gTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTcagcag                                      <  1:4124993/57‑1 (MQ=255)
                                               gTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCagcagag                                    <  1:3257406/59‑1 (MQ=255)
                                               gTTTTTCGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACGccc                        <  1:1589916/71‑1 (MQ=255)
                                               gTTTTTCGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACGccc                        <  1:2442515/71‑1 (MQ=255)
                                                tttttCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGata                                                   >  1:300331/1‑43 (MQ=255)
                                                tttttCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTcag                                         <  1:3682944/53‑1 (MQ=255)
                                                tttttCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTcag                                         <  1:1551770/53‑1 (MQ=255)
                                                tttttCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGAc                            <  1:90109/66‑1 (MQ=255)
                                                  tttCGATCCAGATCCCCTGACGGTTGATCTTCGCttt                                                       >  1:2188643/1‑37 (MQ=255)
                                                  tttCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGat                                                    >  1:4131593/1‑40 (MQ=255)
                                                  tttCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTa                                                >  1:976679/1‑44 (MQ=255)
                                                  tttCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTa                                                >  1:3539007/1‑44 (MQ=255)
                                                  tttCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTcag                                         <  1:3611347/51‑1 (MQ=255)
                                                  tttCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTcag                                         <  1:234736/51‑1 (MQ=255)
                                                  tttCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCAt                      <  1:2566564/70‑1 (MQ=255)
                                                     cGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTcag                                         >  1:144660/1‑48 (MQ=255)
                                                      gATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTCAGCGGGGCAGGAGACGCCCATgccg                  <  1:676450/70‑1 (MQ=255)
                                                       aTCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTcag                                         <  1:1678418/46‑1 (MQ=255)
                                                       aTCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATaccgac                >  1:2805520/1‑71 (MQ=255)
                                                         ccAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATaccgaccg              >  1:1803164/1‑71 (MQ=255)
                                                            gATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCggg            >  1:21813/1‑70 (MQ=255)
                                                             aTCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTcagcag                                      >  1:1853814/1‑43 (MQ=255)
                                                              tCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTcagca                                       >  1:1003515/1‑41 (MQ=255)
                                                              tCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTcagca                                       >  1:3798517/1‑41 (MQ=255)
                                                              tCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAg         <  1:2251038/71‑1 (MQ=255)
                                                                   tGACGGTTGATCTTCGCTTTGATATTACGGTCagcagagcag                                 <  1:2459745/42‑1 (MQ=255)
                                                                     aCGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATgcgc  >  1:2693676/1‑71 (MQ=255)
                                                                     aCGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATgc    <  1:864029/69‑1 (MQ=255)
                                                                     aCGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATgc    <  1:366663/69‑1 (MQ=255)
                                                                     |                                                                      
CGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGC  >  minE/1015010‑1015149

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: