Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,438,207	A→G	Y860H (TAT→CAT)	ydbK ←	putative pyruvate‑flavodoxin oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,438,207	0	A	G	100.0%	45.5 / NA	14	Y860H (TAT→CAT)	ydbK	putative pyruvate‑flavodoxin oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (12/2);  total (12/2)

GCAAATTGATCCAGCAGACGCAGCACGCGGACACGGTGTTGATCGACCGTCAGGCGGAAACCAAGGCCAAATTC‑GGCATTATCTTCAAATAGAGAGTTCGCCCATGCCGGCCCACGACCGTTGGCATCGGTGGTATACGGTGTAGAGGGCAGGTTACCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGCGATCAACATCCGGTCGCCATAGAGCTGAGTCAGTAATTTAATATACGGCGTCTCGCCACAACCGGAGCAAGCACCTG  >  NC_000913/1438071‑1438338                                                                                                                                          |                                                                                                                                    gCAAATTGATCCAGCAGACGCAGCACGCGGACACGGTGTTGATCGACCGTCAGGCGGAAACCAAGGCCAAATTC‑GGCATTATCTTCAAATAGAGAGTTCGCCCATGCCGGCCCACGACCGTTGGCATCGGTGGTATGCgg                                                                                                                                  >  1:166736/1‑140 (MQ=255)     aTTGATCCAGCAGACGCAGCACGCGGACACGGTGTTGATCGACCGTCAGGCGGAAACCAAGGCCAAATTC‑GGCATTATCTTCAAATAGAGAGTTCGCCCATGCCGGCCCACGACCGTTGGCATCGGTGGTATGCGGTGTa                                                                                                                              <  2:166736/140‑1 (MQ=255)     aTTGATCCAGCAGACGCAGCACCCGGACACGTTGTCGATCGACCGTCAAGCGGAAACCAAGGCCAAATTC‑GGCATTATCTTCAAATAGAGAGTTCGCCCATGCCGGCCCACGACCGTTGGCATCGGTGGTATGCGGTGTa                                                                                                                              <  2:233671/140‑1 (MQ=255)          cccAGCGGACGCAAAAAGGGGACCCGGAGTTGATCGGCCGTCAGGCGGAAGCCAAGGCCAAATTC‑GGCATTACCTTCAAAAAGAAAGTACGCCCAAGCGGGCCCACGCCCGTTGGCACCGGTGGTATGCGGTGTaagggg                                                                                                                         <  2:426747/139‑6 (MQ=255)                                               ccGTCAGGCGGAAACCAAGGCCAAATTC‑GGCATTATCTTCAAATAGAGAGTTCGCCCATGCCGGCCCACGACCGTTGGCATCGGTGGTATGCGGTGTAGAGGGCAGGTTACCGCCATAAATTGAAGAACAGCCAGAGGCg                                                                                    >  1:547627/1‑140 (MQ=255)                                               ccGTCAGGCGGAAACCAAGGCCAAATTC‑GGCATTATCTTCAAATAGAGAGTTCGCACATGCCGGCCCACGACCGTTGGCATCGGTGGTATGCGGTGTAGAGGTCAGGTTACCGCAATCAATTGAAGAACAGCCAGTGGCg                                                                                    >  1:364073/1‑140 (MQ=255)                                                 gTCAGGCGGAAACCAAGGCCAAATTC‑GGCATTATCTTCAAATAGAGAGTTCGCCCATGCCGGCCCACGACCGTTGGCATCGGTGGTATGCGGTGTAGAGGGCCGGTTACCGCCATAAATTGAAGAACAGCCAGTGGCGtt                                                                                  >  1:402281/1‑140 (MQ=255)                                                 gTCAGGCGGAAACCAAGGCCAAATTC‑GGCATTATCTTCAAATAGAGAGTTCGCCCATGCCGGCCCACGACCGTTGGCATCGGTGGTATGCGGTGTAGAGGGCAGGTTACCGCCATAAGTTGAAGTACAGCCAGTGGCGtt                                                                                  >  1:376574/1‑140 (MQ=255)                                                         gAAACCAAGGCCAAATTC‑GGCATTATCTTCAAATAGAGAGTTCGCCCATGCCGGCCCACGACCGTTGGCATCGGTGGTATGCGGTGTAGAGGGCAGGTTACCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGCGATCa                                                                          >  1:687614/1‑140 (MQ=255)                                                                aGGCCAAATTC‑GGCATTATCTTCAAATAGAGAGTTCGCCCATGCCGGCCCACGACCGTTGGCATCGGTGGTATGCGGTGTAGAGGGCAGGTTACCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGCGATCAACATCCg                                                                   >  1:488344/1‑140 (MQ=255)                                                                 ggCCAAATTC‑GGCATTATCTTCAAATAGAGAGTTCGCCCATGCCGGCCCACGACCGTTGGCATCGGTGGTATGCGGTGTAGAGGGCAGGTTACCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGCGATCAACATCCgg                                                                  >  1:574624/1‑140 (MQ=255)                                                                  cccAAATTC‑GGCATTATCTTAAAATAGAGAGTCCCCCCAGGCCGCCCCACGCCCGTGGGCATCGGTGGTATGGGTTGTAGAGGGCGGGTTACCGCCATAAATTGAAGAACAGCCAGAGGCGTTAGCGATCAACATCCGGt                                                                 <  2:547627/139‑1 (MQ=255)                                                                   ccAAATTTGGGCTTTATCTCAAAAGAGGGAGTCCGCCCCTGCCGGCCCACGCCGGTTGGCATCGGTGGTATGGGGTGTAGAGGGCGGTTTTCCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGCGATCAACATCCGGt                                                                 <  2:679825/140‑1 (MQ=255)                                                                             gCATTATCTTCAAATAGAGAGTTCGCCCATGCCGGCCCACGACCGTTGGCATCGGTGGTATGCGGTGTAGAGGGCAGGTTACCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGCGATCAACATCCGGTCGCCATagag                                                       >  1:193139/1‑140 (MQ=255)                                                                                              agagTTCGCCCATGCCGGCCCACGACCGTTGGCATCGGTGGTATGCGGTGTAGAGGGCAGGTTACCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGCGATCAACATCCGGTCGCCATAGAGCTGAGTCAGTAATTTaa                                      >  1:688042/1‑140 (MQ=255)                                                                                                  ttCGCCCATGCCGGCCCACGACCGTTGGCATCGGTGGTATGCGGTGTAGAGGGCAGGTTACCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGCGATCAACATCCGGTCGCCATAGAGCTGAGTCAGTAATTTAatatg                                  >  1:716977/1‑139 (MQ=255)                                                                                                          tGCCGGCCCACGACCGTTGGCATCGGTGGTATGCGGTGTAGAGGGCAGGTTACCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGCGATCAACATCCGGTCGCCATAGAGCTGAGTCAGTAATTTAATATGCGGCGtct                          >  1:653999/1‑140 (MQ=255)                                                                                                                   ccgcccGTTGGCATGGGGGGTATGGGGTGAAGGGGCCGGTTTCCCGCCGTAAATTGAAGAACAGCCAGTGGCGTTAGCGATCAACATCCGGTCGCCATAGAGCTGAGTCAGTAATTTAATATGCGGCGTCTCGCCACAAc                 <  2:364073/136‑1 (MQ=255)                                                                                                                     gccggtgggCATGGGTGGTATGGGGTGTAGGGGGCAGGTTACCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGCGATCAACATCCGGTCGCCATAGAGCTGAGTCAGTAATTTAATATGCGGCGTCTCGCCACAACCg               <  2:687614/133‑1 (MQ=255)                                                                                                                            ggcttgggtggtAGGGGGTGTGGGGGGCGGGTTCCCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGCGATCAACATCCGGTCGCCATAGAGCTGAGTCAGTAATTTAATATGCGGCGTCTCGCCACAACCGGagcaag        <  2:716977/134‑1 (MQ=255)                                                                                                                                  ggtggtAGGGGGTGTAGGGGGAAGTTTACCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGCGATCAACATCCGGTCGCCATAGAGCTGAGTCAGTAATTTAATATGCGGCGTCTCGCCACAACCGGAGCAAGCACCTg  <  2:574624/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                    GCAAATTGATCCAGCAGACGCAGCACGCGGACACGGTGTTGATCGACCGTCAGGCGGAAACCAAGGCCAAATTC‑GGCATTATCTTCAAATAGAGAGTTCGCCCATGCCGGCCCACGACCGTTGGCATCGGTGGTATACGGTGTAGAGGGCAGGTTACCGCCATAAATTGAAGAACAGCCAGTGGCGTTAGCGATCAACATCCGGTCGCCATAGAGCTGAGTCAGTAATTTAATATACGGCGTCTCGCCACAACCGGAGCAAGCACCTG  >  NC_000913/1438071‑1438338

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