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	96.0%	70.5 / NA	25	Q128* (CAG→TAG)	emrY	tripartite efflux pump membrane subunit EmrY
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (19/5);  minor base C (1/0);  total (20/5)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

CTGGTAGATTCATTTTGACCGGTGAAGTCTCAGTTTCTCTTCCTTTAAGTAAGGTTAAGCATAATGTCAGGACGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCT‑GTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATGGCATTTGCTACACCAAACGAGGTGATAACCCAGGTGCCTTCGTCTGTTGATGCTCCCAGAA  >  NC_000913/2481552‑2482037                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   cTGGTAGATTCATTTTGACCGGTGAAGTCTCAGTTTCTCTTCCTTTAAGGAAGGATAAGCATAATGTCAGGACGATAATCCCCCTAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCT‑ATGACAGt                                                                                                                                                                                                                                             >  1:323023/1‑251 (MQ=255)      aGATTCATTTTGACCGGTGAAGTCTCAGTTTCTCTTCCTTTAAGTAAGGTTAAGCATAATGTCAGGACGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCT‑ATGACAGTGGaa                                                                                                                                                                                                                                         <  1:263363/250‑1 (MQ=255)                 gACCGGTGAAGTCTCAGTTTCTCTTCCTTTAAGTAAGGTTAAGCATAATGTCAGGACGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCT‑ATGACAGTGGAATTAACGGcccc                                                                                                                                                                                                                              <  1:347184/250‑1 (MQ=255)                                 gTTTCTCTTCCTTTAAGTAAGGTTAAGCATAATGTCAGGACGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCT‑ATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAAc                                                                                                                                                                                                             >  1:558090/1‑251 (MQ=255)                                   ttCTCTTCCTTTAAGTAAGGTTAAGCATAATGTCAGGACGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCT‑ATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACg                                                                                                                                                                                                            >  1:129113/1‑250 (MQ=255)                                         cccTTTAAGAAAGGTTAAGCATAATGTCAGGACAATAATCCCCAGGGGGAAATTAATAAAAAATCTCCACCCCCAGAAAAAGTTACCAAAATTACACCCCCCAAATATCGCCCCACATAGCGGAGCAATAATCCCGGTCATTGACAATAATGCCAGAGAAAATTTTTTTTTTTCGGGCGAATATTTCCTTAAAAAAAACCTCT‑ATGACAGGGAAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTa                                                                                                                                                                                                      <  1:161925/249‑1 (MQ=255)                                                                          gATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCT‑ATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCCCATCAAGATTGGTCGATAGGc                                                                                                                                                                    >  2:182265/1‑251 (MQ=255)                                                                                                    aaaaTATCCACCCCCAGCTAAAGTTATCAAAAATAAACCCCCCCAATACCGCCCAACAAACCGGAGCGATAATCCCGGTCATGGACAATAATGCAAGAAAAAATGTTTTTTTTTGTGGCGAATATTTCTTAAATAAAAACCTCA‑ATGACGGTGAAATTACCGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGaaaaa                                                                                                                                          <  1:436465/251‑1 (MQ=255)                                                                                                            cAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCT‑ATGACAGTGGAATTACCGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGa                                                                                                                                   >  2:93825/1‑250 (MQ=255)                                                                                                                         agatATCAAAAATATACCAGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGACAGAGCAAATGATCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCT‑ATGACAGTGGAATTAACGGCACCGCAATTAACCCCTGAACGACTCTAAAAAATATCAGCACCACAAGCTTGGTCGATAGGCTCCAAATTACAGAAGACAGCGAACAAAAAGGGACATAAAGaaaaaa                                                                                                                      >  2:494028/4‑250 (MQ=255)                                                                                                                                   aaTATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCT‑ATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAAtt                                                                                                           >  2:131596/1‑251 (MQ=255)                                                                                                                                          ccGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCT‑ATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGGGACTGAAAGTAAAAAAAATCTTAATTAGCCTAt                                                                                                    >  2:123379/1‑251 (MQ=255)                                                                                                                                             cccAATATCCGCCCACATAGTGGAGTGATAATCACGTTTAGGCACCAAAAGGCCAGAAAAAGTTTTTTTTTTTCGGGGGGATAATTCCTTAA‑AAAAAACCCTTAGGACAGGGAAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATGCGCCTATTCt                                                                                                 <  1:463308/251‑1 (MQ=255)                                                                                                                                                                                ggTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCT‑ATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGGGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCAAGTAACAGGGATCGCAATgg                                                              >  1:47115/1‑251 (MQ=255)                                                                                                                                                                                 gTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAACCACT‑CTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGCCTCTAAAACATATCAGCAATTCAAGCTTGCTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGCCTCAAAGGCAAAATAATCTTAATTCTCCTATTCTTTGTGCCAACCTGCCCGTAACAGGGAACGCAATGGc                                                             >  2:205810/1‑251 (MQ=255)                                                                                                                                                                                    aTTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCT‑ATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATGGCAtt                                                          >  1:221333/1‑251 (MQ=255)                                                                                                                                                                                         ccATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCT‑ATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATGGCATTTGCt                                                      >  1:178214/1‑250 (MQ=255)                                                                                                                                                                                                     agCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCT‑ATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCAAACCTGCCAGTAAAAGGGATCGCAATGGCATTTGCTACACCAAACGAgg                                         >  1:320135/1‑251 (MQ=255)                                                                                                                                                                                                      gCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCT‑ATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATGGCATTTGCTACACCAAACGAGGt                                        >  1:377005/1‑251 (MQ=255)                                                                                                                                                                                                                   tttCTGGCGGATAATTCCTTAATAATAAACTCT‑ATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATGGCATTTGCTACACCAAACGAGGTGATAACCCAGGTg                           >  1:496827/1‑251 (MQ=255)                                                                                                                                                                                                                      atGGCGGATAATTCCTTAATAATAAACTCT‑ATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCACGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATGGCATTTGATACACCACACGAGGTGATAACCCAGGTGCCt                        >  2:185828/2‑251 (MQ=255)                                                                                                                                                                                                                             aTAATTCCTTAATATTAAACTCT‑ATGACAGCGGACTTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGATACAAAAGTGACTGAAAGTaa                                                                                                                         >  1:91896/1‑147 (MQ=255)                                                                                                                                                                                                                                   ccTTAATAATAAACTCT‑ATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATGGCATTTGCTACACCAAACGAGGTGATAACCCAGGTGCCTTCGTCTGTTGATg           >  1:477969/1‑251 (MQ=255)                                                                                                                                                                                                                                          aataaACTCT‑ATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGAc                                                                                                                                  >  1:318588/1‑125 (MQ=255)                                                                                                                                                                                                                                              aaCTCT‑ATGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATGGCATTTGCTACACCAAACGAGGTGATAACCCAGGTGCCTTCGTCTGTTGATGCTCCCAGaa  >  1:347768/1‑249 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   CTGGTAGATTCATTTTGACCGGTGAAGTCTCAGTTTCTCTTCCTTTAAGTAAGGTTAAGCATAATGTCAGGACGATAATCCCCATAGGGACATTGATTAAAAATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCT‑GTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATGGCATTTGCTACACCAAACGAGGTGATAACCCAGGTGCCTTCGTCTGTTGATGCTCCCAGAA  >  NC_000913/2481552‑2482037

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