Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_bme_pgi	1,529,203	C→A	pseudogene (1302/2037 nt)	rhsE →	pseudogene, Rhs family

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_bme_pgi	1,529,203	0	C	A	100.0%	17.5 / NA	11	pseudogene (1302/2037 nt)	rhsE	pseudogene, Rhs family
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (8/3);  total (8/3)

CCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGA  >  NC_000913_3_bme_pgi/1529162‑1529343                                          |                                                                                                                                             ccGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCCGGGGCTTATCCGCGGAAACGGCACTACGCCGCGGCGCGGGGAGGAGGATGAATGGGGCAACCAGCTTACTGCGGCGGAGCCGCATCACCCGCACCAGcc                                   >  2:257896/1‑149 (MQ=2)         aaaaGCTCATCTTTATCACGGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACAGCCATACGGCGTGGCGCGGGGAGGATGATGAATGGGGCAACCAGCTTAATGACGAGAACCCGCATCACCTGCACCAGCCGTACCGTc                           >  1:477922/1‑149 (MQ=11)          aaaGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGATGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTTATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCt                          <  2:195614/149‑1 (MQ=11)             gctCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCCTCCCGTCTGcc                       >  1:101223/3‑149 (MQ=11)                    tttATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGgcag                >  1:411860/1‑149 (MQ=255)                       atcACTCCGACCACCGGGGACCGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAACGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGgcagcag             <  1:31963/147‑1 (MQ=14)                        tcACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGGGTATGATGACTGGGGCAACCAGCTTAATGAGGAGAAGCCGCATCACCTGCAACGGCCGTCCCGTCTACCAGGgcggcagc            >  1:54801/2‑149 (MQ=11)                        tcACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCGGCACCAGCCGTACCGTCTGCCAGGgcagcagc            >  1:13047/2‑149 (MQ=255)                           cTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCatg         >  1:259649/1‑149 (MQ=255)                                 ccaccGGGTACTGCCGCTGGCGCTTAGCAGCGACGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCTTCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATgagg   <  1:145186/149‑1 (MQ=11)                                  caccGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCCGGGCACCCGCATGATgagga  >  2:503146/1‑149 (MQ=14)                                          |                                                                                                                                             CCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGA  >  NC_000913_3_bme_pgi/1529162‑1529343

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