Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,528,423	A→G	pseudogene (534/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,528,423	0	A	G	100.0%	27.4 / NA	14	pseudogene (534/2049 nt)	rhsE	RhsE protein in rhs element
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (7/7);  total (7/7)

CCTGAACAGCCTGGTGTATGACCGTGACTACGGGTGGAATGACAACGGCGACCTGGTGCGCATCAGCGGCCCGCGACAGACGCGGGAATACGGCTACAGCGCCACGGGCAGGCTGGAGAGTGTGCGCACCCTCGCACCAGACCTGGACATCCGCATCCCGTATGCCACGGACCCGGC  >  NC_000913/1528288‑1528464                                                                                                                                        |                                          ccTGAACAGCCTGGTGTATGACCGTGACTACGGGTGGAATGACAACGGCGACCTGGTGCGCATCAGCGGCCCGCGACAGACGCGGGAATACGGCTACAGCGCCACGGGCAGGCTGGAGAGTGTGCGCACCCTCGCGCCAg                                       <  2:1439836/140‑1 (MQ=21)             ggTGTATGACCGTGACTACGGGTGGAATGACAACGGCGACCTGGTGCGCATCAGCGGCCCGCGACAGACGCGGGAATACGGCTACAGCGCCACGGGCAGGCTGGAGAGTGTGCGCACCCTCGCGCCAGACCTGGACATcc                           >  1:1058772/1‑140 (MQ=37)                gtATGACCGTGACTACGGGTGGAATGACAACGGCGACCTGGTGCGCATCAGCGGCCCGCGACAGACGCGGGAATACGGCTACAGCGCCACGGGCAGGCTGGAGAGTGTGCGCACCCTCGCGCCAGACCTGGACATCCGCa                        <  1:463294/140‑1 (MQ=18)                 tATGACCGTGACTACGGGTGGAATGACAACGGCGACCTGGTGCGCATCAGCGGCCCGCGACAGACGCGGGAATACGGCTACAGCGCCACGGGCAGGCTGGAGAGTGTGCGCACCCTCGCGCCAGACCTGGACATCCGCAt                       <  1:535967/140‑1 (MQ=17)                     aCCGTGACTACGGGTGGAATGACAACGGCGACCTGGTGCGCATCAGCGGCCCGCGACAGACGCGGGAATACGGCTACAGCGCCACGGGCAGGCTGGAGAGTGTGCGCACCCTCGCGCCAGACCTGGACATCCGCATCCCg                   <  1:125439/140‑1 (MQ=17)                       cGTGACTACGGGTGGAATGACAACGGCGACCTGGTGCGCATCAGCGGCCCGCGACAGACGCGGGAATACGGCTACAGCGCCACGGGCAGGCTGGAGAGTGTGCGCACCCTCGCGCCAGACCTGGACATCCGCATCCCGTa                 <  1:135239/140‑1 (MQ=17)                         tGACTACGGGTGGAATGACAACGGCGACCTGGTGCGCATCAGCGGCCCGCGACAGACGCGGGAATACGGCTACAGCGCCACGGGCAGGCTGGAGAGTGTGCGCACCCTCGCGCCAGACCTGGACATCCGCATCCCGTATg               >  1:1509465/1‑140 (MQ=18)                               cGGGTGGAATGACAACGGCGACCTGGTGCGCATCAGCGGCCCGCGACAGACGCGGGAATACGGCTACAGCGCCACGGGCAGGCTGGAGAGTGTGCGCACCCTCGCGCCAGACCTGGACATCCGCATCCCGTATGCCACgg         >  2:511115/1‑140 (MQ=18)                                gggTGGAATGACAACGGCGACCTGGTGCGCATCAGCGGCCCGCGACAGACGCGGGAATACGGCTACAGCGCCACGGGCAGGCTGGAGAGTGTGCGCACCCTCGCGCCAGACCTGGACATCCGCATCCCGTATGCCACGGa        >  1:172380/1‑140 (MQ=18)                                     gAATGACAACGGCGACCTGGTGCGCATCAGCGGCCCGCGACAGACGCGGGAATACGGCTACAGCGCCACGGGCAGGCTGGAGAGTGTGCGCACCCTCGCGCCAGACCTGGACATCCGCATCCCGTATGCCACGGACCcgg   >  1:426113/1‑140 (MQ=17)                                     gAATGACAACGGCGACCTGGTGCGCATCAGCGGCCCGCGACAGACGCGGGAATACGGCTACAGCGCCACGGGCAGGCTGGAGAGTGTGCGCACCCTCGCGCCAGACCTGGACATCCGCATCCCGTATGCCACGGACCcgg   >  2:1002645/1‑140 (MQ=17)                                      aaTGACAACGGCGACCTGGTGCGCATCAGCGGCCCGCGACAGACGCGGGAATACGGCTACAGCGCCACGGGCAGGCTGGAGAGTGTGCGCCCCCTCGCGCCAGACCTGGACATCCGCATCCCGTATGCCACGGACCcggc  <  2:910574/140‑1 (MQ=14)                                      aaTGACAACGGCGACCTGGTGCGCATCAGCGGCCCGCGACAGACGCGGGAATACGGCTACAGCGCCACGGGCAGGCTGGAGAGTGTGCGCACCCTCGCGCCAGACCTGGACATCCGCATCCCGTATGCCACGGACCcggc  >  1:604524/1‑140 (MQ=14)                                      aaTGACAACGGCGACCTGGTGCGCATCAGCGGCCCGCGACAGACGCGGGAATACGGCTACAGCGCCACGGGCAGGCTGGAGAGTGTGCGCACCCTCGCGCCAGACCTGGACATCCGCATCCCGTATGCCACGGACCcggc  <  2:753814/140‑1 (MQ=14)                                                                                                                                        |                                          CCTGAACAGCCTGGTGTATGACCGTGACTACGGGTGGAATGACAACGGCGACCTGGTGCGCATCAGCGGCCCGCGACAGACGCGGGAATACGGCTACAGCGCCACGGGCAGGCTGGAGAGTGTGCGCACCCTCGCACCAGACCTGGACATCCGCATCCCGTATGCCACGGACCCGGC  >  NC_000913/1528288‑1528464

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: ‑