Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,014,659 C→A 28.8% V470V (GTG→GTT fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,014,6590CA28.8% 91.3 / 28.9 52V470V (GTG→GTTfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base C (18/19);  new base A (8/7);  total (26/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGAT  >  Exported/1014590‑1014723
                                                                     |                                                                
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)
      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)
                              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:836799/49‑1 (MQ=255)
                                     tAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGt                                                  <  1:3664711/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:4858543/1‑71 (MQ=255)
                                            cgccgTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCg                     <  1:4932291/71‑1 (MQ=255)
                                            cgccgTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCg                     <  1:5149872/71‑1 (MQ=255)
                                                 tgtTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCAt                            <  1:4861554/59‑1 (MQ=255)
                                                 tgtTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCAt                            <  1:4868250/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:118167/1‑70 (MQ=255)
                                                        tACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTg          <  1:2412482/70‑1 (MQ=255)
                                                          cAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGc         <  1:3500963/69‑1 (MQ=255)
                                                          cAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCa        <  1:1257073/70‑1 (MQ=255)
                                                              cATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGa   >  1:1866316/1‑71 (MQ=255)
                                                               aTCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGAt  >  1:785037/1‑71 (MQ=255)
                                                                  cGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATAc                          >  1:3708461/1‑44 (MQ=255)
                                                                  cGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATAc                          >  1:3604764/1‑44 (MQ=255)
                                                                   gTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCg                     <  1:5604358/48‑1 (MQ=255)
                                                                   gTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCg                     <  1:5680458/48‑1 (MQ=255)
                                                                   gTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCa        >  1:72032/1‑61 (MQ=255)
                                                                     |                                                                
CGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGAT  >  Exported/1014590‑1014723

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: