Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,014,740 C→A 18.7% R443R (CGG→CGT fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,014,7400CA18.7% 134.5 / 14.6 53R443R (CGG→CGTfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base C (29/14);  new base A (4/6);  total (33/20)
Fisher's exact test for biased strand distribution p-value = 1.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCC  >  Exported/1014671‑1014803
                                                                     |                                                               
gCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAatccatcc                                                                >  1:4892227/1‑71 (MQ=255)
    ttGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCAtccgtcc                                                            >  1:2943055/1‑71 (MQ=255)
     taccTTTCGCCAGCATGATCATGCTGCCGCCGTGGGATTGCAGCAGATCCACGTAGGAGTCCAtacggcct                                                           >  1:2604965/3‑70 (MQ=255)
     tGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCAtccgtccg                                                           >  1:2437893/1‑71 (MQ=255)
      gCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCAtccgtccg                                                           >  1:153873/1‑70 (MQ=255)
      gCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCAtccgt                                                              <  1:4345819/67‑1 (MQ=255)
      accTTTCGCCAGCATGATCATGCTGCCGCCGTGGGATTGCAGCAGATCCACGTAGGAGTCCAtacggcctg                                                          >  1:3270859/2‑69 (MQ=255)
      accTTTCGCCAGCATGATCATGCTGCCGCCGTGGGATTGCAGCAGATCCACGTAGGAGTCCAtacggcctg                                                          >  1:3405713/2‑69 (MQ=255)
      accTTTCGCCAGCATGATCATGCTGCCGCCGTGGGATTGCAGCAGATCCACGTAGGAGTCCAtacggcctg                                                          >  1:5384115/2‑69 (MQ=255)
      accTTTCGCCAGCATGATCATGCCGCCGCCGTGGGATTGCAGCAGATCCACGTAGGAGTCCAtacggcct                                                           >  1:5142767/2‑69 (MQ=255)
        cTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCcggcg                                                        >  1:2273606/1‑71 (MQ=255)
          ttCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGt                                                      >  1:3373890/1‑71 (MQ=255)
           tCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGt                                                      <  1:1482870/70‑1 (MQ=255)
           tCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTc                                                     >  1:4617686/1‑71 (MQ=255)
             gCCAGCATGATCATGCTGCCGCCGTGGGATTGCAGCAGATCCACGTAGGAGTCCATACGGCCTGCggtggt                                                   <  1:208187/71‑1 (MQ=255)
             gCCAGCATGATCATGCTGCCGCCGTGGGATTGCAGCAGATCCACGTAGGAGTCCATACGGCCTGCGGtg                                                     <  1:2335407/69‑2 (MQ=255)
              ccAGCATGATCATGCTGCCGCCGTGGGATTGCAGCAGATCCACGTAGGAGTCCATACGGCCTGCGGTGGtt                                                  <  1:5185377/71‑1 (MQ=255)
              ccAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCAtccg                                                               <  1:5585939/58‑1 (MQ=255)
              ccAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGt                                                   <  1:3198675/70‑1 (MQ=255)
                aGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTgg                                                >  1:1936462/1‑71 (MQ=255)
                aGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTgg                                                <  1:3163240/71‑1 (MQ=255)
                  cATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGt                                                   <  1:5507567/66‑1 (MQ=255)
                  cATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGc                                              >  1:2295277/1‑71 (MQ=255)
                    tgttCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCa                                            <  1:3217070/68‑1 (MQ=255)
                    tGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTg                                                 >  1:1444583/1‑66 (MQ=255)
                    tGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGcc                                             <  1:4603298/70‑1 (MQ=255)
                     gATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGtt                                                  <  1:2871881/64‑1 (MQ=255)
                      aTCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAatccatcc                                                                >  1:3933305/1‑49 (MQ=255)
                      aTCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAatccatcc                                                                >  1:2968097/1‑49 (MQ=255)
                      aTCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGcc                                             >  1:520481/1‑68 (MQ=255)
                       tCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCg                                                    >  1:4616898/1‑60 (MQ=255)
                         aTACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAagag                                        >  1:5288592/1‑70 (MQ=255)
                            ctgccgccgTGGGATTGCAGCAGATCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGGCCAAGTGAAc                                     >  1:1768185/4‑70 (MQ=255)
                              tCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCg                                   >  1:3303521/1‑70 (MQ=255)
                               ccgccgTGGGATTGCAGCAGATCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGGCCAAGTGAAcctg                                  <  1:2496049/70‑3 (MQ=255)
                               ccGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCAGTCCGGCGGTCGTTGGGCCAAGAGAACCg                                   >  1:532879/1‑69 (MQ=255)
                                    cTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAaga                                         >  1:73698/1‑58 (MQ=255)
                                    cTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAaga                                         >  1:3919429/1‑58 (MQ=255)
                                    cTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAaga                                         >  1:5226654/1‑58 (MQ=255)
                                    cTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAaga                                         >  1:2758279/1‑58 (MQ=255)
                                    cTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAaga                                         >  1:3284340/1‑58 (MQ=255)
                                          ttGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGt                                                   <  1:1779207/42‑1 (MQ=255)
                                          ttGCAGCAGATCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGGCCAAGTGAAcctg                                  <  1:2145841/59‑3 (MQ=255)
                                           tGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTc                     >  1:1543561/1‑71 (MQ=255)
                                           tGCAGCAGATCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGGCCAAgtga                                       <  1:5144256/53‑1 (MQ=255)
                                           tGCAGCAGATCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAAccgg                      >  1:3406880/1‑68 (MQ=255)
                                            gCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCg                                                    <  1:3510167/39‑1 (MQ=255)
                                            gCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCg                                                    <  1:2092348/39‑1 (MQ=255)
                                            gCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGt                                                   >  1:1125083/1‑40 (MQ=255)
                                            gCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGt                                                   >  1:61690/1‑40 (MQ=255)
                                            gCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGt                                                   >  1:1098524/1‑40 (MQ=255)
                                            gCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTcc                    >  1:2276558/1‑71 (MQ=255)
                                                ttGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCaa                                           <  1:2636206/44‑1 (MQ=255)
                                                ttGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCaa                                           <  1:3044223/44‑1 (MQ=255)
                                                 agATCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAAccg                       >  1:1166034/2‑62 (MQ=255)
                                                 agATCCACGTAGGAGTCCATACGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAAcc                        >  1:4609802/2‑62 (MQ=255)
                                                  gATCGACATAAGAATCCATCCGTCCGGCGGTCGTTg                                                 >  1:1441623/1‑36 (MQ=255)
                                                  gATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCaa                                           >  1:4746208/1‑42 (MQ=255)
                                                  gATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAaga                                         >  1:5664825/1‑44 (MQ=255)
                                                  gATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAag                                          >  1:63370/1‑43 (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)
                                                                     |                                                               
GCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCC  >  Exported/1014671‑1014803

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: