Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,014,650 G→C 39.7% A473A (GCC→GCG fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,014,6500GC39.7% 39.0 / 65.6 63A473A (GCC→GCGfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base G (15/23);  new base C (11/14);  total (26/37)
Fisher's exact test for biased strand distribution p-value = 7.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-01

CCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCA  >  Exported/1014585‑1014717
                                                                 |                                                                   
cccTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCa                                                                  <  1:5435483/69‑1 (MQ=255)
cccTGCGCCAATACAGCGGCCGGACCACCGATACAGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCa                                                                  <  1:3040774/69‑1 (MQ=255)
  cTGCGCCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGcgt                                                                 <  1:2813926/68‑3 (MQ=255)
  cTGCGCCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTcgg                                                              >  1:5648506/1‑69 (MQ=255)
  cTGCGCCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTcg                                                               >  1:5141058/1‑69 (MQ=255)
  cTGCGCCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTc                                                                <  1:227697/69‑1 (MQ=255)
  cTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCa                                                                  >  1:5498553/1‑67 (MQ=255)
  cTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCg                                                              <  1:406601/71‑1 (MQ=255)
     cgcCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGtaa                                                           >  1:4676832/1‑69 (MQ=255)
     cgcCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGta                                                            <  1:2331422/70‑2 (MQ=255)
     cgcCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGta                                                            <  1:2271464/70‑2 (MQ=255)
     cgcCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGt                                                             <  1:2188830/69‑1 (MQ=255)
     cgcCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGt                                                             <  1:4717182/69‑1 (MQ=255)
     cgcCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGt                                                             <  1:3096114/69‑1 (MQ=255)
     cgcCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTc                                                            <  1:1776041/70‑1 (MQ=255)
        cAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACct                                                        >  1:5625313/1‑71 (MQ=255)
        cAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACct                                                        >  1:2937934/1‑71 (MQ=255)
         agtACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACct                                                        >  1:3344664/3‑70 (MQ=255)
          gtACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACctgc                                                      >  1:401206/2‑71 (MQ=255)
          aTACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCAt                                                                 <  1:5582778/60‑1 (MQ=255)
           tACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATcc                                                               <  1:2153796/61‑1 (MQ=255)
           tACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACctgct                                                     >  1:477448/1‑71 (MQ=255)
               gccgCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTgctgact                                                 >  1:2804282/4‑71 (MQ=255)
                ccgCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTgctgactg                                                >  1:1923016/3‑71 (MQ=255)
                   ccGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCgg                                             >  1:1332908/1‑71 (MQ=255)
                     ggACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGtt                                           <  1:3406118/71‑1 (MQ=255)
                        ccgccGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTAcc                                        >  1:1428085/4‑71 (MQ=255)
                          accGATACTGCCAAGGTAGAAGCCTCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGtt                                           <  1:1485907/66‑1 (MQ=255)
                           ccGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCttt                                     <  1:2124863/71‑1 (MQ=255)
                           ccGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCa                                                           >  1:5358990/1‑49 (MQ=255)
                              aTACTGCCGAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGc                                  >  1:4718648/1‑71 (MQ=255)
                               tACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCg                                              <  1:584932/58‑1 (MQ=255)
                                   gCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATcc                                                               >  1:1468642/1‑37 (MQ=255)
                                   gCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCAc                                                          >  1:1089860/1‑42 (MQ=255)
                                    ccgaGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCTTTCGCCAGCAt                            >  1:4711726/4‑71 (MQ=255)
                                     cAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGt                                            >  1:3420216/1‑54 (MQ=255)
                                     cAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATg                           <  1:3314451/71‑1 (MQ=255)
                                      aaGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGa                          <  1:4642660/71‑1 (MQ=255)
                                         gTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGc                                  >  1:9102/1‑60 (MQ=255)
                                         gTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCa                       <  1:2046162/71‑1 (MQ=255)
                                          tAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGt                                            <  1:3664711/49‑1 (MQ=255)
                                          tAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGt                                            <  1:836799/49‑1 (MQ=255)
                                          tAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCa                       >  1:2626776/1‑70 (MQ=255)
                                          tAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCAt                      <  1:887462/71‑1 (MQ=255)
                                          tAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCTTTCGCCAGCATGATCAt                      <  1:2955299/71‑1 (MQ=255)
                                           aGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCAt                      >  1:2383907/1‑70 (MQ=255)
                                           aGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCTTTCGCCAGCATGATCAtg                     >  1:875607/1‑70 (MQ=255)
                                            gAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCTTTCGCCAGCATGATCAtg                     <  1:2329419/70‑2 (MQ=255)
                                             aaGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCTTTCGCCAGCATGATCAtgc                    <  1:3357362/70‑3 (MQ=255)
                                             aaGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCTTTCGCCAGCATGATCAtg                     <  1:3328683/69‑2 (MQ=255)
                                             aaGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCTTTCGCCAGCATGATCATGCt                   <  1:530231/71‑1 (MQ=255)
                                                 cgccgTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCg               <  1:4932291/71‑1 (MQ=255)
                                                 cgccgTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCg               <  1:5149872/71‑1 (MQ=255)
                                                 cgccgTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCg               >  1:4858543/1‑71 (MQ=255)
                                                      tgtTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCAt                      <  1:4868250/59‑1 (MQ=255)
                                                      tgtTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCAt                      <  1:4861554/59‑1 (MQ=255)
                                                        ttttttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCAt                      <  1:4304837/57‑1 (MQ=255)
                                                         tttttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTgcgc       >  1:3348219/1‑71 (MQ=255)
                                                          ttttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTc                                    <  1:328727/41‑1 (MQ=255)
                                                             tACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTg    <  1:2412482/70‑1 (MQ=255)
                                                             tACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTg    >  1:118167/1‑70 (MQ=255)
                                                               cAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGc   <  1:3500963/69‑1 (MQ=255)
                                                               cAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCa  <  1:1257073/70‑1 (MQ=255)
                                                                 |                                                                   
CCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCA  >  Exported/1014585‑1014717

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: