Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,351,590	C→T	V25I (GTC→ATC)	dcuS ←	sensor histidine kinase DcuS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,351,590	0	C	T	100.0%	206.1 / NA	65	V25I (GTC→ATC)	dcuS	sensor histidine kinase DcuS
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (25/40);  total (25/40)

CGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGCGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGACCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATAT‑AGGGGAAGAGACTCTGGC‑AACGGAAACTGCCAGTGCTGTATGAAGATTCCGGGGCTATGCTTATAGCGATAATCATACTGATGAGAGAG  >  NC_000913/4351343‑4351834                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       cGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCAt                                                                                                                                                                                                                                                     >  2:139147/1‑251 (MQ=255)  gAAGACTGTGCATATCGGTAACGAAAATAAACAGCAGAGCGTTGCGTTGGCGTACGGGTTCCGCGATGGTCGGGATGTCACTCTCCTGTGGTTTTTTCTGCAAGCCCTGACGGAGTGCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGAGTCATATCACTGAGTGGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCAtt                                                                                                                                                                                                                                                    <  2:398306/251‑1 (MQ=255)   aaGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTa                                                                                                                                                                                                                                                   <  1:520846/251‑1 (MQ=255)      gcTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGa                                                                                                                                                                                                                                                <  2:470754/250‑1 (MQ=255)      aCTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGa                                                                                                                                                                                                                                                <  1:485712/251‑1 (MQ=255)             atatCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACtgt                                                                                                                                                                                                                                         >  1:386481/1‑251 (MQ=255)              tatCGGTAACGCCAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCAGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGCCGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCCTATCACTGCTTTGCGTGAAGTAAATCAGCTGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTATGTTCACtgtg                                                                                                                                                                                                                                        >  1:288353/1‑251 (MQ=255)               atCGGTAACGAGAATAAACAGCAGATCGTTGCGCGGGCGTACGGCTTCCGCGATGGCCGGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGTCGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCGCGTGTCATAGCACTGATTTGCGAGAAGTAAATCTGATGCACCACCAATAGCACCGAGAACAGTGCCGCACTGATCATTAAGATCACTGTgg                                                                                                                                                                                                                                       <  2:296458/251‑1 (MQ=255)               atCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTgg                                                                                                                                                                                                                                       >  1:533034/1‑251 (MQ=255)               atCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCAATGTgg                                                                                                                                                                                                                                       >  1:558549/1‑251 (MQ=255)                 cGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGt                                                                                                                                                                                                                                      >  1:516036/1‑250 (MQ=255)                        gACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAAtt                                                                                                                                                                                                                              <  1:374649/251‑1 (MQ=255)                          cAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCATGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGAGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGCTCATTAAGATCACTGTGGTACTCAATTTc                                                                                                                                                                                                                            >  2:441693/1‑251 (MQ=255)                             tAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATTGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGCTCATTAAGATCACTGTGGTACTCAATTTCATc                                                                                                                                                                                                                         >  2:18175/1‑251 (MQ=255)                                   gcagATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGCACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGt                                                                                                                                                                                                                   >  1:85407/1‑251 (MQ=255)                                    cagATGGTTGCGTGTGCGTTCGGCTTCCGCGATGGCCGGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAATTAAATCAGATGCACCACCAATAGCACCGAGAACAGGACCGCACTGATCATTAAGATCACTGTGGTACTCAAGTGCATCGGACGtt                                                                                                                                                                                                                  <  2:305882/251‑1 (MQ=255)                                    cagATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCTCCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGtt                                                                                                                                                                                                                  <  1:61158/251‑1 (MQ=255)                                    cagATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGtt                                                                                                                                                                                                                  <  2:494800/251‑1 (MQ=255)                                                           cttacgcgATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTTAATCAGATGAACCACCAATAGCACCGAGCACAGTACCGCACTGATCATTAAGATCACTGAGGTACTCAATTTTATCGGACGTTTGAGTAACATGCGGTAGGGCaat                                                                                                                                                                                           >  1:233181/5‑251 (MQ=255)                                                           cTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCaat                                                                                                                                                                                           >  1:55388/1‑251 (MQ=255)                                                                          ggATGCCCCTCTCCTGCGGTTTTTTCTGCAAGCTCTGACTGATTTCCGGCGAGTCGGAGAGGGTACGCGCCCCTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCtcatca                                                                                                                                                                              <  2:436532/249‑1 (MQ=255)                                                                                      ccTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCTGCGAGTCGGCGAGGGTACACGCCACTGCCAGTACATGGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCGCCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCAGCGGACGTTTGCGTAACATGCGGTAGGTCAATGAATGTCTCATCAGCTTCCTtgtgt                                                                                                                                                                  <  2:490138/249‑1 (MQ=255)                                                                                      ccTGCGGTTTTTTCGGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTtgtgt                                                                                                                                                                  <  2:255549/249‑1 (MQ=255)                                                                                      ccTGCGGGTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGGGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAAGCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATTATTAAGATCACTGGGGTACTCAAGTTCATCGGACGTTTGCGTAACAGGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGa                                                                                                                                                                <  2:288353/251‑1 (MQ=255)                                                                                         gCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACGGCACTGATCATTAAGATCACTGTGGTACTGAATTTCATCGGGCGTTTGCGTAACATGCGGTAGGGCAATGGATTTCTCATCAGCTGCCTGGTGTGACaa                                                                                                                                                             <  2:573693/251‑1 (MQ=255)                                                                                         gCGGTTTTTTCTGCAAGCCCTGACGGATTGCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACaa                                                                                                                                                             <  1:16470/251‑1 (MQ=255)                                                                                              tttttCTGCAAGCCCTTACGGCTTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGATAGCCTATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTGGCTTAACATGCAGTAGGGCAATGAATGTCTCAGCTGATTCCTTGTGTGACAAAttg                                                                                                                                                         <  2:129739/250‑2 (MQ=255)                                                                                                       aaGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCAtt                                                                                                                                               >  1:504541/1‑251 (MQ=255)                                                                                                               aCGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCACATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTAtga                                                                                                                                       <  2:185501/251‑1 (MQ=255)                                                                                                                 ggATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTGGTTGGCTAGCCCATCGCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATGTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGGGTGACAAATTTCTTAGGCATTATCTCtgatg                                                                                                                                     <  1:481376/251‑1 (MQ=255)                                                                                                                  gATTTCCGGCGAGTCGGCCAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCACATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCATATGCACCACCAATATCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACTTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCtgatga                                                                                                                                    <  1:333015/251‑1 (MQ=255)                                                                                                                         ggCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGt                                                                                                                             <  1:294077/251‑1 (MQ=255)                                                                                                                         ggCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCAGGAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCTGTAGGGCAATGAATGTCTCAGCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGt                                                                                                                             <  1:579067/251‑1 (MQ=255)                                                                                                                                    gAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGCGAAGTAAATCAGATGCACCACCAATAGCACCGCGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTATTTccaatgg                                                                                                                  >  2:518966/1‑245 (MQ=255)                                                                                                                                     aGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGGAATTTAAAgggg                                                                                                                 >  1:524134/1‑250 (MQ=255)                                                                                                                                             cgcCACTGCCAGTGCCTTGTTGTCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTTAATTTCATCGGACGTTTGCGTAACAGGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATTAGGCGGGTAAGTTAAAGGGAGTTAGAAt                                                                                                         <  2:303527/251‑1 (MQ=255)                                                                                                                                             cgcCACTGCCAGTGCCTTGGTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAAt                                                                                                         <  1:180567/251‑1 (MQ=255)                                                                                                                                               ccACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCCTTAAGATCACTGTGGTACTCAATTTCATCGGGCGTTTGCGTATCATGAGGGAGTGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGCGGCGGGTGATTCATATGGAGTaaaacac                                                                                                        >  2:517799/1‑245 (MQ=255)                                                                                                                                                        gTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCtata                                                                                              <  1:224563/251‑1 (MQ=255)                                                                                                                                                               gtaggCTAGCCCCTATCGTGGCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCGTGTGTGACAAATTGCTTAAGCATTAGCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATAT‑Agggta                                                                                      <  2:134924/248‑3 (MQ=255)                                                                                                                                                                 tgGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCtg                                                                                                                                        >  1:114850/1‑200 (MQ=255)                                                                                                                                                                 tgGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCtg                                                                                                                                        <  2:114850/200‑1 (MQ=255)                                                                                                                                                                      agcccagtccgtatCATATCACTGATTTGCGAGAAGTAAATCAGATGCAACACCAAGAGCACTGCTCACAGTACCGCACTTATCAATAAGATCACTGTGGGTCTCAATTTCTACGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGATTCCTTGTGTGACAAATTTCTGAAGCATTATCTCTGATGAGTCGTGTAATTCAAAGGGCGTAAGAAGGATTGTCtatat‑atgtgaatgggct                                                                               <  2:152674/238‑13 (MQ=255)                                                                                                                                                                       gCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCAACGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGCCGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGAGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGACAAATCAAAGGGAGTAAGCAAGTTTGGATATAC‑AGGGGGAGAGActt                                                                              >  1:224607/1‑250 (MQ=255)                                                                                                                                                                           aTCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGGAATTCTAAGGGAGTAAGAATGATTGGGTATAT‑AGGGGAAGAGACTCTGGc                                                                          >  2:165320/1‑251 (MQ=255)                                                                                                                                                                               cgtgTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATAT‑AGGGGAAGAGACTCTGGC‑AAc                                                                      <  1:205147/247‑1 (MQ=255)                                                                                                                                                                               cgtgTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATGATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATAT‑AGGGGAAGAGACTCTGGc                                                                          <  1:577226/244‑1 (MQ=255)                                                                                                                                                                                 tgTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCAATGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAAGGAGGAAGAAAGGATGGCTATAT‑AGGGGGAGAGACTCTGGC‑AACgg                                                                    >  1:237220/2‑250 (MQ=255)                                                                                                                                                                                  gTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACTGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATAT‑AGGGGAAGAGACTCTGGC‑AACGGaa                                                                  <  2:449972/251‑1 (MQ=255)                                                                                                                                                                                            tGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGGAATTCAAAGGGGGTAAGAAAGATTTGCTATAT‑AGGGGGAGAGACTCTGGC‑AAAGGAAAATGcccggg                                                        >  2:111760/1‑247 (MQ=255)                                                                                                                                                                                                 tGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATAT‑AGGGGAAGAGACTCTGGC‑AACGGAAACTGCCAGTGCTGt                                                    <  2:2979/250‑1 (MQ=255)                                                                                                                                                                                                   cGAGAAGTAAATCATATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTGCATCGGATGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCGTCCTTGTGTGACAAATTTCTTAAGCATTATCGCTGATGAGGCTGGTAATTCAAAGGGAGTAAGAATGATTGGCGATAT‑AGGGGAAGAGACTCTGGC‑AACGGAAACTGCCAGTGCTGTATg                                                 <  1:270868/251‑1 (MQ=255)                                                                                                                                                                                                       aaGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAAATCAAAGGGGGGAAGAATCATTTGCTATAT‑AGGGGGAGAGACCCTGGC‑CAAAGCAACTGCCAGGGCTTTATGa                                                >  1:543172/1‑248 (MQ=255)                                                                                                                                                                                                       aaGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGGAATTTAAAGGGGGGAAGAATGATTGGCTATAT‑AGGGGGAGGGAACCTTGC‑AACCGAAACTGCCTGTGCTGTTTGa                                                >  2:465885/1‑248 (MQ=255)                                                                                                                                                                                                           aaaTCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATAT‑AGGGGAAGAGACTCTGGC‑AACGGAAACTGCCAGTGCTGTATGAAGATTcc                                         <  2:13142/251‑1 (MQ=255)                                                                                                                                                                                                                    gCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAAGGATTGGCTATAT‑AGGGGAAGAGACTCTGGC‑AACGGAAACTGCCAGTGCTGTATGAAGATTCCGGGGCTATg                                <  2:55388/251‑1 (MQ=255)                                                                                                                                                                                                                              tAGCACCGAGCACAGTACCGCACTGATCATTAAGCTCACTGTGGTACTCACTTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGGAATTCCAACGGGGTTAGAATGATTGGCTATATCAGGGGA‑GAGACTCTGGC‑GACGGCAACTGCCACTGCTGTATTAAGATTCCGGGGCTagacctataacg                       >  2:247089/1‑239 (MQ=255)                                                                                                                                                                                                                               aGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATAT‑AGGGGAAGAGACTCTGGC‑AACGGAAACTGCCAGTGCTGTATGAAGATTCCGGGGCTATGCTTATAGCGAt                     <  1:106531/251‑1 (MQ=255)                                                                                                                                                                                                                                   ccGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGGACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATATCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTTAAAAGGAGGAAGAATGGTTTGATATAT‑AGGGGGAGAGAATCTGGCCAACGGAAA‑TGCCCGAGCTTTATGAAGATTCCGGGGCTATGCTTATAACGGTAATc                 >  1:47980/1‑251 (MQ=255)                                                                                                                                                                                                                                        aaCAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATAT‑AGGGGAAGAGACTCTGGC‑AACGGAAACTGCCAGTGCTGTATGAAGATTCCGGGGCTATGCTTATAGCGATAATCATACt            <  1:18175/251‑1 (MQ=255)                                                                                                                                                                                                                                              acctcACTGATCAGTAATATCGCTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTTACAAATGTCTTAAGCAATATCTCTCATGAGTCGGGTAATTCAAAGGGAGTAAGATTGATTGGCTATCT‑AGGGGAACAGACTCTGGC‑AACGTAAACTGCCAGTGCTGTATGAAGATTCCGGGGCTATGCTTATTGCGATAATCATACtgatga       <  1:46187/246‑1 (MQ=255)                                                                                                                                                                                                                                              aCCGCACTGATCATTAAGATCACTGTGGTACGCAATTTCAGCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTAGAT‑AGGGGAAGAGACGCTGGC‑AACGTAAACTGCCAGTGCTGGATGAAGATTCCGGGGCTATGCTTATAGCAATAATCATACtgatga       <  2:233181/250‑1 (MQ=255)                                                                                                                                                                                                                                                cGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTTAAAGGGGGTAAGAATGGTTGGCTATAT‑AGGGGAAGAGACTCTGGC‑AACGGAAACTGCCAGTTCTGTATGAAGATTTCGGGGCTATGCTTATAGCGCTAATAATACtggtgat      >  2:582284/1‑248 (MQ=255)                                                                                                                                                                                                                                                 gCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATAT‑AGGGGAAGAGACTCTGGC‑AACGGAAACTGCCAGTGCTGTATGAAGAGTCCGGGGCTATGCTTATAGCGATAATCATACTGACgagaga   <  2:533034/251‑1 (MQ=255)                                                                                                                                                                                                                                                  cACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATAT‑AGGGGAAGAGACTCTGGC‑AACGGAAACTGCCAGTGCTGTATGAAGATTCCGGGGCTATGCTTATAGCGATAATCATACTGATGagagag  <  1:312600/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       CGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGCGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGACCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATAT‑AGGGGAAGAGACTCTGGC‑AACGGAAACTGCCAGTGCTGTATGAAGATTCCGGGGCTATGCTTATAGCGATAATCATACTGATGAGAGAG  >  NC_000913/4351343‑4351834

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