Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,015,076 C→T 9.8% Q331Q (CAG→CAA) ‡ fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,015,0760CT 9.8% 212.9 / 11.0 71Q331Q (CAG→CAA) ‡fumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base C (28/36);  new base T (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

CTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAG  >  minE/1015006‑1015142
                                                                      |                                                                  
ctcCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATcccc                                                                    >  1:2102682/1‑71 (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)
      cTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACgg                                                              <  1:2770773/71‑1 (MQ=255)
      cTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACgg                                                              <  1:3294437/71‑1 (MQ=255)
      cTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACgg                                                              <  1:1286753/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: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:3798517/1‑41 (MQ=255)
                                                                  tCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTcagca                                >  1:1003515/1‑41 (MQ=255)
                                                                  tCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAg  <  1:2251038/71‑1 (MQ=255)
                                                                      |                                                                  
CTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAG  >  minE/1015006‑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: