Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,014,548 T→C 15.5% E507E (GAA→GAG fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,014,5480TC15.5% 161.7 / 20.4 71E507E (GAA→GAGfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base T (33/27);  new base C (5/6);  total (38/33)
Fisher's exact test for biased strand distribution p-value = 7.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.19e-01

TTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATAC  >  Exported/1014480‑1014618
                                                                    |                                                                      
ttgtcatcGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAGCTc                                                                      >  1:1254032/4‑71 (MQ=255)
     atcGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAGCTCCGGa                                                                  >  1:2639127/1‑70 (MQ=255)
     atcGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAGCTCCGGa                                                                  >  1:5086069/1‑70 (MQ=255)
        cacaAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGatat                                                               >  1:1079912/1‑70 (MQ=255)
        cacaAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGatat                                                               >  1:1175426/1‑70 (MQ=255)
        cacaAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATAtt                                                              >  1:453674/1‑71 (MQ=255)
           aaGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTc                                                             >  1:5501114/1‑69 (MQ=255)
                tAAACGCCGGGAAATCTTCTACTTCGATTTTCAAGATAGCTTCCATACCCAGCTCCGGa                                                                  >  1:275311/1‑59 (MQ=255)
                tAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCg                                                                    <  1:1877283/57‑1 (MQ=255)
                tAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACAt                                                      >  1:2076490/1‑71 (MQ=255)
                  aaCGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCaa                                                           <  1:10312/64‑1 (MQ=255)
                    cGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCa                                                  <  1:4026277/71‑1 (MQ=255)
                    cGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCa                                                  <  1:4032262/71‑1 (MQ=255)
                    cGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCa                                                  <  1:358922/71‑1 (MQ=255)
                       cGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAGCTCCGGATAAGCGACGCACTCCAgat                                               <  1:1093234/71‑3 (MQ=255)
                        gggAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGct                                              <  1:2595090/71‑1 (MQ=255)
                         ggAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAGCTCCg                                                                    <  1:3488393/48‑1 (MQ=255)
                         ggAAATCTTCTACTTCGATTTTCCAGATAGATTCCATACCCAGCTCCGGATAAGCGACGCACTCCAgatg                                              >  1:5477308/1‑67 (MQ=255)
                           aaaTCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTc                                                                      <  1:1354722/44‑1 (MQ=255)
                           aaaTCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTc                                                                      <  1:4764458/44‑1 (MQ=255)
                           aaaTCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCtt                                           <  1:1585056/71‑1 (MQ=255)
                               cTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATa                                       >  1:5071271/1‑71 (MQ=255)
                                ttCTACTTCGATTTTCCAGATAGCTTCCATACCCAGCTcc                                                                     <  1:4732056/40‑1 (MQ=255)
                                ttCTACTTCGATTTTCCAGATAGCTTCCATACCCAGCTcc                                                                     <  1:2219373/40‑1 (MQ=255)
                                ttCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTc                                                             <  1:1886663/48‑1 (MQ=255)
                                ttCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTc                                                             <  1:351022/48‑1 (MQ=255)
                                ttCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATa                                       >  1:3467154/1‑70 (MQ=255)
                                   cACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAgg                                                >  1:5163477/1‑58 (MQ=255)
                                   cACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATa                                       >  1:388147/1‑67 (MQ=255)
                                   cACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATa                                       >  1:3412052/1‑67 (MQ=255)
                                   cACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACtt                                    <  1:4335428/70‑1 (MQ=255)
                                   cACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATAATCAACACATTCCAGGCTCTTAATACtt                                    <  1:2827364/70‑1 (MQ=255)
                                   cACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCAGGATATTCAACACATTCCAgg                                                >  1:3877106/1‑58 (MQ=255)
                                   cACTTCAATTTTCCAGATGGCTTCCACTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTc                                   >  1:5353610/1‑71 (MQ=255)
                                    aCTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTGAACACATTCCAGGCTCTTAATACTTcc                                  >  1:5670052/1‑71 (MQ=255)
                                    aCTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTcc                                  >  1:3431314/1‑71 (MQ=255)
                                     cTTCGATTTTCCAGATAGCTTCCATACCCAGCTCCGGat                                                                 <  1:5003720/39‑1 (MQ=25)
                                     cTTCGATTTTCCAGATAGCTTCCATACCCAGCTCCGGat                                                                 <  1:3316729/39‑1 (MQ=25)
                                     cTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGc                                               <  1:499907/57‑1 (MQ=255)
                                      ttCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAgg                                                >  1:1280102/1‑55 (MQ=255)
                                        cAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCaa                                                           >  1:3713759/1‑42 (MQ=255)
                                        cAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAgg                                                >  1:1895583/1‑53 (MQ=255)
                                        cAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAg                                                 >  1:4543578/1‑52 (MQ=255)
                                           ttttCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCa                                                  <  1:3468184/48‑1 (MQ=255)
                                           ttttCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCa                                                  <  1:2060825/48‑1 (MQ=255)
                                           ttttCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATa                                       >  1:4637325/1‑59 (MQ=255)
                                           ttttCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTgcgc                            >  1:4608584/1‑70 (MQ=255)
                                            tttCCAGATGGCTTCCATTCCCAGTTCCGGATATTc                                                             <  1:5380084/36‑1 (MQ=255)
                                             ttCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACtt                                    >  1:5137394/1‑60 (MQ=255)
                                             ttCCAGATGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCaa                         <  1:839897/70‑1 (MQ=255)
                                                cAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATa                       <  1:4875235/70‑1 (MQ=255)
                                                   aTGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTaa                                         >  1:1433669/1‑49 (MQ=255)
                                                      gCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGc                >  1:2413168/1‑71 (MQ=255)
                                                      gCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGc                >  1:838174/1‑71 (MQ=255)
                                                      gCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGc                <  1:3358408/71‑1 (MQ=255)
                                                       cTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGcc               <  1:1576348/71‑1 (MQ=255)
                                                       cTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGcc               >  1:836335/1‑71 (MQ=255)
                                                       cTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGcc               <  1:1383428/71‑1 (MQ=255)
                                                         tCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGc                   >  1:4193030/1‑65 (MQ=255)
                                                          ccATTCCCAGTTCCGGATATTCAACACATTCCAgg                                                >  1:652606/1‑35 (MQ=255)
                                                          ccATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGa            >  1:3342986/1‑71 (MQ=255)
                                                             ttCCCAGTTCCGGATATTCAACACATTCCAGGCTCtt                                           <  1:381481/37‑1 (MQ=255)
                                                             ttCCCAGTTCCGGATATTCAACACATTCCAGGCTCtt                                           <  1:716008/37‑1 (MQ=255)
                                                             ttCCCAGTTCCGGATATTCAACACATTCCAGGCTCtt                                           <  1:1178778/37‑1 (MQ=255)
                                                               cccAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCaa                         >  1:4928335/1‑53 (MQ=255)
                                                               cccAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATa                       >  1:5630077/1‑55 (MQ=255)
                                                               cccAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGaccacc       >  1:5442996/1‑71 (MQ=255)
                                                               cccAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGaccacc       >  1:1637925/1‑71 (MQ=255)
                                                                ccAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGaccacc       <  1:1525995/70‑1 (MQ=255)
                                                                 cAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCg              >  1:1622604/1‑62 (MQ=255)
                                                                   gTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGcc               <  1:5281274/59‑1 (MQ=255)
                                                                    ttCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATAc  <  1:607906/71‑1 (MQ=255)
                                                                    |                                                                      
TTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATAC  >  Exported/1014480‑1014618

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