Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 759,237 A→G 15.8% T178A (ACC→GCC)  sucA → 2‑oxoglutarate decarboxylase, thiamine triphosphate‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913759,2370AG15.8% 161.7 / 25.8 76T178A (ACC→GCC) sucA2‑oxoglutarate decarboxylase, thiamine triphosphate‑binding
Reads supporting (aligned to +/- strand):  ref base A (35/29);  new base G (6/6);  total (41/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.97e-02

GCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACTTTCAATAGCGAAGAGAAAAAACGCTTCTTAAGCGAACTGACCGCCGCT  >  NC_000913/759138‑759335
                                                                                                   |                                                                                                  
gCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCCGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCa                                                                                                    <  2:1614452/100‑1 (MQ=255)
  cAGCGGCAAAGAAACCATGAATCTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCAcc                                                                                                  <  2:2067121/100‑1 (MQ=255)
  cAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCAcc                                                                                                  >  2:778114/1‑100 (MQ=255)
  cAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCAcc                                                                                                  <  2:713313/100‑1 (MQ=255)
   aGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCGCCg                                                                                                 >  2:273682/1‑100 (MQ=255)
    gCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCga                                                                                                <  1:2379028/100‑1 (MQ=255)
     cGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCGCCgaa                                                                                               >  1:1760314/1‑100 (MQ=255)
     cGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCgaa                                                                                               <  2:891425/100‑1 (MQ=255)
      ggCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCgaag                                                                                              >  1:1027973/1‑100 (MQ=255)
        cAAAGAAACCATGAAACTCGGCGAGTTGCTGGAAGCCCTCAAGCAAACATACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCgaagaa                                                                                            >  2:898189/1‑100 (MQ=255)
         aaaGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCACAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGaaa                                                                                           >  2:367333/1‑100 (MQ=255)
             aaaCCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACg                                                                                       <  2:390927/100‑1 (MQ=255)
                ccATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCGCCGAAGAAAAACGCTg                                                                                    <  2:767629/100‑1 (MQ=255)
                 cATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTgg                                                                                   >  1:416026/1‑100 (MQ=255)
                 cATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTgg                                                                                   >  2:880844/1‑100 (MQ=255)
                   tGAAACTCGGCGGGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGAt                                                                                 <  2:158323/100‑1 (MQ=255)
                    gAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATc                                                                                <  1:523579/100‑1 (MQ=255)
                          cGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGCTTGGTGCCGAGTATATGCACATTACCAGCGCCGAAGAAAAACGCTGGATCCAACAg                                                                          >  2:2176938/1‑100 (MQ=255)
                             cGAGCTGCTGGAAGCACTCAAGCAACCCCACAGCGGCCCGATTGATGCCGAGTATATGCACATTACCAGCACCGAAGAAACACGCTAGATCCAACAGCGt                                                                       <  2:705312/100‑1 (MQ=255)
                                gctgctGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATc                                                                    <  2:404186/100‑1 (MQ=255)
                                 ctgctgGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCGCCGAAGAAAAACGCTGGATCCAACAGCGTATCg                                                                   <  1:1892713/100‑1 (MQ=255)
                                 ctgctgGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTAGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTACCg                                                                   >  2:160347/1‑100 (MQ=255)
                                  tgctgGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCGCCGAAGAAAAACGCTGGATCCAACAGCGTATCGa                                                                  <  1:1386151/100‑1 (MQ=255)
                                  tgctgGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGa                                                                  >  1:2155788/1‑100 (MQ=255)
                                   gctgGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAg                                                                 >  2:2337607/1‑100 (MQ=255)
                                   gctgGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAg                                                                 >  2:1025942/1‑100 (MQ=255)
                                   gctgGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAg                                                                 >  1:865717/1‑100 (MQ=255)
                                   gctgGAAGCCCTCAAGCAAACCTACTGCGGCCCGATAGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAg                                                                 >  1:1689908/1‑100 (MQ=255)
                                    ctgGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGt                                                                <  1:1185276/100‑1 (MQ=255)
                                      ggAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCt                                                              <  2:1013812/100‑1 (MQ=255)
                                         aGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGt                                                           <  1:2223780/100‑1 (MQ=255)
                                          gCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTc                                                          >  1:910284/1‑100 (MQ=255)
                                            ccTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGCAGAAAAACGCTGGATCTAACAGCGCATCAAGTCTAGTagc                                                        >  1:760909/1‑100 (MQ=255)
                                            ccTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTcgc                                                        >  1:617737/1‑100 (MQ=255)
                                            ccTCAAGCAAACCTACTGCGGCCCGATAACTGCCGCGTATAAGCAAACTCCCTGCACCGCAGAAAAACGCTGGAACCAACAGCTTAGCGCGTATGGt                                                           >  2:1790493/1‑97 (MQ=255)
                                             cTCAAGCAAACCTACTGCGGCCCGATTGGTGGCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTcgcg                                                       <  1:1548777/100‑1 (MQ=255)
                                             cTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCGCCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTcgcg                                                       >  1:1769307/1‑100 (MQ=255)
                                             cTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTcgcg                                                       <  1:1548952/100‑1 (MQ=255)
                                               cAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCGCCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTcgcgcg                                                     <  2:1760293/100‑1 (MQ=255)
                                               cAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGACAAACGCTGGCTCCAACAGCGTATCGAGTCTGGTagcgcg                                                     >  2:446926/1‑100 (MQ=255)
                                               cAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGATGGATCCAACAGCGTATCGAGTCTGGTcgcgcg                                                     >  1:1196336/1‑100 (MQ=255)
                                                 aGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGAc                                                   >  1:1057170/1‑100 (MQ=255)
                                                   cAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCGCCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACtt                                                 <  1:1499008/100‑1 (MQ=255)
                                                      aCCTACTGCGGCCCGCTTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCGACAGCGTATCGAGTCTGGTCGCGCGACTTTCa                                              <  1:2281151/100‑1 (MQ=255)
                                                        cTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACTTTCAAt                                            <  1:1463137/100‑1 (MQ=255)
                                                           cTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACTTTCAATAGc                                         <  1:642680/100‑1 (MQ=255)
                                                           cTGCGGACCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTATGGTCGAGCGACTTTCAATa                                           >  2:1937381/1‑98 (MQ=255)
                                                            tGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACTTTCAATAGCg                                        >  2:738153/1‑100 (MQ=255)
                                                            tGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAAAAGCGTATCGAGTCTGGTCGCGCGACTTTCAATAGCg                                        >  2:984577/1‑100 (MQ=255)
                                                             gCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACTTTCAATAGCGa                                       <  1:2336110/100‑1 (MQ=255)
                                                               ggCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACTTTCAATAGCGAag                                     >  2:155592/1‑100 (MQ=255)
                                                                gCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACTTTCAATAGCGAaga                                    >  2:43838/1‑100 (MQ=255)
                                                                  ccGATTGGTGCCGAGTATATGCACATTACCATCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACTTTCAATAGCGAAgaga                                  <  1:432490/100‑1 (MQ=255)
                                                                   cGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACTTTCAATAGCGAAGAGaa                                 >  1:1245263/1‑100 (MQ=255)
                                                                     aTTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACTTTCAATAGCGAAGAGaaaa                               <  1:778123/100‑1 (MQ=255)
                                                                       tGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACTTTCAATAGCGAAGAGaaaaaa                             <  1:1473782/100‑1 (MQ=255)
                                                                          tGCCGAGTATATGCACATTACCAGCGCCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACTTTCAATAGCGAAGAGAAAAAACg                           >  1:596254/1‑99 (MQ=255)
                                                                            ccGAGTATCTGCACATTACCAGCACCGAAGAAAACCGCTGGCTCCACCAGCGTATCGAGTCTGGTCGCGCGACTTCTAATTGCGAAGTGAAAAAACGc                          <  2:1394366/98‑1 (MQ=255)
                                                                                gTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACTTTCAATAGCGAAGAGAAAAAACGCTTCTTa                    >  1:1827567/1‑100 (MQ=255)
                                                                                 tataTGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACTTTCAATAGCGAAGAGAAAAAACGCTTCTTa                    <  1:1102816/99‑1 (MQ=255)
                                                                                  atatGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACTTTCAATAGCGAAGAGAAAAAACGCTTCTTAAg                  >  1:1254200/1‑100 (MQ=255)
                                                                                   tatGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACTTTCAATAGCGAAGAGAAAAAACGCTTCTTAAGc                 <  1:288222/100‑1 (MQ=255)
                                                                                   tatGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACTTTCAATAGCGAAGAGAAAAAACGCTTCTTAAGc                 >  1:81951/1‑100 (MQ=255)
                                                                                   tatGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACTTTCAATAGCGAAGAGAAAAAACGCTTCTTAAGc                 <  1:1633800/100‑1 (MQ=255)
                                                                                     tGCACAAGACCAGAACAGAAGAAAAACGCTGGACACAACAGTGTATCGAGTCTGGTCGCGCGACTTTTAAAAAAAAAAAAAAAAAACGCTTCTTAAGCg                <  2:45217/99‑1 (MQ=255)
                                                                                      gCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACTTTCAATAGCGAAGAGAAAAAACGCTTCTTAAGCgga              >  1:1235448/1‑98 (MQ=255)
                                                                                       cacaTTACCAGCGCCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACTTTCAATAGCGAAGAGAAAAAACGCTTCTTAAGCGAAc             <  2:1207084/100‑1 (MQ=255)
                                                                                       cacaTTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACTTTCAATAGCGAAGAGAAAAAACGCTTCTTAAGCGAAc             >  1:2163277/1‑100 (MQ=255)
                                                                                         caTTACCAGCGCCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACTTTCAATAGCGAAGAGAAAAAACGCTTCTTAAGCGAACTg           >  2:414595/1‑100 (MQ=255)
                                                                                           ttACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACTTTCAATAGCGAAGAGAAAAAACGCTTCTTAAGCGAACTGAc         >  1:2249879/1‑100 (MQ=255)
                                                                                              ccAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACTTTCAATAGCGAAGAGAAAAAACGCTTCTTAAGCGAACTGAccgc      >  2:1345083/1‑100 (MQ=255)
                                                                                                 gCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACTTTCAATAGCGAAGAGAAAAAACGCTTCTTAAGCGAACTGACcgccgc   <  2:416023/100‑1 (MQ=255)
                                                                                                 gCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACTTTCAATAGCGAAGAGAAAAAACGCTTCTTAAGCGAACTGACcgccgc   <  1:271923/100‑1 (MQ=255)
                                                                                                 gCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACTTTCAATAGCGAAGAGAAAAAACGCTTCTTAAGCGAACTGACcgccgc   >  2:1933917/1‑100 (MQ=255)
                                                                                                  cACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACTTTCAATAGCGAAGAGAAAAAACGCTTCTTAAGCGAACTGAccgccg    >  2:1892263/1‑98 (MQ=255)
                                                                                                  cACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACTTTCAATAGCGAAGAGAAAAAACGCTTCTTAAGCGAACTGACCGCCGCt  <  1:1933939/100‑1 (MQ=255)
                                                                                                  cACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACTTTCAATAGCGAAGAGAAAAAACGCTTCTTAAGCGAACTGACCGCCGCt  >  1:1025812/1‑100 (MQ=255)
                                                                                                   |                                                                                                  
GCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAACAGCGTATCGAGTCTGGTCGCGCGACTTTCAATAGCGAAGAGAAAAAACGCTTCTTAAGCGAACTGACCGCCGCT  >  NC_000913/759138‑759335

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