Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AE015451	4,714,831	C→A	A21D (GCC→GAC)	PP_4173 →	Sensor histidine kinase/response regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	4,714,831	0	C	A	100.0%	205.9 / NA	61	A21D (GCC→GAC)	PP_4173	Sensor histidine kinase/response regulator
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (32/29);  total (32/29)

CGCCGGCGCCCTGGCCAGTACCTACAGCGACCAGGACGACTGGCTGTGGGTGCTCGACCAGTTGCTCGACCAGCCGCCGACGCGCTGGCCGCGCAGCCTGTACCCGGTGCGTTACGGCGTCAGTGGCAACCAGCAGGTGGCCCGTGCCCTGGGTCTTGAACCAATCGATCAGAACGCCGCCGTCAAGGCCCTGGCCGAGGGAGAACCCACCCCATGAGCAAACGCTTGAGCTGGGACATTCATACCCGCACGCAGATCATCAGTCTCGGCCCCGCCCTGCTACTGACCTTGCTCTTGATCAGCTTCTTTACCTTCGTACGCATCCAGGACCTGCGCCAGGAACTCAACCACACCGGGCAACTGATCGCCAACCAGTTGGCACCGGCGTCGGAATATGGGGTGATCTCGGGCAACAACGAAGTGCTTGAAGGGCTGATGCGTGCCACCCTCAGCATCCCGCACGTGCGCTTTCTGGAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTTCCAGG  >  AE015451/4714557‑4715121                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                                                   cgcCGGCGACCTGGCCAGTACCTACAGCGACCAGGACGACTGGCTGTGGGTGCTCGACCAGTTGCTCGACCAGCCGCCGACGCGCTGGCCGCGCAGCCTGTACCCGGTGCGTTACGGCGTCAGTGGCAACCAGCAGGTGGCCCGTGCCCTGGGTCTTGAACCAATCGATCAGAACGCCGCCGTCAAGGCCCTGGCCGAGGGAGAACCCACCCCATGAGCAAACGCTTGAGCTGGGACATTCATACCCGCACGCAGATCATCAGTCTCGGCCCCGACCTGCTACTGACCTTGCTCTTGATCa                                                                                                                                                                                                                                                                          <  1:416696/301‑1 (MQ=255)      ccgcCCTGGCCAGTACCGACAGCGACCAGGACGACTGGCTGTGGGTGCTCGACCAGTTGCTCGACCAGCCGCCGACGCGCTGGCCGCGCAGCCTGTACCCGGTGCGTTACGGCGTCAGTGGCAACCAGCAGGTGGCCCGTGCCCTGGGTCTTGAACCAATCGATCAGAACGCCGCCGTCAAGGCCCTGGCCGAGGGAGAACCCACCCCATGAGCAAACGCTTGAGCTGGGACATTCATACCCGCACGCAGATCATCAGTCTCGGCCCCGACCTGCTACTGACCTTGCTCTTGATCAGcttc                                                                                                                                                                                                                                                                     <  1:500276/300‑1 (MQ=255)          ccTGGCCAGTACCTACAGCGACCAGGACGACTGGCTGTGGGTGCTCGACCAGTTGCTCGACCAGCCGCCGACGCGCTGGCCGCGCAGCCTGTACCCGGTGCGTTACGGCGTCAGTGGCAACCAGCAGGTGGCCCGTGCCCTGGGTCTTGAACCAATCGATCAGAACGCCGCCGTCAAGGCCCTGGCCGAGGGAGAACCCACCCCATGAGCAAACGCTTGAGCTGGGACATTCATACCCGCACGCAGATCATCAGTCTCGGCCCCGACCTGCTACTGACCTTGCTCTTGATCAGCTTCTTTa                                                                                                                                                                                                                                                                 <  1:12445/301‑1 (MQ=255)             ggCCAGTACCTACAGCGACCAGGACGACTGGCTGTGGGTGCTCGACCAGTTGCTCGACCAGCCGCCGACGCGCTGGCCGCGCAGCCTGTACCCGGTGCGTTACGGCGTCAGTGGCAACCAGCAGGTGGCCCGTGCCCTGGGTCTTGAACCAATCGATCAGAACGCCGCCGTCAAGGCCCTGGCCGAGGGAGAACCCACCCCATGAGCAAACGCTTGAGCTGGGACATTCATACCCGCACGCAGATCATCAGTCTCGGCCCCGACCTGCTACTGACCTTGCTCTTGATCAGCTTCTTTACCt                                                                                                                                                                                                                                                              <  1:542974/301‑1 (MQ=255)             ggCCAGTACCTACAGCGACCAGGACGACTGGCTGTGGGTGCTCGACCAGTTGCTCGACCAGCCGCCGACGCGCTGGCCGCGCAGCCTGTACCCGGTGCGTTACGGCGTCAGTGGCAACCAGCAGGTGGCCCGTGCCCTGGGTCTTGAACCAATCGATCAGAACGCCGCCGTCAAGGCCCTGGCCGAGGGAGAACCCACCCCATGAGCAAACGCTTGAGCTGGGACATTCATACCCGAACGCAGATCATCAGTCTCGGCCCCGACCTGCTACTGACCTTGCTCTTGATCAGCTTCTTTACCt                                                                                                                                                                                                                                                              >  2:416101/1‑301 (MQ=255)                 aagaccgaccgCGACCAAGACGAATGGCTGTGGGTGCTCGACCAGGTGCTCGACCAGCCGCCGACGCGCTGGGCGCACAGCCTGTAACCGGTGCGTTACGGCGTCAGTGGCAACCAGCAGGTGGCCCGTGCCCTGGGTCTTGAACCAATCGATCAGAACGCCGCCGTCAAGGCCCTGGCCGAGGGAGAACCCACCCCATGAGCAAACGCTTGAGCTGGGACATTCATACCCGCACGCAGATCATCAGTCTCGGCCCCGACCTGCTACTGACCTTGCTCTTGATCAGCTTCTTTACCTTCGt                                                                                                                                                                                                                                                          <  2:539584/291‑1 (MQ=255)                         cAGCGACCAGGACGACTGGCTGTGGGTGCTCGACCAGTTGCTCGACCAGCCGCCGACGCGCTGGCCGCGCAGCCTGTACCCGGTGCGTTACGGCGTCAGTGGCAACCAGCAGGTGGCCCGTGCCCTGGGTCTTGAACCAATCGATCAGAACGCCGCCGTCAAGGCCCTGGCCGAGGGAGAACCCACCCCATGAGCAAACGCTTGAGCTGGGACATTCATACCCGCACGCAGATCCTCAGTCTCGGCCCCGACCTGCTACTGACCTTGCCCTTGATCCGCTTCTTTACCTTCGTACGCATcc                                                                                                                                                                                                                                                  >  2:211911/1‑301 (MQ=255)                         cAGCGACCAGGACGACTGGCTGTGGGTGCTCGACCAGTTGCTCGACCAGCCGCCGACGCGCTGGCCGCGCAGCCTGTACCCGGTGCGTTACGGCGTCAGTGGCAACCAGCAGGTGGCCCGTGCCCTGGGTCTTGAACCAATCGATCAGAACGCCGCCGTCAAGGCCCTGGCCGAGGGAGAACCCACCCCATGAGCAAACGCTTGAGCTGGGACATTCATACCCGCACGCAGATCATCAGTCTCGGCCCCGACCTGCTACTGACCTTGCTCTTGATCAGCTGCTTTACCTTCGTACGCATcc                                                                                                                                                                                                                                                  >  2:505410/1‑301 (MQ=255)                            cGACCAGGACGACTGGCTGTGGGTGCTCGACCAGTTGCTCGACCAGCCGCCGACGCGCTGGCCGCGCAGCCTGTACCCGGTGCGTTACGGCGTCAGTGGCAACCAGCAGGTGGCCCGTGCCCTGGGTCTTGAACCAATCGATCAGAACGCCGCCGTCAAGGCCCTGGCCGAGGGAGAACCCACCCCATGAGCAAACGCTTGAGCTGGGACATTCATACCCGCACGCAGATCATCAGTCTCGTCCCCGACCTGCTACTGACCTTGCTCTTGATCAGCTTCTTTACCTTCGTACTCATCCAgg                                                                                                                                                                                                                                               >  2:225359/1‑301 (MQ=255)                                  ggACGACTGGCTGTGGGTGCTCGACCAGTTGCTCGACCAGCCGCCGACGCGCTGGCCGCGCAGCCTGTACCCGGTGCGTTACGGCGTCAGTGGCAACCAGCAGGTGGCCCGTGCCCTGGGTCTTGAACCAATCGATCAGAACGCCGCCGTCAAGGCCCTGGCCGAGGGAGAACCCACCCCATGAGCAAACGCTTGAGCTGGGACATTCATACCCGCACGCAGATCATCAGTCTCGGCCCCGACCTGCTACTGACCTTGCTCTTGATCAGCTTCTTTACCTTCGTACGCATCCAGGACCttc                                                                                                                                                                                                                                         >  2:308571/1‑299 (MQ=255)                                  ggACGACTGGCTGTGGGTGCCCGACCAGTTGCTCGACCAGCCGCCGACGCGCTGGCCGCGCAGCCTGTACCCGGTGCGTTACGGCGTCAGTGGCAACCAGCCGGTGGCCCGTGCCCTGGGTCTTGAACCAATCGATCAGAACGCCGCCGTCAAGGCTCTGGCCGCGGGAGAACCCACCCCATGAGCACACGCTTGAGCTGGGACATTCATAACCGCACGCAGATCATCAGTCTAGGCCGCGACCTGCTACTGACCTTCCTCTTGATCAGCTTCTTCACCTTCGTACTCATCCAGGACCTgc                                                                                                                                                                                                                                         >  2:472900/1‑301 (MQ=255)                                        cTGGCTGTGGGTGCTCGACCAGTTGCTCGACCAGCCGCCGACGCGCTGGCCGCGCAGCCTGTACCCGGTGCGTTACGGCGTCAGTGGCAACCAGCAGGTGGCCCGTGCCCTGGGTCTTGAACCAATCGATCAGAACGCCGCCGTCAAGGCCCTGGCCGAGGGAGAACCCACCCCATGAGCAAACGCTTGAGCTGGGACATTCATACCCGCACGCAGATCATCAGTCTCGGCCCCGACCTGCTACTGACCTTGCTCTTGATCAGCTTCTTTACCTTCGTACGCATCCAGGACCTGCGCCAgg                                                                                                                                                                                                                                   <  2:215866/301‑1 (MQ=255)                                           gCTGTGGGTGCTCGACCAGTTGCTCGACCAGCCGCCGACGCGCTGGCCGCGCAGCCTGTACCCGGTGCGTTACGGCGTCAGTGGCAACCAGCAGGTGGCCCGTGCCCTGGGTCTTGAACCAATCGATCAGAACGCCGCCGTCAAGGCCCTGGCCGAGGGAGAACCCACCCCATGAGCAAACGCTTGAGCTGGGACATTCATACCCGCACGCAGATCATCAGTCTCGGCCCCGACCTGCTACTGACCTTGCTCTTGATCAGCTTCTTTACCTTCGTACGCATCCAGGACCTGCGCCAGGAAc                                                                                                                                                                                                                                >  1:545991/1‑301 (MQ=255)                                                 gtgcgccgacaggtgggtCGACCAGCCGCCGACGCGCTGGCCGCGCAGCCTGTACCCGGTGCGTTACGGCGTCAGTGGCAACCAGCAGGTGGCCCGTGCCCTGGGTCTTGAACCAATCGATCAGAACGCCGCCGTCAAGGCCCTGGCCGAGGGAGAACCCACCCCATGAGCAAACGCTTGAGCTGGGACATTCATACCCGCACGCAGATCATCAGTCTCGGCCCCGACCTGCTACTGACCTTGCTCTTGATCAGCTTCTTTACCTTCGTACGCATCCAGGACCTGCGCCAGGAACTCAAcc                                                                                                                                                                                                                          <  2:271397/284‑1 (MQ=255)                                                     cTCGACCAGTTGCTCGACCAGCCGCCGACGCGCTGGCCGCGCAGCCTGTACCCGGTGCGTTACGGCGTCAGTGGCAACCAGCAGGTGGCCCGTGCCCTGGGTCTTGAACCAATCGATCAGAACGCCGCCGTCAAGGCCCTGGCCGAGGGAGAACCCACCCCATGAGCAAACGCTTGAGCTGGGACATTCATACCCGCACGCAGATCATCAGTCTCGGCCCCGACCTGCTACTGACCTTGCTCTTGATCAGCTTCTTTACCTTCGTACGCATCCAGGACCTGCGCCAGGAACTCAACCacac                                                                                                                                                                                                                      <  1:308690/301‑1 (MQ=255)                                                         aCCAGTTGCTCGACCAGCCGCCGACGCGCTGGCCGCGCAGCCTGTACCCGGTGCGTTACGGCGTCAGTGGCAACCAGCAGGTGGCCCGTGCCCTGGGTCTTGAACCAATCGATCAGAACGCCGCCGTCAAGGCCCTGGCCGAGGGAGAACCCACCCCATGAGCAAACGCTTGAGCTGGGACATTCATACCCGCACGCAGATCATCAGTCTCGGCCCCGACCTGCTACTGACCTTGCTCTTGATCAGCTTCTTTACCTTCGTACGCATCCAGGACCTGCGCCAGGAACTCAACCACACCggg                                                                                                                                                                                                                  >  1:226482/1‑301 (MQ=255)                                                             gTTGCTCGACCAGCCGCCGACGCGCTGGCCGCGCAGCCTGTACCCGGTGCGTTACGGCGTCAGTGGCAACCAGCAGGTGGCCCGTGCCCTGGGTCTTGAACCAATCGATCAGAACGCCGCCGTCAAGGCCCTGGCCGAGGGAGAACCCACCCCATGAGCAAACGCTTGAGCTGGGACATTCATACCCGCACGCAGATCATCAGTCTCGGCCCCGACCTGCTACTGACCTTGCTCTTGATCAGCTTCTTTACCTTCGTACGCATCCAGGACCTGCGCCAGGAACTCAACCACACCGGGCAAc                                                                                                                                                                                                              <  1:505618/301‑1 (MQ=255)                                                                   cgaacAGCCGCCGACGCGCTGGCCGCGCAGGCTGTACCCGGTGCGTTACGGCGTCAGTGGCAACCAGCAGGTGGCCCGTGCCCTGGGTCTTGAACCAATCGATCAGAACGCCGCCGTCAAGGCCCTGGCCGAGGGAGAACCCACCCCATGAGCAAACGCTTGAGCTGGGACATTCATACCCGCACGCAGATCATCAGTCTCGGCCCCGACCTGCTACTGACCTTGCTCTTGATCAGCTTCTTTACCTTCGTACGCATCCAGGACCTGCGCCAGGAACTCAACCACACCGGGCAACTGATCg                                                                                                                                                                                                        <  2:271810/297‑1 (MQ=255)                                                                     aCCAGCCGCCGACGCGCTGGCCGCGCAGCCTGTACCCGGTGCGTTACGGCGTCAGTGGCAACCAGCAGGTGGCCCGTGCCCTGGGTCTTGAACCAATCGATCAGAACGCCGCCGTCAAGGCCCTGGCCGAGGGAGAACCCACCCCATGAGCAAACGCTTGAGCTGGGACATTCATACCCGCACGCAGATCATCAGTCTCGGCCCCGACCTGCTACTGACCTTGCTCTTGATCAGCTTCTTTACCTTCGTACGCATCCAGGACCTGCTCCAGGAACTCAACCACCCCGGGCAACTGATCGcc                                                                                                                                                                                                      >  2:751888/1‑301 (MQ=255)                                                                      ccAGCCGCCGACGCGCTGGCCGCGCAGCCTGTACCCGATGCGTTACGGCGTCAGTGGCAACCAGCAGGTGGCCCGTGCCCTGGGTCTTGAACCAATCGATCAGAACGCCGCCGTCAAGGCCCTGGCCGAGGGAGAACCCACCCCATGAGCAAACGCTTGAGCTGGGACATTCATACCCGCACGCAGATCATCAGTCTCGGCCCCGACCTGCTACTGACCTTGCTCTTGATCAGCTTCTTTACCTTCGTACGCATCCAGGACCTGCGCCAGGAACTCAACCACACCGGGCAACTGATCGCCa                                                                                                                                                                                                     >  2:210640/1‑301 (MQ=255)                                                                        aggcggcgACGCGCTGGCCGCGCAGCATGTACCCGGTGCGTCACGGCGTCAGTGGCAACCAGCAGGTGGCCGGTGCCCTGTGTCTTGAACCACTCGATCAGAACGCCGTCGTCAAGGCGCTGGGCGAAGGAGAACCCACACCATGAGCAAACGCTTGAGCTGGGACATTCATGCCCGCACGCAGATCATCAGTCTCGGCGCCGACCTGCTACTGAACTTGCTCCTGATCAGCTTCTTTACCTTCGTACGCATCCAGGAGCTGCGCCAGGAACTCAACCACACCGGGCAACTGATCGCCAAc                                                                                                                                                                                                   <  2:76550/295‑1 (MQ=255)                                                                          ccgccgTCGCGCTGGCCGCTCAGCCTGTACCCGGTGCGTTACGGCGTCAGTGGCAACCAGCAGGTGGCCCGTGCCCTGGGTCTTGAACCAATCGATCAGAACGCCGCCGTCAAGGCCCTGGCCGAGGGAGAACCCACCCCATGAGCAAACGCTTGAGCTGGGACATTCATACCCGCACGCAGATCATCAGTCTCGGCCCCGACCTGCTACTGACCTTGCTCTTGATCAGCTTCTTTACCTTCGTACGCATCCAGGACCTGCGCCAGGAACTCAACCACACCGGGCAACTGATCGCCAACCa                                                                                                                                                                                                 <  1:473102/301‑1 (MQ=255)                                                                            gccgACGCGCTGGCCGCGCAGCCTGTACCCGGTGCGTTACGGCGTCAGTGGCAACCAGCAGGTGGCCCGTGCCCTGGGTCTTGAACCAATCGATCAGAACGCCGCCGTCAAGGCCCTGGCCGAGGGAGAACCCACCCCATGAGCAAACGCTTGAGCTGGGACATTCATACCCGCACGCAGATCATCAGTCTCGGCCCCGACCTGCTACTGACCTTGCTCTTGATCAGCTTCTTTACCTTCGTACGCATCCAGGACCTGCGCCAGGAACTCAACCCCACCGGGCAACTGATCGCCAACCAGt                                                                                                                                                                                               >  1:217479/1‑301 (MQ=255)                                                                              cgACGCGCTGGCCGCGCAGCCTGTACCCGGTGCGTTACGGCGTCAGTGGCAACCAGCAGGTGGCCCGTGCCCTGGGTCTTGAACCAATCGATCAGAACGCCGCCGTCAAGGCCCTGGCCGAGGGAGAACCCACCCCATGAGCAAACGCTTGAGCTGGGACATTCATACCCGCACGCAGATCATCAGTCTCGGCCCCGACCTGCTACTGACCTTGCTCTTGATCAGCTTCTTTACCTTCGTACGCATCCAGGACCTGCCCCAGGAACGCAACCACAACGGGCAACTGATACCCAACCacgtg                                                                                                                                                                                             >  2:536840/1‑297 (MQ=255)                                                                                       ggCCGCGCAGCCTGTACCCGGTGCGTTACGGCGTCAGTGGCAACCAGCAGGTGGCCCGTGCCCTGGGTCTTGAACCAATCGATCAGAACGCCGCCGTCAAGGCCCTGGCCGAGGGAGAACCCACCCCATGAGCAAACGCTTGAGCTGGGACATTCATACCCGCACGCAGATCATCAGTCTCGGCCCCGACCTGCTACTGACCTTGCTCTTGATCAGCTTCTTTACCTTCGTACGCATCCAGGACCTGCGCCAGGAACTCAACAAAACCGGGAAACTGATCGCCAACCAATTGGCAACGGCg                                                                                                                                                                                    >  2:248191/1‑301 (MQ=255)                                                                                        gCCGCGCAGCCTGTGCCCGGTGCGTTACGGCGTCAGTGGCAACCAGCAGGTGGCCCGTGCCCTGGGTCTTGAACCAATCGATCAGAACGCCGCCGTCAAGGCCCTGGCCGAGGGAGAACCCACCCCATGAGCAAACGCTTGAGCTGGGACATTCATACCCGCACGCAGATCATCAGTCTCGGCCCCGACCTGCTACTGACCTTGCTCTTGATCAGCTTCTTTACCTTCGTACGCATCCAGGACCTGCGCCAGGAACTCAACCACACCGGGCAACTGATCGCCAACCAGTTGGCACCGGCGt                                                                                                                                                                                   <  2:393843/301‑1 (MQ=255)                                                                                                           gtgcctttcGGCGTGAGTGCCAACCAGCAGGTGGCCCGTGCCCTGGGTCTTGAACCAATCGATCAGAACGCCGCCGTCAAGGCCCTGGCCGAGGGCGAACCCACCCCATGAGCAAACGCTTGAGCTGGGACATTCATACCCGCACGCAGATCATCAGTCTCGGCCCCGACCTGCTACTGACCTTGCTCTTGATCAGCTTCTTTACCTTCGTACGCATCCAGGACCTGCGCCAGGAACTCAACCACACCGGGCAACTGATCGCCAACCAGTTGGCACCGGCGTCGGAATATGGGGTGAtctc                                                                                                                                                                <  2:301586/293‑1 (MQ=255)                                                                                                              cGTTACGGCGTCAGTGGCAACCAGCAGGTGCCGCGTGCCCTGGGTCTTGAACCAATCGATCAGAACGCCGCCGTCAAGGCCCTGGCCGAGGGAGAACCCACCCCATGAGCAAACGATTGAGCTGGGACATTCATACCCGCACGCAGATAATCAGTCTCGGCCCCGACCTGCGACTGACCTTGCGCTTGATCAGGTTCGTTACCTTCGTAAGAATCCCGGACGTGCGACAGGAAATCAAACACACGCGGCAACTGATCGCCAAACAGGTGGAATCGGCGTGGGaagaacggggacgagcgag                                                                                                                                                             >  2:518368/1‑284 (MQ=255)                                                                                                                     gCGTCAGTGGCAACCAGCAGGTGGCCCGTGCCCTGGGTCTTGAACCAATCGATCAGAACGCCGCCGTCAAGGCCCTGGCCGAGGGAGAACCCACCCCATGAGCAAACGCTTGAGCTGGGACATTCATACCCGCACGCAGATCATCAGTCTCGGCCCCGACCTGCTACTGACCTTGCTCTTGATCAGCTTCTTTACCTTCGTACGCATCCAGGACCTGCGCCAGGAACTCAACCACACCGGGCAACTGATCGCCAACCAGTTGGCACCGGCGTCGGAATATGGGGTGATCTCGGGCaacaac                                                                                                                                                      <  1:136913/301‑1 (MQ=255)                                                                                                                              gCAACCAGCAGGTGGACCGGTCCCGGGGTCTTGGACCTATCGATCAGAACGCCGCCGTCAAGGCCCTGGCCGAGGGAGAACCCACCCCATGAGCAAACGCTTGAGCTGGGACATTCATACCCGCACGCAGATCATCAGTCTCGGCCCCGACCTGCTACTGACCTTGCTCTTGATCAGCTTCTTTACCTTCGTACGCATCCAGGACCTGCGCCAGGAACTCAACCACACCGGGCAACTGATCGCCAACCAGTTGGCACCGGCGTCGGAATATGGGGTGATCTCGGGCAACAACGAAGTGCtt                                                                                                                                             <  2:531521/301‑1 (MQ=255)                                                                                                                               cAACCAGCAGGTGGCCCGTGCCCTGGGTCTTGAACCAATCGATCAGAACGCCGCCGTCAAGGCCCTGGCCGAGGGAGAACCCACCCCATGAGCAAACGCTTGAGCTGGGACATTCATACCCGCACGCAGATCATCAGTCTCGGCCCCGACCTGCTACTGACCTTGCTCTTGATCAGCTTCTTTACCTTCGTACGCATCCAGGACCTGCGCCAGGAACTCAACCACACCGGGCAACTGATCGCCAACCAGTTGGCACCGGCGTCGGAATATGGGGTGATCTCGGGCAACAACGAAGTGCTTg                                                                                                                                            <  1:197003/301‑1 (MQ=255)                                                                                                                                      cagGTGGCCCGTGCCCTGGGTCTTGAACCAATCGATCAGAACGCCGCCGTCAAGGCCCTGGCCGAGGGAGAACCCACCCCATGAGCAAACGCTTGAGCTGGGACATTCATACCCGCACGCAGATCATCAGTCTCGGCCCCGACCTGCTACTGACCCTGCTCTTGATCAGCTTCTTTACCTTCGTACGCATCCAGGACCTGCGCCAGGAACTCAACCACACCGGGCAACTGATCGCCAACCAGTTGGCACCGGCGTCGGAATATGGGGTGATCGCGGCAAAAAACGAACGGCTGGAAGGGcc                                                                                                                                     >  2:635477/1‑300 (MQ=255)                                                                                                                                                        ggTCTTGAACCAATCGATCAGAACGCCGCCGTCAAGGCCCTGGCCGAGGGAGAACCCACCCCATGAGCAAACGCTTGAGCTGGGACATTCATACCCGCACGCAGATCATCAGTCTCGGCCCCGACCTGCTACTGACCTTGCTCTTGATCAGCTCCTTTACCTTCGTACGCACCCAGGACTTTCGCCAGTAACTCCACCACACCGGTCAACTGCTCGCCAACCAGCTGTCACCGCCGTCGGAATCGGGCGTGATCTCCGGCAACACCGACATGCTTGAGGGCCTGTTGCGTGCCCCCCTCac                                                                                                                   >  2:264793/1‑300 (MQ=255)                                                                                                                                                              gAACCAATCGATCAGAACGCCGCCGTCAAGGCCCTGGCCGAGTGAGAACCCACCCCATGAGCAAACGCTTGAGCTGGGACATTCATACCCGCACGCAGATCATCAGTCTCGGCCCCGACCTGCTACTGACCTTGCTCTTGATCAGCTTCTTTACCTTCGTACGCATCCAGGACCTGCGCCAGGAACTCAACCACACCGGGCAACTGATCGCCAACCAGTTGGCACCGGCGTCGGAATATGGGGTGATCGCGGGCAACAACGAAGTGCTTGAAGGGCTGATGCGTGCCACCCTCAGCATccc                                                                                                             >  2:156340/1‑301 (MQ=255)                                                                                                                                                                  caagcGATCAGAACGCCGCCGTCAAGGCCCTGGCCGAGGGAGAACCCACCACATGAGCAAACGCTTGAGCTGGGACATTCATACCCGCACGCAGATCATCAGTCTCGGCCCCGACCTGCTACTGACCTTGCTCTTGATCAGCTTCTTTACCTTCGTACGCATCCAGGACCTGCGCCAGGAACTCAACCACACCGGGCAACTGATCGCCAACCAGTTGGCACCGGCGTCGGAATATGGGGTGATCTCGGGCAACAACGAAGTGCTTGAAGGGCTGATGCGTGCCACCCTCAGCATCCCGCAc                                                                                                         <  1:264898/297‑1 (MQ=255)                                                                                                                                                                       gATCAGAACGCCGCCGTCAAGGCCCTGGCCGAGGGAGAACCCACCCCATGAGCAAACGCTTGAGCTGGGACATTCATACCCGCACGCAGATCATCAGTCTCGGCCCCGACCTGCTACTGACCTTGCTCTTGATCAGCTTCTTTACCTTCGTACGCATCCAGGACCTGCGCCAGGAACTCAACCACACCGGGCAACTGATCGCCAACCAGTTGGCACCGGCGTCGGAATATGGGGTGATCTCGGGCAACAACGAAGTGCTTGAAGGGCTGATGCGTGCCACCCTCAGCATCCCGCACGTgcg                                                                                                    <  1:225449/301‑1 (MQ=255)                                                                                                                                                                                  aagcagctcagCCCTGGCCGAGGGAGAACCCACCCCATGAGCAAACGCTTGAGCTGGGACATTCATACCCGCACGCAGATCATCAGTCCCGGCCCCGACCTGCTACTGACCTTGCTCTTGATCAGCTTCTTTACCTTCGTACGCATCCAGGACCTGCGCCAGGAACCCAACCACACCGGTCAACTTATCGCTAACCAGTTGGCACCGGCGACGGCATCTGGGGTTAACTCGGTCAACAACCAACGGCTTGAAGGGCTGATGCGTGCAAGCATCAACATCCCCCACgcacatgtctgtgaag                                                                                         >  2:385934/11‑286 (MQ=255)                                                                                                                                                                                   gccctgcaagCCCTGGCCCAAGGGGAAACCACCCCATGAGCAAACGCTTGAGCTTGGAAATTCATACCCGCACGCAGATCATCAATCTCGGCCCCGACCTGCTACTGGCCTTGCTCTTGATCAGCTTCTTTACCTTCGTACGCATCCAGGACCTGCGCCAGGAACTCAACCACACCGGGCAACTGATCGCCAACCAGTTGGCACCGGCGTCGGAATATGGGGTGATCTCGGGCAACAACGAAGTGCTTGAAGGGCTGATGCGTGCCACCCTCAGCATCCCGCACGTGCGCTTTCTGGAAGt                                                                                        <  2:545761/292‑1 (MQ=255)                                                                                                                                                                                         aaGGCCCTGGCCGAGGGAGAACCCACCCCATGAGCAAACGCTTGAGCTGGGACATTCATACCCTCACGCAGATCATCAGTCTCGGCCCCGACCTGCTACTGACCCTGCTCTTGATCAGCTTCTTTACCTTCGTACGCATCCAGGACCTGCGCCAGGAACTCAACCACACCGGGCAACTGATCGCCAACCAGTTGGCACCGGCGTCGGAATATAGGGTGATCTCGGGCAACAACGAAGTGCTTGAAGGGCTGATGCGTGCCACCCTCAGCATCCCGCACGTGCGCTTTCTGGAAGTGCAGGa                                                                                  <  1:635744/301‑1 (MQ=255)                                                                                                                                                                                                      aGGGAGAACCCACCACATGAGCAAACGCTTGAGCTGGGACATTCATACCCGCACGCAGATCATCAGTCTCGGCCCCGACCTGCTACTGACCTTGCTCTTGATCAGCTTCTTTACCTTCGTACGCATCCAGGACCTGCGCCAGGAACTCAACCACACCGGGCAACTGATCGCCAACCAGTTGGCACCGGCGTCGGAATATGGGGTGATCTCGGGCAACAACGAAGTGCTTGAAGGGCTGATGCGTGCCACCCTCAGCATCCCGCACGTGCGCTTTCTGGAAGTGCAGGACAGCCGCAACcac                                                                     <  1:678456/301‑1 (MQ=255)                                                                                                                                                                                                          agaACCCACCCCATGAGCAAACGCTTGAGCTGGGACATTCATACCCGCACGCAGATCATCAGTCTCGGCCCCGACCTGCTACTGACCTTGCTCTTGATCAGCTTCTTTACCTTCGTACGCATCCAGGACCTGCGCCAGGAACTCAACCACACCGGGCAACTGATCGCCAACCCGTTGGCACCGCCGTCGGAATATGGGGTGATCTCGGGCAACAACGACGTGCTTGAAGGGCTGATGCGTGCCACCCTCAGCATCACGCACGTGCGCTCTCTGGAAGTGCAGGACAGCCGCAACCACATcc                                                                 >  1:446258/1‑301 (MQ=255)                                                                                                                                                                                                               ccacccCATGAGCAAACGCTTGAGCTGGGACATTCATACCCGCACGCAGATCATCAGTCTCGGCCCCGACCTGCTACTGACCTTGCTCTTGATCAGCTTCTTTACCTTCGTACGCATCCAGGACCTGCGCCAGGAACTCAACCACACCGGGCAACTGATCGCCAACCAGTTGGCACCTGCGTCGGAATATGGTGTGATCTCGGGTAACAACGACGTCCTTGAAGGGCTGATGCGTGCTACCGTCAGCATCCCGCACGTTCACTTTCTGCAAGTGCAAGACatgtgcaccaacctcctggtg                                                            >  2:201927/1‑281 (MQ=255)                                                                                                                                                                                                                     cATGAGCAAACGCTTTGGCTGTGACATTCATACCCGCACGCAGATCATCAGTCTCGGCCCCGACCTGCTACTGACCTTGCTCTTGATCAGCTTCTTTACCTTCGTACGCATCCAGGACCTGCGCCAGGAACTCAACCACACCGGGCAACTGATCGCCAACCAGTTGGCACCGGCGTCGGAATATGGGGTGATCTCGGGCAACAACGAAGTGCTTGAAGGGCTGATGCGTGCCACCCTCAGCATCCCGCACGTGCGCTTTCTGGAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTc                                                      <  2:217394/301‑1 (MQ=255)                                                                                                                                                                                                                      aTGAGCAAACGCTTGAGCTGGGACATTCATACCCGCACGCAGATCATCAGTCTCGGCCCCGACCTGCTACTGACCTTGCTCTTGATCAGCTTCTTTACCTTCGTACGCATCCAGGACCTGCGCCAGGAACTCAACCACACCGGGCAACTGATCGCCAACCAGTTGGCACCGGCGTCGGAATATGGGGTGATCTCGGGCAACAACGAAGTGCTTGAAGGGCTGATGCGTGCCACCCTCAGCATCCCGCACGTGCGCTTTCTGGAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCg                                                     <  1:210724/301‑1 (MQ=255)                                                                                                                                                                                                                       tGAGCAAACGCTTGAGCTGGGACATTCATACCCGCACGCAGATCATCAGTCTCGGCCCCGACCTGCTACTGACCTTGCTCTTGATCAGCTTCTTTACCTTCGTACGCATCCAGGACCTGCGCCAGGAACTCAACCACACCGGGCAACTGATCGCCAACCAGTTGGCACCGGCGTCGGAATATGGGGTGATCTCGGGCAACAACGAAGTGCTTGAAGGGCTGATGCGTGCCACCCTCAGCATCCCGCACGTGCGCTTTCTGGAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGa                                                    <  1:416266/301‑1 (MQ=255)                                                                                                                                                                                                                                 cTTGAGCTGGGACATTCATACCCGCACGCAGATCATCAGTCTCGGCCCCGACCTGCTACTGACCTTGCTCTTGATCAGCTTCTTTACCTTCGTACGCATCCAGGACCTGCGCCAGGAACTCAACCACACCGGGCAACTGATCGCCAACCAATTGGCACCGGCGTCGGCATATGGGGTGATCTCGGGCAACAACGAAGTGCTTGAAGGGCTGATGCATGCCACCCTCAGCGTCCCCCACGTGCTCTTTCTGGAAGTGCAGGCCAGCCGCAACCACATCCTGTTGTACGTCGAGCAGTCcgac                                          >  1:83807/1‑301 (MQ=255)                                                                                                                                                                                                                                        tGGGACATTCATACCCGCACGCAGATCATCAGTCTCGGCCCCGACCTGCTACTGACCTTGCTCTTGATCAGCTTCTTTACCTTCGTACGCATCCAGGACCTGCGCCAGGAACTCAACCACACCGGGCAACTGATCGCCAACCAGTTGGCACCGGCGTCGGAATATGGGGTGATCTCGGGCAACAACGAAGTGCTTGAAGGGCTGATGCGTGCCACCCTCAGCATCCCGCACGTACGCTTTCTGAAAGTGCAGGACAGCCGCAACCACATCCGGGTTTACCTCGAGCCGTCCGTCGACAGTg                                   >  2:631548/1‑301 (MQ=255)                                                                                                                                                                                                                                             cATTCATACCCGCACGCAGATCATCAGTCTCGGCCCCGACCTGCTACTGACCTTGCTCTTGATCAGCTTCTTTACCTTCGTACGCATCCAGGACCTGCGCCAGGAACTCAACCACACCGGGCAACTGATCGCCAACCAGTTGGCACCGGCGTCGGAATATGGGGTGATCTCGGGCAACAACGAAGTGCTTGAAGGGCTGATGCGTGCCACCCTCAGCATCCCGCACGTGCGCTTTCTGGAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGACAGTGCCAAc                              >  1:64814/1‑301 (MQ=255)                                                                                                                                                                                                                                              aTTCATACCCGCACGCAGATCATCAGTCTCGGCCCCGACCTGCTACTGACCTTGCTCTTGATCAGCTTCTTTACCTTCGTACGCATCCAGGACCTGCGCCAGGAACTCAACCACACCGGGCAACTGATCGCCAACCAGTTGGCACCGGCGTCGGAATATGGGGTGATCTCGGGCAACAACGAAGTGCTTGAAGGGCTGATGCGTGCCACCCTCAGCATCCCGCACGTGCGCTTTCTGGAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAAcc                             >  1:151883/1‑301 (MQ=255)                                                                                                                                                                                                                                              aTTCATACCCGCACGCAGATCATCAGTCTCGGCCCCGACCTGCTACTGACCTTGCTCTTGATCAGCTTCTTTACCTTCGTACGCATCCAGGACCTGCGCCAGGAACTCAACCACACCGGGCAACTGATCGCCAACCAGTTGGCACCGGCGTCGGAATATGGGGTGATCTCGGGCAACAACGAAGTGCTTGAAGGGCTGATGCGTGCCACCCTCAGCATCCCGCACGTGCGCTTTCTGGAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAAcc                             >  1:293808/1‑301 (MQ=255)                                                                                                                                                                                                                                                 cATACCCGCACGCAGATCATCAGTCTCGGCCCCGACCTGCTACTGACCTTGCTCTTGATCAGCTTCTTTACCTTCGTACGCATCCAGGACCTGCGCCAGGAACTCAACCACACCGGGCAACTGATCGCCAACCAGTTGGCACCGGCGTCGGAATATGGGGTGATCTCGGGCAACAACGAAGTGCTTGAAGGGCTGATGCTTGCCACCCTCAGCATCCCGCACGTGCGCTTTCTGGAAGTGCAGGACATGCGCAACCACCTCCTGGTTTCCGTCGatccgtccagcggagtgcccaaccccc                          >  2:621479/1‑275 (MQ=255)                                                                                                                                                                                                                                                 cATACCCGCACGCAGATCATCAGTCTCGGCCCCGACCTGCTACTGACCTTGCTCTTGATCAGCTTCTTTACCTTCGTACGCATCCAGGACCTGCGCCAGGAACTCAACCACACCGGGCAACTGATCGCCAACCAGTTGGCACCGGCGTCGGAATATGGGGTGATCTCGGGCAACAACGAAGTGCTTGAAGGGCTGATGCGTGCCACCCTCAGCATCCCGCACGTGCGCTTTCTGGAAGTGCAGGACAGCCGCAACCCCATCCTGGTGTACGTCGAGCAGTCAGACGAAAGTGCCAACcgcg                          >  2:118821/1‑301 (MQ=255)                                                                                                                                                                                                                                                 cATACCCGCACGCAGATCATCAGTCTCGGCCCCGACCTGCTACTGACCTTGCTCTTGATCAGCTTCTTTACCTTCGTACGCATCCAGGACCTGCGCCAGGAACTCAACCACACCGGGCAACTGATCGCCAACCAGTTGGCACCGGCGTCGGAATATGGGGTGATCTCGGGCAACAACGAAGTGCTTGAAGGGCTGATGCGTGCCACCCTCAGCATCCCGCACGTGCGCTTTCTGGAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGATCAGTCCGACGAAAGTGCCAACcgcc                          >  1:430224/1‑300 (MQ=255)                                                                                                                                                                                                                                                      ccGCACGCAGATCATCAGTCTCGGCCCCGACCTGCTACTGACCTTGCTCTTGATCAGCTTCTTTACCTTCGTACGCATCCAGGACCTGCGCCAGGAACTCAACCACACCGGGCAACTGACCGCCAACCAGCTGGCACCGGCGTCGGCATATGGGGTGATCTCGGGCAACACCGAAGTGCTCGAAGGCCTGATGCGTGCCACCCTCAGCATCCCGCACGTGCGCTTTCTGGAAGTGTAGGACAGCCGCAACCACCTCCTGGCGTACGCCGAGCAGTCCGACGAAAGTGCCAAcgccgcgaaa                     >  2:73689/1‑292 (MQ=255)                                                                                                                                                                                                                                                        gcacgcaGATCATCAGTCTCGGCCCCGACCTGCTACTGACCTTGCTCTTGATCAGCTTCTTTACCTTCGTACGCATCCAGGACCTGCGCCAGGAACTCAACCACACCGGGCAACTGATCGCCAACCAGTTGGCACCGGCGTCGGAATATGGGGTGATCTCGGGCAACAACGAAGTGCTTGAAGGGCTGATGCGTGCCACCCTCAGCATCCCGCACGTGCGCTTTCTGGAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACa                   <  1:211995/301‑1 (MQ=255)                                                                                                                                                                                                                                                                 tcatcaGTCTCGGCCCCGACCTGCTACTGACCTTGCTCTTGATCAGCTTCTTTACCTTCGTACGCATCCAGGACCTGCGCCAGGAACTCAACCACACCGGGCAACTGATCGCCAACCAGTTGGCACCGGCGTCAGAATATGGGGTGATCTCGGGCAACAACGAAGTGCTTGAAGGGCTGATGCGTGCCACCCTCAGCATCCCGCACGTGCGCTTTCTGGAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGgt          <  1:243743/301‑1 (MQ=255)                                                                                                                                                                                                                                                                   atcaGTCTCGGCCCCGACCTGCTACTGACCTTGCTCTTGATCAGCTTCTTTACCTTCGTACGCATCCAGGACCTGCGCCAGGAACTCAACCACACCGGGCAACTGATCGCCAACCAGTTGGCACCGGCGTCGGAATCTGGGGTGATCTCGGGCAACAACGAAGTGCTTGAAGGGCTGATGCGTGCCACCCTCAGCATCCCGCACGTGCGCTTTCTGGAAGTGCCGGACAGCCGCAACCACATCCTGGTGTACGTCGCGCAGTCAGACGCAAGTGCCAACcgcgcgacacagccagaggtgt        >  2:385614/1‑285 (MQ=255)                                                                                                                                                                                                                                                                     cagtgccGGCCCCGACCGCCTACTGACCTGGCTCTTGGTCAACTTCTTTACCTTCGAACGCATCCAGGACCTGCGCCAGGAACGCAACCACGCCGGGCAACTGATCGCCGACCAGTTGGCACCGGCGTCGGAATATGGGGTGATCTCGGACAACAACGAAGTGCTTGAAGGGCTGATGCGTGCCACCCTCAGCGTCCCGCACGTGCGCTTTCTGGAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGGCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTTc      <  1:369695/295‑1 (MQ=255)                                                                                                                                                                                                                                                                     caGTCTCGGCCCCGACCTGCTACTGACCTTGCTCTTGATCAGCTTCTTTACCTTCGTACGCATCCAGGACCTGCGCCAGGAACTCAACCACACCGGGCAACTGATCGCCAACCAGTTGGCACCGGCGCCGGAATATGGGGTGATCTCGGGCAACAACGACGTGCTTGAAGGCATGATGCGTGCCACCCTCAGCATCCCGCACGTGCGCTTTCTGGAAGTGCAGGACAGCCGCAACCACATCTCGCTGTACGTCGAGCAGTCCGACGACAGGGCCAACCGCGCGCACCAGGTCGCGGTGTTc      >  2:369547/1‑301 (MQ=255)                                                                                                                                                                                                                                                                       gtgtcGGACCCGACCTGCTACTGACCTTGCTCTTGATCAGCTTCTTTACCTTCGTACGCATCCAGGACCTGCGCCAGGAACTCAACCACACCGGGCAACTGATCGCCAACCAGTTGGCACCGGCGTCGGAATATGGGGTGATCTCGGGCAACAACGAAGTGCTTGAAGGGCTGATGCGTGCCACCCTCAGCATCCCGCACGTGCGCTTTCTGGAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTTCCa    <  1:690676/298‑1 (MQ=255)                                                                                                                                                                                                                                                                         ctcGGCCCCGACCTGCTACTGACCTTGCTCTTGATCAGCTTCTTTACCTTCGTACGCATCCAGGACCTGCGCCAGGAACTCAACCACACCGGGCAACTGATCGCCAACCAGTTGGCACCGGCGTCGGAATATGGGGTGATCTCGGGCAACAACGAAGTGCTTGAAGGGCTGATGCGTGCCACCCTCAGCATCCCGCACGTGCGCTTTCTGGAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTTCCAgg  >  2:490466/1‑301 (MQ=255)                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                                                   CGCCGGCGCCCTGGCCAGTACCTACAGCGACCAGGACGACTGGCTGTGGGTGCTCGACCAGTTGCTCGACCAGCCGCCGACGCGCTGGCCGCGCAGCCTGTACCCGGTGCGTTACGGCGTCAGTGGCAACCAGCAGGTGGCCCGTGCCCTGGGTCTTGAACCAATCGATCAGAACGCCGCCGTCAAGGCCCTGGCCGAGGGAGAACCCACCCCATGAGCAAACGCTTGAGCTGGGACATTCATACCCGCACGCAGATCATCAGTCTCGGCCCCGCCCTGCTACTGACCTTGCTCTTGATCAGCTTCTTTACCTTCGTACGCATCCAGGACCTGCGCCAGGAACTCAACCACACCGGGCAACTGATCGCCAACCAGTTGGCACCGGCGTCGGAATATGGGGTGATCTCGGGCAACAACGAAGTGCTTGAAGGGCTGATGCGTGCCACCCTCAGCATCCCGCACGTGCGCTTTCTGGAAGTGCAGGACAGCCGCAACCACATCCTGGTGTACGTCGAGCAGTCCGACGAAAGTGCCAACCGCGCGCAACAGGTCGAGGTGTTCCAGG  >  AE015451/4714557‑4715121

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