Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,014,728 A→G 23.1% Y447Y (TAT→TAC fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,014,7280AG23.1% 160.6 / 46.9 91Y447Y (TAT→TACfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base A (40/30);  new base G (13/8);  total (53/38)
Fisher's exact test for biased strand distribution p-value = 8.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCG  >  Exported/1014660‑1014798
                                                                    |                                                                      
aCCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACAt                                                                       >  1:2767763/1‑70 (MQ=255)
    gctggctgCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGaa                                                                  >  1:3249287/1‑71 (MQ=255)
    gctgactgCGGTTACCTTTCGCCAGCATGATCATGCTGCCGCCGTGGGATTGCAGCAGATCCACGTAGGag                                                                  >  1:472056/1‑70 (MQ=255)
     ctggctgCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGa                                                                   <  1:3328560/69‑1 (MQ=255)
     ctggctgCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCAACATAAGa                                                                   <  1:2458362/69‑1 (MQ=255)
     ctgactgCGGTTACCTTTCGCCAGCATGATCATGCTGCCGCCGTGGGATTGCAGCAGATCCACGTa                                                                      >  1:5578970/1‑66 (MQ=255)
       gactgCGGTTACCTTTCGCCAGCATGATCATGCTGCCGCCGTGGGATTGCAGCAGATCCACGTAGGAgtc                                                                >  1:2740101/3‑70 (MQ=255)
         ctgCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAatccat                                                             >  1:919753/1‑71 (MQ=255)
         ctgCGGTTACCTTTCGCCAGCATGATCATGCTGCCGCCGTGGGATTGCAGCAGATCCACGTAGGAgtccat                                                             <  1:635382/71‑1 (MQ=255)
          tgCGGTTACCTTTCGCCAGCATGATCATGCTGCCGCCGTGGGATTGCAGCAGATCCACGTAGGAgtccata                                                            >  1:33357/1‑70 (MQ=255)
           gCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAatccatcc                                                           >  1:4892227/1‑71 (MQ=255)
              gTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAatcca                                                              <  1:3370554/65‑1 (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:3405713/2‑69 (MQ=255)
                 accTTTCGCCAGCATGATCATGCTGCCGCCGTGGGATTGCAGCAGATCCACGTAGGAGTCCAtacggcctg                                                     >  1:3270859/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)
                    tttCGCCAGCATGATCATGCTGCCGCCGTGGGATTGCAGCAGATCCACGTAGGAgtcca                                                              <  1:4624087/59‑1 (MQ=255)
                     ttCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGt                                                 >  1:3373890/1‑71 (MQ=255)
                      tCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACa                                                                        >  1:2347324/1‑47 (MQ=255)
                      tCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACa                                                                        >  1:1007613/1‑47 (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)
                        gCCAGCATGATCATACTTCCGCCCTGCGCTAGCAGTTGATCGACATAAGaa                                                                  >  1:179943/1‑51 (MQ=255)
                         ccAGCATGATCATGCTGCCGCCGTGGGATTGCAGCAGATCCACGTAGGAGTCCATACGGCCTGCGGTGGtt                                             <  1:5185377/71‑1 (MQ=255)
                         ccAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAatc                                                                <  1:3114174/52‑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)
                               tGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATa                                                                      <  1:4895398/40‑1 (MQ=255)
                               tGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATa                                                                      <  1:469015/40‑1 (MQ=255)
                               tGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATa                                                                      <  1:2032453/40‑1 (MQ=255)
                               tGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTg                                            >  1:1444583/1‑66 (MQ=255)
                               tGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGcc                                        <  1:4603298/70‑1 (MQ=255)
                                gATCATACTTCCGCCCTGTGCTTGCAGTTGATCGACATAAGaa                                                                  <  1:1828306/43‑1 (MQ=255)
                                gATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGtt                                             <  1:2871881/64‑1 (MQ=255)
                                 aTCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAatccatcc                                                           >  1:2968097/1‑49 (MQ=255)
                                 aTCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAatccatcc                                                           >  1:3933305/1‑49 (MQ=255)
                                 aTCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAat                                                                 <  1:4334084/43‑1 (MQ=255)
                                 aTCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAat                                                                 <  1:430467/43‑1 (MQ=255)
                                 aTCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAat                                                                 <  1:3905659/43‑1 (MQ=255)
                                 aTCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGcc                                        >  1:520481/1‑68 (MQ=255)
                                  tCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAg                                                                    <  1:1201293/39‑1 (MQ=255)
                                  tCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAg                                                                    <  1:1640576/39‑1 (MQ=255)
                                  tCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAg                                                                    <  1:5372804/39‑1 (MQ=255)
                                  tCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCg                                               >  1:4616898/1‑60 (MQ=255)
                                    aTACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAagag                                   >  1:5288592/1‑70 (MQ=255)
                                      aCTTCCGCCCTGCGCTTGCAGTTGATCGACATAAg                                                                    >  1:917123/1‑35 (MQ=255)
                                      aCTTCCGCCCTGCGCTTGCAGTTGATCGACATAAg                                                                    >  1:276670/1‑35 (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)
                                           cGCCCTGCGCTTGCAGTTGATCGACATAAGAatcca                                                              >  1:1688855/1‑36 (MQ=255)
                                           cGCCCTGCGCTTGCAGTTGATCGACATAAGAatcca                                                              >  1:4922593/1‑36 (MQ=255)
                                               cTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAaga                                    >  1:3284340/1‑58 (MQ=255)
                                               cTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAaga                                    >  1:2758279/1‑58 (MQ=255)
                                               cTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAaga                                    >  1:73698/1‑58 (MQ=255)
                                               cTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAaga                                    >  1:3919429/1‑58 (MQ=255)
                                               cTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAaga                                    >  1:5226654/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:1098524/1‑40 (MQ=255)
                                                       gCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGt                                              >  1:61690/1‑40 (MQ=255)
                                                       gCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGt                                              >  1:1125083/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)
                                                                    |                                                                      
ACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCG  >  Exported/1014660‑1014798

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: