Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,299,682	A→G	V450A (GTG→GCG)	mdtP ←	putative multidrug efflux pump outer membrane channel

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,299,682	0	A	G	100.0%	127.9 / NA	35	V450A (GTG→GCG)	mdtP	putative multidrug efflux pump outer membrane channel
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (21/14);  total (21/14)

ACGCGGCAGACATTTTATTTTTTCTCGACGACGGGACCTGCCTGATACCCGCCGCCCAGCGATTTCATCAACTGAATGCTTTGGATCACCCGGCGGCTGTCCAGCATCAGTAATGACATCTCTTCGGCAAGCACTGGCAACCGGGCTTCGGTGGCCTGTAAGCGGCTGGTTAAGCCGCGCTGATAGGCGGCCTCGGCAGCGCGCTGGGTAAAGCGCGTGGCTTCCACGCGTTCAGCCTGCATTTCTCGCTCGTCGTTGAGCGT  >  NC_000913/4299550‑4299812                                                                                                                                     |                                                                                                                                   aCGCGGCAGACATTTTATTTTTTCTCGACGACGGGACCTGCCTGATACCCGCCGCCCAGCGATTTCATCAACTGAATGCTTTGGATCACCCGGCGGCTGTCCAGCATCAGTAATGACATCTCTTCGGCAAGCGCTGGCaa                                                                                                                             <  1:292659/140‑1 (MQ=255) aCGCGGCAGACATTTTATTTTTTCTCGACGACGGGACCTGCCTGATACCCGCCGCCCAGCGATTTCATCAACTGAATGCTTTGGATCACCCGGCGGCTGTCCAGCATCAGTAATGACATCTCTTCGGCAAGCGCTGGCaa                                                                                                                             <  2:180836/140‑1 (MQ=255)          aCATTTTATTTTTTCTCGACGACGGGACCTGCCTGATACCCGCCGCCCAGCGATTTCATCAACTGAATGCTGTGGATCACCCGGCGGCTGTCCAGCATCAGTAATGACATCTCTTCGGCAAGCGCTGGCAACCGGGCTTc                                                                                                                    >  2:767359/1‑140 (MQ=255)            attttattttTTCTCGACGACGGGACCTGCCTGATACCCGCCGCCCAGCGATTTCATCAACTGAATGCTTTGGATCACCCGGCGGCTGTCCAGCATCAGTAATGACATCTCTTCGGCAAGCGCTGGCAACCGGGCTTCgg                                                                                                                  <  2:861950/140‑1 (MQ=255)            attttattttTTCTCGACGACGGGACCTGCCTGATACCCGCCGCCCAGCGATTTCATCAACTGAATGCTTTGGATCACCCGGCGGCTGTCCAGCATCAGTAATGACATCTCTTCGGCAAGCGCTGGCAACCGGGCTTCgg                                                                                                                  <  2:176668/140‑1 (MQ=255)                tattttTTCTCGACGACGGGACCTGCCTGATACCCGCCGCCCAGCGATTTCATCAACTGAATGCTTTGGATCACCCGGCGGCTGTCCAGCATCAGTAATGACATCTCTTCGGCAAGCGCTGGCAACCGGGCTTCGGTGGc                                                                                                              <  1:767359/140‑1 (MQ=255)                     tttCTCGACGACGGGACCTGCCTGATACCCGCCGCCCAGCGATTTCATCAACTGAATGCTTTGGATCACCCGGCGGCTGTCCAGCATCAGTAATGACATCTCTTCGGCAAGCGCTGGCAACCGGGCTTCGGTGGCCTGTa                                                                                                         >  1:339154/1‑140 (MQ=255)                      ttCTCGACGACGGGACCTGCCTGATACCCGCCGCCCAGCGATTTCATCAACTGAATGCTTTGGATCACCCGGCGGCTGTCCAGCATCAGTAATGACATCTCTTCGGCAAGCGCTGGCAACCGGGCTTCGGTGGCCTGTaa                                                                                                        >  2:833211/1‑140 (MQ=255)                         tcGACGACGGGACCTGCCTGATACCCGCCGCCCAGCGATTTCATCAACTGAATGCTTTGGATCACCCGGCGGCTGTCCAGCATCAGTAATGACATCTCTTCGGCAAGCGCTGGCAACCGGGCTTCGGTGGCCTGTAAGCg                                                                                                     <  2:482132/140‑1 (MQ=255)                          cgacgaCGGGACCTGCCTGATACCCGCCGCCCAGCGATTTCATCAACTGAATGCTTTGGATCACCCGGCGGCTGTCCAGCATCAGTAATGACATCTCTTCGGCAAGCGCTGGCAACCGGGCTTCGGTGGCCTGTAAGCgg                                                                                                    <  1:194912/140‑1 (MQ=255)                          cgacgaCGGGACCTGCCTGATACCCGCCGCCCAGCGATTTCATCAACTGAATGCTTTGGATCACCCGGCGGCTGTCCAGCATCAGTAATGACATCTCTTCGGCAAGCGCTGGCAACCGGGCTTCGGTGGCCTGTAAGCgg                                                                                                    <  1:833211/140‑1 (MQ=255)                                       tgcctgATACCCGCCGCCCAGCGATTTCATCAACTGAATGCTTTGGATCACCCGGCGGCTGTCCAGCATCAGTAATGACATCTCTTCGGCAAGCGCTGGCAACCGGGCTTCGGTGGCCTGTAAGCGGCTGGTTAAGCcgc                                                                                       <  2:446095/140‑1 (MQ=255)                                         cctgATACCCGCCGCCCAGCGATTTCATCAACTGAATGCTTTGGATCACCCGGCGGCTGTCCAGCATCAGTAATGACATCTCTTCGGCAAGCGCTGGCAACCGGGCTTCGGTGGCCTGTAAGCGGCTGGTTAAGCCgcgc                                                                                     >  1:542569/1‑140 (MQ=255)                                                   gccgccCAGCGATTTCATCAACTGAATGCTTTGGATCACCCGGCGGCTGTCCAGCATCAGTAATGACATCTCTTCGGCAAGCGCTGGCAACCGGGCTTCGGTGGCCTGTAAGCGGCTGGTTAAGCCGCGCTGATAggcgg                                                                           >  2:1024109/1‑140 (MQ=255)                                                          aGCGATTTCATCAACTGAATGCTTTGGATCACCCGGCGGCTGTCCAGCATCAGTAATGACATCTCTTCGGCAAGCGCTGGCAACCGGGCTTCGGTGGCCTGTAAGCGGCTGGTTAAGCCGCGCTGATAGGCGGCCTCGGc                                                                    >  1:183836/1‑140 (MQ=255)                                                               tttCATCAACTGAATGCTTTGGATCACCCGGCGGCTGTCCAGCATCAGTAATGACATCTCTTCGGCAAGCGCTGGCAACCGGGCTTCGGTGGCCTGTAAGCGGCTGGTTAAGCCGCGCTGATAGGCGGCCTCGGCAgcgc                                                               <  2:54159/140‑1 (MQ=255)                                                                                     aTCACCCGGCGGCTGTCCAGCATCAGTAATGACATCTCTTCGGCAAGCGCTGGCAACCGGGCTTCGGTGGCCTGTAAGCGGCTGGTTAAGCCGCGCTGATAGGCGGCCTCGGCAGCGCGCTGGGTAAAGCGCGTGGCTTc                                         <  1:1024109/140‑1 (MQ=255)                                                                                         cccGGCGGCTGTCCAGCATCAGTAATGACATCTCTTCGGCAAGCGCTGGCAACCGGGCTTCGGTGGCCTGTAAGCGGCTGGTTAAGCCGCGCTGATAGGCGGCCTCGGCAGCGCGCTGGGTAAAGCGCGTGGCTTCCAcg                                     >  2:726755/1‑140 (MQ=255)                                                                                          ccGGCGGCTGTCCAGCATCAGTAATGACATCTCTTCGGCAAGCGCTGGCAACCGGGCTTCGGTGGCCTGTAAGCGGCTGGTTAAGCCGCGCTGATAGGCGGCCTCGGCAGCGCGCTGGGTAAAGCGCGt                                               >  1:655176/1‑129 (MQ=255)                                                                                          ccGGCGGCTGTCCAGCATCAGTAATGACATCTCTTCGGCAAGCGCTGGCAACCGGGCTTCGGTGGCCTGTAAGCGGCTGGTTAAGCCGCGCTGATAGGCGGCCTCGGCAGCGCGCTGGGTAAAGCGCGt                                               <  2:655176/129‑1 (MQ=255)                                                                                              cggcTGTCCAGCATCAGTAATGACATCTCTTCGGCAAGCGCTGGCAACCGGGCTTCGGTGGCCTGTAAGCGGCTGGTTAAGCCGCGCTGATAGGCGGCCTCGGCAGCGCGCTGGGTAAAGCGCGTGGCTTCCACGCGTTc                                >  1:323934/1‑140 (MQ=255)                                                                                                 cTGTCCAGCATCAGTAATGACATCTCTTCGGCAAGCGCTGGCAACCGGGCTTCGGTGGCCTGTAAGCGGCTGGTTAAGCCGCGCTGATAGGCGGCCTCGGCAGCGCGCTGGGTAAAGCGCGTGGCTTCCACGCGTTCAGc                             <  2:1063132/140‑1 (MQ=255)                                                                                                    tCCAGCATCAGTAATGACATCTCTTCGGCAAGCGCTGGCAACCGGGCTTCGGTGGCCTGTAAGCGGCTGGTTAAGCCGCGCTGATAGGCGGCCTCGGCAGCGCGCTGGGTAAAGCGCGTGGCTTCCACGCGTTCAGCCTg                          >  2:18447/1‑140 (MQ=255)                                                                                                              gTAATGACATCTCTTCGGCAAGCGCTGGCAACCGGGCTTCGGTGGCCTGTAAGCGGCTGGTTAAGCCGCGCTGATAGGCGGCCTCGGCAGCGCGCTGGGTAAAGCGCGTGGCTTCCACGCGTTCAGCCTGCATTTctcgc                >  1:484574/1‑140 (MQ=255)                                                                                                                   gACATCTCTTCGGCAAGCGCTGGCAACCGGGCTTCGGTGGCCTGTAAGCGGCTGGTTAAGCCGCGCTGATTGGCGGCCTCGGCAGCGCGCTGGGTAAAGCGCGTGGCTTCCACGCGTTCAGCCTGCATTTCTCGCtcgtc           >  2:378398/1‑140 (MQ=255)                                                                                                                   gACATCTCTTCGGCAAGCGCTGGCAACCGGGCTTCGGTGGCCTGTAAGCGGCTGGTTAAGCCGCGCTGATAGGCGGCCTCGGCAGCGCGCTGGGTAAAGCGCGTGGCTTCCACGCGTTCAGCCTGCATTTCTCGCtcgtc           >  1:953317/1‑140 (MQ=255)                                                                                                                       tctcTTCGGCAAGCGCTGGCAACCGGGCTTCGGTGGCCTGTAAGCGGCTGGTTAAGCCGCGCTGATAGGCGGCCTCGGCAGCGCGCTGGGTAAAGCGCGTGGCTTCCACGCGTTCAGCCTGCATTTCTCGCTCGTCGTTg       >  1:1921/1‑140 (MQ=255)                                                                                                                       tctcTTCGGCAAGCGCTGGCAACCGGGCTTCGGTGGCCTGTAAGCGGCTGGTTAAGCCGCGCTGATAGGCGGCCTCGGCAGCGCGCTGGGTAAAGCGCGTGGCTTCCACGCGTTCAGCCTGCATTTCTCGCTCGTCGTTg       >  2:833182/1‑140 (MQ=255)                                                                                                                        ctctTCGGCAAGCGCTGGCAACCGGGCTTCGGTGGCCTGTAAGCGGCTGGTTAAGCCGCGCTGATAGGCGGCCTCGGCAGCGCGCTGGGTAAAGCGCGTGGCTTCCACGCGTTCAGCCTGCATTTCTCGCTCGTCGTTGa      >  1:775786/1‑140 (MQ=255)                                                                                                                        ctctTCGGCAAGCGCTGGCAACCGGGCTTCGGTGGCCTGTAAGCGGCTGGTTAAGCCGCGCTGATAGGCGGCCTCGGCAGCGCGCTGGGTAAAGCGCGTGGCTTCCACGCGTTCAGCCTGCATTTCTCGCTCGTCGTTGa      >  2:7598/1‑140 (MQ=255)                                                                                                                         tctTCGGCAAGCGCTGGCAACCGGGCTTCGGTGGCCTGTAAGCGGCTGGTTAAGCCGCGCTGATAGGCGGCCTCGGCAGCGCGCTGGGTAAAGCGCGTGGCTTCCACGCGTTCAGCCTGCATTTCTCGCTCGTCGTTGAg     >  2:41633/1‑140 (MQ=255)                                                                                                                         tctTCGGCAAGCGCTGGCAACCGGGCTTCGGTGGCCTGTAAGCGGCTGGTTAAGCCGCGCTGATAGGCGGCCTCGGCAGCGCGCTGGGTAAAGCGCGTGGCTTCCACGCGTTCAGCCTGCATTTCTCGCTCGTCGTTGAg     >  2:397375/1‑140 (MQ=255)                                                                                                                         tctTCGGCAAGCGCTGGCAACCGGGCTTCGGTGGCCTGTAAGCGGCTGGTTAAGCCGCGCTGATAGGCGGCCTCGGCAGCGCGCTGGGTAAAGCGCGTGGCTTCCACGCGTTCAGCCTGCATTTCTCGCTCGTCGTTGAg     >  1:748580/1‑140 (MQ=255)                                                                                                                          ctTCGGCAAGCGCTGGCAACCGGGCTTCGGTGGCCTGTAAGCGGCTGGTTAAGCCGCGCTGATAGGCGGCCTCGGCAGCGCGCTGGGTAAAGCGCGTGGCTTCCACGCGTTCAGCCTGCATTTCTCGCTCGTCGTTGAGc    <  2:748580/140‑1 (MQ=255)                                                                                                                            tCGGCAAGCGCTGGCAACCGGGCTTCGGTGGCCTGTAAGCGGCTGGTTAAGCCGCGCTGATAGGCGGCCTCGGCAGCGCGCTGGGTAAAGCGCGTGGCTTCCACGCGTTCAGCCTGCATTTCTCGCTCGTCGTTGAGCGt  >  1:366665/1‑140 (MQ=255)                                                                                                                                     |                                                                                                                                   ACGCGGCAGACATTTTATTTTTTCTCGACGACGGGACCTGCCTGATACCCGCCGCCCAGCGATTTCATCAACTGAATGCTTTGGATCACCCGGCGGCTGTCCAGCATCAGTAATGACATCTCTTCGGCAAGCACTGGCAACCGGGCTTCGGTGGCCTGTAAGCGGCTGGTTAAGCCGCGCTGATAGGCGGCCTCGGCAGCGCGCTGGGTAAAGCGCGTGGCTTCCACGCGTTCAGCCTGCATTTCTCGCTCGTCGTTGAGCGT  >  NC_000913/4299550‑4299812

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