Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,015,100 A→G 43.2% N323N (AAT→AAC fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,015,1000AG43.2% 35.4 / 111.5 95N323N (AAT→AACfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base A (30/24);  new base G (19/22);  total (49/46)
Fisher's exact test for biased strand distribution p-value = 4.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCA  >  Exported/1015033‑1015170
                                                                   |                                                                      
tatatTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGata                                                                        <  1:4241631/68‑1 (MQ=255)
tatatTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTa                                                                     <  1:1601331/71‑1 (MQ=255)
  tatTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGata                                                                        >  1:685870/1‑66 (MQ=255)
  tatTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGata                                                                        >  1:3493630/1‑66 (MQ=255)
    tttGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGt                                                                 <  1:562367/71‑1 (MQ=255)
     ttGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGt                                                                 >  1:5620087/1‑70 (MQ=255)
     tggCCTGGGTTGTGTTCCAGTTTTTCGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGt                                                                 >  1:3042698/3‑70 (MQ=255)
      ggCCTGGGTTGTGTTCCAGTTTTTCGATCCAGATACCTTCGCGGTTGATTTTAGCTTTAATGTTACGGt                                                                 >  1:5529572/2‑69 (MQ=255)
       gCCTGGGTTGTGTTCCAGTTTTTCGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTcag                                                              <  1:493611/71‑1 (MQ=255)
       gCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTcag                                                              >  1:1544333/1‑71 (MQ=255)
        ccTGGGTTGTGTTCCAGTTTTTCGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTcag                                                              <  1:2585167/70‑1 (MQ=255)
         cTGGGTTGTGTTCCAGTTTTTCGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTcagcg                                                            <  1:765686/71‑2 (MQ=255)
           ggATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTcagcag                                                           >  1:5459541/1‑70 (MQ=255)
           ggATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTcagca                                                            >  1:2699696/1‑69 (MQ=255)
           ggATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTcagca                                                            >  1:3118813/1‑69 (MQ=255)
           ggATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTAGGTcagcag                                                           >  1:2488653/1‑69 (MQ=255)
            gATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCagcagag                                                         >  1:4853656/1‑71 (MQ=255)
            gATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCagcagag                                                         <  1:4837800/71‑1 (MQ=255)
            gATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCagcagag                                                         >  1:3866169/1‑71 (MQ=255)
            gATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCagcagag                                                         >  1:4189154/1‑71 (MQ=255)
            gATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTAGCTTTGATATTACGGTCagcaga                                                          <  1:489632/70‑1 (MQ=255)
              ttATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCagcagagc                                                        <  1:2001343/70‑1 (MQ=255)
              ttATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCagcagag                                                         >  1:943428/1‑69 (MQ=255)
               tgtGTTCCAGTTTTTCGATGCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTCagcggagc                                                        <  1:4692495/67‑1 (MQ=255)
                gtGTTCCAGTTTTTCGATCCATATACCTTCGCGGTTGATTTTCGCTTTAATGTTACCGTCagcggagcag                                                      <  1:3984026/69‑1 (MQ=255)
                gtGTTCCAGTTTTTCGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTCagcggagcag                                                      <  1:2617848/69‑1 (MQ=255)
                gtGTTCCAGTTTTTCGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTCagcggagcag                                                      <  1:812238/69‑1 (MQ=255)
                 tGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGga                                                    <  1:4598355/71‑1 (MQ=255)
                 tGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGga                                                    <  1:5354398/71‑1 (MQ=255)
                 tGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGga                                                    <  1:3690239/71‑1 (MQ=255)
                 tGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGga                                                    <  1:1212523/71‑1 (MQ=255)
                 tGTTCCAGTTTTTCGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGga                                                    <  1:19175/71‑1 (MQ=255)
                  gTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGga                                                    <  1:2179569/70‑1 (MQ=255)
                  gTTCCAGTTTTTCGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGga                                                    <  1:1146097/70‑1 (MQ=255)
                   ttCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGgaga                                                  >  1:3603061/1‑71 (MQ=255)
                   ttCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGgaga                                                  >  1:1604107/1‑71 (MQ=255)
                   ttCCAGTTTTTCGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGgaga                                                  <  1:1525592/71‑1 (MQ=255)
                   ttCCAGTTTTTCGATCCAGATACCTTAGCGGTTGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGgaga                                                  <  1:4348262/71‑1 (MQ=255)
                    tCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGAc                                                 >  1:3339198/1‑71 (MQ=255)
                     ccAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGAc                                                 <  1:5284528/70‑1 (MQ=255)
                     ccAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACg                                                >  1:55464/1‑71 (MQ=255)
                     ccAGTTTTTCGATCCAGATCCCCTGACGGGTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGAc                                                 <  1:99413/70‑1 (MQ=255)
                     ccAGTTTTTCGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTa                                                                     >  1:5324530/1‑50 (MQ=255)
                     ccAGTTTTTCGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACg                                                >  1:1805971/1‑71 (MQ=255)
                      cAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCagcagagcag                                                      <  1:425225/64‑1 (MQ=255)
                      cAGTTTTTCGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGGCAGCGGAGCAGGAGACGc                                               <  1:4591599/71‑1 (MQ=255)
                       aGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCagcagagca                                                       >  1:2883009/1‑62 (MQ=255)
                        gTTTTTCGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTCagcggagcag                                                      <  1:988283/62‑1 (MQ=255)
                        gTTTTTCGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTCagcggagcag                                                      <  1:1152191/62‑1 (MQ=255)
                        gTTTTTCGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTCagcggagcag                                                      <  1:3605047/62‑1 (MQ=255)
                         tttttCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTcagcag                                                           <  1:1178992/56‑1 (MQ=255)
                         tttttCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTcag                                                              <  1:939032/53‑1 (MQ=255)
                         tttttCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTcag                                                              <  1:1659209/53‑1 (MQ=255)
                         tttttCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTcag                                                              <  1:4554022/53‑1 (MQ=255)
                         tttttCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGga                                                    <  1:3194806/63‑1 (MQ=255)
                           tttCGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTCagcggagca                                                       >  1:195876/1‑58 (MQ=255)
                           tttCGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGgaga                                                  >  1:1768927/1‑63 (MQ=255)
                           tttCGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACGCCCAtg                                          >  1:5106137/1‑70 (MQ=255)
                              cGATCCAGATACCTTCGGGGTTGATTTTCGCTTTAATGTTACGGTCAGCGGAGCTGGAGACGCCCATgcc                                        <  1:959164/70‑1 (MQ=255)
                              cGATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACGCCCATgcc                                        >  1:631225/1‑70 (MQ=255)
                               gATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACGCCCATgccga                                      <  1:2196846/71‑1 (MQ=255)
                                  ccAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTcagca                                                            >  1:1112260/1‑46 (MQ=255)
                                  ccAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTcagca                                                            >  1:1884346/1‑46 (MQ=255)
                                  ccAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACGCCCATgccgaccg                                   >  1:2889007/1‑71 (MQ=255)
                                     gATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTcagca                                                            <  1:4953286/43‑1 (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)
                                             gACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGgag                                                   >  1:5123485/1‑44 (MQ=255)
                                             cgcGGTTGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAgg                                                     >  1:1368071/3‑42 (MQ=25)
                                             cgcGGTTGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAgg                                                     >  1:2607614/3‑42 (MQ=25)
                                              gcGGTTGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACGCCCATGCCGACCGGGCAGGATgcgc                       >  1:678811/2‑71 (MQ=255)
                                              gcGGTTGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACGCCCATGCCGACCGGGCAGGATgc                         <  1:3486527/68‑1 (MQ=255)
                                              aCGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGcc                                              <  1:5279903/48‑1 (MQ=255)
                                                ggTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCg                     <  1:2419461/71‑1 (MQ=255)
                                                 gTTGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACGCCCATGCCGACCGGGCAGGATGCGCCGt                    <  1:1919808/71‑1 (MQ=255)
                                                  ttGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACGCCCa                                            >  1:4669302/1‑46 (MQ=255)
                                                  ttGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACGCCCa                                            >  1:2768811/1‑46 (MQ=255)
                                                  ttGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACGCCCa                                            >  1:492989/1‑46 (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:1126322/1‑71 (MQ=255)
                                                   tGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGa                  >  1:4244530/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)
                                                          cGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCa                                            >  1:226671/1‑38 (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)
                                                                   |                                                                      
TATATTTGCCTGGATTATGTTCCAGTTTTTCGATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCA  >  Exported/1015033‑1015170

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: