Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	506,730	(G)5→4	coding (353/480 nt)	ybaK ←	Cys‑tRNA(Pro) and Cys‑tRNA(Cys) deacylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	506,726	0	G	.	100.0%	127.0 / NA	31	coding (357/480 nt)	ybaK	Cys‑tRNA(Pro) and Cys‑tRNA(Cys) deacylase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (20/11);  total (20/11)

CGGCACAATGCTTAATCGCGGCGGGCGATATCAGCAAATTTGGCATCGAGGATCTTTGCCAGATCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCGCCGGAAACATAAATAGTGGCAAATTCTTGTGCGGGGGCGTCGATAATCGTTGGCAGACGTTTTTTCTGCCCCAGTGGGCTAATCCCCCCAACCAGGTATCCCGTCGAACGCTGCGCGACCATCGGATCGGCCATCTCAACTTTCTTGGCACCCAGCGCT  >  NC_000913/506592‑506852                                                                                                                                       |                                                                                                                               cGGCACAATGCTTAATCGCGGCGGGCGATATCAGCAAATTTGGCATCGAGGATCTTTGCCAGATCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCGCCGGAAACATAAATAGTGGCAAATTCTTGTGCggggcg                                                                                                                           <  2:1257189/140‑3 (MQ=255)  ggCACAATGCTTAATCGCGGCGGGCGATATCAGCAAATTTGGCATCGAGGATCTTTGCCAGATCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCGCCGGAAACATAAATAGTGGCAAATTCTTGTGCggggcgg                                                                                                                          >  2:224688/1‑137 (MQ=255)         tGCTTAATCGCGGCGGGCGATATCAGCAAATTTGGCATCGAGGATCTTTGCCAGATCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCGCCGGAAACATAAATAGTGGCAAATTCTTGTGCgggg                                                                                                                             <  2:995493/130‑1 (MQ=255)         tGCTTAATCGCGGCGGGCGATATCAGCAAATTTGGCATCGAGGATCTTTGCCAGATCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCGCCGGAAACATAAATAGTGGCAAATTCTTGTGCgggg                                                                                                                             >  1:995493/1‑130 (MQ=255)                   cggcggGCGATATCAGCAAATTTGGCATCGAGGATCTTTGCCAGATCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCGCCGGAAACATAAATAGTGTCAAATTCATGTGC‑GAGGCGTCGATAATCGTTGGCAGa                                                                                                        >  1:447412/1‑140 (MQ=255)                   cggcggGCGATATCAGCAAATTTGGCATCGAGGATCTTTGCCAGATCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCGCCGGAAACATAAATAGTGGCAAATTCTTGTGC‑GGGGCGTCGATAATCGTTGGCAGa                                                                                                        >  1:927848/1‑140 (MQ=255)                              atCAGCAAATTTGGCATCGAGGATCTTTGCCAGATCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCGCCGGAAACATAAATAGTGGCAAATTCTTGTGC‑GGGGCGTCGATAATCGTTGGCAGACGTTTTTTCTg                                                                                             <  2:927848/140‑1 (MQ=255)                                        ttGGCATCGAGGATCTTTGCCAGATCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCGCCGGAAACATAAATAGTGGCAAATTCTTGTGC‑GGGGCGTCGATAATCGTTGACAGACGTTTTTTCTGCCCCAGTggg                                                                                   <  1:896968/140‑1 (MQ=255)                                          ggCATCGAGGATCTTTGCCAGATCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCGCCGGAAACATAAATAGTGGCAAATTCTTGTGC‑GGGGCGTCGATAATCGTTGGCAGACGTTTTTTCTGCCCCAGTGGGCt                                                                                 >  1:192636/1‑140 (MQ=255)                                                  ggATCTTTGCCAGATCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCGCCGGAAACATAAATAGTGGCAAATTCTTGTGC‑GGGGCGTCGATAATCGTTGGCAGACGTTTTTTCTGCCCCAGTGGGCTAATccccc                                                                         <  2:192636/140‑1 (MQ=255)                                                         tGCCAGATCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCGCCGGAAACATAAATAGTGGCAAATTCTTGTGC‑GGGGCGTCGATAATCGTTGGCAGACGTTTTTTATGCCCCAGTGGGCTAATCCCCCCAACCAg                                                                  >  1:537559/1‑140 (MQ=255)                                                           ccAGATCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCGCCGGAAACATAAATAGTGGCAAATTCTTGTGC‑GGGGCGTCGATAATCGTTGGCAGACGTTTTTTCTGCCCCAGTGGGCTAATCCCCCCAACCAGGt                                                                >  1:293214/1‑140 (MQ=255)                                                              gATCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCGCCGGAAACATAAATAGTGGCAAATTCTTGTGC‑GGGGCGTCGATAATCGTTGGCAGACGTTTTTTCTGCCCCAGTGGGCTAATCCCCCCAACCAGGTATc                                                             <  2:379766/140‑1 (MQ=255)                                                               aTCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCGCCGGAAACATAAATAGTGGCAAATTCTTGTGC‑GGGGCGTCGATAATCGTTGGCAGACGTTTTTTCTGCCCCAGTGGGCTAATCCCCCCa                                                                       >  1:551621/1‑129 (MQ=255)                                                               aTCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCGCCGGAAACATAAATAGTGGCAAATTCTTGTGC‑GGGGCGTCGATAATCGTTGGCAGACGTTTTTTCTGCCCCAGTGGGCTAATCCCCCCAACCAGGTATcc                                                            >  1:1081018/1‑140 (MQ=255)                                                               aTCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCGCCGGAAACATAAATAGTGGCAAATTCTTGTGC‑GGGGCGTCGATAATCGTTGGCAGACGTTTTTTCTGCCCCAGTGGGCTAATCCCCCCAACCAGGTATcc                                                            >  1:134524/1‑140 (MQ=255)                                                               aTCGCCTGCCGCCAGTTAGATATCCAGTCCGCGCTTGCCGCCGGAAACATAAATAGTGGCAAATTCTTGTGC‑GGGGCGTCGATAATCGTTGGCAGACGTTTTTTCTGCCCCAGTGGGCTAATCCCCCCAAc                                                                     >  2:1065715/1‑131 (MQ=255)                                                               aTCGCCTGCCGCCAGTTAGATATCCAGTCCGCGCTTGCCGCCGGAAACATAAATAGTGGCAAATTCTTGTGC‑GGGGCGTCGATAATCGTTGGCAGACGTTTTTTCTGCCCCAGTGGGCTAATCCCCCCAAc                                                                     <  1:1065715/131‑1 (MQ=255)                                                               aTCGCCTGCCGACAGTTCGATATCCAGTCCGCGCTTGCCGCCGGAAACATAAATAGTGGCAAATTCTTGTGC‑GGGGCGTCGATAATCGTTGGCAGACGTTTTTTCTGCCCCAGTGGGCTAATCCCCCCa                                                                       <  2:551621/129‑1 (MQ=255)                                                                tCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCGCCGGAAACATAAATAGTGGCAAATTCTTGTGC‑GGGGCGTCGATAATCGTTGGCAGACGTTTTTTCTGCCCCAGTGGGCTAATCCCCCCAACCAGGTATccc                                                           >  2:161044/1‑140 (MQ=255)                                                                        cgccAGTTCGATATCCAGTCCGCGCTTGCCGCCGGAAACATAAATAGTGGCAAATTCTTGTGC‑GGGGCGTCGATAATCGTTGGCAGACGTTTTTTCTGCCCCAGTGGGCTAATCCCCCCAACCAGGTATCCCGTCGAACg                                                   <  1:1506033/140‑1 (MQ=255)                                                                           cAGTTCGATATCCAGTCTGCGCTTGCCGCCGGAAACATAAATAGTGGCTAATTCTTGTGC‑GGGGCGTCGATAATCGTTGGCAGACGTTTTTTCTGCCCCAGTGGGCTAATCCCCCCAACCAGGTATCCCGTCGAACGCTg                                                <  2:447412/140‑1 (MQ=255)                                                                              ttCGATATCCAGTCCGCGCTTGCCGCCGGAAACATAAATAGTGGCAAATTCTTGTGC‑GGGGCGTCGATAATCGTTGGCAGACGTTTTTTCTGCCCCAGTGGGCTAATCCCCCCAACCAGGTATCCCGTCGAACGCTgcgc                                             <  2:1081018/140‑1 (MQ=255)                                                                                         gTCCGCGCTTGCCGCCGGAAACATAAATAGTGGCAAATTCTTGTGC‑GGGGCGTCGATAATCGTTGGCAGACGTTTTTTCTGCCCCAGTGGGCTAATCCCCCCAACCAGGTATCCCGTCGAACGCTGCGCGACCatcggat                                  >  2:51581/1‑140 (MQ=255)                                                                                          tCCGCGCTTGCCGCCGGAAACATAAATAGTGGCAAATTCTTGTGC‑GGGGCGTCGATAATCGTTGGCAGACGTTTTTTCTGCCCCAGTGGGCTAATCCCCCCACCCAGGTACCCCATCAAACGCTGCGCGACCatcggatc                                 >  1:589040/1‑140 (MQ=255)                                                                                          tCCGCGCTTGCCGCCGGAAACATAAATAGTGGCAAATTCTTGTGC‑GGGGCGTCGATAATCGTTGGCAGACGTTTTTTCTGCCCCAGTGGGCTAATCCCCCCAACCAGGTAGCCCATCCATCGCTGCGCGACCatcggatc                                 >  1:560190/1‑140 (MQ=255)                                                                                           ccGCGCTTGCCGCCGGAAACATAAATAGTGGCAAATTCTTGTGC‑GGGGCGTCGATAATCGTTGGCAGACGTTTTTTCTGCCCCAGTGGGCTAATCCCCCCAACCAGGTATCCCGTCGAACGCTGCGCGACCatcggatcg                                >  2:338983/1‑140 (MQ=255)                                                                                            cgcgCTTGCCGCCGGAAACATAAATAGTGGCAAATTCTTGTGC‑GGGGCGTCGATAATCGTTGGCAGACGTTTTTTCTGCCCCAGTGGGCTAATCCCCCCAACCAGGTATCCCGTCGAACGCTGCGCGACCatcggatcgg                               >  2:256371/1‑140 (MQ=255)                                                                                              cgcTTGCCGCCGGAAACATAAATAGTGGCAAATTCTTGTGC‑GGGGCGTCGATAATCGTTGGCAGACGTTTTTTCTGCCCCAGTGGGCTAATCCCCCCAACCAGGTATCCCGTCGAACGCTGCGCGACCATCGGATCGGcc                             >  2:313743/1‑140 (MQ=255)                                                                                               gcTTGCCGCCGGAAACATAAATAGTGGCAAATTCTTGTGC‑GGGGCGTCGATAATCGTTGGCAGACGTTTTTTCTGCCCCAGTGGGCTAATCCCCCCAACCAGGTATCCCGTCGAACGCTGCGCGACCATCGGATCGGCCa                            >  1:793558/1‑140 (MQ=255)                                                                                               gcTTGCCGCCGGAAACATAAATAGTGGCAAATTCTTGTGC‑GGGGCGTCGATAATCGTTGGCAGACGTTTTTTCTGCCCCAGTGGGCTAATCCCCCCAACCAGGTATCCCGTCGAACGCTGCGCGACCATCGGATCGGCCa                            >  2:1398186/1‑140 (MQ=255)                                                                                                cTTGCCGCCGGAAACATAAATAGTGGCAAATTCTTGTGC‑GGGGCGTCGATAATCGTTGGCAGACGTTTTTTCTGCCCCAGTGGGCTAATCCCCCCAACCAGGTATCCCGTCGAACGCTGCGCGACCATCGGATCGGCCAt                           <  1:161044/140‑1 (MQ=255)                                                                                                      gccgGAAACATAAATAGTGGCAAATTCTTGTGC‑GGGGCGTCGATAATCGTTGGCAGACGTTTTTTCTGCCCCAGTGGGCTAATCCCCCCAACCAGGTATCCCGTCGAACGCTGCGCGACCATCGGATCGGCCATCTCAAc                     <  2:793558/140‑1 (MQ=255)                                                                                                                         gCAAATTCTTGTGC‑GGGGCGTCGATAATCGTTGGCAGACGTTTTTTCTGCCCCAGTGGGCTAATCCCCCCAACCAGGTATCCCGTCGAACGTTGCGCGACCATCAGATGGGCCATCTCAACTTTCTTGGCACCCAGCGCt  >  2:519104/1‑140 (MQ=255)                                                                                                                         gCAAATTCTTGTGC‑GGGGCGTCGATAATCGTTGGCAGACGTTTTTTCTGCCCCAGTGGGCTAATCCCCCCAACCAGGTATCCCGTCGAACGCTGCGCGACGATCGGGTGAGCCATCTAAACTGACTTGGCAACCAGCGCt  >  2:718770/1‑140 (MQ=255)                                                                                                                                       |                                                                                                                               CGGCACAATGCTTAATCGCGGCGGGCGATATCAGCAAATTTGGCATCGAGGATCTTTGCCAGATCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCGCCGGAAACATAAATAGTGGCAAATTCTTGTGCGGGGGCGTCGATAATCGTTGGCAGACGTTTTTTCTGCCCCAGTGGGCTAATCCCCCCAACCAGGTATCCCGTCGAACGCTGCGCGACCATCGGATCGGCCATCTCAACTTTCTTGGCACCCAGCGCT  >  NC_000913/506592‑506852

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