Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,015,130 A→G 40.3% G313G (GGT→GGC fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,015,1300AG40.3% 42.7 / 81.4 77G313G (GGT→GGCfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base A (26/20);  new base G (10/21);  total (36/41)
Fisher's exact test for biased strand distribution p-value = 6.17e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

CGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCA  >  Exported/1015063‑1015199
                                                                   |                                                                     
cGATCCAGATACCTTCGGGGTTGATTTTCGCTTTAATGTTACGGTCAGCGGAGCTGGAGACGCCCATgcc                                                                     <  1:959164/70‑1 (MQ=255)
cGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACGCCCATgcc                                                                     >  1:631225/1‑70 (MQ=255)
 gATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACGCCCATgccga                                                                   <  1:2196846/71‑1 (MQ=255)
    ccAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACGCCCATgccgaccg                                                                >  1:2889007/1‑71 (MQ=255)
          ccccTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAg                                                           >  1:2236504/1‑70 (MQ=255)
          ccccTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCACATACCGACCGGGCAg                                                           >  1:2294138/1‑70 (MQ=255)
           cccTCACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGa                                                         <  1:3037613/71‑1 (MQ=255)
             cTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGCCCGGGCAGGAt                                                        >  1:5521129/1‑70 (MQ=255)
             cTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATg                                                       >  1:4567470/1‑71 (MQ=255)
             cTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATg                                                       >  1:3316979/1‑71 (MQ=255)
              tGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATgc                                                      >  1:2567914/1‑71 (MQ=255)
                gcGGTTGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACGCCCATGCCGACCGGGCAGGATgcgc                                                    >  1:678811/2‑71 (MQ=255)
                gcGGTTGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACGCCCATGCCGACCGGGCAGGATgc                                                      <  1:3486527/68‑1 (MQ=255)
                  ggTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCg                                                  <  1:2419461/71‑1 (MQ=255)
                   gTTGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACGCCCATGCCGACCGGGCAGGATGCGCCGt                                                 <  1:1919808/71‑1 (MQ=255)
                    ttGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACGCCCATGCCGACCGGGCAGGATGCGCCGt                                                 >  1:938789/1‑70 (MQ=255)
                    ttGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGt                                                 >  1:4734037/1‑70 (MQ=255)
                     tGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGa                                               >  1:2493085/1‑71 (MQ=255)
                     tGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGa                                               >  1:4244530/1‑71 (MQ=255)
                     tGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGa                                               >  1:1126322/1‑71 (MQ=255)
                      gATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACGCCCATGCCGACCGGGCAGGATGCGCCGTGAc                                              >  1:5690256/1‑71 (MQ=255)
                      gATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACGCCCATGCCGACCGGGCAGGATGCGCCGTGAc                                              >  1:5197910/1‑71 (MQ=255)
                       aTCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATaccgaccg                                                                >  1:4283364/1‑52 (MQ=255)
                        ttttCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACGCCCATGCCGACCGGGCAGGa                                                         <  1:1891375/56‑1 (MQ=255)
                        ttttCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACGCCCATGCCGACCGGGCAGGATGCGCCGTGACGt                                            <  1:519775/69‑1 (MQ=255)
                            cGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGt                                            <  1:2708153/67‑1 (MQ=255)
                             gCTTTAATGTTACGGTCAGCGGAGCAGGAGACGCCCATGCCGACCGGGCAGGATGCGCCGTGAc                                              >  1:1503710/1‑64 (MQ=255)
                             gCTTTAATGTTACGGTCAGCGGAGCAGGAGACGCCCATGCCGACCGGGCAGGATGCGCCGTGACGTGGCa                                        <  1:820874/70‑1 (MQ=255)
                                  gATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGa                                   <  1:5118496/70‑1 (MQ=255)
                                   atatTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAAt                                 >  1:5249545/1‑71 (MQ=255)
                                     atTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCa                               >  1:4396330/1‑71 (MQ=255)
                                      ttACGGTCAGCGGAGCAGGAGACGCCCATGCCGACCGGGCAGGATGCGCCGTGACGTGGCAGACGGATAAc                              >  1:4556881/1‑71 (MQ=255)
                                        aCGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGAt                                                        >  1:2811572/1‑43 (MQ=255)
                                        aCGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGa                                               >  1:2761736/1‑52 (MQ=255)
                                        aCGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGa                                               >  1:146559/1‑52 (MQ=255)
                                        aCGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGa                                               >  1:3645640/1‑52 (MQ=255)
                                        aCGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCAcgc                            >  1:4825718/1‑71 (MQ=255)
                                         cGGTCAGCGGAGCAGGAGACGCCCATGCCGACCGGGCAGGATGCGCCGTGACGTGGCAGACGGATAAcgcg                           <  1:4992893/71‑1 (MQ=255)
                                          ggTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCAcgc                            <  1:124909/69‑1 (MQ=255)
                                          ggTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCgg                          <  1:5372228/71‑1 (MQ=255)
                                              agcggagcagGAGACGCCCATGCCGACCGGGCAGGATGCGCCGTGACGTGGCAGACGGATAACGCGAATGt                      >  1:5579161/1‑71 (MQ=255)
                                              agcagagcagGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATg                       <  1:4931657/70‑1 (MQ=255)
                                                cggagcagGAGACGCCCATGCCGACCGGGCAGGATGCGCCGt                                                 <  1:4277022/40‑1 (MQ=38)
                                                cggagcagGAGACGCCCATGCCGACCGGGCAGGATGCGCCGTGa                                               <  1:2728780/42‑1 (MQ=38)
                                                cagagcagGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCAcgcg                           <  1:746093/64‑1 (MQ=255)
                                                cagagcagGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCg                    <  1:3972203/71‑1 (MQ=255)
                                                 ggagcagGAGACGCCCATGCCGACCGGGCAGGATGCGCCGTGACGTGGCAGACGGATAACGCGAATGTCg                    <  1:4468816/69‑1 (MQ=255)
                                                 ggagcagGAGACGCCCATGCCGACCGGGCAGGATGCGCCGTGACGTGGCAGACGGATAACGCGAATGTCGt                   <  1:1066112/70‑1 (MQ=255)
                                                 agagcagGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCg                    <  1:1080131/70‑1 (MQ=255)
                                                 agagcagGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCg                    <  1:4041384/70‑1 (MQ=255)
                                                   agcagGAGACGCCCATGCCGACCGGGCAGGATGCGCCGTGACGTGGCAGACGGATAAcgcg                           <  1:2555907/61‑1 (MQ=255)
                                                     cagGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGc               <  1:44878/71‑1 (MQ=255)
                                                      agGAGACGCCCATGCCGACCGGGCAGGATGCGCCGt                                                 >  1:2062798/1‑36 (MQ=255)
                                                      agGAGACGCCCATGCCGACCGGGCAGGATGCGCCGTGACGTGGCAGACGGATAACGCGAATGTCGTGcgcg              <  1:1667475/71‑1 (MQ=255)
                                                       ggAGACGCCCATGCCGACCGGGCAGGATGCGCCGTGACGTGGCAGACGGATAACGCGAATGTCGTGcgcg              >  1:5622372/1‑70 (MQ=255)
                                                       ggAGACGCCCATGCCGACCGGGCAGGATGCGCCGTGACGTGGCAGACGGATAACGCGAATGTCGTGCGCGa             <  1:4194267/71‑1 (MQ=255)
                                                       ggAGACGCCCATACCGACCGGGCAGGATGCGCCGTGAc                                              >  1:43725/1‑38 (MQ=255)
                                                       ggAGACGCCCATACCGACCGGGCAGGATGCGCCGTGAc                                              >  1:3562757/1‑38 (MQ=255)
                                                         agaCGCCCATGCCGACCGGGCAGGATGCGCCGTGACGTGGCAGACGGATAACGCGAATGTCGTGCGCGAAg           <  1:1877428/71‑1 (MQ=255)
                                                         agaCGCCCATGCCGACCGGGCAGGATGCGCCGTGACGTGGCAGACGGATAACGCGAATGTCGTGCGCGAAg           <  1:2178286/71‑1 (MQ=255)
                                                         agaCGCCCATGCCGACCGGGCAGGATGCGCCGTGACGTGGCAGACGGATAACGCGAATGTCGTGCGCGAAg           <  1:550342/71‑1 (MQ=255)
                                                         agaCGCCCATGCCGACCGGGCAGGATGCGCCGTGACGTGGCAGACGGATAACGCGAATGTCGTGCGCGAAg           <  1:3884826/71‑1 (MQ=255)
                                                          gaCGCCCATGCCGACCGGGCAGGATGCGCCGTGACGTGGCAGACGGATAACGCGAATGTCGTGCGCGAAg           >  1:3118423/1‑70 (MQ=255)
                                                            cGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTAt        >  1:515647/1‑71 (MQ=255)
                                                             gCCCATGCCGACCGGGCAGGATGCGCCGTGACGTGGCAGACGGATAACGCGAATGTCGTGCGCGAAGTAtt       <  1:2853639/71‑1 (MQ=255)
                                                              cccATGCCGACCGGGCAGGATGCGCCGTGACGTGGCAGACGGATAACGCGAATGTCGTgcgc               <  1:2705980/62‑1 (MQ=255)
                                                              cccATGCCGACCGGGCAGGATGCGCCGTGACGTGGCAGACGGATAACGCGAATGTCGTgcgc               <  1:2241003/62‑1 (MQ=255)
                                                               ccATGCCGACCGGGCAGGATGCGCCGTGACGTGGCAGACGGATAACGCGAATGTCGTGCGCGAAGTATTTa     <  1:3622182/71‑1 (MQ=255)
                                                                cATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTac    >  1:3480242/1‑71 (MQ=255)
                                                                cATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTac    >  1:1605888/1‑71 (MQ=255)
                                                                 aTACCGACCGGGCAGTATGCGCCGTGACGTGGCAGGCGAAt                                 <  1:4568092/41‑1 (MQ=255)
                                                                 aTACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAAt                                 <  1:3047222/41‑1 (MQ=255)
                                                                 aTACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAAt                                 <  1:1725854/41‑1 (MQ=255)
                                                                 aTACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAAt                                 <  1:1732405/41‑1 (MQ=255)
                                                                 aTACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAAt                                 <  1:4156888/41‑1 (MQ=255)
                                                                  tACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTa         <  1:1456498/64‑1 (MQ=255)
                                                                  tACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTAttt      <  1:2331947/67‑1 (MQ=255)
                                                                  tACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTAttt      <  1:5133979/67‑1 (MQ=255)
                                                                  tACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTacca  >  1:5632792/1‑71 (MQ=255)
                                                                  tACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTacca  >  1:761644/1‑71 (MQ=255)
                                                                  tACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTacca  >  1:3721887/1‑71 (MQ=255)
                                                                  tACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTacca  >  1:4369366/1‑71 (MQ=255)
                                                                   |                                                                     
CGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCA  >  Exported/1015063‑1015199

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: