Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,015,766 A→C 16.4% V101V (GTT→GTG fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,015,7660AC16.4% 102.1 / 19.9 55V101V (GTT→GTGfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base A (26/20);  new base C (5/4);  total (31/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

CGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCC  >  minE/1015699‑1015835
                                                                   |                                                                     
cGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACg                                                                    >  1:290796/1‑71 (MQ=255)
 ggTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGc                                                                   >  1:3045087/1‑71 (MQ=255)
   tCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCt                                                                 <  1:33834/71‑1 (MQ=255)
   tCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCt                                                                 >  1:1415530/1‑71 (MQ=255)
   tCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCt                                                                 >  1:1480523/1‑71 (MQ=255)
     cATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGcc                                                                  <  1:949855/68‑1 (MQ=255)
      aTACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTc                                                              <  1:1151863/71‑1 (MQ=255)
      aTACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTc                                                              >  1:3650727/1‑71 (MQ=255)
       tACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCtt                                                                <  1:2996289/68‑1 (MQ=255)
        acgcGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTggc                                                            >  1:1523251/4‑71 (MQ=255)
        acacGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTcgc                                                            >  1:599650/1‑71 (MQ=255)
               ggccTTTTTTACCGACGATAATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCAc                                                                     >  1:2483132/3‑55 (MQ=255)
               gCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTCGCAGAACGcc                                                                  >  1:2014258/1‑58 (MQ=255)
                cccTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATg                                                     <  1:2857399/70‑1 (MQ=255)
                  cTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATg                                                     <  1:2847556/68‑1 (MQ=255)
                  cTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCg                                                  >  1:1683019/1‑71 (MQ=255)
                   ttttttACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTcg                                                             >  1:415978/1‑59 (MQ=255)
                    tttttACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCtt                                                                <  1:1053496/55‑1 (MQ=255)
                    tttttACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTcgccgc                                                         >  1:3572142/1‑62 (MQ=255)
                     ttttACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAg                                               >  1:1662617/1‑71 (MQ=255)
                     ttttACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAg                                               >  1:3476318/1‑71 (MQ=255)
                      tttACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGt                                              >  1:3416392/1‑71 (MQ=255)
                       ttACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGt                                              <  1:1854993/70‑1 (MQ=255)
                        tACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGtt                                             <  1:2008371/70‑1 (MQ=255)
                        tACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCgg                                                 <  1:1229421/66‑1 (MQ=255)
                                 taatCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAAt                                   >  1:3900731/1‑71 (MQ=255)
                                   atCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTg                                 >  1:3847882/1‑71 (MQ=255)
                                    tCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTa                                            >  1:2662165/1‑59 (MQ=255)
                                     cgcgGTGCCGGTATCCTGACAGGTTGGCAGAACGCCt                                                                 <  1:2241640/37‑1 (MQ=255)
                                     cgcgGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAg                                               <  1:1979499/55‑1 (MQ=255)
                                      gcgGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTg                                 >  1:79115/1‑68 (MQ=255)
                                      gcgGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCa                               <  1:361761/70‑1 (MQ=255)
                                        ggTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAgcg                            <  1:800964/71‑1 (MQ=255)
                                            ccGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCca                         >  1:294181/1‑70 (MQ=255)
                                            ccGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCca                         >  1:2654439/1‑70 (MQ=255)
                                            ccGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCca                         >  1:3525072/1‑70 (MQ=255)
                                                tATCCTTACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTg                                 <  1:2764797/58‑1 (MQ=255)
                                                 aTCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACatat                    >  1:1576625/1‑70 (MQ=255)
                                                       gcAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTcgt             <  1:1186741/70‑3 (MQ=255)
                                                         aGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGaa                                    >  1:923204/1‑46 (MQ=255)
                                                         aGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCAtt            >  1:1658305/1‑70 (MQ=255)
                                                            ttGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCg        <  1:476506/71‑1 (MQ=255)
                                                            ttGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCg        <  1:2832034/71‑1 (MQ=255)
                                                            ttGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCg        <  1:3139360/71‑1 (MQ=255)
                                                            ttGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCg        <  1:2836493/71‑1 (MQ=255)
                                                            ttGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCg        <  1:2719413/71‑1 (MQ=255)
                                                            ttGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCg        <  1:2065601/71‑1 (MQ=255)
                                                            tcggCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCg        <  1:2486585/69‑1 (MQ=255)
                                                            tcggCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCg        >  1:1036365/3‑71 (MQ=255)
                                                             tGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGaa                                    >  1:3794479/1‑42 (MQ=255)
                                                              ggCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTcgt             >  1:2607250/1‑62 (MQ=255)
                                                              ggCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGc       <  1:1344389/70‑1 (MQ=255)
                                                              ggCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCt      >  1:3631046/1‑71 (MQ=255)
                                                               gCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTg     >  1:529851/1‑71 (MQ=255)
                                                                  gAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGcc  >  1:100230/1‑71 (MQ=255)
                                                                   aaCGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCAttt           <  1:1878438/61‑1 (MQ=255)
                                                                   |                                                                     
CGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCC  >  minE/1015699‑1015835

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: