Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,015,097 G→A 5.8% I324I (ATC→ATT fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,015,0970GA 5.8% 224.5 / 3.1 69I324I (ATC→ATTfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base G (43/22);  new base A (2/2);  total (45/24)
Fisher's exact test for biased strand distribution p-value = 6.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

ATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCA  >  Exported/1015032‑1015161
                                                                 |                                                                
atatatTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGatat                                                              <  1:3536867/70‑1 (MQ=255)
atatatTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATAtt                                                             <  1:951932/71‑1 (MQ=255)
     tttGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGt                                                        >  1:2682209/1‑71 (MQ=255)
           tGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTcag                                                     >  1:2317259/1‑68 (MQ=255)
             gATTATGTTCCAGTTTTTGGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCagcagag                                                >  1:1421823/1‑71 (MQ=255)
             gATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCagcagag                                                <  1:2602334/71‑1 (MQ=255)
              aTTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCagcagag                                                >  1:463143/1‑70 (MQ=255)
              aTTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCagcagag                                                >  1:3226318/1‑70 (MQ=255)
              aTTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCagcagag                                                <  1:2092071/70‑1 (MQ=255)
                tgtGTTCCAGTTTTTCGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTCagcggagcag                                             >  1:519123/3‑71 (MQ=255)
                tATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGa                                                                 >  1:1293695/1‑51 (MQ=255)
                  tGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCagcagagcag                                             >  1:3240928/1‑69 (MQ=255)
                  tGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCagcagagcag                                             >  1:788041/1‑69 (MQ=255)
                  tGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCagcagagcag                                             >  1:2833994/1‑69 (MQ=255)
                  tGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGga                                           <  1:3470646/71‑1 (MQ=255)
                     tCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTcagca                                                   >  1:2387993/1‑60 (MQ=255)
                     tCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTcag                                                     >  1:2123218/1‑58 (MQ=255)
                       cAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGga                                           >  1:1528039/1‑66 (MQ=255)
                       cAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGga                                           >  1:2236462/1‑66 (MQ=255)
                        aGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATAtt                                                             <  1:1435428/47‑1 (MQ=255)
                        aGTTTTTCGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAatg                                                               >  1:1878441/1‑44 (MQ=255)
                          tttttCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCa                                   <  1:2873524/71‑1 (MQ=255)
                           ttttCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGAc                                        <  1:3178995/65‑1 (MQ=255)
                            tttCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATa                                 >  1:485722/1‑71 (MQ=255)
                            tttCGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACGCCCAt                                  <  1:1496480/70‑1 (MQ=255)
                              tCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGgaga                                         >  1:2360816/1‑61 (MQ=255)
                              tCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATacc                               >  1:3622661/1‑71 (MQ=255)
                               cGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACg                                       <  1:2870580/62‑1 (MQ=255)
                               cGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATacc                               >  1:2623169/1‑70 (MQ=255)
                                gATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATac                                >  1:2664655/1‑68 (MQ=255)
                                 aTCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGcc                                     >  1:2033016/1‑62 (MQ=255)
                                 aTCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATaccgac                            >  1:790025/1‑71 (MQ=255)
                                 aTCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATaccgac                            >  1:1527639/1‑71 (MQ=255)
                                 aTCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATaccgac                            >  1:858511/1‑71 (MQ=255)
                                 aTCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATaccgac                            >  1:3142814/1‑71 (MQ=255)
                                  tCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATaccgacc                           >  1:916533/1‑71 (MQ=255)
                                   ccAGATCCCCTGACGGTTGATCTTCGCTTTGATATTa                                                            >  1:3629084/1‑37 (MQ=255)
                                   ccAGATCCCCTGACGGTTGATCTTCGCTTTGATATTa                                                            >  1:3496109/1‑37 (MQ=255)
                                   ccAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATaccgaccg                          >  1:1643423/1‑71 (MQ=255)
                                     aGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCggg                        >  1:1920539/1‑71 (MQ=255)
                                     aGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCggg                        >  1:304036/1‑71 (MQ=255)
                                        tCCCCTGACGGTTGATCTTCGCTTTGATATTACgg                                                         >  1:519633/1‑35 (MQ=255)
                                        tCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTcag                                                     >  1:1458843/1‑39 (MQ=255)
                                        tCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATa                                 >  1:2240153/1‑59 (MQ=255)
                                         ccccTGACGGTTGATCTTCGCTTTGATATTACGGTcagca                                                   >  1:251179/1‑40 (MQ=255)
                                         ccccTGACGGTTGATCTTCGCTTTGATATTACGGTcagca                                                   >  1:2147126/1‑40 (MQ=255)
                                             tGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATaccga                             >  1:2903237/1‑58 (MQ=255)
                                             tGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATaccga                             >  1:1248502/1‑58 (MQ=255)
                                             tGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATgc                >  1:943717/1‑71 (MQ=255)
                                                 ggTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATg                 >  1:653342/1‑66 (MQ=255)
                                                 ggTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCg            >  1:1964882/1‑71 (MQ=255)
                                                  gTTGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACGCCCATGCCGACCGGGCAGGATGCGCCGt           <  1:1594562/71‑1 (MQ=255)
                                                  gTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATaccga                             <  1:1819366/53‑1 (MQ=255)
                                                  gTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATaccga                             <  1:2799486/53‑1 (MQ=255)
                                                  gTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATaccga                             <  1:1562135/53‑1 (MQ=255)
                                                  gTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATgc                <  1:1409692/66‑1 (MQ=255)
                                                  gTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATgc                <  1:1999935/66‑1 (MQ=255)
                                                  gTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGcc             <  1:973160/69‑1 (MQ=255)
                                                  gTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGcc             <  1:3526707/69‑1 (MQ=255)
                                                   ttGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGgaga                                         >  1:3225770/1‑40 (MQ=255)
                                                   ttGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGt           >  1:1798943/1‑70 (MQ=255)
                                                   ttGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGt           >  1:1145489/1‑70 (MQ=255)
                                                    tGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGt           <  1:1945554/69‑1 (MQ=255)
                                                        cTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATaccgac                            <  1:3016320/48‑1 (MQ=255)
                                                        cTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTg     >  1:2465474/1‑71 (MQ=255)
                                                         ttCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAgg                    <  1:1187418/55‑1 (MQ=255)
                                                           cGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTgg    <  1:3726445/69‑1 (MQ=255)
                                                           cGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTgg    <  1:356847/69‑1 (MQ=255)
                                                            gCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCa  <  1:2406213/70‑1 (MQ=255)
                                                                 |                                                                
ATATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCA  >  Exported/1015032‑1015161

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: