Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,673,441	A→G	S21P (TCA→CCA)	mdtJ ←	multidrug efflux system transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,673,441	0	A	G	100.0%	240.5 / NA	68	S21P (TCA→CCA)	mdtJ	multidrug efflux system transporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (31/37);  total (31/37)

TCCCACAGCGCATAAGCTACGCCTAAGGCGATTTTTTTAACGGCGAAAGAGAGAAATATATACGACAGAGAAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCC  >  NC_000913/1673308‑1673575                                                                                                                                      |                                                                                                                                       tCCCACAGCGCATAAGCTACGCCTAAGGCGATTTTTTTAACGGCGAAAGAGAGAAATATATACGACAGAGAAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGGCAGCg                                                                                                                                   >  1:218757/1‑139 (MQ=255) tCCCACAGCGCATAAGCTACGCCTAAGGCGATTTTTTTAACGGCGAAAGAGAGAAATATATACGACAGAGAAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGGCAGCg                                                                                                                                   >  2:1244599/1‑139 (MQ=255) tCCCACAGCGCATAAGCTACGCCTAAGGCGATTTTTTTAACGGCGAAAGAGAGAAATATATACGACAGAGAAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGGCAGCg                                                                                                                                   <  1:1397931/139‑1 (MQ=255)  cccACAGCGCATAAGCTACGCCTAAGGCGATTTTTTTAACGGCGAAAGAGAGAAATATATACGACAGAGAAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGGCAGCGt                                                                                                                                  >  1:1404583/1‑139 (MQ=255)   ccACAGCGCATAAGCTACGCCTAAGGCGATTTTTTTAACGGCGAAAGAGAGAAATATATACGACAGAGAAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGGCAGCGTa                                                                                                                                 <  1:219526/139‑1 (MQ=255)     acaGCGCATAAGCTACGCCTAAGGCGATTTTTTTAACGGCGAAAGAGAGAAATATATACGACAGAGAAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGGCAGCGTAcc                                                                                                                               <  1:1244599/139‑1 (MQ=255)         cgcATAAGCTACGCCTAAGGCGATTTTTTTAACGGCGAAAGAGAGAAATATATACGACAGAGAAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGGCAGCGTACCGGTa                                                                                                                           >  1:438501/1‑139 (MQ=255)         cgcATAAGCTACGCCTAAGGCGATTTTTTTAACGGCGAAAGAGAGAAATATATACGACAGAGAAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGGCAGCGTACCGGTa                                                                                                                           >  1:798647/1‑139 (MQ=255)         cgcATAAGCTACGCCTAAGGCGATTTTTTTAACGGCGAAAGAGAGAAATATATACGACAGAGAAATCATCACCAGCATTAAAATAAAGCCGCCATTT‑CCTCACTGACGCTCGCCCATTTCATTGGCAGCGTACCGGTaa                                                                                                                          >  1:355175/1‑139 (MQ=255)          gcATAAGCTACGCCTAAGGCGATTTTTTTAACGGCGAAAGAGAGAAATATATACGACAGAGAAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGGCAGCGTACCGGTaa                                                                                                                          >  1:1127614/1‑139 (MQ=255)                gCTACGCCTAAGGCGATTTTTTTAACGGCGAAAGAGAGAAATATATACGACAGAGAAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGGCAGCGTACCGGTAATTTCTg                                                                                                                    <  2:438501/139‑1 (MQ=255)                          aGGCGATTTTTTTAACGGCGAAAGAGAGAAATATATACGACAGAGAAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGGCAGCGTACCGGTAATTTCTGTAGCAATAGc                                                                                                          <  2:218757/139‑1 (MQ=255)                           ggCGATTTTTTTAACGGCGAAAGAGAGAAATATATACGACAGAGAAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGGCAGCGTACCGGTAATTTCTGTAGCAATAGcc                                                                                                         <  2:1404583/139‑1 (MQ=255)                           ggCGATTTTTTTAACGGCGAAAGAGAGAAATATATACGACAGAGAAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGGCAGCGTACCGGTAATTTCTGTAGCAATAGcc                                                                                                         <  2:533200/139‑1 (MQ=255)                               aTTTTTTTAACGGCGAAAGAGAGAAATATATACGACAGAGAAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGGCAGCGTACCGGTAATTTCTGTAGCAATAGCCAGAc                                                                                                     >  2:598359/1‑139 (MQ=255)                                tttttttAACGGCGAAAGAGAGAAATATATACGACAGAGAAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGGCAGCGTACCGGTAATTTCTGTAGCAATAGCCAGAcc                                                                                                    <  2:798647/139‑1 (MQ=255)                                    tttAACGGCGAAAGAGAGAAATATATACGACAGAGAAATCATCACCAGCATTAAAATAAAGCCGCCATTT‑CCTCACTGACGCTCGCCCATTTCATTGGCAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCtaata                                                                                               <  2:355175/139‑1 (MQ=255)                                           gcaaaaGAGAGAAATATATACGACAGAGAAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGGCAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATcc                                                                                         <  2:1116367/136‑1 (MQ=255)                                           gCGAAAGAGAGAAATATATACGACAGAGAAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGGCAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATcc                                                                                         <  2:159750/139‑1 (MQ=255)                                             gAAAGAGAGAAATATATACGACAGAGAAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGGCAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCaa                                                                                       <  2:491055/139‑1 (MQ=255)                                                     gaAATATATACGACAGAGAAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGGCAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAAtata                                                                               <  1:598359/139‑1 (MQ=255)                                                       aaTATATACGACAGAGAAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGGCAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATTCAATAAATATACa                                                                             <  1:1080298/139‑1 (MQ=255)                                                                gACAGAGAAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGGCAGCGTACCGGTAATTTCTGTAGCaa                                                                                                              >  1:1268001/1‑97 (MQ=255)                                                                gACAGAGAAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGGCAGCGTACCGGTAATTTCTGTAGCaa                                                                                                              <  2:1268001/97‑1 (MQ=255)                                                                gACAGAGAAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGGCAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTtc                                                                    <  2:810578/139‑1 (MQ=255)                                                                  cAGAGAAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGGCAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGt                                                                         <  1:1552404/132‑1 (MQ=255)                                                                  cAGAGAAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGGCAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGt                                                                         >  2:1552404/1‑132 (MQ=255)                                                                   agagAAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGGCAGCGTACCGGTAATTTCTg                                                                                                                    <  1:1278869/88‑1 (MQ=255)                                                                   agagAAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGGCAGCGTACCGGTAATTTCTg                                                                                                                    >  2:1278869/1‑88 (MQ=255)                                                                   agagAAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGGCAGCGTACCGGTAATTTCTGTAGCa                                                                                                               <  1:1406756/93‑1 (MQ=255)                                                                   agagAAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGGCAGCGTACCGGTAATTTCTGTAGCa                                                                                                               >  2:1406756/1‑93 (MQ=255)                                                                   agagAAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGGCAGCGTACCGGTAATTTCTGTAGCAATAGCCAGAc                                                                                                     >  2:319057/1‑103 (MQ=255)                                                                   agagAAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGGCAGCGTACCGGTAATTTCTGTAGCAATAGCCAGAc                                                                                                     <  1:319057/103‑1 (MQ=255)                                                                     agaAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGGCAGCGTa                                                                                                                                 >  2:864938/1‑73 (MQ=255)                                                                     agaAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGGCAGCGTa                                                                                                                                 <  1:864938/73‑1 (MQ=255)                                                                     agaAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGGCAGCGTACCGGTAATTTCTGTAGCAATAGCCAGa                                                                                                      >  1:807371/1‑100 (MQ=255)                                                                     agaAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGGCAGCGTACCGGTAATTTCTGTAGCAATAGCCAGa                                                                                                      <  2:807371/100‑1 (MQ=255)                                                                     agaAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGGCAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCaa                                                                                       <  2:799066/115‑1 (MQ=255)                                                                     agaAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGGCAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCaa                                                                                       >  1:799066/1‑115 (MQ=255)                                                                     agaAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGGCAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACAt                                                                            >  1:691197/1‑126 (MQ=255)                                                                     agaAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGGCAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACAt                                                                            <  2:691197/126‑1 (MQ=255)                                                                     agaAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGGCAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTc                                                                        >  2:1374604/1‑130 (MQ=255)                                                                     agaAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGGCAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTc                                                                        <  1:1374604/130‑1 (MQ=255)                                                                     agaAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGGCAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTg                                                               >  1:1589291/1‑139 (MQ=255)                                                                            atcaCCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGGCAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAgaga                                                        >  1:419767/1‑139 (MQ=255)                                                                               aCCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGGCAGCGTACCGGTAATTTCTGTAGCAATAGcc                                                                                                         <  2:469521/87‑1 (MQ=255)                                                                               aCCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGGCAGCGTACCGGTAATTTCTGTAGCAATAGcc                                                                                                         >  1:469521/1‑87 (MQ=255)                                                                                 cAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGGCAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACa                                                                             >  2:482522/1‑113 (MQ=255)                                                                                 cAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGGCAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACa                                                                             <  1:482522/113‑1 (MQ=255)                                                                                      ttAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGGCAGCGTACCGGTAATTTCTGTAGc                                                                                                                >  2:685340/1‑73 (MQ=255)                                                                                      ttAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGGCAGCGTACCGGTAATTTCTGTAGc                                                                                                                <  1:685340/73‑1 (MQ=255)                                                                                      ttAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGGCAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACAt                                                                            <  1:139540/109‑1 (MQ=255)                                                                                      ttAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGGCAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACAt                                                                            >  2:139540/1‑109 (MQ=255)                                                                                              aaGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGGCAGCGTACCGGTAAt                                                                                                                         <  1:1012638/56‑1 (MQ=255)                                                                                              aaGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGGCAGCGTACCGGTAAt                                                                                                                         >  2:1012638/1‑56 (MQ=255)                                                                                              aaGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGGCAGCGTACCGGTAATTTCTGTAGc                                                                                                                <  1:41013/65‑1 (MQ=255)                                                                                              aaGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGGCAGCGTACCGGTAATTTCTGTAGc                                                                                                                >  2:41013/1‑65 (MQ=255)                                                                                              aaGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGGCAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCaa                                                                                       <  2:720200/90‑1 (MQ=255)                                                                                              aaGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGGCAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCaa                                                                                       >  1:720200/1‑90 (MQ=255)                                                                                              aaGCCGCCATCTCCCTCACTGACGCTCGCCCATTTCATTGGCAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCt                                      >  1:532379/1‑139 (MQ=255)                                                                                                    ccATTTCCCTCACTGACGCTCGCCCATTTCATTGGCAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTc                                >  2:660651/1‑139 (MQ=255)                                                                                                           ccTCACTGACGCTCGCCCATTTCATTGGCAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGaaa                         <  1:660651/139‑1 (MQ=255)                                                                                                                      cTCGCCCATTTCATTGGCAGCGTACCGGTAATTTCTGTAGCAATAGcc                                                                                                         <  1:1442752/48‑1 (MQ=255)                                                                                                                      cTCGCCCATTTCATTGGCAGCGTACCGGTAATTTCTGTAGCAATAGcc                                                                                                         >  2:1442752/1‑48 (MQ=255)                                                                                                                         gCCCATTTCATTGGCAGCGTACCGGTAATTTCTGTAGCAATAGc                                                                                                          >  1:1267337/1‑44 (MQ=255)                                                                                                                         gCCCATTTCATTGGCAGCGTACCGGTAATTTCTGTAGCAATAGc                                                                                                          <  2:1267337/44‑1 (MQ=255)                                                                                                                                  aTTGGCAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGcccc  <  2:532379/139‑1 (MQ=255)                                                                                                                                  aTTGGCAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGcccc  <  2:419767/139‑1 (MQ=255)                                                                                                                                      |                                                                                                                                       TCCCACAGCGCATAAGCTACGCCTAAGGCGATTTTTTTAACGGCGAAAGAGAGAAATATATACGACAGAGAAATCATCACCAGCATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCC  >  NC_000913/1673308‑1673575

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