Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP009273	1,191,676	C→T	H366H (CAC→CAT)	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,676	0	C	T	100.0%	79.7 / NA	36	H366H (CAC→CAT)	icd	e14 prophage; isocitrate dehydrogenase, specific for NADP+
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (22/14);  total (22/14)

GGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCA  >  CP009273/1191533‑1191796                                                                                                                                                |                                                                                                                         ggtatcggtatcGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATggg                                                                                                                     >  1:369571/1‑149 (MQ=255)      cggtatcgCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGa                                                                                                                <  1:530253/149‑1 (MQ=255)       ggtatcgCCCCTGGTGCAAACATCGGTGACGAATGCGCCATGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGAc                                                                                                               >  2:116125/1‑149 (MQ=255)        gtatcgCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGAct                                                                                                              >  1:460181/1‑148 (MQ=255)               cccTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCgg                                                                                                       <  2:638640/149‑1 (MQ=255)                ccTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCCTATGGGTTGTACTGAAGCGGc                                                                                                      >  2:280851/1‑149 (MQ=255)                       cAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGAcctg                                                                                               >  1:227453/1‑146 (MQ=255)                       cAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGGTGGACTGAAGCGGCTGAcctg                                                                                               >  1:673234/1‑146 (MQ=255)                        aaaCATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGCGATGATGCTGCGCCATATGGGTTGGACGGAAGCGGCTGAcctga                                                                                              >  1:190277/1‑145 (MQ=255)                            aTCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGt                                                                                          >  1:588813/1‑149 (MQ=255)                            aTCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGt                                                                                          <  2:35001/149‑1 (MQ=255)                            aTCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAACCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGt                                                                                          <  1:268610/149‑1 (MQ=38)                            aTCGGGGACGAATGCGCCCCGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGt                                                                                          <  1:1349/149‑1 (MQ=38)                                             ccTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATAATGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGGGGGATTGAGGGGGCTGGGCTGAGTGTTAAAGGTATGGacgacg                                                                         >  1:232504/1‑144 (MQ=255)                                                            ccacGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCaatgccaaga                                                          <  1:593119/149‑9 (MQ=21)                                                              acGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGCAGGGCCAATccatggcaggccc                                                        >  2:659651/1‑137 (MQ=25)                                                              acGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGAcc                                                        >  1:590768/1‑149 (MQ=21)                                                              acGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGAcc                                                        >  1:387406/1‑149 (MQ=21)                                                              acGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGGTAAAGGTATGGGAGGCGCAATCattgccacgagc                                                        >  2:659655/1‑138 (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)                                                                                                                                                |                                                                                                                         GGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCA  >  CP009273/1191533‑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: ‑