Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	250,930	A→G	I425T (ATC→ACC)	pepD ←	cytosol nonspecific dipeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	250,930	0	A	G	100.0%	116.0 / NA	30	I425T (ATC→ACC)	pepD	cytosol nonspecific dipeptidase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (22/8);  total (22/8)

CGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGTTTACCCAGC  >  NZ_CP009273/250802‑251061                                                                                                                                 |                                                                                                                                    cGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGGTGTTCGGCg                                                                                                                            >  2:519919/1‑138 (MQ=255) cGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGGTGTTCGGCg                                                                                                                            <  1:519919/138‑1 (MQ=255) cGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGGTGTTCGGCGt                                                                                                                           >  2:178397/1‑139 (MQ=255)                               agTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGGTGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACg                                                                                             >  1:148808/1‑139 (MQ=255)                                 tgtgGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGGTGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCa                                                                                               >  1:263839/1‑135 (MQ=255)                                 tgtgGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGGTGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCa                                                                                               <  2:263839/135‑1 (MQ=255)                                 tgtgGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGGTGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTa                                                                                           >  1:347682/1‑139 (MQ=255)                                   tgGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGGTGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTAcc                                                                                         <  2:148808/139‑1 (MQ=255)                                   tgGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGGTGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTAc                                                                                          >  2:280268/1‑138 (MQ=255)                                   tgGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGGTGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTAc                                                                                          >  1:31872/1‑138 (MQ=255)                                   tgGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGGTGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTAc                                                                                          >  2:116307/1‑138 (MQ=255)                                           tGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGGTGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACg                                                                                             >  1:250584/1‑127 (MQ=255)                                           tGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGGTGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACg                                                                                             <  2:250584/127‑1 (MQ=255)                                           tGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGGTGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCAt                                                                                 >  1:193233/1‑139 (MQ=255)                                           tGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGGTGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCAt                                                                                 >  1:84515/1‑139 (MQ=255)                                           tGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGGTGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCAt                                                                                 >  1:392954/1‑139 (MQ=255)                                                    gCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGGTGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCg                                                                            >  1:71043/1‑135 (MQ=255)                                                    gCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGGTGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCg                                                                            <  2:71043/135‑1 (MQ=255)                                                                    gTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGGTGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTgc                                                        >  2:343654/1‑139 (MQ=255)                                                                    gTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGGTGTTCGGCGTCTCGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTgc                                                        >  2:313326/1‑139 (MQ=255)                                                                         tttCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGGTGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTgccagcc                                                   <  2:193233/139‑1 (MQ=255)                                                                         tttCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGGTGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTgccagc                                                    >  2:318175/1‑138 (MQ=255)                                                                                                  aCCACATTCCAGGCCCGCGTGGATAATCTGGGTGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGtt                          >  1:408365/1‑139 (MQ=255)                                                                                                           cAGGCCCGCGTGGATAATCTGGGTGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCgcc                  <  1:343654/138‑1 (MQ=255)                                                                                                            aGGCCCGCGTGGATAATCTGGGTGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCgccag                >  1:30564/1‑139 (MQ=255)                                                                                                            aGGCCCGCGTGGATAATCTGGGTGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCgccag                >  1:223477/1‑139 (MQ=255)                                                                                                                   cgTGGATAATCTGGGTGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGtt          >  1:470866/1‑138 (MQ=255)                                                                                                                   cgTGGATAATCTGGGTGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGtt          >  2:234659/1‑138 (MQ=255)                                                                                                                   cgTGGATAATCTGGGTGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGtt          >  2:353006/1‑138 (MQ=255)                                                                                                                           aTCTGGGTGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGTTTACCCAGc  <  1:280268/138‑1 (MQ=255)                                                                                                                                 |                                                                                                                                    CGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGTTTACCCAGC  >  NZ_CP009273/250802‑251061

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