Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	404,238	2 bp→TT	coding (536‑537/1116 nt)	dgcC →	diguanylate cyclase DgcC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	404,238	0	C	T	100.0%	198.9 / NA	55	P179L (CCC→CTT)	dgcC	diguanylate cyclase DgcC
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (31/24);  total (31/24)
*	NC_000913	404,239	0	C	T	100.0%	198.3 / NA	55	P179L (CCC→CTT)	dgcC	diguanylate cyclase DgcC
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (31/24);  total (31/24)

GCCTGGCAGATAGCGAGCAGGGCCGTCGATCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAACATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGACGGGCGTGTATAACCGACGTCATTGGGAAACTATGTTACGCAATGAATTTGATAACTGTCGGCGGCATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACA  >  NC_000913/403952‑404523                                                                                                                                                                                                                                                                                               ||                                                                                                                                                                                                                                                                                             gCCTGGCAGATAGCGAGCAGGGCCGTCGATCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCTTATTATTGTCAtt                                                                                                                                                                                                                                                                                  >  1:466269/1‑300 (MQ=255)     ggCAGATAGCGAGCAGGGCCGTCGATCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACGCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCTTATTATTGTCCTTTATc                                                                                                                                                                                                                                                                              >  2:701355/1‑300 (MQ=255)     ggCAGATAGCGAGCAGGGCCGTCGATCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCTTATTATTGTCATTTATc                                                                                                                                                                                                                                                                              >  1:161722/1‑300 (MQ=255)      gcaggtggCGAGCAGGGCCGTCAATCCGCTTGGCCGGGAAATTTACAACTTAGAAACCGATTCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCTTATTATTGTCATTTATcc                                                                                                                                                                                                                                                                             <  2:76812/293‑1 (MQ=255)      gCAGATAGCGAGCAGGGCCGTCGATCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCTTATTATTGTCATTTATcc                                                                                                                                                                                                                                                                             <  1:297611/300‑1 (MQ=255)          aTAGCGAGCAGGGCCGTCGATCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGCATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCTTATTATTGTCATTTATCCTctg                                                                                                                                                                                                                                                                         >  1:268806/1‑300 (MQ=255)           tAGCGAGCAGGGCCGTCGATCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGTTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCTTATTATTGTCATTTATCCTctgc                                                                                                                                                                                                                                                                        <  1:643184/300‑1 (MQ=255)                  cagtgCCGTTGATCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAATATTAGCGGGAATGTGGGTAGGCGGAATGGGCGTAAATGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGCCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACTGTGTCGTTCAATAGTGCACCGCTGGAATGGTGGCTCTCCCTTCTTATTATTGTCATTTATCCTCTGCTGTTTgg                                                                                                                                                                                                                                                                 <  2:350088/296‑1 (MQ=255)                    gggCCGTCGATCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGNCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGNTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCTTATTATTGTCATTTATCCTCTGCTGTTTGGCt                                                                                                                                                                                                                                                               >  2:645298/1‑300 (MQ=255)                                            aaaTTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCTTATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCa                                                                                                                                                                                                                                       <  2:806425/300‑1 (MQ=255)                                                                      gTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGNGCCGCTGGAATGGTGGCTCTCCCTTCTTATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAACATAAACGCAGGTTg                                                                                                                                                                                                             >  1:559864/1‑300 (MQ=255)                                                                      gTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCTTATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAACATAAACGCAGGTTg                                                                                                                                                                                                             >  1:318293/1‑300 (MQ=255)                                                                      gTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCTTATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTAGCGGAACATAAACGCAGGTTg                                                                                                                                                                                                             >  1:127455/1‑300 (MQ=255)                                                                            gCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCTTATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTCCCAGACGGCAACCCAACTGGCGGAACATAAACGCAGGTTGCAGGTc                                                                                                                                                                                                       >  2:228176/1‑300 (MQ=255)                                                                                    gtgGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCTTATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAACATAAACGCAGGTTGCAGGTCATGAGTAc                                                                                                                                                                                               >  1:113177/1‑300 (MQ=255)                                                                                           ggCGTAGTGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATCGNGCGCCGCTGGAATGGTGGCTCTCCCTTCTTATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAACATAAACGCAGGTTGCAGGTCGTGAGTAACCGCGAc                                                                                                                                                                                        >  2:195731/1‑300 (MQ=255)                                                                                              gTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCTTATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAACATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGc                                                                                                                                                                                     <  1:795559/300‑1 (MQ=255)                                                                                               tAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCTTATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAACATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCa                                                                                                                                                                                    >  1:614201/1‑300 (MQ=255)                                                                                                                  gcCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGTTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGTTTGCTCTCCCTTCTTATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAACATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGACGGGCGTGTATAACCg                                                                                                                                                                 <  2:264632/300‑1 (MQ=255)                                                                                                                  gcCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCATTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCTTATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAAACAAACTGGCGGAACATAAACGCAGGTTGCAGGTCATGAGTAACCGCGACGGCATGACGGGCGTGTATAACCg                                                                                                                                                                 >  2:433210/1‑300 (MQ=255)                                                                                                                    cTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCAACCTCGAGCTGACGGGCATTACTGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCTTATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAACATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGACGGGCGTGTATAACCGAc                                                                                                                                                               >  1:219776/1‑300 (MQ=255)                                                                                                                                 gTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCTTATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAACATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGACGGGCGTGTATAACCGACGTCATTGGGAAAc                                                                                                                                                  >  2:11106/1‑300 (MQ=255)                                                                                                                                 gTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCTTATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAACATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGACGGGCGTGTATAACCGACGTCATTGGGAAAc                                                                                                                                                  >  1:786660/1‑300 (MQ=255)                                                                                                                                              tgtCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCTTATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAACATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGACGGGCGTGTATAACCGACGTCATTGGGAAACTATGTTACGCAAt                                                                                                                                     >  2:25958/1‑300 (MQ=255)                                                                                                                                                   gAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCTTATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAACATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGAAGAGCGTGTATAACCGACGTCATTGGGAAGCTATGTTACGCACTGAAtt                                                                                                                                >  2:295122/1‑300 (MQ=255)                                                                                                                                                          attgtggCAGGCGTCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCTTATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAACATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGACGGGCGTGTATAACCGACGTCATTGGGAAACTATGTTACGCAATGAATTTGATAAc                                                                                                                         <  2:109521/295‑1 (MQ=255)                                                                                                                                                            gggggCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCTTATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAACATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGACGGGCGTGTATAACCGACGTCATTGGGAAACTATGTTACGCAATGAATTTGATAACTg                                                                                                                       <  2:128196/300‑1 (MQ=255)                                                                                                                                                              gggCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCTTATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAACATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGACGGGCGTGTATAACCGACGTCATTGGGAAACTATGTTACGCAATGAATTTGATAACTGTc                                                                                                                     >  1:538552/1‑300 (MQ=255)                                                                                                                                                                           ccGTCTGTTTTTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCTTATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAACATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGACGGGCGTGTATAACCGACGTCATTGGGAAACTATGTTACGCAATGAATTTGATAACTGTCGGCGGCATAATcg                                                                                                        <  1:528863/300‑1 (MQ=255)                                                                                                                                                                                     gTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCTTATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAACATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGACGGGCGTGTATAACCGACGTCATTGGGAAACTATGTTACGCAATGAATTTGATAACTGTCGGCGGCATAATCGCGATGCAACg                                                                                              >  1:283929/1‑300 (MQ=255)                                                                                                                                                                                      tCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCTTATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAACATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGACGGGCGTGTATAACCGACGTCATTGGGAAACTATGTTACGCAATGAATTTGATAACTGTCGGCGGCATAATCGCGATGCAACGt                                                                                             <  2:289065/300‑1 (MQ=255)                                                                                                                                                                                           ggTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGCGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTATTATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAACATAGACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGAAGGGCGTGTATAACCGAAGTCATTGGGAACCTATGTTACGCAATGAATTTGATAACTGTCGGAAGCGTAATCGCGATGCAACGTTACTg                                                                                        >  2:241247/1‑300 (MQ=255)                                                                                                                                                                                            gTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCTTATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAACATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGACGGGCGTGTATAACCGACGTCATTGGGAAACTATGTTACGCAATGAATTTGATAACTGTCGGCGGCATAATCGCGATGCAACGTTACTGa                                                                                       >  1:742503/1‑300 (MQ=255)                                                                                                                                                                                                 gtgtTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCAGTACCGTGTCGTTCAATAGTGCGCCGCTTGAATGGTGGCTCTCCCTTCTTATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAACATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGACGGGCGTGTATAACCGACGTCATTGGGAAACTATGTTACGCAATGAATTTGATAACTGTCGGCGGCATAATCGCGATGCAACGTTACTGATtatc                                                                                  <  2:409229/300‑1 (MQ=255)                                                                                                                                                                                                   gtTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCTTATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAACATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGACGGGCGTGTATAACCGACGTCATTGGGAAACTATGTTACGCAATGAATTTGATAACTGTCGGCGGCCTAATCGCGATGCAACGTTAATGATtatcga                                                                                >  2:398226/1‑300 (MQ=255)                                                                                                                                                                                                    ttGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCTTATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAACATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGACGGGCGTGTATAACCGACGTCATTGGGAAACTATGTTACGCAATGAATTTGATAACTGTCGGCGGCATAATCGCGATGCAACGTTACTGATtatcgat                                                                               <  1:613342/300‑1 (MQ=255)                                                                                                                                                                                                           tggtTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCTTATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAACATAAACGCAGGTTGCAGGTCATGAGTACCCTCGACGGCATGACGGGCGTGTATAACCGACGTCATTGGGAAACTATGTTACGCAATGAATTTGATAACTGTCGGCGGCATAATCGCGATGCAACGTTACTGATTATCGATATCGAcc                                                                        <  2:362972/300‑1 (MQ=255)                                                                                                                                                                                                               ttCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCTTATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAACATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGACGGGCGTGTATAACCGACGTCATTGGGAAACTATGTTACGCAATGAATTTGATAACTGTCGGCGGCATAATCGCGATGCAACGTTACTGATTATCGATATCGACCAttt                                                                    >  1:772886/1‑300 (MQ=255)                                                                                                                                                                                                                     ccTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCTTATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAACATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGACGGGCGTGTATAACCGACGTCATTGGGAAACTATGTTACGCAATGAATTTGATAACTGTCGGCGGCATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAagag                                                              <  1:768966/300‑1 (MQ=255)                                                                                                                                                                                                                      cTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCTTATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAACATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGACGGGCGTGTATAACCGACGTCATTGGGAAACTATGTTACGCAATGAATTTGATAACTGTCGGCGGCATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGc                                                             >  2:747372/1‑300 (MQ=255)                                                                                                                                                                                                                              ccTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCTTATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAACATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGACGGGCGTGTATAACCGACGTCATTGGGAAACTATGTTACGCAATGAATTTGATAACTGTCGGCGGCATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGa                                                     >  2:410752/1‑300 (MQ=255)                                                                                                                                                                                                                                 cGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTTGCTCTCCCTTCTTATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAACATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGACGGGCGTGTATAACCGACGTCATTGGGAAACTATGTTACGCAATGAATTTGATAACTGTCGGCGGCATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATAc                                                  <  2:324335/300‑1 (MQ=255)                                                                                                                                                                                                                                  gAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCTTATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGANGGCAACCAAACTGGCGGAACATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGACGGGCGTGTATAACCGACGTCATTGGGAAACTATGTTACGCAATGAATTTGATAACTGTCGGCGGCATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATAcc                                                 <  1:580697/300‑1 (MQ=255)                                                                                                                                                                                                                                   aGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCTTATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAACATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGACGGGCGTGTATAACCGACGTCATTGGGAAACTATGTTACGCAATGAATTTGATAACTGTCGGCGGCATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCt                                                <  1:67733/300‑1 (MQ=255)                                                                                                                                                                                                                                        aCGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCTTATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAACATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGACGGGCGTGTATAACCGACGTCATTGGGAAACTATGTTACGCAATGAATTTGATAACTGTCGGCGGCATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGc                                           <  1:428139/300‑1 (MQ=255)                                                                                                                                                                                                                                          gggCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCTTATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAACATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGACGGGCGTGTATAACCGACGTCATTGGGAAACTATGTTACGCAATGAATTTGATAACTGTCGGCGGCATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCa                                         >  2:365851/1‑300 (MQ=255)                                                                                                                                                                                                                                            gCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCTTATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAACATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGACGGGTNTGTATAACCGACGTCATTGGGAAACTATGTTACGCAATGAATTTGATCACTGTCGGCGGCATAATCGCGATGCAACGTTAATGCTTATCGATATCGACCATTTCAGGCGCATCAACGATGCCTAGGGCCata                                       >  2:188979/1‑299 (MQ=255)                                                                                                                                                                                                                                            gCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCTTATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAACATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGACGGGCGTGTATAACCGACGTCATTGGGAAACTATGTTACGCAATGAATTTGATAACTGTCGGCGGCATAATCGCGATGCAACGTTACTGCTTATTGATATCGACCATTTAAAGAGCATCGACGATACCTGGGGCCatg                                       >  2:309210/1‑300 (MQ=255)                                                                                                                                                                                                                                            gCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCTTATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAACATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGACGGGCGTGTATAACCGACGTCATTGGGAAACTATGTTACGCAATGAATTTGATAACTGTCGGCGGCATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTAAAGAGCATCAACGATACCTGGGGCCatg                                       >  1:293972/1‑300 (MQ=255)                                                                                                                                                                                                                                               ttACCGTGTCGTTCAATAGTGCGCCGCTTGAATGGTTGCTCTCCCTTCTTATTACTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAACATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGACGGGCGTGTATAACCGACGTCATTGGGAAACTATGTTACGCAATGAATTTGATAACTGTCGGCGGCATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCatgatg                                    <  2:780901/300‑1 (MQ=255)                                                                                                                                                                                                                                               ttACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCTTATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAACATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGACGGGCGTGTATAACCGACGTCATTGGGAAACTATGTTACGCAATGAATTTGATAACTGTCGGCGGCATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCatgatg                                    <  1:249668/300‑1 (MQ=255)                                                                                                                                                                                                                                                                aGTGCGCCGCTGGAATGGTGGCTCTCCCTTCTTATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAACATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGACGGGCGTGTATAACCGACGTCATTGGGAAACTATGTTACGCAATGAATTTGATAACTGTCGGCGGCATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGAtt                   <  1:175102/300‑1 (MQ=255)                                                                                                                                                                                                                                                                 gTGCGCCGCTGGAATGGTGGCTCTCCCTTCTTATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAACATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGACGGGCGTGTATAACCGACGTCATTGGGAAACTATNTTACGCAATGAATTTGATAACTGTCGGCGGCATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACNTGGGGCCATGATGTGGGCGATGAAGCGATtg                  <  2:592397/300‑1 (MQ=255)                                                                                                                                                                                                                                                                      ccgcTGGAATGGTGGCTCTCCCTTCTTATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAACATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGACGGGCGTGTATAACCGACGTCATTGGGAAACTATGTTACGCAATGAATTTGATAACTGTCGGCGGCATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGgcg             <  1:664975/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                 gtggCTCTCCCTTCTTATTATTGTCATTTATNCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAACATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGNCGGGCGTGTATAACCGACGTCATTGGGAAACTATGTTACGCAATGAATTTGATAACTGTCGGCGGCATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGaaa  >  2:572372/1‑298 (MQ=255)                                                                                                                                                                                                                                                                                               ||                                                                                                                                                                                                                                                                                             GCCTGGCAGATAGCGAGCAGGGCCGTCGATCCGCTTAGCCGGGAAATTTACAACTTAAAAACCGATGCAGTATTAGCGGGAATGTGGGTAGGCGTAATGGGCGTAAACGTGCTGCCTTCCACCGCGATGTTGATGATTATGTGTCTGAATTTGATGGGGGCAGGCGGCCCCCGTCTGTTTGTCGCGGGTCTGGTGTTGATGGTGGTTTCCTGCCTTGTCACCCTCGAGCTGACGGGCATTACCGTGTCGTTCAATAGTGCGCCGCTGGAATGGTGGCTCTCCCTTCCCATTATTGTCATTTATCCTCTGCTGTTTGGCTGGGTCAGCTACCAGACGGCAACCAAACTGGCGGAACATAAACGCAGGTTGCAGGTCATGAGTACCCGCGACGGCATGACGGGCGTGTATAACCGACGTCATTGGGAAACTATGTTACGCAATGAATTTGATAACTGTCGGCGGCATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACA  >  NC_000913/403952‑404523

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