Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,015,463 T→C 17.0% L202L (TTA→TTG) ‡ fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,015,4630TC17.0% 191.2 / 29.9 88L202L (TTA→TTG) ‡fumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base T (37/36);  new base C (6/9);  total (43/45)
Fisher's exact test for biased strand distribution p-value = 5.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAAT  >  Exported/1015395‑1015532
                                                                    |                                                                     
ggaCAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTaa                                                                     >  1:666597/1‑71 (MQ=255)
   cAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGc                                                                  >  1:4211564/1‑71 (MQ=255)
          ccgTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGtt                                                           >  1:5261989/1‑71 (MQ=255)
          ccgTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGtt                                                           <  1:3999617/71‑1 (MQ=255)
            gTACCGAGGGTACACATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGATTTGGTTTc                                                         >  1:5678796/1‑71 (MQ=255)
            gTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTc                                                         >  1:3458471/1‑71 (MQ=255)
            gTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTc                                                         >  1:3451956/1‑71 (MQ=255)
                 cAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGa                                                     >  1:2520491/1‑70 (MQ=255)
                 cAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGAt                                                    <  1:2068252/71‑1 (MQ=255)
                 cAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGAt                                                    <  1:2061669/71‑1 (MQ=255)
                  aGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGa                                                     >  1:1410421/1‑69 (MQ=255)
                     gTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATagag                                                <  1:3478245/71‑1 (MQ=255)
                     gTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATa                                                   >  1:1316827/1‑68 (MQ=255)
                      tGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTgg                                                             <  1:2885404/57‑1 (MQ=255)
                      tGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAgaga                                               <  1:1084939/71‑1 (MQ=255)
                      tGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAgaga                                               <  1:3181154/71‑1 (MQ=255)
                        cgcATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGAt                                              <  1:4789657/70‑1 (MQ=255)
                        cgcATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATa                                             >  1:5430970/1‑71 (MQ=255)
                          cATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATAc                                            <  1:2189701/70‑1 (MQ=255)
                           aTCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGt                                          >  1:5617337/1‑71 (MQ=255)
                           aTCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGt                                          >  1:377443/1‑71 (MQ=255)
                            ttttctCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTg                                                              <  1:3261203/48‑1 (MQ=255)
                            ttttctCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGt                                                            <  1:4397292/50‑1 (MQ=255)
                              ttctCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGtttt                                       >  1:399374/1‑71 (MQ=255)
                              ttctCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTAGTAGAGATACGttt                                        >  1:3581867/1‑68 (MQ=255)
                              ttctCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGa                                                     >  1:5182584/1‑57 (MQ=255)
                              ttctCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGa                                                     >  1:5276588/1‑57 (MQ=255)
                              ttctCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGa                                                     >  1:466936/1‑57 (MQ=255)
                              ttctCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCt                                        >  1:3010719/1‑70 (MQ=255)
                               tctcGACGAGGAAGTTTTTCAGTTTGCCGGGCGTCAGCATGGCTTTGGTTTCCTGGTAGAGAt                                              <  1:4334658/63‑1 (MQ=255)
                               tctcGACGAGGAAGTTTTTCAGTTTGCCGGGCGTCAGCAGGGCTTTGGTTTCCTGGTAGAGAt                                              <  1:5227161/63‑1 (MQ=255)
                               tctcGACGAGGAAGTTTTTCAGTTTGCCGGGCGTCAGCAGGGCTTTGGTTTCCTGGTAGAGAt                                              <  1:460578/63‑1 (MQ=255)
                               tctcGACGAGGAAGTTTTTCAGTTTGCCGGGCGTCAGCAGGGCTTTGGTTTCCTGGTAGAGAt                                              <  1:2717369/63‑1 (MQ=255)
                               tctcGACGAGGAAGTTTTTCAGTTTGCCGGGCGTCAGCAGGGCTTTGGTTTCCTGGTAGAGAt                                              <  1:2530501/63‑1 (MQ=255)
                               tctcAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTaa                                                                     >  1:3294336/1‑40 (MQ=255)
                               tctcAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATAc                                            <  1:4694561/65‑1 (MQ=255)
                               tctcAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATAc                                            <  1:5426159/65‑1 (MQ=255)
                                ctcAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCt                                        <  1:4258569/68‑1 (MQ=255)
                                 tcAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATagag                                                >  1:4370668/1‑59 (MQ=255)
                                  cAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCttt                                                               >  1:105406/1‑43 (MQ=255)
                                   aaCCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTaa                                                                     >  1:2427887/1‑36 (MQ=255)
                                   aaCCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGt                                     <  1:3157312/68‑1 (MQ=255)
                                   aaCCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGt                                     <  1:3086951/68‑1 (MQ=255)
                                   aaCCAGGCAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGt                                     <  1:5590001/68‑1 (MQ=255)
                                        ggTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGa                                                     <  1:1424968/47‑1 (MQ=255)
                                         gTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAac                            <  1:1417246/71‑1 (MQ=255)
                                         gTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAac                            <  1:4787093/71‑1 (MQ=255)
                                           aaTTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCt                                                       >  1:2734046/1‑42 (MQ=255)
                                           aaTTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCt                                                       >  1:2897564/1‑42 (MQ=255)
                                           aaTTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGcc                                >  1:625148/1‑65 (MQ=255)
                                           aaTTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAacca                          >  1:3940664/1‑71 (MQ=255)
                                            aTTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTg                                      <  1:4766448/58‑1 (MQ=255)
                                            aTTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAaccac                         <  1:3980521/71‑1 (MQ=255)
                                            aTTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAaccac                         <  1:4047212/71‑1 (MQ=255)
                                             tttttCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATa                                             >  1:1575859/1‑50 (MQ=255)
                                             tttttCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATa                                             >  1:3440364/1‑50 (MQ=255)
                                             tttttCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATa                                             >  1:5466074/1‑50 (MQ=255)
                                              ttttCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATAAGTTTTGttgg                                  >  1:2718277/1‑58 (MQ=255)
                                              ttttCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAaccacc                        <  1:3273283/70‑1 (MQ=255)
                                               tttCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAgccgccgc                      >  1:5472185/1‑69 (MQ=255)
                                               tttCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTc                                         <  1:4008520/52‑1 (MQ=255)
                                                 tCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGca                                >  1:2195172/1‑58 (MQ=255)
                                                 tCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGc                                 <  1:4831137/58‑1 (MQ=255)
                                                  cAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCt                    <  1:3567305/70‑1 (MQ=255)
                                                   aGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCtt                   <  1:1191864/70‑1 (MQ=255)
                                                   aGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACg                                           <  1:3921417/46‑1 (MQ=255)
                                                   aGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCg                               >  1:3458132/1‑58 (MQ=255)
                                                   aGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCg                               >  1:4799489/1‑58 (MQ=255)
                                                    gTTTCCCCAGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTg                 >  1:1338241/1‑71 (MQ=255)
                                                       tCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGa                              <  1:5494009/55‑1 (MQ=255)
                                                       tCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGc               >  1:1573830/1‑70 (MQ=255)
                                                       tCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGc               <  1:3608082/70‑1 (MQ=255)
                                                       tCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGc               <  1:3885495/70‑1 (MQ=255)
                                                         cccGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGa             <  1:1886545/70‑1 (MQ=255)
                                                         cccGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGa             <  1:862705/70‑1 (MQ=255)
                                                            ggCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACa         <  1:668128/71‑1 (MQ=255)
                                                             gCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCg                               >  1:1130160/1‑48 (MQ=255)
                                                             gCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCg                               >  1:3268327/1‑48 (MQ=255)
                                                               gTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCtt                                       >  1:2286384/1‑38 (MQ=255)
                                                                tCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCg                               <  1:5152596/45‑1 (MQ=255)
                                                                tCAGTAACGCTTTGGTTTCCTAATAGAGATACGTCTTg                                      <  1:594511/38‑1 (MQ=255)
                                                                 cAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCg                               <  1:1611392/44‑1 (MQ=255)
                                                                  aGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTgg                >  1:247607/1‑58 (MQ=255)
                                                                  aGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGaa   >  1:3129290/1‑71 (MQ=255)
                                                                  aGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGaa   >  1:1553698/1‑71 (MQ=255)
                                                                  aGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGaa   >  1:5405961/1‑71 (MQ=255)
                                                                   gTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGttt                                   <  1:2255237/38‑1 (MQ=255)
                                                                   gTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAAt  <  1:3288018/71‑1 (MQ=255)
                                                                    tAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCg                               <  1:1165737/41‑1 (MQ=255)
                                                                    tAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAaccaccac                      <  1:1364505/50‑1 (MQ=255)
                                                                    tAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAaccaccac                      <  1:629610/50‑1 (MQ=255)
                                                                    tAACGCTTTGGGTTCCTGATAGAGATACGTCTTGTTTGCCGAaccaccac                      <  1:1187300/50‑1 (MQ=255)
                                                                    |                                                                     
GGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAAT  >  Exported/1015395‑1015532

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: