Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,744,751	Δ1 bp	coding (1295/1374 nt)	mdtK →	multidrug efflux pump MdtK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,744,751	0	C	.	100.0%	128.1 / NA	29	coding (1295/1374 nt)	mdtK	multidrug efflux pump MdtK
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (15/14);  total (15/14)

GGTATATCAGATTTCTGACTCAATCCAGGTGATTGGCAGTGGGATTTTGCGTGGTTATAAAGATACGCGTTCCATTTTCTATATTACCTTTACGGCTTACTGGGTGCTGGGCTTGCCAAGCGGCTATATTCTGGCACTGACCGATCTGGTCGTTGAACCTATGGGGCCAGCAGGCTTCTGGATAGGCTTTATTATTGGCCTGACGTCGGCAGCCATTATGATGATGTTGCGTATGCGGTTCCTGCAACGTCTGCCGTCAGCCATCATTCTGCAACGAGCATCCCGCTAATAAAGACAAGGCGCAACCTTCACGGGTTGCGCCTGTATTTTTACGCAGGCTGGAGCGTTGCGCCAATCCCGTCTTCGTCTGGCTGTAATTTCAGAGCGTTACACAGAGTTTGCCGACTGGCATCTGTTATCTCAACGTGTATATCCCCGGTAACTTCCCCTT  >  NC_000913/1744542‑1744992                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                  ggTATATCAGATTTCTGACTCAATCCAGGTGATTGGCAGTGGGATTTTGCGTGGTTATAAAGATACGCGTTCCATTTTCTATATTACCTTTACGGCTTACTGGGTGCTGGGCTTGCCAAGCGGCTATATTCTGGCACTGACCGATCTGGTCGTTGAACCTATGGGGCCAGCAGGCTTCTGGATAGGCTTTATTATTGGCCTGACGTCGG‑AGCCATTATGATGATGTTGCGTATGCGGTTCCTGCAACGTCt                                                                                                                                                                                                         <  1:208912/251‑1 (MQ=255)                tGACTCAATCCAGGTGATTGGCAGTGGGATTTTGCGTGGTTATAAAGATACGCGTTCCATTTTCTATATTACCTTTACGGCTTACTGGGTGCTGGGCTTGCCAAGCGGCTATATTCTGGCACTGACCGATCTGGTCGTTGAACCTATGGGGCCAGCAGGCTTCTGGATAGGCTTTATTATTGGCCTGACGTCGG‑AGCCATTATGATGATGTTGCGTATGCGGTTCCTGCAACGTCTGCCGTCAGCcatcat                                                                                                                                                                                          >  1:234714/1‑251 (MQ=255)                                aTTGGCAGTGGGATTTTGCGTGGTTATAAAGATACGCGTTCCATTTTCTATATTACCTTTACGGCTTACTGGGTGCTGGGCTTGCCAAGCGGCTATATTCTGGCACTGACCGATCTGGTCGTTGAACCTATGGGGCCAGCAGGCTTCTGGATAGGCTTTATTATTGGCCTGACGTCGG‑AGCCATTATGATGATGGTGCGTATGCGGTTCCTGCAACGTCTGCCGTCAGCCATCATTCTGCAACGAGCATcc                                                                                                                                                                          <  1:63911/251‑1 (MQ=255)                                aTTGGCAGAGGGATTTTGCGTGGTTATAAAGATACGCGTTCCATTTTCTATATTACCTTTACGGCTTACTGGGTGCTGGGCTTGCCAAGCGGCTATATTCTGGCACTGACCGATCTGGTCGTTGAACCTATGGGGCCAGCAGGCTTCTGGATAGGCTTTATTATTGGCCTGACGTCGG‑AGCCATTATGATGATGTTGCGTATGCGGTTCCTGCAACGTCTGCCGTCAGCCATCATTCTGCAACGAGCATcc                                                                                                                                                                          >  1:73979/1‑251 (MQ=255)                                  tGNCAGTGGGATTTTGCGTGGTTATAAAGATACGCGTTCCATTTTCTATATTACCTTTACGGCTTACTGGGTGCTGGGCTTGCCAAGCGGCTATATTCTGGCACTGACCGATCTGGTCGTTGAACCTATGGGGCCAGCAGGCTTCTGGATAGGCTTTATTATTGGCCTGACGTCGG‑AGCCATTATGATGATGTTGCGTATGCGGTTCCTGCAACGTCTGCCGTCAGCCATCATTCTGCAACGAGCATccc                                                                                                                                                                         <  2:353708/250‑1 (MQ=255)                                       gTGGGATTTTGCGTGGTTATAAAGATACGCGTTCCATTTTCTATATTACCTTTACGGCTTACTGGGTGCTGGGCTTGCCAAGCGGCTATATTCTGGCACTGACCGATCTGGTCGTTGAACCTATGGGGCCAGCAGGCTTCTGGATAGGCTTTATTATTGGCCTGACGTCGG‑AGCCATTATGATGATGTTGCGTATGCGGTTCCTGCAACGTCTGCCGTCAGCCATCATTCTGCAACGAGCATCCCGCtaat                                                                                                                                                                   >  2:332157/1‑251 (MQ=255)                                                         ataAAGCTACGCGTTCCATTTTCTATATTACCTTTACGGCTTACTGGGTGCTGGGCTTGCCAAGCGGCTATATTCTGGCACTGACCGATCTGGTCGTTGAACCTATGGGGCCAGCAGGCTTCTGGATAGGCTTTATTATTGGCCTGACGTCGG‑AGCCATTATGATGATGTTGCGTATGCGGTTCCTGCAACGTCTGCCGTCAGCcatcat                                                                                                                                                                                          >  2:51490/1‑210 (MQ=255)                                                         ataAAGATACGCGTTCCATTTTCTATATTACCTTTACGGCTTACTGGGTGCTGGGCTTGCCAAGCGGCTATATTCTGGCACTGACCGATCTGGTCGTTGAACCTATGGGGCCAGCAGGCTTCTGGATAGGCTTTATTATTGGCCTGACGTCGG‑AGCCATTATGATGATGTTGCGTATGCGGTTCCTGCAACGTCTGCCGTCAGCcatcat                                                                                                                                                                                          <  1:51490/210‑1 (MQ=255)                                                           aaaGATACGCGTTCCATTTTCTATATTACCTTTACGGCTTACTGGGTGCTGGGCTTGCCAAGCGGCTATATTCTGGCACTGACCGATCTGGTCGTTGAACCTATGGGGCCAGCAGGCTTCTGGATAGGCTTTATTATTGGCCTGACGTCGG‑AGCCATTATGATGATGTTGCGTATGCGGTTCCTGCAACGTCTGCCGTCAGCCATCATTCTGCAACGAGCATCCCGCTAATAAAGACAAGGCGCAACCTTc                                                                                                                                               >  1:139711/1‑251 (MQ=255)                                                              aaTACGCGTTCCATTTTCTATATTACCTTTACGGCTTACTGGGTGCTGGGCTTGCCAAGCGGCTATATTCTGGCACTGACCGATCTGGTCGTTGAACCTATGGGGCCAGCAGGCTTCTGGATAGGCTTTATTATTGGCCTGACGTCGG‑AGCCATTATGATGATGTTGCGTATGCGGTTCCTGCAACGTCTGCCGTCAGCCATCATTCGGCAACGAGCATCCCGCTAATAAAGACAAGGCGCAACCTTCACg                                                                                                                                            <  1:105954/250‑1 (MQ=255)                                                               aTACGCGTTCCATTTTCTATATTACCTTTACGGCTTACTGGGTGCTGGGCTTGCCAAGCGGCTATATTCTGGCACTGACCGATCTGGTCGTTGAACCTATGGGGCCAGCAGGCTTCTGGATAGGCTTTATTATTGGCCTGACGTCGG‑AGCCATTATGATGATGTTGCGTATGCGGTTCCTGCAACGTCTGCCGTCAGCCATCATTCTGCAACGAGCATCCCGCTAATAAAGACAAGGCGCAACCTTCACgg                                                                                                                                           >  1:365412/1‑251 (MQ=255)                                                                     gTTCCATTTTCTATATTACCTTTACGGATTACTGGGTGCTGGGCTTGCCAAGCGGCTATATTCTGGCACTGACCGATCTGGTCGTTGAACCTATGGGGCCAGCAGGCTTCTGGATAGGCTTTAATATTGGCCTGACGTCGG‑AGCCATTATGATGATGTTGCGTATGCGGTTCCTGCAACGTCTGCCGTCAACCATCATTCTGCAACGAGCATCCCGCTAATAAAGAAAAGCCGCAACCTTCACGGGTTgcg                                                                                                                                     >  1:236533/1‑251 (MQ=255)                                                                       tCCATTTTCTATATTACCTTTACGGCTTACTGGGTGCTGGGCTTGCCAAGCGGCTATATTCTGGCACTGACCGGTCTGGTCGTTGAACCTATGGGGCCAGCAGGCTTCTGGATAGGCTTTATTATTGGCCTGACGTCGG‑AGCCATTATGATGATGTTGCGTATGCGGTTCCTGCAACGTCTGCCGTCAGCCATCATTCTGCAACGAGCATCCCGCTAATAAAGACAAGGCGCAACCTTCACGGGTTGCGcc                                                                                                                                   >  1:68094/1‑251 (MQ=255)                                                                                 atntTACCTTTACGGCTTACTGGGTGCTGGGCTTGCCAAGCGGCTATATTCTGGCACTGACCGATCTGGTCGTTGAACCTATGGGGCCAGCAGGCTTCTGGATAGGCTTTATTATTGGCCTGACGTCGG‑AGCCATTATGATGATGTTGCGTATGCGGTTCCTGCAACGTCTGCCGTCAGCCATCATTCTGCAACGAGCATCCCGCTAATAAAGACAAGGCGCAACCTTCACGGGTTGCGCCTGTAttttt                                                                                                                          <  2:365412/250‑1 (MQ=255)                                                                                      cccTTTACGGCTTACTGGGTGCTGGGCTTGCCAAGCGGCTATATTCTGGCACTGACCGATCTGGTCGTTGAACCTATGGGGCCAGCAGGCTTCTGGATAGGCTTTATTATTGGCCTGACGTCGG‑AGCCATTATGATGATGTTGCGTATGCGGTTCCTGCAACGTCTGCCGTCAGCCATCATTCTGCAACGTGCATCCCGCTAATAAAGACAAGGCGCAACCTTCACGGGTTGCGCCTGTATTTTTACGCa                                                                                                                     <  1:162467/249‑1 (MQ=255)                                                                                                cTTACTGGGTGCTGGGCTTGCCAAGCGGCTATATTCTGGCACTGACCGATCTGGTCGTTGAACCTATGGGGCCAGCAGGCTTCTGGATAGGCTTTATTATTGGCCTGACGTCGG‑AGCCATTATGATGATGTTGCGTATGCGGTTCCTGCAACGTCTGCCGTCAGCCATCATTCTGCAACGAGCATCCCGCTAATAAAGACAAGGCGCAACCTTCACGGGTTGCGCCTGTATTTTTACGCAGGCTGGAGCGt                                                                                                          >  1:14421/1‑251 (MQ=255)                                                                                                     tGGGTGCTGGGCTTGCCAAGCGGCTATATTCTGGCACTGACCGATCTGGTCGTTGAACCTATGGGGCCAGCAGGCTTCTGGATAGGCTTTATTATTGGCCTGACGTCGG‑AGCCATTATGATGATGTTGCGTATGCGGTTCCTGCAACGTCTGCCGTCAGCCATCATTCTGCAACGAGCATCCCGCTAATAAAGACAAGGCGCAACCTTCACGGGTTGCGCCTGTATTTTTACGCAGGCTGGAGCGTTgcgc                                                                                                     >  1:400631/1‑251 (MQ=255)                                                                                                                                        ccTGACCGATCTGGTCGTTGAACCTATGGGGCCAGCAGGCTTCTGGATAGGCTTTATTATTGGCCTGACGTCGG‑AGCCATTATGATGATGTTGCGTATGCGGTTCCTGCAACGTCTGCCGTCAGCCATCATTCTGCAACGAGCATCCCGCTAATAAAGACAAGGCGCAACCTTCACGGGTTGCGCCTGTATTTTTACGCATGCTGGAGCGTTGCGCCAATCCCGTCTTCGTCTGGCTGTAATTTCAGAGCg                                                                  <  1:175903/250‑1 (MQ=255)                                                                                                                                           gACCGATCTGGTCGTTGAGCCTATGGGGCCAGCAGGCTTATGGATAGGCTTTATTAGTGGCCTGACGTCGG‑AGCCATTATGATGATGTTGCGTATGCGGTTCCTGCAACGTCTGCCGTCAGCCATCATTCTGCAACGAGCATCCCGCTAATAAAGACAAGGCGCAACCTTCACGGGTTGCGCCTGTATTTTTACGCAGGCTGGAGCGTTGCGCCAATCCCGTCTTCGTCTGGCTGTAATTTCAGAGCGt                                                                 <  2:177820/249‑1 (MQ=255)                                                                                                                                           gACCGATCTGGTCGTTGAACCTATGGGGCCAGCAGGCTTCTGGATAGGCTTTATTATTGGCCTGACGTCGG‑AGCCATTATGATGATGTTGCGTATGCGGTTCCTGCAACGTCTGCCGTCAGCCATCATTCTGCAACGAGCATCCCGCTAATAAAGACAAGGCGCAACCTTCACGGGTTGCGCCTGTATTTTTACGCAGGCTGGAGCGTTGCGCCAATCCCGTCTTCGTCTGGCTGTAATTTCAGAGCGtt                                                                >  1:3271/1‑250 (MQ=255)                                                                                                                                                           gAACCTATGGGGCCAGCATGCTTCTGGATAGGCTTTATTATTGGCCTGACGTCGG‑AGCCATTATGATGATGTTGCGTATGCGGTTCCTGCAACGTCTGCCGTCAGCCATCATTCTGCAACGAGCATCCCGCTAATAAAGACAAGGCGCAACCTTCACGGGTTGCGCCTGTATTTTTACGCAGGCTGGAGCGTTGCGCCAATCCCGTCTTCGTCTGGCTGTAATTTCAGAGCGTTACACAGAGTTTGccta                                                <  2:301952/250‑3 (MQ=255)                                                                                                                                                                    gggCCAGCAGGCTTCTGGATAGGCTTTATTATTGGCCTGACGTCGG‑AGCCATTATGATGATGTTGCGTATGCGGTTCCTGCAACGTCTGCCGTCAGCCATCATTCTGCAACGAGCATCCCGCTAATAAAGACAAGGCGCAACCTTCACGGGTTGCGCCTGTATTTTTACGCAGGCTGGAGCGTTGCGCCAATCCCGTCTTCGTCTGGCTGTAATTTCAGAGCGTTACACAGAGTTTGCCGACTGGCATCTg                                      >  1:365325/1‑251 (MQ=255)                                                                                                                                                                       ccAGCAGGCTTCTGGATAGGCTTTATTATTGGCCTGACGTCGG‑AGCCATTATGATGATGTTGCGTATGCGGTTCCTGCAACGTCTGCCGTCAGCCATCATTCTGCAACGAGCATCCCGCTAATAAAGACAAGGCGCAACCTTCACGGGTTGCGCCTGTATTTTTACGCAGGCTGGAGCGTTGCGCCAATCCCGTCTTCGTCTGGCTGTCATTTCAGAGCGTTACACAGAGTTTGCCGACTGGCATCTGtt                                    >  1:282496/1‑250 (MQ=255)                                                                                                                                                                           cagGCTTCTGGATAGGCTTTATTATTGGCCTGACGTCGG‑AGCCATTATGATGATGTTGCGTATGCGGTTCCTGCAACGTCTGCCGTCAGCCATCATTCTGCAACGAGCATCCCGCTAATAAAGACAAGGCGCAACCTTCACGGGTTGCGCCTGTATTTTTACGCAGGCTGGAGCGTTGCGCCAATCCCGTCTTCGTCTGGCTGTAATTTCAGAGCGTTACACAGAGTTTGCCGACTGGCATCTGTTAtctc                               <  2:234714/251‑1 (MQ=255)                                                                                                                                                                             ggcttcgtgttAGGCTTTATTATTGGCCTGACGTCGG‑AGCCATTATGATGATGTTGCGTATGCGGTTCCTGCACCGTCTGCCGTCAGCCATCATTCTGCAACGAGCATCCCGCTAATAAAGACAAGGCGCAACCTTCACGGGTTGCGCCTGTATTTTTACGCAGGCTGGAGCGTTGCGCCAATCCCGTCTTCGTCTGGCTGTAATTCCAGAGCGTTACACAGAGTTTGCCGACTGGCATCTGTTATCTCaa                             <  2:108975/241‑1 (MQ=255)                                                                                                                                                                                       tAGGCTTTATTATTGGCCTGACGTCGG‑AGCCATTATGATGATGTTGCGTATGCGGTTCCTGCAACGTCTGCCGTCAGCCATCATTCTGCAACGAGCATCCCGCTAATAAAGACAAgg                                                                                                                                                         <  1:245664/117‑1 (MQ=255)                                                                                                                                                                                       tAGGCTTTATTATTGGCCTGACGTCGG‑AGCCATTATGATGATGTTGCGTATGCGGTTCCTGCAACGTCTGCCGTCAGCCATCATTCTGCAACGAGCATCCCGCTAATAAAGACAAgg                                                                                                                                                         >  2:245664/1‑117 (MQ=255)                                                                                                                                                                                                 tattGGCCTGACGGCGG‑AGCCATTATGATGATGTTGCGTATGCGGTTCCTGCAACGTCTGCCGTCAGCCATCATTCTGCAACGAGCATCCCGCTAATAAAGACAAGGCGCAACCTTCACGGGTTGCGCCTGTATTTTTACGCAGGCTGGAGCGTTGCGGCAATCCCGTCTTCGTCTGGCTGTAATTTCAGAGCGTTACACAGAGTTTGCCGACTGGCATCTGTTAACTCAACGTGTATATCCCCGGTAAc          <  1:332157/250‑1 (MQ=255)                                                                                                                                                                                                        cTGACGTCGG‑AGCCATTATGATGATGTTGCGTATGCGGTTCCTGCAACGTCTGCCGTCAGCCATCATTCTGCAACGAGCATCCCGCTAATAAAGACAAGGCGCAACCTTCACGGGTTGCGCCTGTATTTTTACGCAGGCTGGAGCGTTGCGCCAATCCCGTCTTCGTCTGGCTGTAATTTCAGAGCGTTACACAGAGTTTGCCGACTGGCATCTGTTATCTCAACGTGTATATCCCCGGTAACTTCCCCtt  >  1:382672/1‑251 (MQ=255)                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                  GGTATATCAGATTTCTGACTCAATCCAGGTGATTGGCAGTGGGATTTTGCGTGGTTATAAAGATACGCGTTCCATTTTCTATATTACCTTTACGGCTTACTGGGTGCTGGGCTTGCCAAGCGGCTATATTCTGGCACTGACCGATCTGGTCGTTGAACCTATGGGGCCAGCAGGCTTCTGGATAGGCTTTATTATTGGCCTGACGTCGGCAGCCATTATGATGATGTTGCGTATGCGGTTCCTGCAACGTCTGCCGTCAGCCATCATTCTGCAACGAGCATCCCGCTAATAAAGACAAGGCGCAACCTTCACGGGTTGCGCCTGTATTTTTACGCAGGCTGGAGCGTTGCGCCAATCCCGTCTTCGTCTGGCTGTAATTTCAGAGCGTTACACAGAGTTTGCCGACTGGCATCTGTTATCTCAACGTGTATATCCCCGGTAACTTCCCCTT  >  NC_000913/1744542‑1744992

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