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%	123.1 / NA	37	Q128* (CAG→TAG)	emrY	tripartite efflux pump membrane subunit EmrY
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (27/10);  total (27/10)

CCTGGTAGATTCATTTTGACCGGTGAAGTCTCAGTTTCTCTTCCTTTAAGTAAGGTTAAGCATAATGTCAGGACGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCC‑AGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATG‑CCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATGGCATTTGCTACACCAAACGAGGTGATAACCCAGGTGCCTTCGTCTGTTGATGCT  >  NC_000913/2481551‑2482030                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                            ccTGGTAGATTCATTTTGACCGGTGAAGTCTCAGTTTCTCTTCCTTTAAGTAAGGTTAAGCATAATGTCAGGACGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCC‑AGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATG‑CCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTATGACa                                                                                                                                                                                                                                        <  1:647557/250‑1 (MQ=255)    ggTAGATTCAATTTGACCGGTGAAGTCTCAGTTTCTCTTCCTTTAAGTAAGGTTAAGCATAATGTCAGGACGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCC‑AGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATG‑CCAGAGCAAATGTTCTTTTTTCTGGCAGATAATTCCTTAACAATAAACTCTATGACAGt                                                                                                                                                                                                                                      <  2:600680/249‑1 (MQ=255)            cATTTTGACCGGTGAAGTCTCAGTTTCTCTTCCTTTAAGTAAGGTTAAGCATAATGTCAGGACGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCC‑AGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATG‑CCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTATGACAGTGGAATTAAc                                                                                                                                                                                                                             <  1:219198/250‑1 (MQ=255)                  gACCGGTGAAGTCTCAGTTTCTCTTCCTTTAAGTAAGGTTAAGCATAATGTCAGGACGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCC‑AGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATG‑CCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTATGACAGTGGAATTAACGGCCCCg                                                                                                                                                                                                                      <  1:162343/251‑1 (MQ=255)                    ccGGTGAAGTCTCAGTTTCTCTTCCTTTAAGNAAGGTTAAGCATAATGTCAGGACGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCC‑AGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATG‑CCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTATGACAGTGGAATTAACGGCCCCGcc                                                                                                                                                                                                                    <  2:337977/251‑1 (MQ=255)                      ggTGAAGTCTCAGTTTCTCTTCCTTTAAGTAAGGTTAAGCATAATGTCAGGACGATAATCCCCATAGGGAAATTGATTAAAAATATCCAACCCC‑AGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATG‑CCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTATGACAGTGGAATTACCGGCACCGcctt                                                                                                                                                                                                                  >  2:690388/1‑249 (MQ=255)                                      ctctTCCTTTAAGTAAGGTTAAGCATAATGTCAGGACGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCC‑AGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATG‑CCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAAAAATAAACTATATAACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGAcc                                                                                                                                                                                                  >  2:549595/1‑250 (MQ=255)                                      ctctTCCTTTAAGTAAGGTTAAGCATAATGTCAGGACGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCC‑AGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGACATTGACCATAATG‑CCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTATGACAGTGGAATTAACGGCCCCGACATTAACCCCTAAACAAct                                                                                                                                                                                                  >  2:454366/1‑251 (MQ=255)                                           cctgtaagtaagGTTAAGAAAAATGTAAGGACAAAAACCCCCATAGGAAAATTGATTAAAAA‑ACCCAACCCCCAGCAAAAGTTACCACAAAAATAACCCCCCAATTTCGGCCCACATACCGGGACGATAATCACGGTCATTCACAATAAGG‑CCAGAGAAATTGTTTTTTTTTCTGGCGGATATTTCCTAAATAAAAACCCCTAGGACAGTGGATTTAGCGGCCCCGCCATTAACCCCTGAACGACTCTaaa                                                                                                                                                                                             <  1:582332/247‑1 (MQ=255)                                              ttAAGTAAGGTTAAGCATAATGTCAGGACGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCC‑AGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATG‑CCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTATGACAGTGGAATTAACGACCCCGCCATTAACCACTGAACGACTCTaaaaa                                                                                                                                                                                           >  1:463940/1‑250 (MQ=255)                                                        ttAAGCATAATGTCAGGACGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCC‑AGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATG‑CCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAAAGACTCTAAAAAATATCAGcac                                                                                                                                                                                 >  1:26302/1‑250 (MQ=255)                                                         aaagaaTAATCTCAGGACAAAAACCCCCAAAGGGAAATTAATTAAAAATCTCCCCCCCC‑AGCTAAATTTATCAAAAAAAAACCCCCCCAATACCGCCCCACATATCGGAGATATATTCACGGTTATTGACCA‑AATGACAAGAGCAAGTTTTTTTTTTTCTGGCGAAAAATCCTTTAATAATAAACTCTAGGACAGAGAAATGAACGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACAt                                                                                                                                                                               <  2:410751/246‑1 (MQ=255)                                                              aTAATGTCAGGACGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCC‑AGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATG‑CCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAACCTCTATGACa                                                                                                                                                                                                                                        >  2:435375/1‑189 (MQ=255)                                                              aTAATGTCAGGACGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCC‑AGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATG‑CCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTATGACa                                                                                                                                                                                                                                        <  1:435375/189‑1 (MQ=255)                                                                  tGTCAGGACGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCC‑AGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATG‑CCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTgg                                                                                                                                                                      >  1:128316/1‑251 (MQ=255)                                                                                    cATAGGGACATTGATTAAAAATCTCCAACCCC‑CGCTAACGTTATCACAAATCTAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATG‑CCAGAGCAACTGTTCTTTTTTCTGGCGGATAATTCCTTAATAATACACTCTATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACCTCAAGATTGGTCATTAGGCTACACAtt                                                                                                                                                     >  2:239569/1‑250 (MQ=255)                                                                                                       aaTATCCAACCCC‑AGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATG‑CCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTATGACAGTGGaa                                                                                                                                                                                                                                  >  2:352258/1‑154 (MQ=255)                                                                                                       aaTATCCAACCCC‑AGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATG‑CCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTATGACAGTGGaa                                                                                                                                                                                                                                  <  1:352258/154‑1 (MQ=255)                                                                                                                cccc‑AGCTAAAGTTATCACAAATATAACCGCACAATATCCGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATG‑CCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGCCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGGGACTg                                                                                                                         >  2:409910/1‑250 (MQ=255)                                                                                                                      gCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATG‑CCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGt                                                                                                                    >  1:416785/1‑250 (MQ=255)                                                                                                                         aaaGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATG‑CCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCCAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTaaa                                                                                                                 >  2:500190/1‑250 (MQ=255)                                                                                                                           aGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATG‑CCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGGGACTGAAAGTAAAaat                                                                                                              >  1:260827/1‑251 (MQ=255)                                                                                                                               aTCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATG‑CCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATc                                                                                                          >  1:298080/1‑251 (MQ=255)                                                                                                                                 caaaaaTATAACCGCCTAATAGCGGCCAAAATATGGGGGCAAAAACCCGGGTTATGGACCATAAGG‑CCAGAGAAATTTTTCTTTTTTTGGGCGGATATTTCTTTAAAAAAAACCACTATGAGAGGGAAATAAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCtt                                                                                                        <  1:540111/248‑1 (MQ=255)                                                                                                                                                ccAATATAGGCCTACATATCGAAGAGATAATCACGGGCATTGAATATAATG‑CCAGAGCCCATGTTTTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTACGACAGTGGAATTAACGGCCCACCAATTAACCCCTGCAAGACTCTAAAAAATATCAGCCCATCAAGATTGGTCGATCGGCTAACCATTAATGAAGACAGCGAACATACAGAGACTGAAAGTAAAAATAATATGAATTCGCCTAttatt                                                                                         >  2:69478/1‑248 (MQ=255)                                                                                                                                                            cacaTATCGGAGCGATAATCACGGTCATTGACCATAATG‑CCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGCACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCTCCTATTCTTTGTGCCAACCTg                                                                             >  2:10895/1‑251 (MQ=255)                                                                                                                                                                     gAGCGATAATCACGGTCATTGACCATAATG‑CCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTNGCCTATTCTTTGTGCCAACCTGCCAGTAAc                                                                     >  2:357475/1‑250 (MQ=255)                                                                                                                                                                      aGCGATAATCACGGTCATTGACCATAATG‑CCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGAGCCAACCTGCCAGTAAc                                                                     >  2:420269/1‑249 (MQ=255)                                                                                                                                                                                aCGGTCATTGACCATAATG‑CCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAAt                                                         >  1:518241/1‑251 (MQ=255)                                                                                                                                                                                      cttGACGATAATG‑CCAGAGCAAATGTTCTTTTTTCTGGCGAATAATTCCTTACTAATAAACTATATGAAAATGTAATTAACGGCCCCGCCATTTAAACCTGAACGACTCTAAAAAATATCAGCACATCACGATTGGTCGCTAGGCTCCACATTAATGCAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTCGAGCCCACCTGCCAGTATAAGGGATAGCAATAGCAtt                                                   >  2:88489/2‑251 (MQ=255)                                                                                                                                                                                            cATAATG‑CCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAACACATCAAAATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTAACAACCTGCCAGTAACAGGGATCGCAATGTCAGTTGCtcc                                             >  1:1496/1‑249 (MQ=255)                                                                                                                                                                                                           aaaTGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATGGCATTTGCTACACCAAACGAGGTg                                >  1:475020/1‑250 (MQ=255)                                                                                                                                                                                                                   ttttttCTGGCGGATAATTCCTTAATAATAAACTCTATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATGGCATTTGCTACACCAAACGAGGTGATAAccc                         >  1:228545/1‑249 (MQ=255)                                                                                                                                                                                                                        tCTGGCGGATAATTCCTTAATAATAAACTCTATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATGGCATTTGCCACACCAAACGAGGTGATAACCCAGGTGcc                  >  1:77588/1‑251 (MQ=255)                                                                                                                                                                                                                                aTAATTCCTTAATAATAAACTCTATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATGGCATTTGCTACACCAAACGAGGTGATAACCCAGGTGCCTTCGTCTg          >  1:230333/1‑251 (MQ=255)                                                                                                                                                                                                                                       cTTAATAATAAACTATATGACAGTGGAATTCACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGAAAGCGAAAAAAAAGTGACTGAAAGTAAACATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATGGCATTTGCTACACCAAACGAGGTGATAACCCAGGTGCCTTCGTCTGTTGATGc   >  2:88880/1‑251 (MQ=255)                                                                                                                                                                                                                                        ttAATAATAAACTCTATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGAACCAGGGATCGCAATGGCATTTGCTACACCAAACGAGGTGATAACCCAGGTGCCTTCGTCTGTTGTTGcc  >  2:662826/1‑250 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                            CCTGGTAGATTCATTTTGACCGGTGAAGTCTCAGTTTCTCTTCCTTTAAGTAAGGTTAAGCATAATGTCAGGACGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCC‑AGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATG‑CCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATGGCATTTGCTACACCAAACGAGGTGATAACCCAGGTGCCTTCGTCTGTTGATGCT  >  NC_000913/2481551‑2482030

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: ‑