Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,015,078 G→C 9.8% Q331E (CAG→GAG) ‡ fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,015,0780GC 9.8% 206.8 / 10.3 71Q331E (CAG→GAG) ‡fumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base G (28/36);  new base C (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 = 9.52e-01

CCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAG  >  minE/1015009‑1015142
                                                                     |                                                                
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)
   cTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACgg                                                              <  1:2770773/71‑1 (MQ=255)
   cTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACgg                                                              <  1:1286753/71‑1 (MQ=255)
   cTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACgg                                                              <  1:3294437/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:1176771/71‑1 (MQ=255)
                        tatatTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTa                                         <  1:1990288/71‑1 (MQ=255)
                        tatatTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTa                                         <  1:2056677/71‑1 (MQ=255)
                         atatTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTa                                         <  1:452900/70‑1 (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)
                            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:398089/71‑1 (MQ=255)
                                         tGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGga                        <  1:4049330/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:1551770/53‑1 (MQ=255)
                                                 tttttCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTcag                                  <  1:3682944/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:3798517/1‑41 (MQ=255)
                                                               tCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTcagca                                >  1:1003515/1‑41 (MQ=255)
                                                               tCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAg  <  1:2251038/71‑1 (MQ=255)
                                                                    tGACGGTTGATCTTCGCTTTGATATTACGGTCagcagagcag                          <  1:2459745/42‑1 (MQ=255)
                                                                     |                                                                
CCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAG  >  minE/1015009‑1015142

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: