Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP006584	2,265,885	C→A	A155A (GCC→GCA)	LY180_10860 →	terminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006584	2,265,885	0	C	A	100.0%	16.7 / NA	12	A155A (GCC→GCA)	LY180_10860	terminase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (8/4);  total (8/4)

TTAATCCGAACGTCGCTAACCGCAACAAAGGCGGGCGTCGCAAACCGAAAAAGAATTTTTTCAGTGACGAGGCCATCGAAAAGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCCGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCGCCTTTGCCCGGCTGGTTGACGTTGACCTGACCGGTGACCCAATTGTCCTGGGCAATAACGGCGCAAAACTGATTTTTCTCGGCACCAACTCCAACACCG  >  CP006584/2265746‑2266147                                                                                                                                            |                                                                                                                                                                                                                                                                       ttAATCCGAACGTCGCTAACCGCAACAAAGGCGGGCGTCGCAAACCGAAAAAGAATTTTTTCAGTGACGAGGCCATCGAAAAGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGCATAcct                                                                                                          >  1:265652/1‑296 (MQ=14)                    ccGCAACAAAGGCGGGCGTCGCAAACCGAAAAAGAATTTTTTCAGTGACGAGGCCATCGAAAAGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCGCCTTTGCCCGGCTg                                                                                       >  1:211871/1‑298 (MQ=14)                           aaaGGCGGGCGTCGCAAACCGAAAAAGAATTTTTTCAGTGACGAGGCCATCGAAAAGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGAAGGCGTATGTGTTCCGCGAATACATCATCGCCTTTGCCCGGCTGgttgacg                                                                                >  1:242125/1‑298 (MQ=14)                                               gAAAAAGAATTTTTTCAGTGACGAGGCCATCGAAAAGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCGCCTTTGCCCGGCTGGTTGACGTTGACCTGACCGGTGACCCa                                                            >  1:294556/1‑298 (MQ=14)                                                              cAGTGACGAGGCCATCGCAAAGCTGGCGCAGATTTTCTTTGAGCAGACTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATCTCCTGAAATCCCGCAAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAACACAGTACATAACCAGATTTGTCTGTAGTCACTTACACCGCAGGCGTATGGGTCCACCGATCACACACTCGCCGTTGCCCGGCTGGTTTAAGTTGACCATCACGGGCACCCCATTGGTCTGGGCAAt                                             >  2:106912/1‑298 (MQ=2)                                                                           aTCGAAAAGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCGCCTTTGCCCGGCTGGTTGACGTTGACCTGACCGGTGACCCAATTGTCCTGGGCAATAACGGCGCAAAAc                                >  1:214117/1‑298 (MQ=17)                                                                               aaaaGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCGCCTTTGCCCGGCTGGTTGACGTTGACCTGACCGGTGACCCAATTGTCCTGGGCAATAACGGCGCAAAACTGAt                            >  1:173082/1‑298 (MQ=17)                                                                                       gAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCGCCTTTGCCCGGCTGGTTGACGTTGACCTGACCGGTGACCCAATTGTCCTGGGCAATAACGGCGCAAAACTGATTTTTCTCg                    <  1:331861/298‑1 (MQ=17)                                                                                        aGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCGCCTTTGCCCGGCTGGTTGACGTTGACCTGACCGGTGACCCAATTGTCCTGGGCAATAACGGCGCAAAACTGATTTTTCTCgg                   >  1:152094/1‑298 (MQ=17)                                                                                            gATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCAGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCGCCTTTGCCCGGCTGGTTGACGTTGACCTGACCGGTGACCCAATTGTCCTGGGCAATAACGGCGCAAAACTGATTTTTCTCGGCAcc               <  1:106912/298‑1 (MQ=17)                                                                                               tttCTGTAAGAATTCTTCCGAATATCAGGTCGACTGGGAGCGCGAAGGGCTTGAGCACCGCATCCGAGATAGCCTGAAATCCCGCCAGATGGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCGCCTTTGCCCGGCTGGTTGACGTTGACCTGACCGGTGACCCAATTGTCCTGGGCAATAACGGCGCAAAACTGATTTTTCTCGGCACCAAc            <  2:101556/298‑1 (MQ=11)                                                                                                         cAGTCTTTCGCATTCCCTTTTGACTGGATTCGCGCAGGGCTGGGGCACCGCATGCGCGTTCTCCTGAAACCCCGCCAGATTGGCTCCACGTTTTATTTTTCCCGCGAGGCGCTGCTGAGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGTCGTATGTGTTCCGCGAATACATCATCGCCTTTGCCCGGCTGGTTGACGTTGACCTGACCGGTGACCCAATTGTCCTGGGCAATAACGGCGCAAAACTGATTTTTCTCGGCACCAACTCCAACACcg  <  2:173082/298‑1 (MQ=11)                                                                                                                                            |                                                                                                                                                                                                                                                                       TTAATCCGAACGTCGCTAACCGCAACAAAGGCGGGCGTCGCAAACCGAAAAAGAATTTTTTCAGTGACGAGGCCATCGAAAAGCTGGAGCAGATTTTCTTTGAGCAGTCTTTCGAATATCAGTTGCACTGGTATCGCGCCGGGCTTGAGCACCGCATCCGCGATATCCTGAAATCCCGCCAGATTGGCGCGACGTTTTATTTTTCCCGCGAGGCGCTGCTGCGCGCCCTGAAAACCGGTCATAACCAGATTTTTCTGTCGGCCAGTAAAACGCAGGCGTATGTGTTCCGCGAATACATCATCGCCTTTGCCCGGCTGGTTGACGTTGACCTGACCGGTGACCCAATTGTCCTGGGCAATAACGGCGCAAAACTGATTTTTCTCGGCACCAACTCCAACACCG  >  CP006584/2265746‑2266147

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑