Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,220,493	C→T	intergenic (+1292/+1812)	ymgF → / ← ymgD	inner membrane protein YmgF/PF16456 family protein YmgD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,220,493	0	C	T	100.0%	33.0 / NA	11	intergenic (+1292/+1812)	ymgF/ymgD	inner membrane protein YmgF/PF16456 family protein YmgD
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (6/5);  total (6/5)

AAAAATGAATGGTGGTTCCATAACAACTAATGGCATTAATAGCTATGGTGCTTATGCTAATGGGAAAAAAGCATATATTAATTTAGATTATGTGGCACTTGAAACTGTGGCTGATGGAAGTTATGCAGTTGCTATTCGACAAGGTAACATTGATATAAAAAATAGTTCTATTACAACAACAGGCACTAAAGCCCCCATTGCAAAAATATACAATGGTGGAGAGTTATTTTTTTCCAATGTCACCGCGGTATCAAAACAAGATAAAGGAATATCAATTGATGCATCAAATATCGATTCTCAAGCCAAAATAGCACTATTAAGTGTTGAACTTTCAAGTGCTTTGGATAGTATTGATGTTAACAAAACTACAACGGATGTAAGTATCCTTAATCGAAGTATTATCACACCTGGTAATAATGTTCTGGTTAATAATACTGGAGGTGACTTA  >  NC_000913/1220260‑1220707                                                                                                                                                                                                                                          |                                                                                                                                                                                                                       aaaaaTGGATGGTGGTTCCATAACAAATAATGGCATTAATAGCTATGGTGCTTATGCTAATGGGAAAAAAGAATATATTAATTTAGATTATGTGGCACTTGAAACTGTGGCTGATGGAAGTTATGCAGTTGCTATTCGACAAGGTAACATTGATATAAAAAATAGTTCTATTACAACAACAGGCACTAAAGCCCCCATTGCAAAAATATACAATGGTGGAGAGTTATTTTTTTTCAATGTCACAGCGGTAt                                                                                                                                                                                                       >  1:192177/1‑251 (MQ=255)   aaaTGAATGGTGGTTCCATAACAACTAATGGCATTAATAGCTATGGTGCTTATGCTAATGGGAAAAAAGCATATATTAATTTAGATTATGTGGCACTTGAAACTGTGGCTGATGGAAGTTATGCAGTTGCTATTCGACAAGGTAACATTGATATAAAAAATAGTTCTATTACAACAACAGGCACTAAAGCCCCCATTGCAAAAATATACAATGGTGGAGAGTTATTTTTTTTCAATGTCACCGCGGTATc                                                                                                                                                                                                      >  2:27981/1‑250 (MQ=255)                      aacaacTAATGGCATTAATAGCTATGGTGCTTATGCTAATGGGAAAAAAGCATATATTAATTTAGATTATGTGGCACTTGAAACTGTGGCTGATGGAAGTTATGCAGTTGCTATTCGACAAGGTAACATTGATATAAAAAATAGTTCTATTACAACAACAGGCACTAAAGCCCCCATTGCAAAAATATACAATGGTGGAGAGTTATTTTTTTTCAATGTCACCGCGGTATCAAAACAAGATAAAGGAatat                                                                                                                                                                                  >  1:137906/1‑251 (MQ=255)                                             aTGGTGCTTATGCTAATGGGAAAAAAGCATATATTAATTTAGATTATGTGGCACTTGAAACTGTGGCTGATGGAAGTTATGCAGTTGCTATTCGACAAGGTAACATTGATATAAAAAATAGTTCTATTACAACAACAGGCACTAAAGCCCCCATTGCAAAAATATACAATGGTGGAGAGTTATTTTTTTTCAATGTCACCGCGGTATCAAAACAAGATAAAGGAATATCAATTGATGCATCAAATATCGAt                                                                                                                                                           >  1:158021/1‑251 (MQ=255)                                                                             atTAATTTAGATTATGTGGCACTTGAAACTGTGGCTGATGGAAGTTATGCAGTTGCTATTCCACAAGGTAACATTGATATAAAAAATAGTTCTATTACAACAACAGGCACTAAAGCCCCCATTGCAAAAATATACAATGGTGGAGAGTTATTTTTTTTCAATGTCACCGCGGTATCAAAACAAGATAAAGGAATATCAATTGATGCATCAAATATCGATTCTCAAGCCAAAATAGCACTATTAAGTGTTg                                                                                                                            >  1:155155/1‑250 (MQ=255)                                                                                               gCACTTGAAACTGTGGCTGATTGAAGTTATGCAGTTGCTATTCGGTAAGGTAACATTGATATAAAAAAGATTTCTATTACAACAACAGGCACTAAAGACCCCATTGCAAAAATATACAATGGTGGAGAGTTATTTTTTTTCAATGTCACCGCGGTATCAAAACAAGATAAAGGAATATCAATTGATGCATCAAATATCGATTCGCAAGCCAAAATAGCACTATTAAGTGTTGAACTTTCAAGTGCTTTGGa                                                                                                         <  2:149463/251‑1 (MQ=255)                                                                                                                                              aGGTAACATTGATATAAAAAATAGTTCTATTACAACAACAGGCACTAAAGCCCCCATGGCAAAAATATACAATGGTGGAGAGTTATTTTTTTTCAATGTCACCGCGGTATCAAAACAAGATAAAGGAATATCAATTGATGCATCAAATATCGATTCTCAAGCCAAAATAGCACTATTAAGTGGTGATCTTTCAAGTGCTTTGGATAGTATTGATGTTAACAAAACTACAACGGATGTAAGTATCCTTAATc                                                          <  2:137906/251‑1 (MQ=255)                                                                                                                                                                                  aacaGGCACTACAGCCCCCATTGCCAACATATACAATGGTGGAGAGTTATTTTTTTTCAATGTCACCGCGGTATCAAAACAAGCTAAAGGAATATCAATTGATGCATCAAATATCGCTTCTCAAGCCAAAATAGCACTATTAAGTGTTGAACTTTCAAGTGCTTTGGATAGTATTGa                                                                                                >  2:294003/1‑177 (MQ=255)                                                                                                                                                                                            aaaGCCCCCATTGCAAAAATATACAATGGTTGAGAGTTATTTTTTTTCAATGTCACCGCGGTATCAAAACAAGATAAAGGAATATCAATTGATGCATCAAATATCGATTCTCAAGCCAAAATAGCACTATTAAGTGTTGAACTTTCAAGTGCTTTGGATAGTATTGATGTTAACAAAACTACAACGGATGTAAGTATCCTTAATCGAAGTATTATCACACCTGTTAATAATGTTCTTGTTAATAATCCTgg            <  1:266563/251‑1 (MQ=255)                                                                                                                                                                                                    cATTGCAAAAATATACAATGGTGGAGAGTTATTTTTTTTCAATGTCACCGCGGTATCAAAACAAGATAAAGGAATATCAATTGATGCATCAAATATCGATTCTCAAGCCAAAATAGCACTATTAAGTGTTGAACTTTCAAGTGCTTTGGATAGTATTGATGTTAACAAAACTACAACGGATGTAAGTATCCTTAATCGAAGTATTATCACACCTGGTAATAATGTTCTGGTTAATAATACTGGAGGTGACt    <  1:27981/251‑1 (MQ=255)                                                                                                                                                                                                       tGCAAAAATATACAATGGTGGAGAGTTATTTTTTTTCAATGTCACCGCGGTATCAAAACAAGATAAAGGAATATCAATTGATGCATCAAATATCGATTCTCAAGCCAAAATAGCACTATTAAGTGTTGAACTTTCAAGTGCTTTGGATAGTATTGATGTTAACAAAACTACAACGGATGTAAGTATCCTTAATCGAAGTATTATCACACCTGGTAATAATGTTCTGGTTAATAATACTGGAGGTGACTTa  <  1:9209/250‑1 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                       AAAAATGAATGGTGGTTCCATAACAACTAATGGCATTAATAGCTATGGTGCTTATGCTAATGGGAAAAAAGCATATATTAATTTAGATTATGTGGCACTTGAAACTGTGGCTGATGGAAGTTATGCAGTTGCTATTCGACAAGGTAACATTGATATAAAAAATAGTTCTATTACAACAACAGGCACTAAAGCCCCCATTGCAAAAATATACAATGGTGGAGAGTTATTTTTTTCCAATGTCACCGCGGTATCAAAACAAGATAAAGGAATATCAATTGATGCATCAAATATCGATTCTCAAGCCAAAATAGCACTATTAAGTGTTGAACTTTCAAGTGCTTTGGATAGTATTGATGTTAACAAAACTACAACGGATGTAAGTATCCTTAATCGAAGTATTATCACACCTGGTAATAATGTTCTGGTTAATAATACTGGAGGTGACTTA  >  NC_000913/1220260‑1220707

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