Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_bme_pgi	1,529,218	C→T	pseudogene (1317/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,218	0	C	T	100.0%	28.0 / NA	14	pseudogene (1317/2037 nt)	rhsE	pseudogene, Rhs family
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (10/4);  total (10/4)

CCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTG  >  NC_000913_3_bme_pgi/1529162‑1529353                                                         |                                                                                                                                        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)                                         gacTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCgggg     <  2:75447/147‑1 (MQ=255)                                            tGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCCCCTGTACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTg  >  1:206590/1‑149 (MQ=255)                                            tGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGGCAGGGCAGCAGCATGATGAGGAGTCGGGGcgg  >  1:225599/1‑147 (MQ=255)                                                         |                                                                                                                                        CCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTG  >  NC_000913_3_bme_pgi/1529162‑1529353

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