Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,529,865	T→C	pseudogene (1976/2049 nt)	rhsE →	RhsE protein in rhs element

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,529,865	0	T	C	100.0%	23.6 / NA	8	pseudogene (1976/2049 nt)	rhsE	RhsE protein in rhs element
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (2/6);  total (2/6)

TTTTCTGCCAAATGTAGCGATGCATGTATTGAAAAATATAACACTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATT  >  NC_000913/1529810‑1529933                                                        |                                                                     ttttCTGCCAAATGTAGCGATGCATGTATTGAAAAATATAACACTGGGCAACTCTCCGGGTTAGATGGTATAAAAGCAGa                                              >  2:285157/1‑80 (MQ=255)  tttCTGCCAAATGTAGCGATGCATGTATTGAAAAATATAACACGGGGCAACTCGCCGGGTTAGATGGTATAAAAGCAGAt                                             <  1:285145/80‑1 (MQ=255)   ttCTGCCAAATGTAGCGATGCATGTATTGAAAAATATAACACTGGGCAACTCTCCGGGTTAGATGGTATAAAAGCAGATa                                            <  1:285146/80‑1 (MQ=255)         cAAATGTAGCGATGCATGTATTGAAAAATATAACACTGGGCAACTCTCCGGGTTAGATGGTATAAAAGCAGATAATCCAAt                                     <  1:285147/81‑1 (MQ=255)               tAGCGATGCATGTATTGAAAAATATAACACTGGGCAACTCTCCGGGTTAGATGGTATAAAAGCAGATAATCCAATAaaagc                               <  1:285148/81‑1 (MQ=255)                         tGTATTGAAAAATATAACACTGGGCAACTCTCCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAggg                     <  2:285149/81‑1 (MQ=255)                                           acTGGGCAACTCTCCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAat   >  1:285163/1‑81 (MQ=255)                                            cTGGGCAACTCTCCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAatt  <  2:285150/81‑1 (MQ=255)                                                        |                                                                     TTTTCTGCCAAATGTAGCGATGCATGTATTGAAAAATATAACACTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATT  >  NC_000913/1529810‑1529933

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