Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,014,602 G→C 18.4% A489A (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,6020GC18.4% 141.1 / 25.0 65A489A (GCC→GCGfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base G (27/26);  new base C (5/7);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 7.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

GCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCT  >  Exported/1014534‑1014666
                                                                    |                                                                
gCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGc                                                                >  1:838174/1‑71 (MQ=255)
gCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGc                                                                >  1:2413168/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)
    ccATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGa                                                            >  1:3342986/1‑71 (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)
                 cGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCg                                                              >  1:5482691/1‑56 (MQ=255)
                 cGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCg                                                              >  1:4276901/1‑56 (MQ=255)
                 cGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCg                                                              >  1:2233633/1‑56 (MQ=255)
                 cGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCg                                                              >  1:1588932/1‑56 (MQ=255)
                  ggATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGcc                                              <  1:3217333/71‑1 (MQ=255)
                       ttCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCgg                                                             >  1:4067478/1‑51 (MQ=255)
                          ctcacaTTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAg                                       >  1:1360593/3‑70 (MQ=255)
                          aaCACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAg                                       <  1:2520645/70‑1 (MQ=255)
                          aaCACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAg                                       <  1:2641148/70‑1 (MQ=255)
                          aaCACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAg                                       <  1:1980413/70‑1 (MQ=255)
                          aaCACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAg                                       <  1:4134195/70‑1 (MQ=255)
                           acacATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGac                                                           <  1:65243/49‑1 (MQ=255)
                                 tCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCaa                                            >  1:5568517/1‑58 (MQ=255)
                                 tCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAgccgcc                               >  1:924183/1‑71 (MQ=255)
                                    aGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCgtg                            <  1:3329862/71‑1 (MQ=255)
                                    aGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCgt                             <  1:2919898/70‑1 (MQ=255)
                                      gCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAgccgcc                               <  1:1879519/66‑1 (MQ=255)
                                        tctTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGa                                                     <  1:2680070/42‑1 (MQ=255)
                                         ctTAATACTTCCCTGCGCCAATACAGCGGCCGGac                                                           <  1:4384918/35‑1 (MQ=255)
                                         ctTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGttttt                       >  1:4765620/1‑71 (MQ=255)
                                          ttAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCa                                             <  1:462459/48‑1 (MQ=255)
                                           aaaTACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATa                                                   >  1:4431457/2‑41 (MQ=255)
                                            aaTACTTCCCTGCGCCAATACAGCGGCCGGACCACCGAt                                                    >  1:5015040/1‑39 (MQ=255)
                                            aaTACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATa                                                   >  1:1266360/1‑40 (MQ=255)
                                            aaTACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATa                                                   >  1:4257964/1‑40 (MQ=255)
                                            aaTACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGa                                      >  1:1823041/1‑53 (MQ=255)
                                              tACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAgg                                          <  1:3616730/47‑1 (MQ=255)
                                                 ttCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAg                                       <  1:162245/47‑1 (MQ=255)
                                                   cccTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCa               <  1:5435483/69‑1 (MQ=255)
                                                   cccTGCGCCAATACAGCGGCCGGACCACCGATACAGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCa               <  1:3040774/69‑1 (MQ=255)
                                                     cTGCGCCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGcgt              <  1:2813926/68‑3 (MQ=255)
                                                     cTGCGCCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTcgg           >  1:5648506/1‑69 (MQ=255)
                                                     cTGCGCCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTcg            >  1:5141058/1‑69 (MQ=255)
                                                     cTGCGCCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTc             <  1:227697/69‑1 (MQ=255)
                                                     cTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCa               >  1:5498553/1‑67 (MQ=255)
                                                     cTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCg           <  1:406601/71‑1 (MQ=255)
                                                        cgcCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGtaa        >  1:4676832/1‑69 (MQ=255)
                                                        cgcCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGta         <  1:2331422/70‑2 (MQ=255)
                                                        cgcCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGta         <  1:2271464/70‑2 (MQ=255)
                                                        cgcCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGt          <  1:4717182/69‑1 (MQ=255)
                                                        cgcCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGt          <  1:2188830/69‑1 (MQ=255)
                                                        cgcCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGt          <  1:3096114/69‑1 (MQ=255)
                                                        cgcCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTc         <  1:1776041/70‑1 (MQ=255)
                                                           cAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACct     >  1:2937934/1‑71 (MQ=255)
                                                           cAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACct     >  1:5625313/1‑71 (MQ=255)
                                                            agtACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACct     >  1:3344664/3‑70 (MQ=255)
                                                             gtACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACctgc   >  1:401206/2‑71 (MQ=255)
                                                             aTACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGccgccg                              >  1:3263390/1‑44 (MQ=255)
                                                             aTACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCAt              <  1:5582778/60‑1 (MQ=255)
                                                              tACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATcc            <  1:2153796/61‑1 (MQ=255)
                                                              tACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACctgct  >  1:477448/1‑71 (MQ=255)
                                                               aCAGCGGCCGGACCACCGATACTGACAAGGTAGAAGCCGCCGTGTTTTTTAc                    >  1:2034702/1‑52 (MQ=255)
                                                                    ggCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTAc                    <  1:3324028/47‑1 (MQ=255)
                                                                    |                                                                
GCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCT  >  Exported/1014534‑1014666

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: