Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,015,112 A→G 21.3% S319S (TCT→TCC fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,015,1120AG21.3% 123.0 / 30.7 61S319S (TCT→TCCfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base A (22/26);  new base G (5/8);  total (27/34)
Fisher's exact test for biased strand distribution p-value = 7.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTC  >  minE/1015045‑1015180
                                                                   |                                                                    
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)
           aGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGcc                                                        <  1:3148770/71‑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)
             tttttCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGAc                                                           <  1:90109/66‑1 (MQ=255)
               tttCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCAt                                                     <  1:2566564/70‑1 (MQ=255)
                   gATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTCAGCGGGGCAGGAGACGCCCATgccg                                                 <  1:676450/70‑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:2306050/71‑1 (MQ=255)
                                     gTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGt                              <  1:3655157/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:2700594/1‑71 (MQ=255)
                                       tGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGa                            >  1:2511486/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:2779622/1‑71 (MQ=255)
                                           cTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTg                        >  1:2658490/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:3364572/66‑1 (MQ=255)
                                              cGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACg                          <  1:529128/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)
                                                        ttACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAg                                        <  1:3241916/42‑1 (MQ=255)
                                                        ttACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGa                                      >  1:3993592/1‑44 (MQ=255)
                                                        ttACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCa            <  1:3981594/70‑1 (MQ=255)
                                                          aCGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGc                  >  1:3357711/1‑62 (MQ=255)
                                                          aCGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGc                  >  1:1534337/1‑62 (MQ=255)
                                                          aCGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCAcgc         >  1:911659/1‑71 (MQ=255)
                                                           cGGTCAGCGGAGCAGGAGACGCCCATGCCGACCGGGCAGGATGCGCCGTGACGTGGCAGACGGATAAcgcg        <  1:2736628/71‑1 (MQ=255)
                                                            ggTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCAcgcg        <  1:2175118/70‑1 (MQ=255)
                                                                agcagagcagGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTg                        <  1:965196/50‑1 (MQ=255)
                                                                 gcagagcagGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTc  <  1:2543808/71‑1 (MQ=255)
                                                                   |                                                                    
GATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTC  >  minE/1015045‑1015180

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: