Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,015,073 C→A 9.8% G332G (GGG→GGT fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,015,0730CA 9.8% 228.7 / 7.9 71G332G (GGG→GGTfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base C (28/36);  new base A (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.73e-01

TTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAG  >  minE/1015004‑1015142
                                                                     |                                                                     
ttCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATcc                                                                      >  1:274429/1‑71 (MQ=255)
 tctcCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATccc                                                                     >  1:417230/1‑71 (MQ=255)
  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: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: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)
                               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: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:2242600/40‑1 (MQ=255)
                                                     gTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCttt                                                <  1:3742768/40‑1 (MQ=255)
                                                     gTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCttt                                                <  1:2168970/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)
                                                                     |                                                                     
TTCTCCCGCTTTGCGCAGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAG  >  minE/1015004‑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: