Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,015,097 G→A 16.2% I324I (ATC→ATT fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,015,0970GA16.2% 198.3 / 26.3 81I324I (ATC→ATTfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base G (30/37);  new base A (6/7);  total (36/45)
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.36e-01

GATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGA  >  minE/1015031‑1015166
                                                                  |                                                                     
gATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGatat                                                                   <  1:1667321/71‑1 (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)
     atTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACgg                                                              >  1:1734537/1‑71 (MQ=255)
      tttGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGata                                                                    >  1:2838542/1‑64 (MQ=255)
      tttGCCTGGATTATGTTCCAGTTTTTCGATCCAGAACCCCTGACGGTTGATCTTCGCTTTGATATTACgg                                                              >  1:310023/1‑70 (MQ=255)
          ccTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTcag                                                          <  1:1662314/70‑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)
                       ccAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAgg                                                 >  1:2332628/1‑66 (MQ=255)
                        cAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGa                                                                      >  1:3175627/1‑44 (MQ=255)
                         aGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGcc                                          <  1:3148770/71‑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)
                             tttCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGat                                                                     >  1:4131593/1‑40 (MQ=255)
                             tttCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTa                                                                 >  1:3539007/1‑44 (MQ=255)
                             tttCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTa                                                                 >  1:976679/1‑44 (MQ=255)
                             tttCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTcag                                                          <  1:234736/51‑1 (MQ=255)
                             tttCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTcag                                                          <  1:3611347/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)
                                               cgcGGTTGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACGCCCATGCCGACCGGGCAGGATgcg                    >  1:2523653/3‑71 (MQ=255)
                                                aCGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATgcgc                   >  1:2693676/1‑71 (MQ=255)
                                                aCGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATgc                     <  1:864029/69‑1 (MQ=255)
                                                aCGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATgc                     <  1:366663/69‑1 (MQ=255)
                                                   gTTGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACGCCCATGCCGACCGGGCAGGATGCGCCGt                >  1:767302/1‑71 (MQ=255)
                                                   gTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGga                                                <  1:4009767/39‑1 (MQ=255)
                                                   gTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGt                <  1:3655157/71‑1 (MQ=255)
                                                   gTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGt                <  1:2306050/71‑1 (MQ=255)
                                                    ttGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACGCCCATGCCGACCGGGCAGGATGCGCCGt                >  1:163156/1‑70 (MQ=255)
                                                    ttGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACGCCCATGCCGACCGGGCAGGATGCGCCGt                >  1:1126405/1‑70 (MQ=255)
                                                    ttGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACGCCCATGCCGACCGGGCAGGATGCGCCGt                >  1:2901124/1‑70 (MQ=255)
                                                    ttGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGt                >  1:1146590/1‑70 (MQ=255)
                                                     tGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGcc                                          <  1:2929735/43‑1 (MQ=255)
                                                     tGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTg               >  1:626491/1‑70 (MQ=255)
                                                     tGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGa              >  1:2695710/1‑71 (MQ=255)
                                                     tGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGa              >  1:2511486/1‑71 (MQ=255)
                                                     tGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGa              >  1:2700594/1‑71 (MQ=255)
                                                        ttttCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACGCCCATGCCGACCggg                             <  1:3168649/51‑1 (MQ=255)
                                                        ttttCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACGCCCATGCCGACCggg                             <  1:3164221/51‑1 (MQ=255)
                                                        tCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGt           >  1:1859979/1‑71 (MQ=255)
                                                         tttCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACGCCCATGCCGACCggg                             <  1:2070230/51‑1 (MQ=255)
                                                         cTTCGCTTTTATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGa              <  1:173596/67‑1 (MQ=255)
                                                         cTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTg          >  1:2658490/1‑71 (MQ=255)
                                                         cTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTg          >  1:2779622/1‑71 (MQ=255)
                                                          ttCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATgc                     <  1:2952554/59‑1 (MQ=255)
                                                          ttCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTgg         >  1:2984807/1‑71 (MQ=255)
                                                          ttCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTgg         >  1:4109445/1‑71 (MQ=255)
                                                            cGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCa                                        >  1:1369787/1‑38 (MQ=255)
                                                            cGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACg            <  1:529128/66‑1 (MQ=255)
                                                            cGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACg            <  1:3364572/66‑1 (MQ=255)
                                                            cGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTg          <  1:771055/68‑1 (MQ=255)
                                                            cGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTg          <  1:2939146/68‑1 (MQ=255)
                                                            cGCTTTGATATTACGGGCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTg          <  1:2031427/68‑1 (MQ=255)
                                                             gCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCa       <  1:209572/70‑1 (MQ=255)
                                                                 tGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCg   >  1:835920/1‑70 (MQ=255)
                                                                 tGATATTACGGGCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCg   <  1:1151465/70‑1 (MQ=255)
                                                                  gATATTACGGTCAGCAGAGCAGGAGACGCCCATaccgac                                 >  1:1671867/1‑39 (MQ=255)
                                                                  gATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGa  <  1:4014703/70‑1 (MQ=255)
                                                                  |                                                                     
GATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGA  >  minE/1015031‑1015166

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: