Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	5,234	G→A	M1M (GTG→ATG) †	yaaX →	DUF2502 domain‑containing protein YaaX

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	5,234	0	G	A	83.3%	31.9 / 0.6	19	M1M (GTG→ATG) †	yaaX	DUF2502 domain‑containing protein YaaX
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base A (7/8);  total (9/10)
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.15e-01

CAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAG‑GGAAAAAGGAGAAATTCT‑CAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCTACACGGACCGCCGCCACCGCCGCGCCACCATA  >  NC_000913/5015‑5474                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                 caGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GTCTTTTTTATGAAGAAATTATGGAGAAAAACGACAG‑GGAAAAAGTAGAAATTCT‑CAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTCTGAAAAAGATGCAATCTATCGTACTCGCACtt                                                                                                                                                                                                                  <  1:213204/251‑1 (MQ=255)          cTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATCCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAG‑GGAAAAAGGAGAAATTCT‑CAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGt                                                                                                                                                                                                          >  1:22063/1‑251 (MQ=255)                 gtATCCCAATTAGGCCGGGTTTGCCGTTATGTA‑TCCAGTTTTTATATGAAGAAATTATGGAGAAAAACGA‑AGTGTAAAAAGGA‑AAATTTTGCAATAATTGCGGTAACGTCGAGATTAAGACTGAGGATAATAACAACCGGCGTTGTCATTTCGGAATCTCCTGATAGCGTGGCGTATCGTGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATATAGCGTACTCTCACTTTCCCTGGTTCtt                                                                                                                                                                                                      <  2:220562/246‑2 (MQ=255)                         aTCAGGCCGGGTTTGCTTTTGTGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAG‑GGAAAAAGGAGAAATTCT‑CAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATg                                                                                                                                                                                          <  1:10166/251‑1 (MQ=255)                                                              tATGAAGAAATTATGGAGAAAAACGACAG‑GGAAAAAGGAGAAATTCT‑CAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTc                                                                                                                                                       <  1:211156/250‑1 (MQ=255)                                                                     aaaTTATGGAGAAAAACGACAG‑GGAAAAAGGAGAAATTCT‑CAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAAt                                                                                                                                               <  1:85658/251‑1 (MQ=255)                                                                          aTGGAGAAAAATGACAG‑GGAAAAAGGAGAAATTCT‑CAATAAGTGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTACGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCACTGGTACTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATACCGTTAGTCCAGTAAGTAAAATTAc                                                                                                                                            >  2:155065/1‑249 (MQ=255)                                                                               gaAAAACGACAG‑GGAAAAAGGAGAAATTCT‑CAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGGCCCATAACGGGCAATGATAAAAGGAGTAACCTATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAgg                                                                                                                                     <  1:35345/251‑1 (MQ=255)                                                                                                           tct‑cAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTATGAAAAAGATGCAAGCTATCGTTCTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACt                                                                                                           <  1:380565/250‑1 (MQ=255)                                                                                                            ct‑cAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTgg                                                                                                         >  1:14131/1‑251 (MQ=255)                                                                                                                                         ggATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACACATAACGGGCAATGATAAAAGGAGTAACCTATGAAAAAGATGCAATCTATCGTACTCGCACTTTCACTGGTTCTGGTCGCTCCCATGGCAGAACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATAGTGGCTATTACTGGGATGGAGGTCACAGGCGCGACCACGGCt                                                                             >  1:72185/1‑251 (MQ=255)                                                                                                                                          gtgtGCGGAGCATAACAACCGGCGTTCTCATCGCGTAATCTCCAGATATCGACCCATAACGGGCAATCAGAAAAGGAGTAACCTATGAAAAAGATGCAATCGATCGTACTCGCACTTTCCCTGTTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGTCGATCGTTATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCtg                                                                            <  2:197684/248‑1 (MQ=255)                                                                                                                                             tGCGGAGAATAACAACCGTCGTTCTCATCGCGAAATCTCCGGATATCGAACCATAACGGGCAATGATAAAAGGAGTAACCTATGAAAAATATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTAGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTAAGTAAAATTACAGATATGCGATCGTGATAACAGTGGCGATTAATGGGATGGAGTTCACTGGCGCGACCACGGCtgggg                                                                         >  2:274553/1‑249 (MQ=255)                                                                                                                                                       aacaacCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACAt                                                               <  1:88561/251‑1 (MQ=255)                                                                                                                                                                       aTCGCGTAATATCCGGATATCGACACATAACGGGCCATGATCAAAGGAGTAACCTATGAAACAGATGCAATCTATCGTACTCGCACTTTTCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGAAAAATTACAGATAGGCGATCGTGCTCATCGTTGCTATTTCTGTAATGGAGGTCACTGGCGCGACCACGGCTGGTGGGAACAACATTATAAATgactagga                                                >  1:206929/1‑243 (MQ=255)                                                                                                                                                                         cgcgTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAAt                                             >  1:15459/4‑251 (MQ=255)                                                                                                                                                                                   ccGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTgg                                       <  1:374992/247‑1 (MQ=255)                                                                                                                                                                                                   aaCGGGCAATGATAAAAGGAGTAACCTATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCCACACGGACCgccgcc                   >  1:78663/1‑251 (MQ=255)                                                                                                                                                                                                                      aGTAACCTATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATa  >  1:334187/1‑249 (MQ=255)                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                 CAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAG‑GGAAAAAGGAGAAATTCT‑CAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGTCAGTAAAATTACAGATAGGCGATCGTGATAATCGTGGCTATTACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCTACACGGACCGCCGCCACCGCCGCGCCACCATA  >  NC_000913/5015‑5474

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