Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,161,702 C→G A162G (GCA→GGA)  btuB → BtuB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,161,7020CG100.0% 279.9 / NA 105A162G (GCA→GGA) btuBBtuB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (53/52);  total (53/52)

Only 100 of 108 total aligned reads displayed.

ATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGG  >  NC_000913/4161630‑4161773
                                                                        |                                                                       
aTGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCaggagg                                                                      >  1:2578548/1‑76 (MQ=255)
aTGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCaggagg                                                                      <  2:623966/76‑1 (MQ=255)
aTGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCaggagg                                                                      >  2:2272316/1‑76 (MQ=255)
 tGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAggg                                                                     >  1:3135020/1‑76 (MQ=255)
   gTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTg                                                                   <  1:2453530/76‑1 (MQ=255)
   gTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTg                                                                   <  2:2021860/76‑1 (MQ=255)
    ttCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTgg                                                                  >  2:3058219/1‑76 (MQ=255)
    ttCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTgg                                                                  >  1:1777835/1‑76 (MQ=255)
    gtCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTgg                                                                  <  1:483519/75‑1 (MQ=255)
      ccGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTgggg                                                                >  1:380318/1‑76 (MQ=255)
      ccGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTgggg                                                                >  2:2596709/1‑76 (MQ=255)
      ccGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTgggg                                                                <  1:2645727/76‑1 (MQ=255)
       cGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGa                                                               >  1:720763/1‑76 (MQ=255)
       cGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGa                                                               <  2:4106637/76‑1 (MQ=255)
       cGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGa                                                               <  1:2591402/76‑1 (MQ=255)
        gATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGaa                                                              <  1:468604/76‑1 (MQ=255)
        gATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGaa                                                              <  1:1827200/76‑1 (MQ=255)
         gtGCGATAGGCGGGGTGGTGAATATCAGCACGACGCGCGATGAGCCCGGAACGGAAATTTCAGGAGGGTGGGGacg                                                             <  2:1108063/75‑3 (MQ=255)
         gtGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAg                                                             <  2:4466511/75‑1 (MQ=255)
         aTGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAg                                                             <  2:2990924/76‑1 (MQ=255)
         aTGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAg                                                             <  1:4229219/76‑1 (MQ=255)
         aTGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAg                                                             <  2:4223863/76‑1 (MQ=255)
          tGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGc                                                            <  2:3570303/76‑1 (MQ=255)
           gCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCa                                                           >  1:1321625/1‑76 (MQ=255)
           gCAATAGGCGGGGCGGGGAATATCATCACCACGCGCGATCAACCCGGAACGGAAATTTATGGAAGGTGGCGAAGCa                                                           >  2:2396583/1‑76 (MQ=255)
              gtAGGCGGGGTGGGGAATATCATCACGGCGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATa                                                        <  2:1059328/75‑1 (MQ=255)
              aTAGGCGGGGGGGTGAATATCATCACGACGCGCGATGAACCCCGAACGGGAAATTCAGGGGGGGGGGGGAGCAATa                                                        >  1:1953307/1‑76 (MQ=255)
               tAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGACCGGAAATTTCAGGAGGGTGGGGAAGCAATAg                                                       >  2:3731143/1‑76 (MQ=255)
               tAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAg                                                       <  2:217856/76‑1 (MQ=255)
                 ggCGTGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCCGGAGGGTGGGGAAGCAATAGtt                                                     <  2:4140437/76‑1 (MQ=255)
                  gCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTa                                                    <  2:856558/76‑1 (MQ=255)
                   cGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATGTCAGGAGGGTGGGGAACCAATAGTTac                                                   >  2:4073305/1‑75 (MQ=255)
                     gggTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCa                                                 >  1:593489/1‑76 (MQ=255)
                      ggtggtGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAg                                                <  1:3092537/76‑1 (MQ=255)
                          gtgAATATCATCACGGCGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACt                                            <  2:2241480/76‑1 (MQ=255)
                          gtgAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACt                                            <  2:3987677/76‑1 (MQ=255)
                          gtgAATATCATCACGACGCGCGATGAACCCGGAACGGAAACTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACt                                            <  1:3801186/76‑1 (MQ=255)
                          gtgAATATCATACCGACGCGCGATGAACCCGGACCGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTCTCAGAACt                                            <  2:2701002/76‑1 (MQ=255)
                           tgAATATCATCACGACGCGCGATGAGCCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTa                                           <  1:1720996/76‑1 (MQ=255)
                             aaTATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTttg                                         >  2:2737605/1‑76 (MQ=255)
                             aaTATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTatg                                         <  2:4482596/76‑1 (MQ=255)
                                atcatcACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCCGAACTatgatg                                      >  2:56531/1‑76 (MQ=255)
                                atcatcACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTatgatg                                      >  1:1140803/1‑76 (MQ=255)
                                atcatcACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTatgatg                                      >  1:2672023/1‑76 (MQ=255)
                                 tcatcaCGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGt                                     <  1:551617/76‑1 (MQ=255)
                                 tcatcaCGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGt                                     <  2:1920548/76‑1 (MQ=255)
                                 tcatcaCGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGt                                     <  1:59449/76‑1 (MQ=255)
                                   atcaCGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGtct                                   <  1:410803/76‑1 (MQ=255)
                                    tcaCGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGtctc                                  <  2:233315/76‑1 (MQ=255)
                                       cgacgCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTAc                               >  2:4237730/1‑76 (MQ=255)
                                        gacgCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACg                              <  1:2647318/76‑1 (MQ=255)
                                        gacgCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACg                              >  1:989908/1‑76 (MQ=255)
                                        gacgCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACg                              <  1:237155/76‑1 (MQ=255)
                                         acgCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACgc                             >  2:3525915/1‑76 (MQ=255)
                                          cgcgcgATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTATCTACgaa                            >  1:1546795/1‑74 (MQ=255)
                                           gcgcgATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACgcag                           >  2:4024447/1‑76 (MQ=255)
                                           gcgcgATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACgcag                           >  1:1173263/1‑76 (MQ=255)
                                           gcgcgATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACgcag                           >  1:4306757/1‑76 (MQ=255)
                                            cgcgATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTatgat                                       >  1:3872934/1‑63 (MQ=255)
                                            cgcgATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATTATGTCTCTACgcagc                          <  2:75149/76‑1 (MQ=255)
                                            cgcgATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACgcagc                          <  2:1108290/76‑1 (MQ=255)
                                             gcgATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACgcagca                         >  2:4142184/1‑76 (MQ=255)
                                             gcgATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACgcagca                         >  1:2273742/1‑76 (MQ=255)
                                             gcgATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACgcagca                         >  1:3450358/1‑76 (MQ=255)
                                             gcgATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACgcagca                         >  2:2885350/1‑76 (MQ=255)
                                             gcgATGAACCCGCAACGGAAATTTCAGGTGGGAGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACgcagca                         >  2:798860/1‑76 (MQ=255)
                                              cgATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGcaa                        >  1:913844/1‑76 (MQ=255)
                                               taTGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGccac                       >  2:4174644/2‑76 (MQ=255)
                                               gATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGcaac                       >  1:828949/1‑76 (MQ=255)
                                               gATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGcaac                       >  1:3641559/1‑76 (MQ=255)
                                                aTGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGcaaca                      >  2:3033743/1‑76 (MQ=255)
                                                aTGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGcaaca                      >  1:3948878/1‑76 (MQ=255)
                                                 tGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGcaacca                     >  2:2075011/1‑74 (MQ=255)
                                                  gAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCaacaac                    <  2:3179801/76‑1 (MQ=255)
                                                  gAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCaacaac                    >  2:2227825/1‑76 (MQ=255)
                                                  gAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCaacaac                    >  2:2685838/1‑76 (MQ=255)
                                                    aCCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTg                  >  1:1962885/1‑76 (MQ=255)
                                                      ccGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTggg                >  1:2071335/1‑76 (MQ=255)
                                                      ccGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTggg                >  2:1101374/1‑76 (MQ=255)
                                                       cggaacggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTgggg               <  2:4192709/76‑1 (MQ=255)
                                                        ggaacggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTggggg              >  2:3022963/1‑76 (MQ=255)
                                                         gaacggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGa             <  1:617636/76‑1 (MQ=255)
                                                          aacggaaATTTCAGGAGGGTGGGTAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTgggggca            >  1:1383038/1‑74 (MQ=255)
                                                          aacggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGac            >  1:1410154/1‑75 (MQ=255)
                                                          aacggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGAt            <  1:2377100/76‑1 (MQ=255)
                                                          aacggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGAt            <  2:811895/76‑1 (MQ=255)
                                                          aacggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGAt            >  2:4248272/1‑76 (MQ=255)
                                                          aacggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGAt            <  2:469908/76‑1 (MQ=255)
                                                            cggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATaa          <  1:3964911/76‑1 (MQ=255)
                                                            cggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATaa          >  1:2403260/1‑76 (MQ=255)
                                                            cggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATaa          <  2:780530/76‑1 (MQ=255)
                                                            cggaaATTTCAGGAGGGGGGGGGAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATaa          <  2:77798/76‑1 (MQ=255)
                                                             ggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAg         >  1:1698569/1‑76 (MQ=255)
                                                              gaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGa        >  2:4381175/1‑76 (MQ=255)
                                                              gaaATATCAGGAGGGAGGGGACGCAATAGTTATCAGAACTATGCTGTCTCTCCGCAGCAACAGCTGGGGGATCAGa        <  1:798860/76‑1 (MQ=255)
                                                               aaaTTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGac       >  1:3513244/1‑76 (MQ=255)
                                                               aaaTTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGac       <  1:2377506/76‑1 (MQ=255)
                                                                 aTTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTCCGCAGCAACAACTGGGGGATAAGacac     <  1:2575487/76‑1 (MQ=255)
                                                                    tCAGGAGGGTGGGGAAGCATTAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACAGGGGGATAAGACACggg  >  2:4204505/1‑76 (MQ=255)
                                                                    tCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACggg  <  2:720763/76‑1 (MQ=255)
                                                                        |                                                                       
ATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGG  >  NC_000913/4161630‑4161773

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