Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,303,755	C→T	V119I (GTT→ATT)	mdtN ←	putative multidrug efflux pump membrane fusion protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,303,755	0	C	T	96.8%	96.4 / ‑5.8	32	V119I (GTT→ATT)	mdtN	putative multidrug efflux pump membrane fusion protein
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (16/14);  total (16/16)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

CCGCTTCGACCGCCGCACGCTGGGCAACTAATGCATCCACGCCGCTGACGGCGCTGGCGGCTGACTGCGCCTGTAACAATACGGCATTAAGATCCGCTTCTGCGGCGCGCTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATG  >  NC_000913/4303518‑4304000                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                      ccGCTTCGACCGCCGCACGCTGGGCAACTAATGCATCCACGCCGCTGACGGCGCTGGCGGCTGACTGCGCCTGTAACAATACGGCATTAAGATCCGCTTCTGCGGCGAGCTGCGCCGTTCTTGCACGGTTAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGGATCTGGGGCCTGTTTCGCTGCGGCAAGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCAcc                                                                                                                                                                                                                                           >  2:192158/1‑250 (MQ=255)                           aCTAAGGCATCCACGCCGCTGACGGCGCTGGCGGCTGACTGCGCGTGTAACAATACGGCATTAAGATCCGCTTCTGCGGCGCTCTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGGAAGGTATCTGTGGCGGGTTGCGCGGCGGCACGGGCTTTTTTTAGCGGTGCAGTAAGCGGGTCGGCACGAAACTGTTGCGTTTCAACGCGACGCTg                                                                                                                                                                                                                <  2:33716/251‑1 (MQ=255)                                   aTCCACGCCGCTGACGGCGCTGGCGGCTGACTGCGCCTGTAACAATACGGCATTAAGATCCGCTTCTGCGGCGCGCTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGc                                                                                                                                                                                                         >  2:104292/1‑250 (MQ=255)                                   aTCCACGCCGCTGACGGCGCTGGCGGCTGACTGCGCCTGTAACAATACGGCATTAAGATCCGCTTCTGCGGCGCGCTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGc                                                                                                                                                                                                         >  1:11963/1‑250 (MQ=255)                                            gcTGACGGCGCTGGCGGCTGACTGCGCCTGTAACAATACGGCATTAAGATCCGCTTCTGCGGCGCGCTGCGCCGTTCTTGAACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCCGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCt                                                                                                                                                                                               >  1:37652/1‑251 (MQ=255)                                                            gcTGACTGCGCCTGTAACAATACGGCATTAAGATCCGCTTCTGCGGCGCGCTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCgag                                                                                                                                                                               <  1:223343/251‑1 (MQ=255)                                                                        tGTAACAATACGGCATTAAGATCCGCTTCTGCGGCGCGCTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGAATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACAAAACTGCTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGATGCTACAGc                                                                                                                                                                   >  1:40935/1‑251 (MQ=255)                                                                           aaCAATACGGCATTAAGATCCGCTTCTGCGGCGCGCTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCtt                                                                                                                                                                 <  1:67262/250‑1 (MQ=255)                                                                                 aCGGCATTAAGATCCGCTTCTGCGGCGCGCTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTGATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAg                                                                                                                                                          <  1:112152/251‑1 (MQ=255)                                                                                    gCATTAAGATCCGCTTCTGCGGCGCGCTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGAtt                                                                                                                                                       <  1:199770/251‑1 (MQ=255)                                                                                          aGATCCGCTTCTGCGGCGCGCTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCCTTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTc                                                                                                                                                 >  2:121079/1‑251 (MQ=255)                                                                                            aTCCGCTTCTGCGGCGCGCTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCgg                                                                                                                                               >  1:245001/1‑250 (MQ=255)                                                                                                           gcgcTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCGTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGGTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGTCAGATTGGCTTCGTACGGGCGCGGGTCGa                                                                                                                                <  2:139359/251‑1 (MQ=255)                                                                                                              cTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTCTCTGTGGCCTGTTTCGCGGCGGCACGGGATTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATgg                                                                                                                             >  2:278712/1‑250 (MQ=255)                                                                                                                                gTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCt                                                                                                           <  1:115381/251‑1 (MQ=255)                                                                                                                                 tCAACATCTTCCGATGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTg                                                                                                          >  2:304237/1‑251 (MQ=255)                                                                                                                                            cGCTGAGACAAAACCTTCTTTCAGTAATTGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCTAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCt                                                                                               <  1:200412/251‑1 (MQ=255)                                                                                                                                                         ccTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACGGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGTTGaa                                                                                  <  2:235027/251‑2 (MQ=255)                                                                                                                                                                    tAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCCGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCt                                                                        >  2:254706/1‑250 (MQ=255)                                                                                                                                                                      aTGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGCGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGAACGCCAGTTCTAc                                                                      >  1:204312/1‑250 (MQ=255)                                                                                                                                                                         gCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTGCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATg                                                                  <  2:91945/251‑1 (MQ=255)                                                                                                                                                                           ccGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCg                                                                <  1:111057/250‑1 (MQ=255)                                                                                                                                                                             ggTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGc                                                              >  2:301056/1‑251 (MQ=255)                                                                                                                                                                              gTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGcc                                                             >  1:37159/1‑251 (MQ=255)                                                                                                                                                                                   gcgTAAAGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACCGCAAATCGCCCTGTTTGCCTGCCTGGTTGTCGGTGACCGCCCGTTCTACAAGGCGGCCGCTgc                                                        >  2:204127/1‑250 (MQ=255)                                                                                                                                                                                        aTGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTcc                                                   <  1:121079/251‑1 (MQ=255)                                                                                                                                                                                               tgtgGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTAcc                                             >  1:280306/1‑250 (MQ=255)                                                                                                                                                                                                    ccTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCa                                       <  1:305531/251‑1 (MQ=255)                                                                                                                                                                                                                   cACGGGCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGAc                        >  1:70561/1‑251 (MQ=255)                                                                                                                                                                                                                      gggCTTTTTCTACCGTGGCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGAcgcg                     <  1:301056/251‑1 (MQ=255)                                                                                                                                                                                                                                      ggCATTAATCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGCCTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGtt      <  1:5678/250‑1 (MQ=255)                                                                                                                                                                                                                                         aTTAACCTAGTCGTCACCAACCTGTTGCGCGTCAACGCTCCGCTTGGTGAGCATAATTTGCTTATCCAGCGCCGCGCGGGAGCCTTCAGCTTTCGCCAGATTTGCTGCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTGGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGTTATCTGCTGACGCGTAAGCGTCATTAGTTGATg  <  1:89620/251‑1 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                      CCGCTTCGACCGCCGCACGCTGGGCAACTAATGCATCCACGCCGCTGACGGCGCTGGCGGCTGACTGCGCCTGTAACAATACGGCATTAAGATCCGCTTCTGCGGCGCGCTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTTGACTGCCTGGTTGTCGGTGACCGCCAGTTCTACAATGCGGCCGCTGACTTCCGGTACCACATCAATGGTATCTGCTGACGCGTAAGCGTCATTAGTTGATG  >  NC_000913/4303518‑4304000

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