Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	729,582	G→A	intergenic (+242/‑1)	ybfA → / → rhsC	DUF2517 domain‑containing protein YbfA/rhs element protein RhsC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	729,582	0	G	A	100.0%	61.5 / NA	20	intergenic (+242/‑1)	ybfA/rhsC	DUF2517 domain‑containing protein YbfA/rhs element protein RhsC
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (15/5);  total (15/5)

ATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGC  >  NC_000913/729450‑729630                                                                                                                                     |                                                 atatTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAG‑AAATGAGCgg                                          >  1:933795/1‑140 (MQ=255)         gAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAG‑AAATGAGCGGAAAACcgg                                  >  2:722097/1‑140 (MQ=255)         gAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAG‑AAATGAGCGGAAAACcgg                                  >  1:219118/1‑140 (MQ=255)         gAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAG‑AAATGAGCGGAAAACcgg                                  >  1:911514/1‑140 (MQ=255)          aaaTATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAG‑AAATGAGCGGAAAACcggc                                 >  2:273530/1‑140 (MQ=255)            atatatatAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAG‑AAATGAGCGGAAAACcggcag                               >  2:323055/1‑138 (MQ=255)                       tatTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAG‑AAATGAGCGGAAAACCGGCAGCGCGTCAGGGc                    >  2:330422/1‑140 (MQ=255)                         ttttCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAG‑AAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGa                  >  1:1341000/1‑140 (MQ=255)                         ttttCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAG‑AAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGa                  >  2:1341704/1‑140 (MQ=255)                           ttCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTGTAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAG‑AAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACa                <  2:382551/140‑1 (MQ=14)                           ttCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAG‑AAATGAGCGGAAAACCGGCAGCGCGTCAggg                     <  1:956879/135‑1 (MQ=255)                           ttCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAG‑AAATGAGCGGAAAACCGGCAGCGCGTCAggg                     >  2:956879/1‑135 (MQ=255)                           ttCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAG‑AAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACa                <  2:1422056/140‑1 (MQ=255)                           ttCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAG‑AAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACa                >  2:529071/1‑140 (MQ=255)                           ttCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAG‑AAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACa                >  2:980993/1‑140 (MQ=255)                            tCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAG‑AAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACAt               >  1:1023860/1‑140 (MQ=255)                            tCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAG‑AAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACAt               >  1:899437/1‑140 (MQ=255)                            tCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAG‑AAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACAt               >  1:805335/1‑140 (MQ=255)                                aaTTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAG‑AAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACg           <  1:980993/140‑1 (MQ=255)                                aaTTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAG‑AAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACg           <  1:529071/140‑1 (MQ=255)                                         tATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAG‑AAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGc  >  1:1058686/1‑140 (MQ=9)                                                                                                                                     |                                                 ATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGCAGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGC  >  NC_000913/729450‑729630

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