Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,009,674	G→A	G304D (GGC→GAC)	ygeY →	putative peptidase YgeY

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,009,674	0	G	A	100.0%	73.8 / NA	23	G304D (GGC→GAC)	ygeY	putative peptidase YgeY
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (12/11);  total (12/11)

CAGTTTCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAAAGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCTGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGTTTGGCAAAGCGCCGGTTGTTGATAAGTGGAC  >  NC_000913/3009539‑3009800                                                                                                                                        |                                                                                                                               cAGTTTCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCCCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAAATCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGACCTg                                                                                                                            >  2:217991/1‑140 (MQ=255)    tttCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAAAGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGACCTGGtt                                                                                                                         >  1:3623/1‑140 (MQ=255)           cGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAAAGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGACCTGGTTTACCCaa                                                                                                                  <  2:3623/140‑1 (MQ=255)              ccGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAAAGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGACCTGGTTTACCCAACCg                                                                                                               >  2:83073/1‑140 (MQ=255)                             gggCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAAAGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGACCTGGTTTACCCAACCGAATGCTACTTc                                                                                                    <  2:116693/136‑1 (MQ=255)                             gggCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAAAGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGACCTGGTTTACCCAACCGAATGCTACTTc                                                                                                    >  1:116693/1‑136 (MQ=255)                              ggCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAAAGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGACCTGGTTTACCCAACCGAATGCTACTTCCCGAc                                                                                               >  1:72301/1‑140 (MQ=255)                                           aaGGAGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAAAGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGACCTGGTTTACCCAACCGAATGCTACTTCCCGACCTGGAAAGTGgaa                                                                                  <  1:145340/140‑1 (MQ=255)                                                gcgcTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAAAGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGACCTGGTTTACCCAACCGAATGCTACTTCCCGACCTGGAAAGTGGAagaaga                                                                             <  1:70112/140‑1 (MQ=255)                                                              ccGCGCCCTGCCTGCAGTACAGAAAGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGACCTGGTTTACCCAACCGAATGCTACTTCCCGACCTGGAAAGTGGAAGAAGATCACTTCACCGTTa                                                               <  1:176371/140‑1 (MQ=255)                                                                 cgcCCTGCCTGCAGTACAGAAAGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGACCTGGTTTACCCAACCGAATGCTACTTCCCGACCTGGAAAGTGGAAGAAGATCACTTCACCGTTAAAg                                                            <  2:164890/140‑1 (MQ=255)                                                                         ctgcAGTACAGAAAGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGACCTGGTTTACCCAACCGAATGCTACTTCCCGACCTGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTg                                                    <  1:83073/140‑1 (MQ=255)                                                                                 cAGAAAGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGACCTGGTTTACCCAACCGAATGCTACTTCCCGACCTGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTa                                            >  2:290382/1‑140 (MQ=255)                                                                                         tAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGACCTGGTTTACCCAACCGAATGCTACTTCCCGACCTGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGc                                    >  1:230626/1‑140 (MQ=255)                                                                                         tAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGACCTGGTTTACCCAACCGAATGCTACTTCCCGACCTGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGc                                    >  2:127221/1‑140 (MQ=255)                                                                                         tAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGACCTGGTTTACCCAACCGAATGCTACTTCCCGACCTGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGc                                    >  1:134052/1‑140 (MQ=255)                                                                                                        tATGTACAACTACGACCGTCCGTCCTGGACTGACCTGGTTTACCCAACCGAATGCTACTTCCCGACCTGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGTTTGGCAAAgcgc                     <  1:340079/140‑1 (MQ=255)                                                                                                                    cGACCGTCCGTCCTGGACTGACCTGGTTTACCCAACCGAATGCTACTTCCCGACCTGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGTTTGGCAAAGCGCCGGTTGTTGATa         >  1:169412/1‑140 (MQ=255)                                                                                                                      aCCGTCCGTCCTGGACTGACCTGGTTTACCCAACCGAATGCTACTTCCCGACCTGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGTTTGGCAAAGCGCCGGTTGTTGATAAg       <  2:72301/140‑1 (MQ=255)                                                                                                                        cgtccgtcCTGGACTGACCTGGTTTACCCAACCGAATGCTACTTCCCGACCTGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGTTTGGCAAAGCGCCGGTTGTTGATAAGTg     >  2:269362/1‑140 (MQ=255)                                                                                                                         gtccgtccTGGACTGACCTGGTTTACCCAACCGAATGCTACTTCCCGACCTGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGTTTGGCAAAGCGCCGGTTGTTGATAAGTgg    <  1:127221/140‑1 (MQ=255)                                                                                                                          tccgtccTGGACTGACCTGGTTTACCCAACCGAATGCTACTTCCCGACCTGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGTTTGGCAAAGCGCCGGTTGTTGATAAGTGGa   <  1:318456/140‑1 (MQ=255)                                                                                                                           ccgtccTGGACTGACCTGGTTTACCCAACCGAATGCTACTTCCCGACCTGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGTTTGGCAAAGCGCCGGTTGTTGATAAGTGGAc  >  2:194887/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                               CAGTTTCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAAAGCTAACGCGGTTGTTTCTATGTACAACTACGACCGTCCGTCCTGGACTGGCCTGGTTTACCCAACCGAATGCTACTTCCCGACCTGGAAAGTGGAAGAAGATCACTTCACCGTTAAAGCACTGGTGAATGCCTACGAAGGGCTGTTTGGCAAAGCGCCGGTTGTTGATAAGTGGAC  >  NC_000913/3009539‑3009800

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