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%	123.7 / NA	37	E14K (GAA→AAA)	emrE	small multidrug resistance protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (25/12);  total (25/12)

ATCTCCCATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATC  >  NC_002947/5610557‑5610822                                                                                                                        |                                                                                                                                                   aTCTCCCATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCgg                                                                                                                      <  2:161987/150‑1 (MQ=255)    tcCCATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGa                                                                                                                   >  2:113818/1‑150 (MQ=255)        ttatCGATAGGTAATGACTTTACCCCATCCAGGCAGCGTCTAGAATGCCTGCTCAGTAACGGAGGATCCTGCCATGAATGCATATCTCTACCTCGCCGTCGCCCTCTGTGCCAAAGTCATCGCCACCGCCTCTATGAAGGCGGGGAGAgg                                                                                                               >  1:96096/2‑150 (MQ=255)             gATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAg                                                                                                                          <  1:415831/134‑1 (MQ=255)             gATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAg                                                                                                                          >  2:415831/1‑134 (MQ=255)              aTAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAg                                                                                                         <  1:113818/150‑1 (MQ=255)                   tAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGATGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGc                                                                                                    >  2:243859/1‑150 (MQ=255)                          tttACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCt                                                                                                            <  2:554613/135‑1 (MQ=255)                          tttACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCt                                                                                                            >  1:554613/1‑135 (MQ=255)                               cccATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGc                                                                                                    <  2:561660/138‑1 (MQ=255)                               cccATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGc                                                                                                    >  1:561660/1‑138 (MQ=255)                                             gTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGtt                                                                          >  1:39997/1‑150 (MQ=255)                                              tCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTc                                                                                                                                 >  2:88066/1‑94 (MQ=255)                                              tCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTc                                                                                                                                 <  1:88066/94‑1 (MQ=255)                                                 aGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCt                                                                      <  1:529402/150‑1 (MQ=255)                                                        ccGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCAACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCAt                                                               >  1:602243/1‑150 (MQ=255)                                                            tCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCg                                                           >  1:377545/1‑150 (MQ=255)                                                              aGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGtt                                                         >  1:70962/1‑150 (MQ=255)                                                                  aCGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTgg                                                     >  2:604524/1‑150 (MQ=255)                                                                    cgattaTCCTGCAATGAATGCCTATACCTACCTCTCTATTGCCATCTGCGCCAAATTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGTCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTTTGGCATTGCGTTCTGGAt                                                   <  2:96096/145‑1 (MQ=255)                                                                                 aTGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCtggtg                                      >  2:247279/1‑150 (MQ=255)                                                                                          tataCCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCgcagc                             >  1:140727/1‑150 (MQ=255)                                                                                               ctacctCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCtg                                         >  2:16939/1‑133 (MQ=255)                                                                                               ctacctCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCtg                                         <  1:16939/133‑1 (MQ=255)                                                                                                 acctCGCCATTGCCATATGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCgg                      >  2:212824/1‑150 (MQ=255)                                                                                                   ctCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTg                    >  1:287037/1‑150 (MQ=255)                                                                                                    tCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTgg                   >  1:47048/1‑150 (MQ=255)                                                                                                    tCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTgg                   >  2:588901/1‑150 (MQ=255)                                                                                                        cATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGaaa                                                                                                                 <  2:230066/52‑1 (MQ=255)                                                                                                        cATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGaaa                                                                                                                 >  1:230066/1‑52 (MQ=255)                                                                                                             ccATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCg                       >  2:486345/1‑137 (MQ=255)                                                                                                             ccATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCg                       <  1:486345/137‑1 (MQ=255)                                                                                                              cATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTa         >  2:598566/1‑150 (MQ=255)                                                                                                                  tGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGcc     >  1:203049/1‑150 (MQ=255)                                                                                                                   gcgcCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAAGGTGGGTGGCGAGGACATGGCGTGCTGGATGCTGACACTGGGGTGCCGCAGGATTCCTGTTGGGATTGCCTACGCCa    >  2:597194/1‑150 (MQ=255)                                                                                                                    cgcCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCAt   <  1:598566/150‑1 (MQ=255)                                                                                                                     gcCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATc  >  1:228953/1‑150 (MQ=255)                                                                                                                        |                                                                                                                                                   ATCTCCCATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATC  >  NC_002947/5610557‑5610822

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