Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,014,644 T→A 16.3% K475N (AAA→AAT) ‡ fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,014,6440TA16.3% 106.3 / 13.8 43K475N (AAA→AAT) ‡fumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base T (13/23);  new base A (4/3);  total (17/26)
Fisher's exact test for biased strand distribution p-value = 4.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

AATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTT  >  Exported/1014578‑1014714
                                                                  |                                                                      
aaTACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGttttt                                                                       <  1:3355991/68‑1 (MQ=255)
aaTACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTa                                                                     >  1:2288056/1‑70 (MQ=255)
    cTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAgg                                                                 >  1:662359/1‑70 (MQ=255)
       cccTGCGCCTATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCa                                                               <  1:2495912/69‑1 (MQ=255)
       cccTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGtttttttc                                                                    >  1:2246774/1‑62 (MQ=255)
       cccTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTa                                                                     >  1:1088155/1‑63 (MQ=255)
         cTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCg                                                           <  1:3492995/71‑1 (MQ=255)
          tGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGt                                                          >  1:275644/1‑71 (MQ=255)
            cgcCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAg                                                                  >  1:156096/1‑61 (MQ=255)
            cgcCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTc                                                         >  1:266624/1‑70 (MQ=255)
            cgcCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCa                                                        >  1:3186191/1‑71 (MQ=255)
                  tACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTaa                                                                    >  1:1909967/1‑52 (MQ=255)
                  tACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACct                                                     >  1:3255236/1‑68 (MQ=255)
                   aCAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAg                                                                  >  1:1803974/1‑54 (MQ=255)
                    cAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTgctgg                                                <  1:398530/71‑1 (MQ=255)
                    cAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTgctgg                                                <  1:495580/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)
                              accaccGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTAc                                                                    <  1:2414400/41‑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:1037621/1‑71 (MQ=255)
                                  ccGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCttt                                  >  1:3542282/1‑71 (MQ=255)
                                  ccGATACTGCCAATGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACct                                                     >  1:705339/1‑52 (MQ=255)
                                  ccGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTa                                                                     <  1:2675090/36‑1 (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:2549758/1‑67 (MQ=255)
                                                                tttttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCt        >  1:1178129/1‑67 (MQ=255)
                                                                 ttttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCt   >  1:2443948/1‑71 (MQ=255)
                                                                 ttttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCt   >  1:117427/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)
                                                                  |                                                                      
AATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTT  >  Exported/1014578‑1014714

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: