Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,014,668 G→A 13.7% S467S (AGC→AGT fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,014,6680GA13.7% 134.3 / 12.3 51S467S (AGC→AGTfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base G (23/21);  new base A (4/3);  total (27/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

CAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCC  >  Exported/1014598‑1014737
                                                                      |                                                                     
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:3619571/71‑1 (MQ=255)
        ggACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGtt                                                               <  1:3375748/71‑1 (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)
               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)
                             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:2549758/1‑67 (MQ=255)
                                            tttttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCt                               >  1:1178129/1‑67 (MQ=255)
                                             ttttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCt                          >  1:1725149/1‑71 (MQ=255)
                                             ttttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCt                          >  1:117427/1‑71 (MQ=255)
                                             ttttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCt                          >  1:2443948/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:1779442/1‑69 (MQ=255)
                                                 aCAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTg                        >  1:3463334/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)
                                                        tCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACtt                                      <  1:842147/48‑1 (MQ=255)
                                                        tCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATc               >  1:3298645/1‑71 (MQ=255)
                                                        tCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATc               >  1:2316786/1‑71 (MQ=255)
                                                          cGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCg              <  1:3010248/70‑1 (MQ=255)
                                                          cGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGa             <  1:96148/71‑1 (MQ=255)
                                                           gTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTc                                     >  1:3073524/1‑46 (MQ=255)
                                                            tCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGtt                   <  1:1729560/63‑1 (MQ=255)
                                                             cACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACa           >  1:1543197/1‑70 (MQ=255)
                                                              aCCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACAt          >  1:3553830/1‑70 (MQ=255)
                                                              aCCTGCTGACTGCGGTTACCTTTCGCCAGCATGATCATGCTGCCGCCGTGGGATTGCAGCAGATCCAcg           >  1:1696151/1‑68 (MQ=255)
                                                                 tgctgGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGa             >  1:3612077/1‑64 (MQ=255)
                                                                 tgctgGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAg       >  1:1654041/1‑70 (MQ=255)
                                                                 tgctgGCTGCGGTTGCCTTTCGCCAGCATGATCAAACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAg       >  1:2290286/1‑70 (MQ=255)
                                                                   ctgactgCGGTTACCTTTCGCCAGCATGATCATGCTGCCGCCGTGGGATTGCAGCAGATCCACGTAGGagt    >  1:3391991/1‑69 (MQ=255)
                                                                   ctgactgCGGTTACCTTTCGCCAGCATGATCATGCTGCCGCCGTGGGATTGCAGCAGATCCACGTAGGagt    >  1:1372482/1‑69 (MQ=255)
                                                                     ggctgCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGc                       >  1:1696675/1‑50 (MQ=255)
                                                                     ggctgCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGc                       >  1:860589/1‑50 (MQ=255)
                                                                      gctgCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGa                 <  1:3376674/55‑1 (MQ=255)
                                                                      gctgCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGa                 <  1:3372384/55‑1 (MQ=255)
                                                                      gctgCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAatcc  <  1:1438703/70‑1 (MQ=255)
                                                                      |                                                                     
CAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCC  >  Exported/1014598‑1014737

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: