Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 760,485 G→A G594S (GGC→AGC)  sucA → subunit of E1(0) component of 2‑oxoglutarate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913760,4850GA96.0% 137.6 / ‑5.9 50G594S (GGC→AGC) sucAsubunit of E1(0) component of 2‑oxoglutarate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (25/23);  minor base C (1/1);  total (26/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

AAGTTGAGATGAAGCGCCTGCAGGAGCTGGCGAAACGCATCAGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGGGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCACCGCCACGCGGTGATCCACAACCAGTCTAACGGTTCCACTTACACGC  >  NC_000913/760348‑760626
                                                                                                                                         |                                                                                                                                             
aggtgcggaTGAAGCGCCTGCAGGAGCTGGCGAAACGCATCAGCACGGTGCCGGAGGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGt                                                                                                                                          <  2:310286/136‑1 (MQ=255)
 agctGAGATGAAGCGCCTGCAGGAGCTGGCGAAACGCAGCAGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTg                                                                                                                                         <  1:503559/140‑1 (MQ=255)
          gAAGCGCCTGCAGGAGCTGGCGAAACGCATCAGCACGGTGCCGGAAGCAGTGGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGGTGGACTGGAGCGGCACGGAAAAcc                                                                                                                                >  2:302343/1‑143 (MQ=255)
            tgcggCTGCAGGAGCTGGCGAAAGGCATCAGCACGGTGCCGGAAGCAGGTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTc                                                                                                                              <  2:68677/142‑1 (MQ=255)
            aGCGCCTGCAGCAGCTGGCGAAACGCATCAGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTc                                                                                                                              <  1:20179/143‑1 (MQ=255)
             gcgcCTGCAGGAGCTGGCGAAACGCATCAGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCg                                                                                                                             >  1:1079108/1‑143 (MQ=255)
               gcCCGCAGGAGCCGGCGAAACGCCTCAGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAACATTTATGGCGATCGCCACGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCt                                                                                                                           <  1:844596/143‑1 (MQ=255)
                   gCAGGAGCTGGCGAAACGCATCAGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACg                                                                                                                       <  2:627957/143‑1 (MQ=255)
                     aGGAGCTGGCGAAACGCATCAGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGcc                                                                                                                     >  1:401584/1‑143 (MQ=255)
                       gAGCTGGCGAAACGCATCAGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGAGTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCCTACGGCAc                                                                                                                   >  2:237883/1‑143 (MQ=255)
                         gCTGGCGAAACGCATCAGCACGGTGCCGGAAGCAGTTGAAATGCTGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGc                                                                                                                 <  2:974132/143‑1 (MQ=255)
                            ggCGAAACGCATCAGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTgg                                                                                                              >  1:290532/1‑143 (MQ=255)
                                    gCATCAGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGaa                                                                                                      <  2:13068/143‑1 (MQ=255)
                                    gCACCAGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGaa                                                                                                      <  2:291667/143‑1 (MQ=255)
                                     cATCAGCACGGTGCCGGAAGCAGTTGCAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCGAGGCGCTGGCTGCCGGTGAGCAACTGTTCGACTGGCGCGGTGCGGAAAACCTCGCTTACGCCCCGCTGGTTGACGAAg                                                                                                     >  2:477858/1‑143 (MQ=255)
                                           cACGGTGCCGGAAGCAGGTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAATACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTc                                                                                               <  1:723694/143‑1 (MQ=255)
                                                  ccGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGAGTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGCTGAAGGCCTTCCGGCTCg                                                                                        >  1:976324/1‑143 (MQ=255)
                                                   cGGACGCAGGTGAAATGCAGTCTCGCGGTGCCAGGATTTAGGGCGAGCGCCAGGCGAGGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGc                                                                                       <  2:358645/143‑1 (MQ=255)
                                                   cGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGc                                                                                       >  1:642412/1‑143 (MQ=255)
                                                    ggAGGCAGTTGAAATGCAGTCTCGCGTGGCCAAGATGTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGcc                                                                                      <  2:638715/143‑1 (MQ=255)
                                                        gCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGACAACCTCGCTTACGCCACGCTGGTTGTTGACGGCATTCCGGTTCGCCTGTc                                                                                  >  1:186608/1‑143 (MQ=255)
                                                          aGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGCCGATGGCTGCCGGTGAGAAACTGTTCGACTGGCGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCgg                                                                                <  1:485853/143‑1 (MQ=255)
                                                           gTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTCGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGCAGGCATTCCGGTTCGCCTGTCggg                                                                               >  1:415420/1‑143 (MQ=255)
                                                           gTTGAAATGCAGTCGCGCTTTGCCGAGATGTGTGGCGATCGCCAGGCGATGGCTGCCGGGGAGAAACTGTTCGACTGGAGCGGTGCGGAAACCGTCTCTTTCGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCggg                                                                               <  2:913456/143‑1 (MQ=255)
                                                            ttGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGCCTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCggat                                                                              >  2:186437/1‑141 (MQ=255)
                                                            ttGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGATGATGAAGGCAGTCCGGTTCGCCTGTCGGGt                                                                              >  2:338381/1‑143 (MQ=255)
                                                               aaaTGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGaa                                                                           <  1:35092/143‑1 (MQ=255)
                                                                aaTGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCACGCGATGGCTGCCGATGAGAAACTGTTCGACTGGAGCGGGGCGGACAAACTCGCTTACGCCACGCTGGTTGATGCAGTCATTCCCGTTCTCCTGTCGGGTGTAg                                                                          >  2:348206/1‑143 (MQ=255)
                                                                 aTGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGa                                                                         >  1:400928/1‑143 (MQ=255)
                                                                      ctctcGCCTTTCCAAGATCTATGCCGATCGCCAGGCGATGGCTGCCGGTCAGAAACTGTTCGACTGGAGCGTTGCGGAAAACCTCGCTTACCCCACGCTGGTTGATTAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCg                                                                    <  1:435548/142‑1 (MQ=255)
                                                                                 ccAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTAcc                                                         >  1:957467/1‑143 (MQ=255)
                                                                                 ccAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTAcc                                                         >  1:806174/1‑143 (MQ=255)
                                                                                        ttATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACATTCTTcc                                                  <  2:553440/143‑1 (MQ=255)
                                                                                           tGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCAcc                                               <  1:1087499/143‑1 (MQ=255)
                                                                                            ggCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTGGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAGTACTCCGGTCGCGGTACCTTCTTCCACCg                                              >  1:564732/1‑143 (MQ=255)
                                                                                              cGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGGTGAAGGCATTCCGGTTCGCCTGGCGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCTCCGcc                                            >  1:961087/1‑143 (MQ=255)
                                                                                                    ccgggCGATGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCACCGCCACGCgg                                      <  1:890368/140‑1 (MQ=255)
                                                                                                           aTGGCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGGAAACCCCGCTGACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCGGGTCGCGGTACGTTGTGCCACCACCGGCCGGTGGTCca                               >  2:87086/1‑143 (MQ=255)
                                                                                                             ggCTGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAACGCCGCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCACCGCCACGCGGTGATCcaca                             <  2:179968/143‑1 (MQ=255)
                                                                                                                tGCCGGTGAGAAACTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCACCGCCACGCGGTGATCCACaaac                          >  1:868208/1‑141 (MQ=255)
                                                                                                                       gagaAGGTGTGCGACTGGAGCGGTGCGGAAAGCCTCGCTCACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCACCGCCACGCGGTGATCCACAACCAGTCTaa                   <  2:221267/143‑1 (MQ=255)
                                                                                                                        aggaaCTGGTCGACTGGAGCGGTGCGGAAAGCCTCGCGTACGCCACGCTGGTGGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCACCGCCACGCGGTGATCCACAACCAGTCTAAc                  <  1:346481/140‑1 (MQ=255)
                                                                                                                           aaCTGTTCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCACCGCCACGCGGTGATCCACAACCAGTCTAACGGt               >  1:594268/1‑143 (MQ=255)
                                                                                                                            aCTGTTCGACTGGAGCGGTGCGGAAAACCTCTCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCACCGCCACGCGGTGATCCACAACCAGTCTAACGGtt              >  1:833098/1‑143 (MQ=255)
                                                                                                                                ttCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCACCGCCACGCGGTGATCCACAACCAGTCTAACGGTTCCAc          >  1:745678/1‑143 (MQ=255)
                                                                                                                                 tCGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCCGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACCTTCGTCCACCGCCACGCGGTGATCCACAACCCGTCCAACGCGTCCACt         >  2:539822/1‑143 (MQ=255)
                                                                                                                                  cGACTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCACCGCCACGCGGTGATCCACAACCAGTCTAACGGTTCCACtt        <  2:290532/143‑1 (MQ=255)
                                                                                                                                     cTGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGACGGCATTCCGGCTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCACCGCCACGCGGTGATCCACAACCAGTCTAACGGTTCCACTTcca     >  1:399244/1‑143 (MQ=255)
                                                                                                                                      tGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCACCg                                              <  1:730958/101‑1 (MQ=255)
                                                                                                                                      tGGAGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCACCg                                              >  2:730958/1‑101 (MQ=255)
                                                                                                                                        cggCGGTGCGCACGTCCTCGCTTACGCCCCGCTGGTTGATGAAGGCGTTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCACCGCCACGCGGTGATCCACAACCAGTCTAACGGTTCCACTTACAcgc  <  2:77470/142‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                             
AAGTTGAGATGAAGCGCCTGCAGGAGCTGGCGAAACGCATCAGCACGGTGCCGGAAGCAGTTGAAATGCAGTCTCGCGTTGCCAAGATTTATGGCGATCGCCAGGCGATGGCTGCCGGTGAGAAACTGTTCGACTGGGGCGGTGCGGAAAACCTCGCTTACGCCACGCTGGTTGATGAAGGCATTCCGGTTCGCCTGTCGGGTGAAGACTCCGGTCGCGGTACCTTCTTCCACCGCCACGCGGTGATCCACAACCAGTCTAACGGTTCCACTTACACGC  >  NC_000913/760348‑760626

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: