Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,014,794 G→C 22.0% A425A (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,7940GC22.0% 110.6 / 28.6 59A425A (GCC→GCGfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base G (31/15);  new base C (8/5);  total (39/20)
Fisher's exact test for biased strand distribution p-value = 7.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

ACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGC  >  Exported/1014726‑1014863
                                                                    |                                                                     
aCATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGc                                                                     <  1:2207995/71‑1 (MQ=255)
  aTAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCg                                                                   >  1:1045288/1‑71 (MQ=255)
    aaGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGa                                                                 >  1:3023545/1‑71 (MQ=255)
        atccatccGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGAccc                                                              >  1:2789702/1‑70 (MQ=255)
               cgtcccGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCgtagta                                                      <  1:3608595/71‑1 (MQ=255)
                   ctgcggTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCgtag                                                        >  1:2669235/3‑65 (MQ=255)
                    tgcggTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCgtag                                                        >  1:3835326/2‑64 (MQ=255)
                    ggcggTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCg                                                 <  1:1963437/71‑1 (MQ=255)
                      cggTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCggg                                               >  1:3979551/1‑70 (MQ=255)
                       ggTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGgat                                              <  1:2927941/71‑1 (MQ=255)
                        gTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGgatg                                             <  1:1281739/71‑1 (MQ=255)
                              gggCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCttt                                       >  1:2622623/1‑71 (MQ=255)
                              gggCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCtt                                        >  1:82283/1‑70 (MQ=255)
                               ggCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCttt                                       <  1:3307703/70‑1 (MQ=255)
                               ggCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCttt                                       <  1:549856/70‑1 (MQ=255)
                                 ccAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTaaa                                                    <  1:830660/55‑1 (MQ=255)
                                 ccAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTaaa                                                    <  1:1415740/55‑1 (MQ=255)
                                  cAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATg                                   <  1:5629680/71‑1 (MQ=255)
                                   aaGTGAACCTGATGGATAACCGGCAGGTGTTTTCGCCGGACCcgc                                                            <  1:1037267/45‑1 (MQ=255)
                                   aaGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGATg                                   >  1:2967144/1‑70 (MQ=255)
                                   aaGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGATg                                   >  1:2722137/1‑70 (MQ=255)
                                   aaGAGAACGGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATc                                          <  1:961661/63‑1 (MQ=255)
                                     gtgaACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCgtag                                                        <  1:1254013/45‑1 (MQ=255)
                                     gtgaACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCgta                                                         >  1:448847/3‑46 (MQ=255)
                                     gagaACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGga                                               >  1:4586304/1‑56 (MQ=255)
                                     gagaACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTAc                                >  1:47173/1‑71 (MQ=255)
                                      agaACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCttt                                       <  1:770335/63‑1 (MQ=255)
                                       gaACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGATGTACTg                              >  1:5638293/1‑71 (MQ=255)
                                        aaCCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTg                                      <  1:4151197/62‑1 (MQ=255)
                                            ggAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGa                                     >  1:590103/1‑59 (MQ=255)
                                            ggAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGa                                     >  1:86848/1‑59 (MQ=255)
                                             gAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAg                        >  1:808093/1‑71 (MQ=255)
                                                  aTAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTaa                                                     >  1:3987559/1‑37 (MQ=255)
                                                    aaCCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGgatga                                            >  1:314085/1‑44 (MQ=255)
                                                    aaCCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCg                            >  1:3455482/1‑60 (MQ=255)
                                                    aaCCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCg                            >  1:5326622/1‑60 (MQ=255)
                                                     aCCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGa                                     >  1:1864890/1‑50 (MQ=255)
                                                     aCCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCa                 >  1:3814459/1‑70 (MQ=255)
                                                     aCCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCa                 >  1:5414016/1‑70 (MQ=255)
                                                     aCCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCAc                >  1:101591/1‑71 (MQ=255)
                                                      ccTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCAcc               >  1:4527079/1‑71 (MQ=255)
                                                         tCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCtt                                        >  1:549322/1‑43 (MQ=255)
                                                         tCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCAcc               >  1:2011655/1‑68 (MQ=255)
                                                         tCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGt             >  1:3026398/1‑70 (MQ=255)
                                                           cGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGccc                     <  1:2412323/60‑1 (MQ=255)
                                                              ggTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATc                                          <  1:2039985/35‑1 (MQ=25)
                                                                ttttGGCCGGACGCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTatccatc     >  1:4626157/1‑71 (MQ=255)
                                                                ttttGGCCGGACCCGCGTAGTAAATCGGATGATCTTTg                                      <  1:235991/38‑1 (MQ=255)
                                                                ttttGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGt             <  1:1013865/63‑1 (MQ=255)
                                                                ttttGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTatccatc     >  1:3414559/1‑71 (MQ=255)
                                                                ttttCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGATGTACTGCg                            <  1:2721111/48‑1 (MQ=37)
                                                                  ttGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGAt                                    >  1:1479560/1‑38 (MQ=255)
                                                                  ttGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGAt                                    >  1:3959589/1‑38 (MQ=255)
                                                                  ttGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTg                              >  1:2857754/1‑44 (MQ=255)
                                                                  ttGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTatcca       >  1:4590974/1‑67 (MQ=255)
                                                                  ttGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCg   >  1:2099960/1‑71 (MQ=255)
                                                                  ttGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCg   >  1:3490139/1‑71 (MQ=255)
                                                                  ttGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCg   >  1:1742759/1‑71 (MQ=255)
                                                                  ttGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCg   >  1:1345685/1‑71 (MQ=255)
                                                                    ggCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGc  <  1:5612886/70‑1 (MQ=255)
                                                                    |                                                                     
ACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGC  >  Exported/1014726‑1014863

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: