Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,015,418 G→A 9.5% R217R (CGC→CGT fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,015,4180GA 9.5% 173.0 / 6.8 63R217R (CGC→CGTfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base G (22/35);  new base A (3/3);  total (25/38)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01

TTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGA  >  Exported/1015351‑1015487
                                                                   |                                                                     
tttCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTgcgc                                                                    <  1:727259/71‑1 (MQ=255)
      cAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCAtc                                                                 <  1:5247040/68‑1 (MQ=255)
        gAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCAtcttct                                                             >  1:3249245/1‑70 (MQ=255)
         aaGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTtctc                                                            >  1:4577754/1‑70 (MQ=255)
           gTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAAc                                                         <  1:2277619/71‑1 (MQ=255)
                accaATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGt                                                    <  1:814153/71‑1 (MQ=255)
                accaATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGt                                                    <  1:1023314/71‑1 (MQ=255)
                     tAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAAttt                                               <  1:3473114/71‑1 (MQ=255)
                       acgaacgCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAAttttt                                             <  1:2967324/71‑1 (MQ=255)
                       acgaacgCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAAttttt                                             <  1:4536778/71‑1 (MQ=255)
                            cgCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCa                                           >  1:794365/1‑68 (MQ=255)
                            cgCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGtt                                        <  1:5365493/71‑1 (MQ=255)
                             gCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAAttt                                               <  1:3748241/63‑1 (MQ=255)
                             gCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAAttt                                               <  1:1857489/63‑1 (MQ=255)
                               gatatGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGtttgc                                     <  1:4848153/70‑3 (MQ=255)
                                 tatGGTACGGCGGGCAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTgccg                                   >  1:631194/1‑70 (MQ=255)
                                 tatGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAAc                                                         <  1:2984148/49‑1 (MQ=255)
                                 tatGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTcccc                                   <  1:4375839/71‑1 (MQ=255)
                                 tatGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTcccc                                   >  1:5045484/1‑71 (MQ=255)
                                  atGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCaa                                                          >  1:551921/1‑47 (MQ=255)
                                  atGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAAc                                                         >  1:1898184/1‑48 (MQ=255)
                                   tGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGttt                                       >  1:1269513/1‑65 (MQ=255)
                                      tACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCa                                           <  1:4972865/58‑1 (MQ=255)
                                      tACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGt                              <  1:142070/71‑1 (MQ=255)
                                      tACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGt                              <  1:4563966/71‑1 (MQ=255)
                                      tACGGAGGACAGGCCGCCGTACCCAGCGTGAGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGt                              <  1:2622271/71‑1 (MQ=255)
                                        cGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTa                                                   >  1:3833740/1‑48 (MQ=255)
                                        cGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGt                              >  1:4328100/1‑69 (MQ=255)
                                        cGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGt                              >  1:1383487/1‑69 (MQ=255)
                                          gaggaCAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAg                                                      >  1:3164576/1‑43 (MQ=255)
                                          gaggaCAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAg                           >  1:2282417/1‑70 (MQ=255)
                                            ggaCAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCa                                           <  1:5363181/52‑1 (MQ=255)
                                            ggaCAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTaa                        >  1:666597/1‑71 (MQ=255)
                                             ggcAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAg                           <  1:4801432/65‑1 (MQ=255)
                                             gaCAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAAttttt                                             <  1:2588765/49‑1 (MQ=255)
                                             gaCAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAAttttt                                             <  1:1109754/49‑1 (MQ=255)
                                             gaCAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAAttttt                                             <  1:799995/49‑1 (MQ=255)
                                             gaCAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAAttttt                                             <  1:4864288/49‑1 (MQ=255)
                                              gcAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGtttg                                      <  1:1713811/54‑2 (MQ=255)
                                               cAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGc                     >  1:4211564/1‑71 (MQ=255)
                                                  gccgccGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTc                             <  1:867798/60‑1 (MQ=255)
                                                    cgccgTACCCAGCGTGCGCATCTTCTCAACCAGGTAAttttt                                             <  1:3312633/42‑1 (MQ=255)
                                                    cgccgTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCa                                           <  1:1375564/44‑1 (MQ=255)
                                                    cgccgTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCa                                           <  1:5326349/44‑1 (MQ=255)
                                                      ccgTACCCAGCGTGCGCATCTTCTCAACCAGGTaa                                                  <  1:1319913/35‑1 (MQ=255)
                                                      ccgTACCCAGCGTGCGCATCTTCTCAACCAGGTAAttttt                                             <  1:2836546/40‑1 (MQ=255)
                                                      ccgTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTc                                            <  1:3432428/41‑1 (MQ=255)
                                                      ccgTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGt                                         >  1:2798167/1‑44 (MQ=255)
                                                      ccgTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGt                                         >  1:2249303/1‑44 (MQ=255)
                                                      ccgTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGtt              <  1:3999617/71‑1 (MQ=255)
                                                      ccgTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGtt              >  1:5261989/1‑71 (MQ=255)
                                                        gTACCGAGGGTACACATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGATTTGGTTTc            >  1:5678796/1‑71 (MQ=255)
                                                        gTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTc            >  1:3458471/1‑71 (MQ=255)
                                                        gTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTc            >  1:3451956/1‑71 (MQ=255)
                                                          aCCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAgca                         >  1:4605434/1‑54 (MQ=255)
                                                             cAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGa        >  1:2520491/1‑70 (MQ=255)
                                                             cAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGAt       <  1:2068252/71‑1 (MQ=255)
                                                             cAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGAt       <  1:2061669/71‑1 (MQ=255)
                                                              aGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGa        >  1:1410421/1‑69 (MQ=255)
                                                                 gTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATagag   <  1:3478245/71‑1 (MQ=255)
                                                                 gTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATa      >  1:1316827/1‑68 (MQ=255)
                                                                  tGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTgg                <  1:2885404/57‑1 (MQ=255)
                                                                  tGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAgaga  <  1:3181154/71‑1 (MQ=255)
                                                                  tGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAgaga  <  1:1084939/71‑1 (MQ=255)
                                                                   |                                                                     
TTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGA  >  Exported/1015351‑1015487

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: