Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	2,135,155	T→C	N181S (AAC→AGC)	dcd ←	dCTP deaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	2,135,155	0	T	C	100.0%	157.7 / NA	42	N181S (AAC→AGC)	dcd	dCTP deaminase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (28/14);  total (28/14)

CATCAGCGTCGTCAGAAATCGTCTCATGGTATCCTCAATGGGCCTGAATTAGTCTTTATCGATTCGGCTGGCTACCGCGCCCTGCTGGTTGCGATATTTCGCATCTTCACGGCGGTTGTAAGGTCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCA  >  NZ_CP009273/2135067‑2135229                                                                                         |                                                                           catcaGCGTCGTCAGAAATCGTCTCATGGTATCCTCAATGGGCCTGAATTAGTCTTTATCGATTCGGCTGGCTACCGCGCCCTgctggctg                                                                          >  2:446073/1‑91 (MQ=255)  atcaGCGTCGTCAGAAATCGTCTCATGGTATCCTCAATGGGCCTGAATTAGTCTTTATCGATTCGGCTGGCTACCGCGCCCTGCTGGCTGc                                                                         >  1:436897/1‑91 (MQ=255)    caGCGTCGTCAGAAATCGTCTCATGGTATCCTCAATGGGCCTGAATTAGTCTTTATCGATTCGGCTGGCTACCGCGCCCTGCTGGCTGCGa                                                                       >  1:473292/1‑91 (MQ=255)        gtcgtcAGAAATCGTCTCATGGTATCCTCAATGGGCCTGAATTAGTCTTTATCGATTCGGCTGGCTACCGCGCCCTGCTGGCTGCGATAtt                                                                   >  1:174908/1‑91 (MQ=255)         tcgtcAGAAATCGTCTCATGGTATCCTCAATGGGCCTGAATTAGTCTTTATCGATTCGGCTGGCTACCGCGCCCTGCTGGCTGCGATAttt                                                                  >  1:528375/1‑91 (MQ=255)          cgtcAGAAATCGTCTCATGGTATCCTCAATGGGCCTGAATTAGTCTTTATCGATTCGGCTGGCTACCGCGCCCTGCTGGCTGCGATATTTc                                                                 <  2:184767/91‑1 (MQ=255)          cgtcAGAAATCGTCTCATGGTATCCTCAATGGGCCTGAATTAGTCTTTATCGATTCGGCTGGCTACCGCGCCCTGCTGGCTGCGATATTTc                                                                 >  1:301787/1‑91 (MQ=255)               gAAATCGTCTCATGGTATCCTCAATGGGCCTGAATTAGTCTTTATCGATTCGGCTGGCTACCGCGCCCTGCTGGCTGCGATATTTCGCATc                                                            <  1:320176/91‑1 (MQ=255)                     gTCTCATGGTATCCTCAATGGGCCTGAATTAGTCTTTATCGATTCGGCTGGCTACCGCGCCCTGCTGGCTGCGATATTTCGCATCTTCAcg                                                      >  2:47908/1‑91 (MQ=255)                        tcATGGTATCCTCAATGGGCCTGAATTAGTCTTTATCGATTCGGCTGGCTACCGCGCCCTGCTGGCTGCGATATTTCGCATCTTCAcggcg                                                   <  2:168961/91‑1 (MQ=255)                             gTATCCTCAATGGGCCTGAATTAGTCTTTATCGATTCGGCTGGCTACCGCGCCCTGCTGGCTGCGATATTTCGCATCTTCACGGCGGTTGt                                              <  2:473292/91‑1 (MQ=255)                               aTCCTCAATGGGCCTGAATTAGTCTTTATCGATTCGGCTGGCTACCGCGCCCTGCTGGCTGCGATATTTCGCATCTTCACGGCGGTTGTaa                                            >  2:139035/1‑91 (MQ=255)                                tCCTCAATGGGCCTGAATTAGTCTTTATCGATTCGGCTGGCTACCGCGCCCTGCTGGCTGCGATATTTCGCATCTTCACGGCGGTTGTAAg                                           >  1:282827/1‑91 (MQ=255)                                   tCAATGGGCCTGAATTAGTCTTTATCGATTCGGCTGGCTACCGCGCCCTGCTGGCTGCGATATTTCGCATCTTCACGGCGGTTGTAAGGTc                                        >  1:697545/1‑91 (MQ=255)                                   tCAATGGGCCTGAATTAGTCTTTATCGATTCGGCTGGCTACCGCGCCCTGCTGGCTGCGATATTTCGCATCTTCACGGCGGTTGTAAGGTc                                        >  1:110426/1‑91 (MQ=255)                                    cAATGGGCCTGAATTAGTCTTTATCGATTCGGCTGGCTACCGCGCCCTGCTGGCTGCGATATTTCGCATCTTCACGGCGGTTGTAAGGTCg                                       >  1:314983/1‑91 (MQ=255)                                    cAATGGGCCTGAATTAGTCTTTATCGATTCGGCTGGCTACCGCGCCCTGCTGGCTGCGATATTTCGCATCTTCACGGCGGTTGTAAGGTCg                                       >  2:541405/1‑91 (MQ=255)                                      aTGGGCCTGAATTAGTCTTTATCGATTCGGCTGGCTACCGCGCCCTGCTGGCTGCGATATTTCGCATCTTCACGGCGGTTGTAAGGTCGCa                                     >  1:11691/1‑91 (MQ=255)                                      aTGGGCCTGAATTAGTCTTTATCGATTCGGCTGGCTACCGCGCCCTGCTGGCTGCGATATTTCGCATCTTCACGGCGGTTGTAAGGTCGCa                                     >  2:855356/1‑91 (MQ=255)                                      aTGGGCCTGAATTAGTCTTTATCGATTCGGCTGGCTACCGCGCCCTGCTGGCTGCGATATTTCGCATCTTCACGGCGGTTGTAAGGTCGCa                                     >  1:258327/1‑91 (MQ=255)                                      aTGGGCCTGAATTAGTCTTTATCGATTCGGCTGGCTACCGCGCCCTGCTGGCTGCGATATTTCGCATCTTCACGGCGGTTGTAAGGTCGCa                                     >  1:464517/1‑91 (MQ=255)                                         ggCCTGAATTAGTCTTTATCGATTCGGCTGGCTACCGCGCCCTGCTGGCTGCGATATTTCGCATCTTCACGGCGGTTGTAAGGTTGCAccg                                  >  1:597279/1‑91 (MQ=255)                                          gCCTGAATTAGTCTTTATCGATTCGGCTGGCTACCGCGCCCTGCTGGCTGCGATATTTCGCATCTTCACGGCGGTTGTAAGGTCGCAccgc                                 >  1:615627/1‑91 (MQ=255)                                          gCCTGAATTAGTCTTTATCGATTCGGCTGGCTACCGCGCCCTGCTGGCTGCGATATTTCGCATCTTCACGGCGGTTGTAAGGTCGCAccgc                                 <  1:366882/91‑1 (MQ=255)                                           ccTGAATTAGTCTTTATCGATTCGGCTGGCTACCGCGCCCTGCTGGCTGCGATATTTCGCATCTTCACGGCGGTTGTAAGGTCGCAccgcc                                >  2:127436/1‑91 (MQ=255)                                           ccTGAATTAGTCTTTATCGATTCGGCTGGCTACCGCGCCCTGCTGGCTGCGATATTTCGCATCTTCACGGCGGTTGTAAGGTCGCAccgcc                                >  1:27129/1‑91 (MQ=255)                                                   aGTCTTTATCGATTCGGCTGGCTACCGCGCCCTGCTGGCTGCGATATTTCGCATCTTCACGGCGGTTGTAAGGTCGCACCGCCGGGCCGGa                        >  1:552598/1‑91 (MQ=255)                                                        ttATCGATTCGGCTGGCTACCGCGCCCTGCTGGCTGCGATATTTCGCATCTTCACGGCGGTTGTAAGGTCGCACCGCCGGGCCGGAAAGCg                   >  1:15869/1‑91 (MQ=255)                                                        ttATCGATTCGGCTGGCTACCGCGCCCTGCTGGCTGCGATATTTCGCATCTTCACGGCGGTTGTAAGGTCGCACCGCCGGGCCGGAAAGCg                   >  1:267608/1‑91 (MQ=255)                                                        ttATCGATTCGGCTGGCTACCGCGCCCTGCTGGCTGCGATATTTCGCATCTTCACGGCGGTTGTAAGGTCGCACCGCCGGGCCGGAAAGCg                   >  1:410828/1‑91 (MQ=255)                                                        ttATCGATTCGGCTGGCTACCGCGCCCTGCTGGCTGCGATATTTCGCATCTTCACGGCGGTTGTAAGGTCGCACCGCCGGGCCGGAAAGCg                   >  1:624520/1‑91 (MQ=255)                                                            cGATTCGGCTGGCTACCGCGCCCTGCTGGCTGCGATATTTCGCATCTTCACGGCGGTTGTAAGGTCGCACCGCCGGGCCGGAAAGCGGCTc               <  2:248173/91‑1 (MQ=255)                                                            cGATTCGGCTGGCTACCGCGCCCTGCTGGCTGCGATATTTCGCATCTTCACGGCGGTTGTAAGGTCGCACCGCCGGGCCGGAAAGCGGCTc               <  2:358984/91‑1 (MQ=255)                                                              aTTCGGCTGGCTACCGCGCCCTGCTGGCTGCGATATTTCGCATCTTCACGGCGGTTGTAAGGTCGCACCGCCGGGCCGGAAAGCGGCTCaa             <  2:565029/91‑1 (MQ=255)                                                                   gctggctACCGCGCCCTGCTGGCTGCGATATTTCGCATCTTCACGGCGGTTGTAAGGTCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTc        <  2:436897/91‑1 (MQ=255)                                                                      ggctACCGCGCCCTGCTGGCTGCGATATTTCGCATCTTCACGGCGGTTGTAAGGTCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAgc     <  2:27129/91‑1 (MQ=255)                                                                      ggctACCGCGCCCTGCTGGCTGCGATATTTCGCATCTTCACGGCGGTTGTAAGGTCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAgc     <  2:615627/91‑1 (MQ=255)                                                                        ctACCGCGCCCTGCTGGCTGCGATATTTCGCATCTTCACGGCGGTTGTAAGGTCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAgcgc   >  1:257145/1‑91 (MQ=255)                                                                         tACCGCGCCCTGCTGGCTGCGATATTTCGCATCTTCACGGCGGTTGTAAGGTCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCa  <  1:17833/91‑1 (MQ=255)                                                                             gcgcCCTGCTGGCTGCGATATTTCGCATCTTCACGGCGGTTGTAAGGTCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCa  <  1:842708/87‑1 (MQ=255)                                                                             gcgcCCTGCTGGCTGCGATATTTCGCATCTTCACGGCGGTTGTAAGGTCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCa  >  2:842708/1‑87 (MQ=255)                                                                                  ctgctgGCTGCGATATTTCGCATCTTCACGGCGGTTGTAAGGTCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAgcgc   >  2:541011/1‑81 (MQ=255)                                                                                  ctgctgGCTGCGATATTTCGCATCTTCACGGCGGTTGTAAGGTCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAgcgc   <  1:541011/81‑1 (MQ=255)                                                                                         |                                                                           CATCAGCGTCGTCAGAAATCGTCTCATGGTATCCTCAATGGGCCTGAATTAGTCTTTATCGATTCGGCTGGCTACCGCGCCCTGCTGGTTGCGATATTTCGCATCTTCACGGCGGTTGTAAGGTCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCA  >  NZ_CP009273/2135067‑2135229

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