Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,014,653 A→G 15.2% D472D (GAT→GAC fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,014,6530AG15.2% 112.4 / 12.2 46D472D (GAT→GACfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base A (17/22);  new base G (4/3);  total (21/25)
Fisher's exact test for biased strand distribution p-value = 6.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGAT  >  Exported/1014585‑1014723
                                                                    |                                                                      
cccTGCGCCTATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCa                                                                        <  1:2495912/69‑1 (MQ=255)
  cTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCg                                                                    <  1:3492995/71‑1 (MQ=255)
   tGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGt                                                                   >  1:275644/1‑71 (MQ=255)
     cgcCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTc                                                                  >  1:266624/1‑70 (MQ=255)
     cgcCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCa                                                                 >  1:3186191/1‑71 (MQ=255)
           tACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACct                                                              >  1:3255236/1‑68 (MQ=255)
             cAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTgctgg                                                         <  1:495580/71‑1 (MQ=255)
             cAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTgctgg                                                         <  1:398530/71‑1 (MQ=255)
                  gCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTgctggc                                                        <  1:1049621/67‑1 (MQ=255)
                  gCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCg                                                    <  1:3287757/71‑1 (MQ=255)
                    cGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCg                                                    <  1:2535967/69‑1 (MQ=255)
                     ggACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCg                                                    <  1:2467512/68‑1 (MQ=255)
                     ggACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGtt                                                 <  1:3375748/71‑1 (MQ=255)
                     ggACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGtt                                                 <  1:3619571/71‑1 (MQ=255)
                          gccGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACctgc                                                            >  1:3591259/2‑55 (MQ=255)
                          gccGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCt                                             >  1:320690/2‑70 (MQ=255)
                          accGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCtt                                            <  1:3607258/71‑1 (MQ=255)
                           ccGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCttt                                           >  1:3542282/1‑71 (MQ=255)
                           ccGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCttt                                           >  1:1037621/1‑71 (MQ=255)
                           ccGATACTGCCAATGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACct                                                              >  1:705339/1‑52 (MQ=255)
                           ccGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACct                                                              >  1:1518011/1‑52 (MQ=255)
                            cGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTc                                          <  1:2295753/71‑1 (MQ=255)
                              aTACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCt                                             <  1:1736244/66‑1 (MQ=255)
                              aTACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGc                                        <  1:1788725/71‑1 (MQ=255)
                               tACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGcc                                       <  1:1873513/71‑1 (MQ=255)
                                 cTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAg                                     <  1:448507/71‑1 (MQ=255)
                                     cAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGcc                                       <  1:1702727/65‑1 (MQ=255)
                                      aaGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTg                                                <  1:1357149/55‑1 (MQ=255)
                                       aGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGAt                               >  1:1908248/1‑71 (MQ=255)
                                        ggTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATc                              <  1:1931665/71‑1 (MQ=255)
                                          tAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACctg                                                             <  1:2254137/38‑1 (MQ=255)
                                          gaGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCgg                                                   <  1:1323498/47‑1 (MQ=255)
                                            gAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCTTTCGCCAGCATGATCAtg                           <  1:2703189/70‑2 (MQ=255)
                                                    cgTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCTTTCGCCAGCATGATCATGCTgccgccg                  <  1:3049470/71‑2 (MQ=255)
                                                     gtgtTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCTTTCGCCAGCATGATCATGCTgccgccg                  <  1:1691351/70‑2 (MQ=255)
                                                        ttttttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTgcg              >  1:3358151/1‑71 (MQ=255)
                                                         tttttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCt                 >  1:1178129/1‑67 (MQ=255)
                                                         tttttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCt                 >  1:2549758/1‑67 (MQ=255)
                                                          ttttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCt            >  1:117427/1‑71 (MQ=255)
                                                          ttttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCt            >  1:2443948/1‑71 (MQ=255)
                                                          ttttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCt            >  1:1725149/1‑71 (MQ=255)
                                                           tttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTgc               >  1:66361/1‑67 (MQ=255)
                                                           tttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCtt           >  1:3504993/1‑71 (MQ=255)
                                                            ttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCtt           <  1:63543/70‑1 (MQ=255)
                                                              aCAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTg          >  1:3463334/1‑69 (MQ=255)
                                                              aCAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTg          >  1:1779442/1‑69 (MQ=255)
                                                                  gCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTg    <  1:1115746/71‑1 (MQ=255)
                                                                    aTCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGAt  >  1:256349/1‑71 (MQ=255)
                                                                    aTCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGAt  >  1:688963/1‑71 (MQ=255)
                                                                    aTCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCACCATGATCATACTTCCGCCCTGCGCTTGCAGTTGa   >  1:1528744/1‑70 (MQ=255)
                                                                    aTCCGTCACCTGCTGGCTGCGGTTGCCTTTAGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGAt  >  1:1630231/1‑71 (MQ=255)
                                                                    |                                                                      
CCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGAT  >  Exported/1014585‑1014723

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: