Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,247,134	A→C	Y467D (TAC→GAC)	lysP ←	lysine:H(+) symporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,247,134	0	A	C	100.0%	45.1 / NA	16	Y467D (TAC→GAC)	lysP	lysine:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (15/1);  total (15/1)

CGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAAATAATCAGGAACAGCGGGATACCAATATACGTTGCCGCTACGCCGCCCCAGTCAATAGTATCTTTCAGGAACGCTTCGTAGTTCTGGCCCAAAGTGATAATCAGACACAGAATGAATGCGAAGATCGGCCCCAGTGGGAAGAAACCTGAACGGTACGGCAGATCGTTAATGTCGTGTCCCTGCAATACGTAACCGCGACGGAAGCGATAGTGGCTAATGGCAATCCCCAGCCAG  >  NC_000913/2246939‑2247378                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     cGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTCGCCGAACCAAATAATCAGGAACAGCGGGATACAAATATACGTTGCCGCTAcgccg                                                                                                                                                                                               >  2:80548/1‑251 (MQ=255)                        aTCAATCCAAATGTTAAATTTTATATCCGCCTGATTTGTCAGGGAAATGGAATTGAAAAAGCACTCTCGGTTGAGAGGGCTTAGCACGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTACTCAGCTTGTCGCCGAACCAAATAATCAGGAACAGCGGGATACCAATATCCGTTGCCGCTCCGCCGCCCCAGTCAATa                                                                                                                                                                                   >  2:262029/1‑240 (MQ=255)                        aTCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTCGCCGAACCAAATAATCAGGAACAGCGGGATACCAATATACGTTGCCGCTACGCCGCCCCAGTCAATa                                                                                                                                                                                   <  1:262029/240‑1 (MQ=255)                            aTCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTCGCCGAACCAAATAATCAGGAACAGCGGGATACCAATATACGTTGCCGCTACGCCGCCCCAGTCAATAGTATCTTTCAGGAAc                                                                                                                                                                    >  1:234973/1‑251 (MQ=255)                                              atatCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTCGCCGAACCAAATAATCAGGAACAGCGGGATACCAATATACGTTGCCGCTACGCCGCCCCAGTCAATAGTATCTTTCAGGAACGCTTCGTAGTTCTGGccc                                                                                                                                                  >  1:290285/1‑251 (MQ=255)                                                           ttGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTCGCCGAACCAAATAATCAGGAACAGCGGGATACCAATATACGTTGCCGCTACGCCGCCCCAGTCAATAGTATCTTTCAGGAACGCTTCGTAGTTCTGGCCCAAAGTGATAATc                                                                                                                                      >  1:81806/1‑250 (MQ=255)                                                             gTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTCGCCGAACCAAATAATCAGGAACAGCGGGATACCAATATACGTTGCCGCTACGCCGCCCCAGTCAATAGTATCTTTCAGGAACGCTTCGTAGTTCTGGCCCAAAGTGATAATc                                                                                                                                      >  1:98255/1‑248 (MQ=255)                                                               cAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTCGCCGAACCAACTAATCAGGAACAGCTGGATAACAATATACTTTGCCGCTACGCCGCCCCAGACAATAGTACCTTTCAGGAACGCTTCGTAGTTCTGGCACTAAGTGATAATCAgtc                                                                                                                                  >  2:103417/1‑248 (MQ=255)                                                                                                gagaGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTCGCCGAACCAAATAATCAGGAACAGCGGGATACCAATATACGTTGCCGCt                                                                                                                                                                                                     >  2:107500/1‑150 (MQ=255)                                                                                                 agagGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCGTTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTCGCCGAACCAAATAATCAGGAACAGCGGGATACCAATATACGTTGCCGCTACGCCGCCCCAGTCAATAGTATCTTTCAGGAACGCTTCGTAGTTCTGGCCCAAAGTGATAATCAGACACAGAATGAATGCGAAGATCGGCCCCAGTGGga                                                                                                 >  1:214601/1‑249 (MQ=255)                                                                                                                     gaagaAACTTATATCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTAAGAAGAGAGTGCCTTTAATCAGCTTGTCGCCGAACCAAATAATCAGGAAAAGCGGGATACCAATAGACGTTGCCGCTACGCCGCAACAGTCAATAGTATCTTTCAGGAACGCTTCGTAGTTCTGGCCCAAAGAGATAATCAGACACAGAATGAATGCGAAGATCGGCCCAAGTGGGAAGAAACCTGAAAGGTACggacg                                                                           >  2:68972/4‑248 (MQ=255)                                                                                                                                  tCTTATCGTTCTGCGGGCACTTCATTTCGCTGTAGCATACGAAATGATTTCCTTTAATCAGCTTGTCGCCGAAACAAATAATCAGGAACAGCGGGGTACCAATATACGTTGACGCTACGACGCCCCAGTCAATAGTATCTTTCAGGAACGCTTAGTAGTTCTGGCCCAAAGTGATAATCAAACACAGAATGAATGCGAAGATCGGCCCAAGTGGGAAGAAACCTGAACGGTACGGCAGATAGTTAATGTCg                                                              >  2:290679/1‑251 (MQ=255)                                                                                                                                         gTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTCGCCGAACCAAATAATCAGGAACAGCGGGATACCAATATACGTTGCCGCTACGCCGCCCCAGTCAATAGTATCTTTCAGGAACGCTTCGTAGTTCTGGCCCAAAGTGATAATCAGACACAGAATGAATGCGAAGATCGGCCCCAGTGGGAAGAAACCTGAACGGTACGGCAGATCGTTAATGTCGTGTCCCt                                                       >  1:272273/1‑251 (MQ=255)                                                                                                                                                               cTGTAGCGTACGAAGAGAGTTCCTTTAATCAGATTGTCGCCGAACCAAGTAATCAGGAACAGCGGGATACCAATATACGTTGCCGCTACGCCGCCCCAGTCAATAGGATCTTTCAGGAACGATTAGTCGTTCTGGCCCACAGTGATAATAAGACACCGAATGAATGCGAAGATCGGCCCCAGTGGTAAGAAACCTGACAGGAACGACCGAACGTTAATGTCGTGTCACTGCaaactgtacccgcgacag                                   >  2:81620/1‑233 (MQ=255)                                                                                                                                                                             gTGAGATCCTTCAATCAGCATGTCGCCGCACCAAATAATCCGGAACAGCGGGATACCAATAAACGTTGCCGCTACGCCGCCACAGTCACTAGTATCTTTCAGGACCGCTTCGTAGTTCTGGCCCAAAGCGATAATCAGACACAGAATGACTGCGACGATCGGCCCCAGTGGGAAGACACCTGAACGGAACGGCAGATGGATAATGTCGTGTCCCTGCAATACGTAACCGCGACGGAACCGagaaaggcta                    >  2:288881/1‑241 (MQ=255)                                                                                                                                                                                              gCTTGTCGCCGAACCAAATAATCAGGAACAGCGTGATACCAATATACGTTGCCGCTACACCGCCCCAGTCAATAGTATCTTTCAGGAACGCTTCGTAGTTCTGGCCCAACGTGATAATCAGCCACAGAATGAATGCGAAGATCGGCCCCAGTGGGAAGAAACCTGAACGGTACGGCCGATAGTTAATATCGTGTCCATGCAATACGTAACCTAGACGGAAGCGATAGGGGccaaaggacatcccccagccg  >  2:28942/1‑231 (MQ=255)                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     CGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCGGTTGAGAGGGCTTAGCAAGGAAGGGAGGAAACTTATTTCTTATCGTTCTGCGGGAACTTCATTTCGCTGTAGCGTACGAAGTGAGTTCCTTTAATCAGCTTGTAGCCGAACCAAATAATCAGGAACAGCGGGATACCAATATACGTTGCCGCTACGCCGCCCCAGTCAATAGTATCTTTCAGGAACGCTTCGTAGTTCTGGCCCAAAGTGATAATCAGACACAGAATGAATGCGAAGATCGGCCCCAGTGGGAAGAAACCTGAACGGTACGGCAGATCGTTAATGTCGTGTCCCTGCAATACGTAACCGCGACGGAAGCGATAGTGGCTAATGGCAATCCCCAGCCAG  >  NC_000913/2246939‑2247378

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