Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 469,429 A→G K187E (AAG→GAG)  mdlA → putative multidrug ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913469,4290AG100.0% 247.7 / NA 68K187E (AAG→GAG) mdlAputative multidrug ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (35/33);  total (35/33)

CGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGC  >  NC_000913/469308‑469563
                                                                                                                         |                                                                                                                                      
cgcTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGgcggcg                                                                                                                       <  2:1432490/139‑1 (MQ=255)
     tgtTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCTGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTc                                                                                                                  <  1:1007383/139‑1 (MQ=255)
               gatgTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAAt                                                                                                        >  2:241610/1‑139 (MQ=255)
                     tACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGc                                                                                                  >  1:818793/1‑139 (MQ=255)
                     tACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGc                                                                                                  >  1:390754/1‑139 (MQ=255)
                       cGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCAc                                                                                                >  1:307097/1‑139 (MQ=255)
                       cGCAAATTAGCTGGCAGTTGACCGTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCAc                                                                                                <  1:1067596/139‑1 (MQ=255)
                            aTTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAgg                                                                                           <  2:390754/139‑1 (MQ=255)
                             ttAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGgcggcg                                                                                                                       <  2:1481162/110‑1 (MQ=255)
                             ttAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGgcggcg                                                                                                                       >  1:1481162/1‑110 (MQ=255)
                               aGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGaaa                                                                                        >  1:42336/1‑139 (MQ=255)
                                    gCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTc                                                                                   >  2:1172803/1‑139 (MQ=255)
                                      aGTTGACCTTATTTTCCCAGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGTAAGCCTCAc                                                                                 <  1:1012845/139‑1 (MQ=255)
                                        ttGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCa                                                                               >  1:560111/1‑139 (MQ=255)
                                          gACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGcc                                                                                     >  2:1305533/1‑131 (MQ=255)
                                          gACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGcc                                                                                     <  1:1305533/131‑1 (MQ=255)
                                          gACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGt                                                                             <  2:307097/139‑1 (MQ=255)
                                           aCCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTa                                                                            <  1:241610/139‑1 (MQ=255)
                                               tATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGAc                                                                                                     >  2:373617/1‑110 (MQ=255)
                                               tATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGAc                                                                                                     <  1:373617/110‑1 (MQ=255)
                                                aTTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGc                                                                       <  2:818793/139‑1 (MQ=255)
                                                      ccTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCt                                                                                    >  2:1053702/1‑120 (MQ=255)
                                                      ccTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCt                                                                                    <  1:1053702/120‑1 (MQ=255)
                                                        tGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCaa                                                               <  2:42336/139‑1 (MQ=255)
                                                                       tGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTa                                                                            >  1:362299/1‑111 (MQ=255)
                                                                       tGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTa                                                                            <  2:362299/111‑1 (MQ=255)
                                                                            gCGATCATGATTAAGCGCAACGGCGGTGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGc                                                                                      >  2:1187238/1‑96 (MQ=255)
                                                                            gCGATCATGATTAAGCGCAACGGCGGTGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGc                                                                                      <  1:1187238/96‑1 (MQ=255)
                                                                                 cATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCg                                                                                                   <  2:906281/78‑1 (MQ=255)
                                                                                 cATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCg                                                                                                   >  1:906281/1‑78 (MQ=255)
                                                                                 cATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGcc                                                                                     <  1:442941/92‑1 (MQ=255)
                                                                                 cATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGcc                                                                                     >  2:442941/1‑92 (MQ=255)
                                                                                         aGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGTCGTTTTCCAGTCTTACTGACCGCACCCAGGAAAGCCTCATCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGtt                              >  1:960970/1‑139 (MQ=255)
                                                                                         aGCGCAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGtt                              >  1:1252653/1‑139 (MQ=255)
                                                                                           cgcAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCAc                                                                                                <  2:1568/71‑1 (MQ=255)
                                                                                           cgcAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCAc                                                                                                >  1:1568/1‑71 (MQ=255)
                                                                                           cgcAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGcc                                                                                     >  1:1379483/1‑82 (MQ=255)
                                                                                           cgcAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGcc                                                                                     <  2:1379483/82‑1 (MQ=255)
                                                                                           cgcAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCAc                                                                                 >  2:1395237/1‑86 (MQ=255)
                                                                                           cgcAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCAc                                                                                 <  1:1395237/86‑1 (MQ=255)
                                                                                           cgcAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCt                                                   >  1:1304806/1‑116 (MQ=255)
                                                                                           cgcAACGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCt                                                   <  2:1304806/116‑1 (MQ=255)
                                                                                              aaCGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTaa                                                                                                         >  2:1128369/1‑59 (MQ=255)
                                                                                              aaCGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTaa                                                                                                         <  1:1128369/59‑1 (MQ=255)
                                                                                              aaCGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGaa                                                                                         >  1:575392/1‑75 (MQ=255)
                                                                                              aaCGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGaa                                                                                         <  2:575392/75‑1 (MQ=255)
                                                                                              aaCGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGc                                                                                      <  1:740843/78‑1 (MQ=255)
                                                                                              aaCGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGc                                                                                      >  2:740843/1‑78 (MQ=255)
                                                                                              aaCGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCttt                                                        <  1:1132105/108‑1 (MQ=255)
                                                                                              aaCGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCttt                                                        >  2:1132105/1‑108 (MQ=255)
                                                                                              aaCGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTg                         >  1:829124/1‑139 (MQ=255)
                                                                                              aaCGGCGATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTg                         >  1:989139/1‑139 (MQ=255)
                                                                                                    gATGCTTTGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCgg                   <  2:560111/139‑1 (MQ=255)
                                                                                                           tGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGc                                                                                      >  2:357343/1‑65 (MQ=255)
                                                                                                           tGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGc                                                                                      <  1:357343/65‑1 (MQ=255)
                                                                                                           tGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCcgcg                    >  1:308904/1‑131 (MQ=255)
                                                                                                           tGCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCcgcg                    <  2:308904/131‑1 (MQ=255)
                                                                                                            gCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGt                                     >  2:1154200/1‑113 (MQ=255)
                                                                                                            gCATGAACGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGt                                     <  1:1154200/113‑1 (MQ=255)
                                                                                                                  aCGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGAcc                                                                                                    <  2:903667/44‑1 (MQ=255)
                                                                                                                  aCGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGAcc                                                                                                    >  1:903667/1‑44 (MQ=255)
                                                                                                                  aCGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAgg                                                                                           >  1:377064/1‑53 (MQ=255)
                                                                                                                  aCGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAgg                                                                                           <  2:377064/53‑1 (MQ=255)
                                                                                                                  aCGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGaa                                                                                         >  2:167457/1‑55 (MQ=255)
                                                                                                                  aCGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGaa                                                                                         <  1:167457/55‑1 (MQ=255)
                                                                                                                   cGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCAc                                                                                                <  2:830469/47‑1 (MQ=255)
                                                                                                                   cGCTTTGAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCAc                                                                                                >  1:830469/1‑47 (MQ=255)
                                                                                                                     cTTTGAGCTGGCACAGGCGCCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGTCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGc  >  2:747443/1‑139 (MQ=255)
                                                                                                                         |                                                                                                                                      
CGCTGTGTTGATTATGATGTCTACGCAAATTAGCTGGCAGTTGACCTTATTTTCCCTGTTGCCGATGCCAGTGATGGCGATCATGATTAAGCGCAACGGCGATGCTTTGCATGAACGCTTTAAGCTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGC  >  NC_000913/469308‑469563

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