Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	731,919	A→G	G779G (GGA→GGG)	rhsC →	Rhs protein with putative toxin domain, putative neighboring cell growth inhibitor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	731,919	0	A	G	100.0%	35.8 / NA	11	G779G (GGA→GGG)	rhsC	Rhs protein with putative toxin domain, putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (5/6);  total (5/6)

CTCTGGCAGCATGAGACCAGACATGCGTACAACGCACAGGGACTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGCTGGCTGTCAGGCATGAAACTCGGCGACACACCGCTGGTGGAGTACACCCGCGACCGCCTGCACCGGGA  >  NC_000913/731878‑732050                                          |                                                                                                                                    ctctGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGCTGGCTGTCAGGCATGAAACTCGGCGACACACCGc                                    <  2:1027747/139‑1 (MQ=255)                        tGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGCTGGCTGTCAGGCATGAAACTCGGCGACACACCGCTGGTGGAGTACACCCGCGACCGc             <  1:97680/139‑1 (MQ=255)                          cGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGCTGGCTGTCAGGCATGAAACTCgg                                               >  1:668515/1‑103 (MQ=255)                          cGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGCTGGCTGTCAGGCATGAAACTCgg                                               <  2:668515/103‑1 (MQ=255)                             aCAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGCTGGCTGTCAGGCATGAAACTCGGCGACACACCGCTGGTGGAGTACACCCGCGACCGCCTGCa        <  2:599675/139‑1 (MQ=19)                               aaCGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGCTGGCTGTCAGGCATGAAACTCGGCGACACACCGCTGGTGGAGTACACCCGCGACCGCCTGCAcc      >  2:401125/1‑139 (MQ=255)                                   cgcaGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGCTGGCTGTCAGGCATGAAACTCGGCGacacac                                       >  1:314648/3‑102 (MQ=255)                                   cgcaGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGCTGGCTGTCAGGCATGAAACTCGGCGacacac                                       >  1:865916/3‑102 (MQ=255)                                   cgcaGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGCTGGCTGTCAGGCATGAAACTCGGCGacacac                                       <  2:314648/100‑1 (MQ=255)                                   cgcaGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGCTGGCTGTCAGGCATGAAACTCGGCGacacac                                       <  2:865916/100‑1 (MQ=255)                                   cgcaGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGCTGGCTGTCAGGCATGAAACTCGGCGACACACCGCTGGTGGAGTACACCCGCGACCGCCTGCACCGGGa  >  1:637500/3‑139 (MQ=255)                                          |                                                                                                                                    CTCTGGCAGCATGAGACCAGACATGCGTACAACGCACAGGGACTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGCTGGCTGTCAGGCATGAAACTCGGCGACACACCGCTGGTGGAGTACACCCGCGACCGCCTGCACCGGGA  >  NC_000913/731878‑732050

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