Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	731,472	T→G	G630G (GGT→GGG)	rhsC →	rhs element protein RhsC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	731,472	0	T	G	100.0%	50.6 / NA	18	G630G (GGT→GGG)	rhsC	rhs element protein RhsC
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (13/5);  total (13/5)

ACGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGTCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTG  >  NC_000913/731341‑731520                                                                                                                                    |                                                 acGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACAGGACACAGTACGATGCGTGGGGAAAAGCCATCCGTACCACGCAGGGCGGGCTGAcgcg                                          <  2:361939/140‑1 (MQ=255)   gCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCCGTACCACGCAGGGCGGGCTGAcgcg                                          <  1:247255/138‑1 (MQ=255)   gCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCCGTACCACGCAGGGCGGGCTGAcgcg                                          >  2:5982/1‑138 (MQ=255)   gCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCCGTACCACGCAGGGCGGGCTGAcgcg                                          <  1:5982/138‑1 (MQ=255)   gCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCCGTACCACGCAGGGCGGGCTGAcgcg                                          >  2:247255/1‑138 (MQ=255)        gCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCCGTACCACGCAGGGCGGGCTGACGCGCAGTATg                                   >  2:402388/1‑140 (MQ=255)                 cgcgGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCCGTACCACGCAGGGCGGGCTGAcgcg                                          >  2:772443/1‑124 (MQ=255)                 cgcgGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCCGTACCACGCAGGGCGGGCTGAcgcg                                          <  1:772443/124‑1 (MQ=255)                 cgcgGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGAt                          >  2:666633/1‑140 (MQ=255)                   cgGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGAtgc                        >  1:178937/1‑140 (MQ=255)                                aCGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTc           >  2:113862/1‑140 (MQ=255)                                aCGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTc           >  1:355730/1‑140 (MQ=255)                                    cgccgGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATcc       >  2:213203/1‑140 (MQ=18)                                    cgccgGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATcc       >  1:561794/1‑140 (MQ=18)                                    cgccgGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATcc       >  2:622732/1‑140 (MQ=17)                                         gTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTg  >  1:751322/1‑140 (MQ=18)                                                 accaccGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGcc    >  1:421935/1‑130 (MQ=18)                                                 accaccGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGcc    <  2:421935/130‑1 (MQ=17)                                                                                                                                    |                                                 ACGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGTCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTG  >  NC_000913/731341‑731520

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