Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP009273	1,191,739	A→G	K387K (AAA→AAG)	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,739	0	A	G	100.0%	28.8 / NA	19	K387K (AAA→AAG)	icd	e14 prophage; isocitrate dehydrogenase, specific for NADP+
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (11/8);  total (11/8)

ACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCA  >  CP009273/1191594‑1191796                                                                                                                                                  |                                                          acGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGAcc                                                        >  1:387406/1‑149 (MQ=21) acGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGAcc                                                        >  1:590768/1‑149 (MQ=21)           gcCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATg                                              >  1:441241/1‑149 (MQ=17)            ccGAAATACGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGa                                             <  2:560706/149‑1 (MQ=17)                   atGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGaa                                      >  1:293756/1‑148 (MQ=17)                   atGCCGGTCAGGACAAAGTAAATCCCGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGaa                                      >  1:346693/1‑148 (MQ=14)                    tGCCGGACAGGACAAAGTAAATCCTGGCTCTATTACTCTCTCCGCGGAGATGATGCTGCGCCATATGTGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGaac                                     <  2:477272/149‑3 (MQ=14)                     gCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGCTTGTTAAAGGTATGGAAGGAGCAATCAATGCCAAGACCGTAACCTATGACTTCGAACg                                    >  2:145131/1‑149 (MQ=17)                          tCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTga                               <  2:347270/149‑1 (MQ=14)                              gACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTCTGGAAGGCGCAATCAATGCCGAGACCGAAACTTATGACTTCGAACGTCTgatgga                           >  2:33271/1‑149 (MQ=14)                              gACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTgatgga                           <  1:405388/149‑1 (MQ=11)                               aCAAAGTAAACCCTGGCCCTATTATTCCCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTgatggaa                          <  1:565237/149‑2 (MQ=11)                                    gTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGCACGTCTGATGGAAGgcgc                     >  1:291335/1‑149 (MQ=1)                                         tCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGCGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTaagc                >  2:633527/1‑147 (MQ=11)                                            tGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGctgc             >  1:220339/1‑149 (MQ=11)                                                    ttattCTGTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGt     <  1:353962/149‑1 (MQ=17)                                                     tatcccctcCCCTGAGATGAGGCACCGCCATATCGGTTGCACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCCCAATCAAGGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGtt    <  2:308358/143‑1 (MQ=17)                                                       tttgctcCGCTGTGATGATGCTGGGCCATATGGGTTAGACTGAGGCAGCTGACGTGATTGTTAAAGGTATGGAAGGCGTAATCTATGTGAAGTCCGTAACCTATGACTTCGTACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCa  <  2:597269/145‑1 (MQ=11)                                                       ttCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGCTGGAAGGCGCTAAGCTGCTGAAATGTTCa  >  1:568616/1‑149 (MQ=14)                                                                                                                                                  |                                                          ACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCA  >  CP009273/1191594‑1191796

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