Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,014,775 C→G 34.1% A432P (GCC→CCC) ‡ fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,014,7750CG34.1% 40.9 / 33.0 41A432P (GCC→CCC) ‡fumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base C (13/14);  new base G (11/3);  total (24/17)
Fisher's exact test for biased strand distribution p-value = 9.58e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAG  >  Exported/1014714‑1014841
                                                             |                                                                  
tGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTc                                                           >  1:1543561/1‑71 (MQ=255)
tGCAGCAGATCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAAccgg                                                            >  1:3406880/1‑68 (MQ=255)
 gCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTcc                                                          >  1:2276558/1‑71 (MQ=255)
      agATCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAAccg                                                             >  1:1166034/2‑62 (MQ=255)
      agATCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAAcc                                                              >  1:4609802/2‑62 (MQ=255)
       gATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGt                                                     <  1:2935683/70‑1 (MQ=255)
         tCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGttt                                                   >  1:116620/1‑70 (MQ=255)
            aCATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGc                                               <  1:2207995/71‑1 (MQ=255)
              aTAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCg                                             >  1:1045288/1‑71 (MQ=255)
                aaGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGa                                           >  1:3023545/1‑71 (MQ=255)
                  gAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCAt                                                                  >  1:429418/1‑46 (MQ=255)
                    atccatccGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGAccc                                        >  1:2789702/1‑70 (MQ=255)
                         tccgtccgGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGttt                                                   <  1:3929613/54‑1 (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:549856/70‑1 (MQ=255)
                                           ggCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCttt                 <  1:3307703/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)
                                                             |                                                                  
TGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAG  >  Exported/1014714‑1014841

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: