Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,481,795	G→A	Q128* (CAG→TAG)	emrY ←	tripartite efflux pump membrane subunit EmrY

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,481,795	0	G	A	100.0%	84.0 / NA	26	Q128* (CAG→TAG)	emrY	tripartite efflux pump membrane subunit EmrY
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (23/3);  total (23/3)

CCTGGTAGATTCATTTTGACCGGTGAAGTCTCAGTTTCTCTTCCTTTAAGTAAGGTTAAGCATAATGTCAGGACGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTC‑TGGCGGATAATTCCTTAATAATAAAC‑TCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATGGCATTTGCTACACCAAACGAGGTGATAACCCAGGTGCCTTCGTCTGTTGATGCTCCCAGAA  >  NC_000913/2481551‑2482037                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   ccTGGTAGATTCATTTTGACCGGTGAAGTCTCAGTTTCTCTTCCTTTAAGTAAGGTTAAGCATAATGTCAGGACGATAATCCCCATAGGGACATTGATGAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTC‑TGGCGGATAATTCCTTAATAATAAAC‑TCTATGACa                                                                                                                                                                                                                                               <  2:94665/250‑1 (MQ=255)              ttttGACCGGTGAAGTCCCAGTTTCTCTTCCTTTAAGTAAGGTTAAGCATAATGTCAGGACGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTC‑TGGCGGATAATTCCTTAATAATAAAA‑TCTATGACAATGGAATTAACGGc                                                                                                                                                                                                                                 >  2:190931/1‑251 (MQ=255)                         gAAGTCTCAGTTTCTCTTCCTTTAAGTAAGGTTAAGCATAATGTCAGGACGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTC‑TGGCGGATAATTCCTTAATAATAAAC‑TCTATGACAGTGGAATTAACGGCCCCGCCATTa                                                                                                                                                                                                                       >  1:406983/1‑250 (MQ=255)                                        ctTCCTTTAAGTAAGGTTAAGCATAATGTCAGGACGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTC‑TGGCGGATAATTCCTTAATAATAAAC‑TCTATGACAGTGGAATTAACGGCCCCGCCCTTAACCCCTGAACGAcact                                                                                                                                                                                                       >  2:444821/1‑248 (MQ=255)                                          tCCTTTAAGTAAGGTTAAGCATAATGTCAGGACGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGCGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTC‑TGGCGGATAATTCCTTAATAATAAAC‑TCTATGACAGTGGAATTAACGGCCCCGCCAGTAACCCCTGAACGACTCTa                                                                                                                                                                                                      >  1:433381/1‑250 (MQ=255)                                                                          cGATAATCCACATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTC‑TGGCGGATAATTCCTTAATAATAAAC‑TCTATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGCTTGGTCGATAgg                                                                                                                                                                     >  2:192029/1‑251 (MQ=255)                                                                           gATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTC‑TGGCGGATAATTCCTTAATAATAAAC‑TCTATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGCCGATAGGc                                                                                                                                                                    >  2:314990/1‑251 (MQ=255)                                                                             tAATCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAATTATCACAAATATAACGGCCCAATATCGGCCCACATATCGGAGCGATAATAAAGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTC‑TGGCGGATAATTCCTTAATAATAAAC‑TCTATGACAGTGGAATTAACGGCCCCGCCCTTAACCCCTGCACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCt                                                                                                                                                                   >  1:73803/1‑250 (MQ=255)                                                                                       ggggacataaataaaaaatccccacccccAGCAAAAGTTACCAAAAATATACCCCCCAAATCCCGGCCCACATATGGGGGCGATAACCACGTTTATAGACAATAATGCCAGAG‑AAATTTTTTTTTTTCCTGGCGGATATTTCCTAAATAATAATC‑TCTAGGAAAGTGAAATAACCGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATg                                                                                                                                                        <  1:197990/226‑1 (MQ=255)                                                                                                             aaaCCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTC‑TGGCGGATAATTCCTTAATAATAAAC‑TCTATGACAGTGGAATTAACGGCCCCTCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTATTGAAGACAGCGAAAAAAACGTGAc                                                                                                                                  >  2:389095/2‑251 (MQ=255)                                                                                                                ccccAGCAAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACCTATCGGAGCGATAATCACAGTCATTGAACATAATGCCAGAGCAAATGATCTTTTTTC‑TGGAGGATACTTCATTAATAATAAAC‑TCTATGGCAGTGGAATTAAAGGCCCCGCAATTAACACCTGAACGACTCTAAAAAAAATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAACAAAAAGTGACTg                                                                                                                                >  2:369113/1‑250 (MQ=255)                                                                                                                                                      cGGCCCACATATCGGAGCGATAATCCCGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTC‑TGGCGGATAATTCCTTAATAATAAAC‑TCTATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCg                                                                                                                                               >  1:140182/1‑197 (MQ=255)                                                                                                                                                       ggCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTC‑TGGCGGATAATTCCTTAATAATAAAC‑TCTATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCa                                                                                         >  1:178812/1‑250 (MQ=255)                                                                                                                                                          caccaTATCGTAGCGATAATCACGTTCATGGCCCTAAAGGCCAGAAAAAATGTTTTTTTTTG‑TGGAGG‑TTATTCCCTAAAAAAAAACCTCAATAACGGTGAAATAAACGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTTCCAACCt                                                                                     <  1:369360/248‑1 (MQ=255)                                                                                                                                                                  cGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTCTTTC‑TGGCGGATAATTCCTTAATAATAAAC‑TCTATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTa                                                                              >  1:266522/1‑250 (MQ=255)                                                                                                                                                                                               aTGCCAGAGCAAATGTTCTTTTTTC‑TGGCGGATAATTCCTTAATAATAAAC‑TCTATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACCGGGATCGCAATGGCATTTGCTACAcc                                                 >  2:327770/1‑250 (MQ=255)                                                                                                                                                                                                     gagCAAATGTTCTTTTTTC‑TGGCGGATACTTCCTTAATAATAAAC‑TCTATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAA‑GTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAAACTGCCCATAACAGGNATCGCCATGGCATTTGCTACACCAAACgggg                                         >  2:230925/1‑248 (MQ=255)                                                                                                                                                                                                             gTTCTTTTTTC‑TGGCGGATAATTCCTTAATAATAAAC‑TCTATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATGGCATTTGCTACACCAAACGAGGTGATAAc                                  >  1:352366/1‑251 (MQ=255)                                                                                                                                                                                                                 ttttttC‑TGGCGGATAATTCCTTAATAATAAAC‑TCTATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAAc                                                                            >  1:155026/1‑205 (MQ=255)                                                                                                                                                                                                                    tttC‑TGGAGGATACTTCCTTACTAATAAAC‑TCTATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTTAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATAAATGAAGACAGCGAAAAAAAAGTGACTGACCGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCACCCTGCCAGTAACAGGGATCGCAATGGCATTAGCTACACCAAACAAGGTGATAACCCAGGtc                           >  2:155689/1‑250 (MQ=255)                                                                                                                                                                                                                       c‑TGGCGGATAATTCCTTAATAATAAAC‑TCTATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATGGCATTTGCTACACCAAACGAGGTGATAACCCAGGTGCCt                        >  1:51693/1‑251 (MQ=255)                                                                                                                                                                                                                            cGGATAATTCCTTAATAATAAAC‑TCTATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATGGCATTTGCTACACCAAACGAGGTGATAACCCAGGTGCCTTCGt                    >  2:10573/1‑251 (MQ=255)                                                                                                                                                                                                                                        taataataaAC‑TCTATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATGGCATTTGCTACACCAAACGAGGTGATAACCCAGGCGCCTTCGTCTGTTGATGCTc        >  2:37105/1‑251 (MQ=255)                                                                                                                                                                                                                                          ataataaAC‑TCTATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATGGCATTTGCTACACCAAACGAGGAGATAACCCAGGTGCCTTCGTCTGTTGATGCTccc      >  1:277296/1‑251 (MQ=255)                                                                                                                                                                                                                                            aataaAC‑TCTATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACAtt                                                                                                                                                            >  2:104723/1‑99 (MQ=255)                                                                                                                                                                                                                                               aaaC‑TCTATGACAGTGGAATTAACGGCCCAGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGCTAGGCTACACATTAATGAAGACAGCGAAAAAACAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATGGCATTTGCTACACCAAACGAGGGGAAAACCAAGGTGCCTTCGTCTGTTGATGCTCCCAGaa  >  2:305737/1‑250 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   CCTGGTAGATTCATTTTGACCGGTGAAGTCTCAGTTTCTCTTCCTTTAAGTAAGGTTAAGCATAATGTCAGGACGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTC‑TGGCGGATAATTCCTTAATAATAAAC‑TCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATGGCATTTGCTACACCAAACGAGGTGATAACCCAGGTGCCTTCGTCTGTTGATGCTCCCAGAA  >  NC_000913/2481551‑2482037

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