Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP009273	1,191,736	G→C	A386A (GCG→GCC)	icd →	e14 prophage; isocitrate dehydrogenase, specific for NADP+

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	1,191,736	0	G	C	100.0%	39.1 / NA	22	A386A (GCG→GCC)	icd	e14 prophage; isocitrate dehydrogenase, specific for NADP+
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (13/9);  total (13/9)

ACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGG‑CTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGA  >  CP009273/1191590‑1191798                                                                                                                                                    |                                                               acccacGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGG‑CTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCaa                                                              >  2:543614/1‑149 (MQ=21)        gTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGG‑CTGACCTGATTGGTAAAGGTATGGGAGGCGCAATCAATGCCAAGACCGTa                                                       >  1:273059/1‑149 (MQ=18)             gcgcCGAAAGAGGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGG‑CTGACCGGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTa                                                  <  2:46792/149‑1 (MQ=17)              cgcCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGG‑CTGACCTGATTGTTAAAGGTATGGGAGGCGCAATCAGTGCCAAGAACGTAActttt                                                 >  1:171861/1‑145 (MQ=18)              cgcCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACCGACGAGGAATGGC‑TGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTAt                                                 <  1:506064/149‑1 (MQ=18)                ccGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCCCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGG‑CTTACTTGAATGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGa                                               <  1:31977/149‑1 (MQ=18)                 cGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGG‑CTGACCTGATTGGTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGAc                                              >  1:729438/1‑149 (MQ=17)                   aaaTATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGG‑CTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACtt                                            >  1:14217/1‑149 (MQ=17)                   aaaTAGGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGG‑CTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACtt                                            <  2:41421/149‑1 (MQ=17)                         gCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGG‑CTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACg                                      >  1:813580/1‑149 (MQ=14)                           cGGTCAGGACAAAGTAAATCCTGGCTCTTTTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGG‑CTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTc                                    <  2:16573/149‑1 (MQ=14)                           cGGGCAGGACATAGTGAATCCTGGCTCTTTTATGCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGG‑CTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTc                                    <  2:468857/149‑1 (MQ=11)                             gTCAGGACAAAGTAAATCCTGTCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGG‑CTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTg                                  >  2:801840/1‑149 (MQ=14)                             gTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGGGATGATGCTGCGCCATATGGGTTGGACTGAATCGG‑CTGGCCTGGCTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTg                                  <  2:751202/149‑1 (MQ=14)                                  gACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGG‑CTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTgatgga                             >  2:16514/1‑149 (MQ=14)                                   aCAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGG‑CTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTgatggga                            >  1:400441/1‑147 (MQ=11)                                                  gCTCTATTATTCTCGCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGG‑CTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGctgctg             <  2:544749/149‑1 (MQ=14)                                                      tattatTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGG‑CTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCATGACCGTAACTTATGACTTTGAACGCCTGATGGAAGGCGCTAAGCTGCTGAAAt         <  1:525884/149‑1 (MQ=14)                                                      tattatTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGG‑CTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAAt         >  1:927107/1‑149 (MQ=14)                                                      tattatTCTCTCAGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGG‑CTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAAt         >  2:800875/1‑149 (MQ=14)                                                       attattCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGG‑CTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACAGTAACTTATGAATTCGAAGGTCTGATGGAAGGCCCTAAGCTGCTGAAATg        >  2:849119/1‑149 (MQ=14)                                                             ctctcCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGG‑CTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGACGCTAAGCTGCTGAACTGTTCaga  >  2:120041/1‑149 (MQ=14)                                                                                                                                                    |                                                               ACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGG‑CTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGA  >  CP009273/1191590‑1191798

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