Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	5,610,676	G→A	E14K (GAA→AAA)	emrE →	small multidrug resistance protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,610,676	0	G	A	100.0%	132.2 / NA	41	E14K (GAA→AAA)	emrE	small multidrug resistance protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (22/19);  total (22/19)

GTTCCACTGCTCGGCCATCAACCTATCTCCCATATCGATAGGTAATG‑ACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTAC  >  NC_002947/5610533‑5610816                                                                                                                                                 |                                                                                                                                             gTTCCACTGCTCGGCCATCAACCTATCTCCCATATCGATAGGTAATG‑ACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCa                                                                                                                                        >  2:386244/1‑150 (MQ=255)   tCCACTGCTCGGCCATCAACTTATCTCCAATGTCGAGACGAAATCAACTTTA‑AGGTTCCAGGCAACGTCTATAACCCCCGCTCAGTAACGTAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAACTCATc                                                                                                                                      <  2:31833/150‑1 (MQ=255)   cttcgatcTCGGCCATCAACCTATCTCCCATATCGATAGGTAATG‑ACTCTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGATTATCCTATCTGGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATc                                                                                                                                      <  1:630668/143‑1 (MQ=255)      aCTGCTCGGCCATCAACCTATCTCCCATATCGATAGGTAATG‑ACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGcc                                                                                                                                   >  2:87284/1‑150 (MQ=255)      aCTGCTCGGCCATCAACCTATCTCCCATATCGATAGGTAATG‑ACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGcc                                                                                                                                   >  2:747943/1‑150 (MQ=255)                            tcCCATATCGATAGGTAATG‑ACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGGAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGa                                                                                                             <  1:1096525/150‑1 (MQ=255)                               cATATCGATAGGTAATG‑ACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAg                                                                                                          <  1:747943/150‑1 (MQ=255)                                     gATAGGTAATG‑ACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGa                                                                                                    <  1:1135146/150‑1 (MQ=255)                                          gTAATG‑ACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAg                                                                                                                                           <  2:530619/106‑1 (MQ=255)                                          gTAATG‑ACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAg                                                                                                                                           >  1:530619/1‑106 (MQ=255)                                          gTAATG‑ACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACg                                                                                               >  2:1011542/1‑150 (MQ=255)                                                  cTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGc                                                                                        <  1:396297/150‑1 (MQ=255)                                                        cccATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCActg                                                                                   >  1:1064488/1‑149 (MQ=255)                                                        cccATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCActg                                                                                   <  2:1064488/149‑1 (MQ=255)                                                               aGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAt                                                                           >  1:778046/1‑150 (MQ=255)                                                                ggCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAtg                                                                          <  1:1011542/150‑1 (MQ=255)                                                                  cAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAgggt                                                                        >  1:551918/1‑150 (MQ=255)                                                                    aCGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAtggtgg                                                                      >  2:439718/1‑150 (MQ=255)                                                                    aCGTCTAGAATGCCCGCTCAGTAACGGAGGACCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAtggtgg                                                                      >  2:1118830/1‑150 (MQ=255)                                                                      gTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGtt                                                                    >  2:910272/1‑150 (MQ=255)                                                                         tAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGc                                                                                                        <  1:747517/111‑1 (MQ=255)                                                                         tAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGc                                                                                                        >  2:747517/1‑111 (MQ=255)                                                                             aTGCCCGCTCAGTAACGGATGATCCTGCCATGAATGCATATACCTACCTCGTCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATg                                                             <  1:263440/150‑1 (MQ=255)                                                                                cccGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCa                                                          >  2:328818/1‑150 (MQ=255)                                                                                  cGCTCAGTAACGGAGGACCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCAtt                                                        <  1:1118830/150‑1 (MQ=255)                                                                                             ggaggaTCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCAAAGCCACCGCATAGACCAACATCTCGTATGCCGTCTTCTGCTTGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCCTTGCGTTCTGGAt                                             >  1:273154/1‑150 (MQ=255)                                                                                                   tCCTGCCATGAATGCCTATACCTACCTCGCCACTTCGGTCTGCGCCAAAGTCATCGCCACCGCCTCGATGACGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGAc                                       >  1:994760/1‑150 (MQ=255)                                                                                                       gCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCtg                                   <  1:451112/150‑1 (MQ=255)                                                                                                         cATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCtggt                                 >  2:126449/1‑150 (MQ=255)                                                                                                          aTGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCtggtg                                >  1:151622/1‑150 (MQ=255)                                                                                                               tGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCg                           >  2:356418/1‑150 (MQ=255)                                                                                                                  cTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCgcag                        <  1:910272/150‑1 (MQ=255)                                                                                                                  cTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCgcag                        >  1:772921/1‑150 (MQ=255)                                                                                                                            ctCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGAt                                                                                                                        <  2:1031067/44‑1 (MQ=255)                                                                                                                            ctCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGAt                                                                                                                        >  1:1031067/1‑44 (MQ=255)                                                                                                                            ctCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATgg                                                            >  2:456681/1‑104 (MQ=255)                                                                                                                            ctCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATgg                                                            <  1:456681/104‑1 (MQ=255)                                                                                                                                 cATTGCCATTTGTGACAACGTCATCGCACCCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGAt         <  2:551918/150‑1 (MQ=255)                                                                                                                                  aTTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGAtt        >  2:857794/1‑150 (MQ=255)                                                                                                                                       cATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTa   <  1:386244/150‑1 (MQ=255)                                                                                                                                        aTCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCGTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTAc  <  1:130618/150‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                             GTTCCACTGCTCGGCCATCAACCTATCTCCCATATCGATAGGTAATG‑ACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTAC  >  NC_002947/5610533‑5610816

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: ‑