Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,015,088 C→T 14.7% K327K (AAG→AAA fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,015,0880CT14.7% 197.4 / 22.2 75K327K (AAG→AAAfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base C (29/35);  new base T (6/5);  total (35/40)
Fisher's exact test for biased strand distribution p-value = 7.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.90e-01

AGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTG  >  minE/1015020‑1015158
                                                                    |                                                                      
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:1176771/71‑1 (MQ=255)
             tatatTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTa                                                         <  1:2056677/71‑1 (MQ=255)
              atatTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTa                                                         <  1:452900/70‑1 (MQ=255)
               tatTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGc                                                                   <  1:1472120/59‑1 (MQ=255)
                atTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACgg                                                      >  1:1734537/1‑71 (MQ=255)
                 tttGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCttt                                                                <  1:1879492/60‑1 (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: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:2168970/40‑1 (MQ=255)
                                     gTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCttt                                                                <  1:3742768/40‑1 (MQ=255)
                                     gTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTcagcag                                               <  1:4124993/57‑1 (MQ=255)
                                     gTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCagcagag                                             <  1:3257406/59‑1 (MQ=255)
                                     gTTTTTCGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACGccc                                 <  1:2442515/71‑1 (MQ=255)
                                     gTTTTTCGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACGccc                                 <  1:1589916/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)
                                                          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:1126405/1‑70 (MQ=255)
                                                               ttGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACGCCCATGCCGACCGGGCAGGATGCGCCGt        >  1:163156/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:2700594/1‑71 (MQ=255)
                                                                tGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGa      >  1:2511486/1‑71 (MQ=255)
                                                                tGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGa      >  1:2695710/1‑71 (MQ=255)
                                                                   tCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGt   >  1:1859979/1‑71 (MQ=255)
                                                                    cTTCGCTTTTATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGa      <  1:173596/67‑1 (MQ=255)
                                                                    cTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTg  >  1:2658490/1‑71 (MQ=255)
                                                                    cTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTg  >  1:2779622/1‑71 (MQ=255)
                                                                    |                                                                      
AGCTCTTCCGGGATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTG  >  minE/1015020‑1015158

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: