Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,301,528	+C	coding (1551/2052 nt)	mdtO ←	putative multidrug efflux pump subunit MdtO

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,301,525	1	.	C	100.0%	119.0 / NA	35	G518G (GGT→GGG)	mdtO	putative multidrug efflux pump subunit MdtO
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (15/20);  total (15/20)

AGACCGATACGAATTTGCAGATAAGTGCGCAGAGCCGTGACTTCCTGCTGGCGTGGGATCCGCATTACTTTACTTAACATA‑CCCAGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCTTCACTTTCAGGCCAGACAAAGGTGTAAATCACCGC  >  NC_000913/4301445‑4301599                                                                                  |                                                                           aGACCGATACGAATTTGCAGATAAGTGCGCAGAGCCGTGACTTCCTGCTGGCGTGGGATCCGCATTACTTTACTTAACATa                                                                             >  1:882263/1‑81 (MQ=255)       aTACGAATTTGCAGATAAGTGCGCAGAGCCGTGACTTCCTGCTGGCGTGGGATCCGCATTACTTTACTTAACATA‑ccccag                                                                      <  1:882219/81‑4 (MQ=255)       aTACGAATTTGCAGATAAGTGCGCAGAGCCGTGACTTCCTGCTGGCGTGGGATCCGCATTACTTTACTTAACATA‑ccccag                                                                      <  2:882218/81‑4 (MQ=255)       aTACGAATTTGCAGATAAGTGCGCAGAGCCGTGACTTCCTGCTGGCGTGGGATCCGCATTACTTTACTTAACATA‑ccccag                                                                      <  2:882217/81‑4 (MQ=255)        tACGAATTTGCAGATAAGTGCGCAGAGCCGTGACTTCCTGCTGGCGTGGGATCCGCATTACTTTACTTAACATA‑ccccagc                                                                     <  2:882220/81‑5 (MQ=255)             aTTTGCAGATAAGTGCGCAGAGCCGTGACTTCCTGCTGGCGTGGGATCCGCATTACTTTACTTAACATACCCCAGCGTgcc                                                                 <  1:882222/81‑1 (MQ=255)              tttGCAGATAAGTGCGCAGAGCCGTGACTTCCTGCTGGCGTGGGATCCGCATTACTTTACTTAACATACCCCAGCGTgc                                                                  >  2:882282/1‑79 (MQ=255)                tGCAGATAAGTGCGCAGAGCCGTGACTTCCTGCTGGCGTGGGATCCGCATTACTTTACTTAACATACCCCAGCGTgccagc                                                              >  2:882278/1‑81 (MQ=255)                  cAGATAAGTGCGCAGAGCCGTGACTTCCTGCTGGCGTGGGATCCGCATTACTTTACTTAACATACCCCAGCGTgccagcc                                                             <  2:882223/80‑1 (MQ=255)                    gATAAGTGCGCAGAGCCGTGACTTCCTGCTGGCGTGGGATCCGCATTACTTTACTTAACATACCCCAGCGTGCCAGCCAGt                                                          >  1:882287/1‑81 (MQ=255)                     aTAAGTGCGCAGAGCCGTGACTTCCTGCTGGCGTGGGATCCGCATTACTTTACTTAACATACCCCAGCGTGCCAGCCAGtt                                                         <  1:882224/81‑1 (MQ=255)                      tAAGTGCGCAGAGCCGTGACTTCCTGCTGGCGTGGGATCCGCATTACTTTACTTAACATACCCCAGCGTGCCAGCCAGttt                                                        <  1:882225/81‑1 (MQ=255)                        aGTGCGCAGAGCCGTGACTTCCTGCTGGCGTGGGATCCGCATTACTTTATTTAACATACCCCAGCGTGCCAGCCAGttttt                                                      >  2:882277/1‑81 (MQ=255)                        aGTGCGCAGAGCCGTGACTTCCTGCTGGCGTGGGATCCGCATTACTTTACTTAACATACCCCAGCGTGCCAGCCAGttttt                                                      <  2:882226/81‑1 (MQ=255)                         gTGCGCAGAGCCGTGACTTCCTGCTGGCGTGGGATCCGCATTACTTTACTTAACATACCCCAGCGTGCCAGCCAGTTTTTg                                                     >  1:882301/1‑81 (MQ=255)                          tGCGCAGAGCCGTGACTTCCTGCTGGCGTGGGATCCGCATTACTTTACTTAACATACCCCAGCGTGCCAGCCAGTTTTTGc                                                    <  1:882227/81‑1 (MQ=255)                           gcgcAGAGCCGTGACTTCCTGCTGGCGTGGGATCCGCATTACTTTACTTAACATACCCCAGCGTGCCAGCCAGTTTTTGCg                                                   >  1:882283/1‑81 (MQ=255)                            cgcAGAGCCGTGACTTCCTGCTGGCGTGGGATCCGCATTACTTTACTTAACATACCCCAGCGTGCCAGCCAGTTTTTGCgg                                                  >  1:882280/1‑81 (MQ=255)                               agagCCGTGACTTCCTGCTGGCGTGGGATCCGCATTACTTTACTTAACATACCCCAGCGTGCCAGCCAGTTTTTGCGGCAg                                               <  2:882228/81‑1 (MQ=255)                               agagCCGTGACTTCCTGCTGGCGTGGGATCCGCATTACTTTACTTAACATACCCCAGCGTGCCAGCCAGTTTTTGCGGCAg                                               <  2:882229/81‑1 (MQ=255)                                gagCCGTGACTTCCTGCTGGCGTGGGATCCGCATTACTTTACTTAACATACCCCAGCGTGCCAGCCAGTTTTTGCGGCAgt                                              >  2:882289/1‑81 (MQ=255)                                gagCCGTGACTTCCTGCTGGCGTGGGATCCGCATTACTTTACTTAACATACCCCAGCGTGCCAGCCAGTTTTTGCGGCAgt                                              <  1:882230/81‑1 (MQ=255)                                          ttCCTGCTGGCGTGGGATCCGCATTACTTTACTTAACATACCCCAGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCtt                                    >  2:882294/1‑81 (MQ=255)                                            ccTGCTGGCGTGGGATCCGCATTACTTTACTTAACATACCCCAGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCtt                                    <  2:882231/79‑1 (MQ=255)                                                ctgGCGTGGGATCCGCATTACTTTACTTAACATACCCCAGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCTTCACttt                              >  2:882275/1‑81 (MQ=255)                                                ctgGCGTGGGATCCGCATTACTTTACTTAACATACCCCAGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCTTCACttt                              <  2:882233/81‑1 (MQ=255)                                                ctgGCGTGGGATCCGCATTACTTTACTTAACATACCCCAGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCTTCACttt                              <  1:882232/81‑1 (MQ=255)                                                       gggATCCGCATTACTTTACTTAACATACCCCAGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCTTCACTTTCAGGc                         >  2:882296/1‑79 (MQ=255)                                                         gATCCGCATTACTTTACTTAACATACCCCAGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCTTCACTTTCAGGCCAGa                     >  2:882291/1‑81 (MQ=255)                                                           tCCGCATTACTTTACTTAACATACCCCAGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCTTCACTTTCAGGCCAGACa                   <  1:882234/81‑1 (MQ=255)                                                            ccGCATTACTTTACTTAACATACCCCAGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCTTCACTTTCAGGCCAGACa                   >  2:882286/1‑80 (MQ=255)                                                               cATTACTTTACTTAACATACCCCAGCGTGCCGGCCAGTTTTTGCGGCAGTGTGCGCGCTTCACTTTCAGGCCAGACAAAgg               <  2:882235/81‑1 (MQ=255)                                                                aTTACTTTACTTAACATACCCCAGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCTTCACTTTCAGGCCAGACAAAGgt              <  2:882237/81‑1 (MQ=255)                                                                aTTACTTTACTTAACATACCCCAGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCTTCACTTTCAGGCCAGACAAAGgt              <  2:882236/81‑1 (MQ=255)                                                                aTTACTTTACTTAACATACCCCAGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCTTCACTTTCAGGCCAGACAAAGgt              >  2:882284/1‑81 (MQ=255)                                                                   actttacttAACATACCCCAGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCTTCACTTTCAGGCCAGACAAAGgtgt            <  2:882238/80‑1 (MQ=255)                                                                       tacttAACATACCCCAGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCTTCACTTTCAGGCCAGACAAAGGTGTAAAt        <  1:882239/80‑1 (MQ=255)                                                                        acttAACATACCCCAGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCTTCACTTTCAGGCCAGACAAAGGTGTAAATc       <  1:882240/80‑1 (MQ=255)                                                                           tAACATACCCCAGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCTTCACTTTCAGGCCAGACAAAGGTGTAAATCACcg   >  1:882281/1‑81 (MQ=255)                                                                             aCATACCCCAGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCTTCACTTTCAGGCCAGACAAAGGTGTAAATCACcgc  <  1:882241/80‑1 (MQ=255)                                                                                  |                                                                           AGACCGATACGAATTTGCAGATAAGTGCGCAGAGCCGTGACTTCCTGCTGGCGTGGGATCCGCATTACTTTACTTAACATA‑CCCAGCGTGCCAGCCAGTTTTTGCGGCAGTGTGCGCGCTTCACTTTCAGGCCAGACAAAGGTGTAAATCACCGC  >  NC_000913/4301445‑4301599

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