Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,951,168 +G coding (37/1818 nt) ilvD → dihydroxyacid dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,951,1681.G96.0% 192.0 / 0.5 76coding (37/1818 nt)ilvDdihydroxyacid dehydratase
Reads supporting (aligned to +/- strand):  ref base . (3/0);  new base G (24/49);  total (27/49)
Fisher's exact test for biased strand distribution p-value = 4.16e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGT‑‑GG‑‑‑GCC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCC  >  NC_000913/3951103‑3951231
                                                                        |                                                                 
gCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGT‑‑GGC‑‑GCGC‑‑CACCGGa                                                         <  2:2487804/76‑1 (MQ=255)
gCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGT‑‑GGC‑‑GCGC‑‑CACCGGa                                                         <  2:3088888/76‑1 (MQ=255)
   tAAGTACCGTCCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGT‑‑GGC‑‑GCGC‑‑CACCGGAATg                                                      <  1:278772/76‑1 (MQ=255)
      gtcccGTTCGGCCCCCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGT‑‑GGC‑‑GCGC‑‑CACCGGAATGAcc                                                   <  2:439035/73‑1 (MQ=255)
      gtcccGTGCCGCCCCCACCACTCAGGGTCGTAATATGGCGGGTGCTCGTGCGCTGT‑‑GGC‑‑GCGC‑‑CACCGGAATGAcc                                                   <  2:425527/73‑1 (MQ=255)
       gcccGTCCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGT‑‑GGC‑‑GCGC‑‑CACCGGAATGAcca                                                  <  2:4702002/74‑2 (MQ=255)
        cccGTTCCGCCCCCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGT‑‑GGC‑‑GCGC‑‑CACCGGAATGACCGa                                                 <  1:364163/75‑1 (MQ=255)
         ccGTCCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGT‑‑GGC‑‑GCGC‑‑CACCGGAATGACCGAc                                                <  1:2225700/76‑1 (MQ=255)
         ccGTCCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGT‑‑GGC‑‑GCGC‑‑CACCGGAATGACCGAc                                                <  1:2620511/76‑1 (MQ=255)
          cGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGT‑‑GGC‑‑GCGC‑‑CACCGGAATGACCGACg                                               <  1:1446818/76‑1 (MQ=255)
          cGTCCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGT‑‑GGC‑‑GCGC‑‑CACCGGAATGACCGACg                                               <  2:2228722/76‑1 (MQ=255)
           gTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGCTCGTGGGCCGT‑‑GGC‑‑GCGC‑‑CCCCCCAATGACCGACGc                                              >  1:1145169/1‑76 (MQ=255)
              ccGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGT‑‑GGC‑‑GCGC‑‑CACCGGAATGACCGACGCCGa                                           <  1:2691822/76‑1 (MQ=255)
               cGCCACCACCACTCATGGTCGTAATATGGCGGGGGCTCGTGCGCTGTGGG‑‑‑GCGC‑‑CCCCGGGATGACCGACGCCGAt                                          >  1:4225566/1‑76 (MQ=255)
               cGCCACCACCACTCATGGTCGTAATATGGCGGGGGCTCGTGCGCTGTGGG‑‑‑GCGC‑‑CCCCGGGATGACCGACGCCGAt                                          >  1:3388269/1‑76 (MQ=255)
                gCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGT‑‑GGC‑‑GCGC‑‑CACCGGAATGACCGACGCCGAtt                                         <  1:1390018/76‑1 (MQ=255)
                  caccaccacTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGT‑‑GGC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTc                                       <  2:660418/76‑1 (MQ=255)
                  caccaccacTCATGGTCGTAATATGGCGGGGGGTCGTGCGCTGT‑‑GGC‑‑GCGC‑‑CCCCGGGATGACCGACGCCGCtttt                                       >  1:3097094/1‑75 (MQ=255)
                   cccaccacTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGT‑‑GGC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCg                                      <  2:918951/75‑1 (MQ=255)
                   cccaccacGAATGGTCGTAATATGGCTGGTGCTCGTGCGCTGT‑‑GGC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCg                                      <  1:379762/75‑1 (MQ=255)
                   accaccacTCATGGTCGTAATATGGCGGGGGGTCGTGCGCTGT‑‑GGCGCGCC‑‑CCCCGGAAGGACGACGCCGCTTTCg                                      >  1:1237892/1‑76 (MQ=255)
                    ccaccacTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGT‑‑GGC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCgg                                     <  2:2155821/76‑1 (MQ=255)
                     caccacTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGT‑‑GGC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGt                                    <  1:1852386/76‑1 (MQ=255)
                      accacTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGT‑‑GGC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTa                                   <  1:1102812/76‑1 (MQ=255)
                        cacTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGT‑‑GGC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAg                                 <  1:1168097/76‑1 (MQ=255)
                        cacTCATGGTCGTAATATGGCGGGGGGTCGTGCGGTGT‑‑GG‑‑‑GGCGGCCCCGGGATGGCCGGCGCCGATTTCGGGAAg                                 >  1:3042391/1‑76 (MQ=255)
                         ccTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGT‑‑GGC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGc                                <  1:206583/75‑1 (MQ=255)
                         acTCATGGTCGTAATATGGCGGGGGGTCGTGCGCTGTGGGG‑‑‑GGC‑‑CCCCCGGATGACCGACGCCGGTTTTGGTAAGc                                >  1:4497931/1‑76 (MQ=255)
                          cTCCTGGGCGTGATATGGCGGGTGCTCGTGCGCTGT‑‑GGC‑‑GCGC‑‑CAGCGGAATGCCCGACGCCGATTTCGGTAAGcc                               <  2:4497983/76‑1 (MQ=255)
                          cTCATGGTCGTAATATGGCGGGGGCTCGTGCGCTGT‑‑GGC‑‑GCGC‑‑CCCCCGAAAGACCGCCCCCCCTTTCGGTACGcc                               >  2:962362/1‑76 (MQ=255)
                           tCATGGTCGTAATATGGCGGGGGGTCGGGGGCTGTGGG‑‑‑GCGC‑‑CCCCGGGATGACCGACGCCGGTTTCTGTTAGGCg                              >  1:1388271/1‑76 (MQ=255)
                             aTGGTCGTAATCTGGCGGGTGCTCGTGCGCTGT‑‑GGC‑‑GCGC‑‑CGCCGGAATGACCGACGCCGATTTCGGTAAGCCGAt                            <  1:1999049/76‑1 (MQ=255)
                              tGGTCGTAATATGGCGGGGGCTCGTGCGCTGT‑‑GGC‑‑GCGC‑‑CACCCGAATGAACGACCCCGACTTCTGTAAGCCGGtt                           >  1:1488430/1‑76 (MQ=255)
                               ggTCGTAATATGGCGGGTGCTCGTGCGCTGT‑‑GGC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTa                          <  2:4601387/76‑1 (MQ=255)
                                gTCGTAATATGGCGGGGGCTCGTGCGCTGT‑‑GGC‑‑GCGC‑‑CCCCGGGATGACCGACGCCGCTTTCGGTGAGCCGATTAt                         >  1:1364857/1‑76 (MQ=255)
                                    tAATATGGCGGGTGCTCGTGCGCTGT‑‑GGC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATcgcg                     <  2:3266334/76‑1 (MQ=255)
                                    tAATATGGCGGGTGCTCGTGCGCTGT‑‑GGC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATcgcg                     <  1:4623945/76‑1 (MQ=255)
                                     aaTATGGCGGGTGCTCGTGCGCTGT‑‑GGC‑‑GCGC‑‑CCCCCGAATGACCGACGCCGCTTTTTGTAAGACCAATATCGCgg                    >  1:1292451/1‑76 (MQ=255)
                                      atatGGCGGGTGCTCGTGCGCTGT‑‑GGC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGt                   <  1:3013357/76‑1 (MQ=255)
                                      atatGGCGGGTGCTCGTGCGCTGT‑‑GGC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGt                   >  2:2229924/1‑76 (MQ=255)
                                      atatGGCGGGTGCTCGTGCGCTGT‑‑GGC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGt                   <  1:3504359/76‑1 (MQ=255)
                                       tatGGCGGGTGCTCGTGCGCTGT‑‑GGC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGtt                  <  1:2572252/76‑1 (MQ=255)
                                       tatGGCGGGTGCGCGTGCGCTGT‑‑GGC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGtt                  <  2:2941710/76‑1 (MQ=255)
                                        atGGCGGGTGCTCGTGCGCTGT‑‑GGC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTtg                 <  2:2456225/76‑1 (MQ=255)
                                        atGGCGGGTGCTCGTGCGCTGT‑‑GGC‑‑GCGC‑‑CACCGGAAAGCACGACCCCGACTTCGGGAAGCTGGTTATCGCGGtct                 >  2:2223444/1‑74 (MQ=255)
                                         tGGCGGGTGCTCGTGCGCTGT‑‑GGC‑‑GCGC‑‑CACCGGAATGACCGCCCCCGATTTCGGTAAGCCGATTATCGCGGtcgg                >  2:4832029/1‑73 (MQ=255)
                                         tGGCGGGTGCTCGTGCGCTGT‑‑GGC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTtgt                <  1:1490144/76‑1 (MQ=255)
                                         tGGCGGGTGCTCGTGCGCTGT‑‑GGC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTtgt                <  2:689925/76‑1 (MQ=255)
                                          cgCGTGTGCTCGTGCGCTGT‑‑GTC‑‑‑GCGC‑‑CACGGGAACGACCGTCGCCGATTTCGGTAAGCCGATTATCGCGGGtgtg               <  1:3831356/75‑1 (MQ=255)
                                           gCGGGTGCTCGTGCGCTGT‑‑GGC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGa              <  1:1864660/76‑1 (MQ=255)
                                              ggTGCTCGTGCGCTGT‑‑GGC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACt           <  2:2113660/76‑1 (MQ=255)
                                              ggTGCTCGTGCGCTGT‑‑GGC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACt           <  2:3709500/76‑1 (MQ=255)
                                              ggTGCTCGTGCGCTGT‑‑GGC‑‑GCGC‑‑AACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACt           <  2:1978114/76‑1 (MQ=255)
                                               gTGCTCGTGCGCTGT‑‑GGC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTc          >  2:2214138/1‑76 (MQ=255)
                                               gTGCTCGTGCGCTGT‑‑GGC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTc          >  1:3556354/1‑76 (MQ=255)
                                               gTGCTCGTGCGCTGT‑‑GGC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTc          <  1:2906302/76‑1 (MQ=255)
                                               gTGCTCGTGCGCTGT‑‑GGC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTc          <  1:586044/76‑1 (MQ=255)
                                               gTGCTCGTGCGCTGT‑‑GGC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTc          <  1:1662362/76‑1 (MQ=255)
                                               gTGCTCGTGCGCTGT‑‑GGC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTc          <  2:2615073/76‑1 (MQ=255)
                                               gTGCTCGTGCGCTGT‑‑GGC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTc          <  1:1845784/76‑1 (MQ=255)
                                               gTGCTCGTGCGCTGT‑‑GGC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTc          >  2:1463843/1‑76 (MQ=255)
                                                tGCTCGTGCGCTGT‑‑GGC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCg         <  1:4029669/76‑1 (MQ=255)
                                                 gCTCGTGCGCTGT‑‑GGC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGt        >  2:3501717/1‑76 (MQ=255)
                                                 gCTCGTGCGCTGT‑‑GGC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGt        <  1:4591927/76‑1 (MQ=255)
                                                  cTCGTGCGCTGG‑‑GGC‑‑GCGC‑‑CACCGGAATGACCGACGCCGACTTCGGTAAGCCGATTATCGCGACTGTGGACTCgct       >  1:289112/1‑74 (MQ=255)
                                                   tCGTGCGCTGT‑‑GGC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTc      <  2:4190382/76‑1 (MQ=255)
                                                   tCGTGCGCTGT‑‑GGC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTc      >  1:2691480/1‑76 (MQ=255)
                                                   tCGTGCGCTGT‑‑GGC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTc      >  1:4015101/1‑76 (MQ=255)
                                                   gcGTGCGCTGT‑‑GGC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTc      <  2:3349316/75‑1 (MQ=255)
                                                     gTGCGCTGT‑‑GGC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCAc    >  2:656481/1‑76 (MQ=255)
                                                     gTGCGCTGT‑‑GGC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCAc    >  1:3459815/1‑76 (MQ=255)
                                                     gTGCGCTGT‑‑GGC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATGGCGGTAAGCCGATTATGGAGGTGGTGAACTCTGTCAc    >  2:1232035/1‑76 (MQ=255)
                                                      tGCGCTGT‑‑GGC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCAcc   <  1:1327567/76‑1 (MQ=255)
                                                      tGCGCTGT‑‑GGC‑‑GCGC‑‑CACCGGAAGGACCGACGCCGATTTCGGTAAGCCGCTTATCGCGCTCGTGAACTCCTTcccc   >  2:1436147/1‑76 (MQ=255)
                                                       gcgcTGT‑‑GGC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCAccc  >  1:3121807/1‑76 (MQ=255)
                                                       gcgcTGT‑‑GGC‑‑GCGC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCAccc  <  1:1176353/76‑1 (MQ=255)
                                                                        |                                                                 
GCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGT‑‑GG‑‑‑GCC‑‑CACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCC  >  NC_000913/3951103‑3951231

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: