Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,567,548	A→G	G324G (GGT→GGC)	mdtM ←	multidrug efflux system protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,567,548	0	A	G	100.0%	157.7 / NA	41	G324G (GGT→GGC)	mdtM	multidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (24/17);  total (24/17)

GCCGTTAAACCATAGCCAGCGGCCGATTTCGACCGAGACCGACATCACCATCAGGATCACCATATTTAGCGATGCGGAGACGGTCCCTTTCGGTAACTTATTGGAAAACAGCGTAAAGCGGAATAAGGTCGGGAAAATCAAACCAATCCCGAAAGCATACAGACTGGTGCCCAGCACCGACCACAGCCAGACGTGCGGCGACAGCAGATTGCCGACAATCAACAGCGAGAGGCCGACCAGTTGAATGGGTACGGCACGCCAGATAAACCGCGGTTCGGTCGGATCTT  >  NC_000913/4567407‑4567693                                                                                                                                              |                                                                                                                                                  gccgTTAAACCATAGCCAGCGGCCGATTTCGACCGAGACCGACATCACCATCAGGATCACCATATTTAGCGATGCGGAGACGGTCCCTTTCGGTAACTTATTGGAAAACAGCGTAAAGCGGAATAAGGTCGGGAAAATCAAGCCAATCCCg                                                                                                                                          <  2:548185/151‑1 (MQ=255)       aaaCCATAGCCAGCGGCCGATTTCGACCGAGACCGACATCACCATCAGGATCACCATATTTAGCGATGCGGAGACGGTCCCTTTCGGTAACTTATTGGAAAACAGCGTAAAGCGGAATAAGGTCGGGAAAATCAAGCCAATCCCGAAAGc                                                                                                                                     <  1:210364/150‑1 (MQ=255)             tAGCCAGCGGCCGATTTCGACCGAGACCGACATCACCATCAGGATCACCATATTTAGCGATGCGGAGACGGTCCCTTTCGGTAACTTATTGGAAAACAGCGTAAAGCGGAATAAGGTCGGGAAAATCAAGCCAATCCCGAAAGCATAc                                                                                                                                 >  2:449680/1‑148 (MQ=255)              aGCCAGCGGCCGATTTCGACCGAGACCGACATCACCATCAGGATCACCATATTTAGCGATGCGGAGACGGTCCCTTTCGGTAACTTATTGGAAAACAGCGTAAAGCGGAATAAGGTCGGGAAAATCAAGCCAATCCCGAAAGCATACAGAc                                                                                                                             <  2:112971/151‑1 (MQ=255)               gCCAGCGGCCGATTTCGACCGAGACCGACATCACCATCAGGATCACCATATTTAGCGATGCGGAGACGGTCCCTTTCGGTAACTTATTGGAAAACAGCGTAAAGCGGAATAAGGTCGGGAAAATCAAGCCAATCCCGAAAGCATACAGACt                                                                                                                            >  1:63483/1‑151 (MQ=255)                   gCGGCCGATTTCGACCGAGACCGACATCACCATCAGGATCACCATATTTAGCGATGCGGAGACGGTCCCTTTCGGTAACTTATTGGAAAACAGCGTAAAGCGGAATAAGGTCGGGAAAATCAAGCCAATCCCGAAAGCATACAGACTGGTg                                                                                                                        <  1:356967/151‑1 (MQ=255)                        cGATTTCGACCGAGACCGACATCACCATCAGGATCACCATATTTAGCGATGCGGAGACGGTCCCTTTCGGTAACTTATTGGAAAACAGCGTAAAGCGGAATAAGGTCGGGAAAATCAAGCCAATCCCGAAAGCATACAGACTGGTGccc                                                                                                                     >  2:9648/1‑149 (MQ=255)                          aTTTCGACCGAGACCGACATCACCATCAGGATCACCATATTTAGCGATGCGGAGACGGTCCCTTTCGGTAACTTATTGGAAAACAGCGTAAAGCGGAATAAGGTCGGGAAAATCAAGCCAATCCCGAAAGCATACAGACTGGTGCCCAGc                                                                                                                  >  1:446047/1‑150 (MQ=255)                             tCGACCGAGACCGACATCACCATCAGGATCACCATATTTAGCGATGCGGAGACGGTCCCTTTCGGTAACTTATTGGAAAACAGCGTAAAGCGGAATAAGGTCGGGAAAATCAAGCCAATCCCGAAAGCATACAGACTGGTGCCCAGCAcc                                                                                                               <  2:533871/150‑1 (MQ=255)                               gaccgagaccgaCATCACCATCAGGATCACCATATTTAGCGATGCGGAGACGGTCCCTTTCGGTAACTTATTGGAAAACAGCGTAAAGCGGAATAAGGTCGGGAAAATCAAGCCAATCCCGAAAGCATACAGACTGGTGCCCAGCACCGAc                                                                                                            <  1:375455/151‑1 (MQ=255)                                accgagaccgaCATCACCATCAGGATCACCATATTTAGCGATGCGGAGACGGTCCCTTTCGGTAACTTATTGGAAAACAGCGTAAAGCGGAATAAGGTCGGGAAAATCAAGCCAATCCCGAAAGCATACAGACTGGTGCCCAGCACCGAcc                                                                                                           >  2:370067/1‑151 (MQ=255)                                      accgaCATCACCATCAGGATCACCATATTTAGCGATGCGGAGACGGTCCCTTTCGGTAACTTATTGGAAAACAGCGTAAAGCGGAATAAGGTCGGGAAAATCAAGCCAATCCCGAAAGCATACAGACTGGTGCCCAGCACCGACCACAGcc                                                                                                     >  2:565849/1‑151 (MQ=255)                                           cATCACCATCAGGATCACCATATTTAGCGATGCGGAGACGGTCCCTTTCGGTAACTTATTGGAAAACAGCGTAAAGCGGAATAAGGTCGGGAAAATCAAGCCAATCCCGAAAGCATACAGACTGGTGCCCAGCACCGACCACAGCCAGACg                                                                                                <  1:6372/151‑1 (MQ=255)                                            aTCACCATCAGGATCACCATATTTAGCGATGCGGAGACGGTCCCTTTCGGTAACTTATTGGAAAACAGCGTAAAGCGGAATAAGGTCGGGAAAATCAAGCCAATCCCGAAAGCATACAGACTGGTGCCCAGCACCGACCACAGCCAGACGt                                                                                               <  2:456307/151‑1 (MQ=255)                                             tCACCATCAGGATCACCATATTTAGCGATGCGGAGACGGTCCCTTTCGGTAACTTATTGGAAAACAGCGTAAAGCGGAATAAGGTCGGGAAAATCAAGCCAATCCCGAAAGCATACAGACTGGTGCCCAGCACCGACCACAGCCAGACGTg                                                                                              >  1:388721/1‑151 (MQ=255)                                                     aGGATCACCATATTTAGCGATGCGGAGACGGTCCCTTTCGGTAACTTATTGGAAAACAGCGTAAAGCGGAATAAGGTCGGGAAAATCAAGCCAATCCCGAAAGCATACAGACTGGTGCCCAGCACCGACCACAGCCAGACGTGCGGCGa                                                                                        >  2:320025/1‑149 (MQ=255)                                                            ccATATTTAGCGATGCGGAGACGGTCCCTTTCGGTAACTTATTGGAAAACAGCGTAAAGCGGAATAAGGTCGGGAAAATCAAGCCAATCCCGAAAGCATACAGACTGGTGCCCAGCACCGACCACAGCCAGACGTGCGGCGACAGCAGAtt                                                                               >  1:101702/1‑151 (MQ=255)                                                            ccATATTTAGCGATGCGGAGACGGTCCCTTTCGGTAACTTATTGGAAAACAGCGTAAAGCGGAATAAGGTCGGGAAAATCAAGCCAATCCCGAAAGCATACAGACTGGTGCCCAGCACCGACCACAGCCAGACGTGCGGCGACAGCAGAtt                                                                               <  1:601405/151‑1 (MQ=255)                                                            ccATATTTAGCGATGCGGAGACGGTCCCTTTCGGTAACTTATTGGAAAACAGCGTAAAGCGGAATAAGGTCGGGAAAATCAAGCCAATCCCGAAAGCATACAGACTGGTGCCCAGCACCGACCACAGCCAGACGTGCGGCGACAGCAGAtt                                                                               >  1:184189/1‑151 (MQ=255)                                                             cATATTTAGCGATGCGGAGACGGTCCCTTTCGGTAACTTATTGGAAAACAGCGTAAAGCGGAATAAGGTCGGGAAAATCAAGCCAATCCCGAAAGCATACAGACTGGTGCCCAGCACCGACCACAGCCAGACGTGCGGCGACAGCAGATTg                                                                              >  2:227678/1‑151 (MQ=255)                                                                  ttAGCGATGCGGAGACGGTCCCTTTCGGTAACTTATTGGAAAACAGCGTAAAGCGGAATAAGGTCGGGAAAATCAAGCCAATCCCGAAAGCATACAGACTGGTGCCCAGCACCGACCACAGCCAGACGTGCGGCGACAGCAGATTGCCGAc                                                                         >  2:228028/1‑151 (MQ=255)                                                                        aTGCGGAGACGGTCCCTTTCGGTAACTTATTGGAAAACAGCGTAAAGCGGAATAAGGTCGGGAAAATCAAGCCAATCCCGAAAGCATACAGACTGGTGCCCAGCACCGACCACAGCCAGACGTGCGGCGACAGCAGATTGCCGACAATCa                                                                    >  1:361499/1‑150 (MQ=255)                                                                         tGCGGAGACGGTCCCTTTCGGTAACTTATTGGAAAACAGCGTAAAGCGGAATAAGGTCGGGAAAATCAAGCCAATCCCGAAAGCATACAGACTGGTGCCCAGCACCGACCACAGCCAGACGTGCGGCGACAGCAGATTGCCGACAATCAAc                                                                  >  1:505584/1‑151 (MQ=255)                                                                             gagaCGGTCCCTTTCGGTAACTTATTGGAAAACAGCGTAAAGCGGAATAAGGTCGGGAAAATCAAGCCAATCCCGAAAGCATACAGACTGGTGCCCAGCACCGACCACAGCCAGACGTGCGGCGACAGCAGATTGCCGACAATCAACAGCg                                                              <  1:398181/151‑1 (MQ=255)                                                                               gaCGGTCCCTTTCGGTAACTTATTGGAAAACAGCGTAAAGCGGAATAAGGTCGGGAAAATCAAGCCAATCCCGAAAGCATACAGACTGGTGCCCAGCACCGACCACAGCCAGACGTGCGGCGACAGCAGATTGCCGACAATCAACAGCgag                                                            <  1:228033/151‑1 (MQ=255)                                                                                   gTCCCTTTCGGTAACTTATTGGAAAACAGCGTAAAGCGGAATAAGGTCGGGAAAATCAAGCCAATCCCGAAAGCATACAGACTGGTGCCCAGCACCGACCACAGCCAGACGTGCGGCGACAGCAGATTGCCGACAATCAACAGCGAGAGGc                                                        >  1:121154/1‑151 (MQ=255)                                                                                         ttCGGTAACTTATTGGAAAACAGCGTAAAGCGGAATAAGGTCGGGAAAATCAAGCCAATCCCGAAAGCATACAGACTGGTGCCCAGCACCGACCACAGCCAGACGTGCGGCGACAGCAGATTGCCGACAATCAACAGCGAGAGGCCGAcc                                                   >  2:145423/1‑150 (MQ=255)                                                                                         ttCGGTAACTTATTGGAAAACAGCGTAAAGCGGAATAAGGTCGGGAAAATCAAGCCAATCCCGAAAGCATACAGACTGGTGCCCAGCACCGACCACAGCCAGACGTGCGGCGACAGCAGATTGCCGACAATCAACAGCGAGAGGCCGAcc                                                   <  1:426722/150‑1 (MQ=255)                                                                                         ttCGGTAACTTATTGGAAAACAGCGTAAAGCGGAATAAGGTCGGGAAAATCAAGCCAATCCCGAAAGCATACAGACTGGTGCCCAGCACCGACCACAGCCAGACGTGCGGCGACAGCAGATTGCCGACAATCAACAGCGAGAGGCCGAc                                                    <  2:265076/149‑1 (MQ=255)                                                                                              tAACTTATTGGAAAACAGCGTAAAGCGGAATAAGGTCGGGAAAATCAAGCCAATCCCGAAAGCATACAGACTGGTGCCCAGCACCGACCACAGCCAGACGTGCGGCGACAGCAGATTGCCGACAATCAACAGCGAGAGGCCGACCAGTTg                                              >  2:15382/1‑150 (MQ=255)                                                                                              tAACTTATTGGAAAACAGCGTAAAGCGGAATAAGGTCGGGAAAATCAAGCCAATCCCGAAAGCATACAGACTGGTGCCCAGCACCGACCACAGCCAGACGTGCGGCGACAGCAGATTGCCGACAATCAACAGCGAGAGGCCGACCAGTTg                                              >  2:339795/1‑150 (MQ=255)                                                                                               aaCTTATTGGAAAACAGCGTAAAGCGGAATAAGGTCGGGAAAATCAAGCCAATCCCGAAAGCATACAGACTGGTGCCCAGCACCGACCACAGCCAGACGTGCGGCGACAGCAGATTGCCGACAATCAACAGCGAGAGGCCGACCAGTTGa                                             <  2:24784/150‑1 (MQ=255)                                                                                                aCTTATTGGAAAACAGCGTAAAGCGGAATAAGGTCGGGAAAATCAAGCCAATCCCGAAAGCATACAGACTGGTGCCCAGCACCGACCACAGCCAGACGTGCGGCGACAGCAGATTGCCGACAATCAACAGCGAGAGGCCGACCAGTTGAAt                                           >  2:198738/1‑151 (MQ=255)                                                                                                aCTTATTGGAAAACAGCGTAAAGCGGAATAAGGTCGGGAAAATCAAGCCAATCCCGAAAGCATACAGACTGGTGCCCAGCACCGACCACAGCCAGACGTGCGGCGACAGCAGATTGCCGACAATCAACAGCGAGAGGCCGACCAGTTGAAt                                           <  1:145426/151‑1 (MQ=255)                                                                                                   tATTGGAAAACAGCGTAAAGCGGAATAAGGTCGGGAAAATCAAGCCAATCCCGAAAGCATACAGACTGGTGCCCAGCACCGACCACAGCCAGACGTGCGGCGACAGCAGATTGCCGACAATCAACAGCGAGAGGCCGACCAGTTGAATggg                                        >  1:75505/1‑151 (MQ=255)                                                                                                                cGTAAAGCGGAATAAGGTCGGGAAAATCAAGCCAATCCCGAAAGCATACAGACTGGTGCCCAGCACCGACCACAGCCAGACGTGCGGCGACAGCAGATTGCCGACAATCAACAGCGAGAGGCCGACCAGTTGAATGGGTACGGCACGcc                             >  2:385682/1‑149 (MQ=255)                                                                                                                              aGGTCGGGAAAATCAAGCCAATCCCGAAAGCATACAGACTGGTGCCCAGCACCGACCACAGCCAGACGTGCGGCGACAGCAGATTGCCGACAATCAACAGCGAGAGGCCGACCAGTTGAATGGGTACGGCACGCCAGATAAACCGCGGTtc             <  1:613068/151‑1 (MQ=255)                                                                                                                                 tCGGGAAAATCAAGCCAATCCCGAAAGCATACAGACTGGTGCCCAGCACCGACCACAGCCAGACGTGCGGCGACAGCAGATTGCCGACAATCAACAGCGAGAGGCCGACCAGTTGAATGGGTACGGCACGCCAGATAAACCGCGGTtcggt          >  2:12108/1‑151 (MQ=255)                                                                                                                                         aTCAAGCCAATCCCGAAAGCATACAGACTGGTGCCCAGCACCGACCACAGCCAGACGTGCGGCGACAGCAGATTGCCGACAATCAACAGCGAGAGGCCGACCAGTTGAATGGGTACGGCACGCCAGATAAACCGCGGTTCGGTCGGATCtt  <  2:534911/151‑1 (MQ=255)                                                                                                                                         aTCAAGCCAATCCCGAAAGCATACAGACTGGTGCCCAGCACCGACCACAGCCAGACGTGCGGCGACAGCAGATTGCCGACAATCAACAGCGAGAGGCCGACCAGTTGAATGGGTACGGCACGCCAGATAAACCGCGGTTCGGTCGGATCtt  >  2:488670/1‑151 (MQ=255)                                                                                                                                         aTCAAGCCAATCCCGAAAGCATACAGACTGGTGCCCAGCACCGACCACAGCCAGACGTGCGGCGACAGCAGATTGCCGACAATCAACAGCGAGAGGCCGACCAGTTGAATGGGTACGGCACGCCAGATAAACCGCGGTTCGGTCGGATCtt  >  1:299918/1‑151 (MQ=255)                                                                                                                                              |                                                                                                                                                  GCCGTTAAACCATAGCCAGCGGCCGATTTCGACCGAGACCGACATCACCATCAGGATCACCATATTTAGCGATGCGGAGACGGTCCCTTTCGGTAACTTATTGGAAAACAGCGTAAAGCGGAATAAGGTCGGGAAAATCAAACCAATCCCGAAAGCATACAGACTGGTGCCCAGCACCGACCACAGCCAGACGTGCGGCGACAGCAGATTGCCGACAATCAACAGCGAGAGGCCGACCAGTTGAATGGGTACGGCACGCCAGATAAACCGCGGTTCGGTCGGATCTT  >  NC_000913/4567407‑4567693

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