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%	189.4 / NA	58	V25I (GTC→ATC)	dcuS	sensor histidine kinase DcuS
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/37);  total (21/37)

AGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGCGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGACCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGG‑AGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAACTGCCA‑GTGCTGTATGAAGATTCCGGGGCTATGCTTATAGCGATAATCATACTGATGAGAGAGGG  >  NC_000913/4351346‑4351836                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                         agtctgtgcaatccggaaaagACAATAAACATCAGATCGTTGCGGTGGCGTTAGGCTTCCGCGATGGCCTGGATGCCACCTTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGCCTGCGAGGGTACGCGCCACGGCCAGTGCCTTGTTGGCTAGCCCAGCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCAtgaa                                                                                                                                                                                                                                                    <  2:133319/231‑4 (MQ=255)             tCGGTAACGACAAGAAACAGCAGGGCGTTGCGTGTGCTTACGTCTTCCGCGAGGGCCTGGATGCCATTCTCCTGCGGTTTTTTCTGCAAGCCCTGGCGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAAGAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTgg                                                                                                                                                                                                                                         <  2:498526/250‑1 (MQ=255)                     gACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCCCCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAAtt                                                                                                                                                                                                                                >  2:183139/1‑251 (MQ=255)                      aCAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAAttt                                                                                                                                                                                                                               <  1:286972/251‑1 (MQ=255)                             aCAGCAGAGCGGTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGAGGCCGCTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGa                                                                                                                                                                                                                        <  2:8126/251‑1 (MQ=255)                                  agATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGttt                                                                                                                                                                                                                   <  1:183139/251‑1 (MQ=255)                                            gtttgcggACGGCTTCCGCGATGGCCTGGTTGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGCCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCTCCACCAAGAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTTTGGTACTCAATTTAATCGGACGTTTGCGTAACATg                                                                                                                                                                                                         <  1:292725/251‑1 (MQ=255)                                                gcgtACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGt                                                                                                                                                                                                     >  1:170873/1‑251 (MQ=255)                                                  gtACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGt                                                                                                                                                                                                     >  1:150928/1‑249 (MQ=255)                                                    gcGGCTTCCGCGAGGGCTTGGATGCCACTCTCCGGCGGTTTTGTCTGCAAGCCCGTACTGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCTGACGTTTGCGTAACATGCGGTAggg                                                                                                                                                                                                 <  2:333737/250‑1 (MQ=255)                                                    aCGGCTTCCGCGATGGCCGGGATGCCCATCTGCTTCGGTTTTTGCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACGCAATGTCATCGGACGTTTGCGTAACATGCGGTAggg                                                                                                                                                                                                 <  2:107989/251‑1 (MQ=255)                                                       gCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCaa                                                                                                                                                                                              <  1:431721/251‑1 (MQ=255)                                                        cTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCaat                                                                                                                                                                                             >  1:324267/1‑251 (MQ=255)                                                          tCCGGGATGGCCGGGATGCCACTCTCCTGGGGTTTTTTCTGCAAGCCCTGCCGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAAGTTCATCGGACGTTTGCGTAACATGCGGTAGGGCaat                                                                                                                                                                                             <  2:149699/249‑1 (MQ=255)                                                                aTGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGtct                                                                                                                                                                                     >  1:474446/1‑251 (MQ=255)                                                                aTGGCCTGGATGCCACTCTCCTGCGGCTTTTTCTGCAAGCCCTGACGGATTTCCTGCGAGTCCGCGAGGGTATGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGGGTCATATCACTACTTTGCGCGAAGTAAATCAGCTGCACCCCCACTAGCACCGTGAACAGTACCGGCCCAATCATTATGATCAATGCGATACCCAATTTAATCGGAAGTATGCGTAACATGCGGTAgagaaaggcaggtct                                                                                                                                                                                     >  2:184333/1‑237 (MQ=255)                                                                 tGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGtctc                                                                                                                                                                                    >  1:314894/1‑251 (MQ=255)                                                                  ggCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCTAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGtctc                                                                                                                                                                                    <  1:108640/250‑1 (MQ=255)                                                                   gCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTGCGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCtcat                                                                                                                                                                                  <  1:468084/251‑1 (MQ=255)                                                                    ccTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATAGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCtcatc                                                                                                                                                                                 >  1:63387/1‑251 (MQ=255)                                                                                 ctcCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTtgt                                                                                                                                                                      <  1:354469/249‑1 (MQ=255)                                                                                       cGGTTTTTTCTGCAAGCCCTGCCGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACaaa                                                                                                                                                              <  1:310532/251‑1 (MQ=255)                                                                                             tttCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCt                                                                                                                                                         <  1:73035/250‑1 (MQ=255)                                                                                                         cTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTAtctc                                                                                                                                            >  1:237243/1‑251 (MQ=255)                                                                                                             gggttttCGGGCGAGTCGTGGAGGGTACTCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAAGCAGATGCTCCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCtgat                                                                                                                                        <  2:155981/247‑1 (MQ=255)                                                                                                               gATTTCCGTCGCGTCGGCGGGGGTACCCGCCACTGCCAGTGACTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCCAGAAGTAAATCAGATGCTCCACCCAGAGCACCGAGTACAGTACCGCGCTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGTCGTTTGCGTAACATGCTGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCAGTATCTCtgatga                                                                                                                                      <  2:334459/251‑1 (MQ=255)                                                                                                               gATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCtgatg                                                                                                                                       >  1:566807/1‑250 (MQ=255)                                                                                                                aTTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAg                                                                                                                                     <  1:44955/251‑1 (MQ=255)                                                                                                                 gttCCGGCGGGGCGGCGAGGGTCCGCGCCACTGCCAGTGCCTTGGTGGCTAGCCCAACTCGTGTCATATCACTGATTTGCTAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAg                                                                                                                                     <  2:448436/249‑1 (MQ=255)                                                                                                                            gcGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCa                                                                                                                         <  1:184333/250‑1 (MQ=255)                                                                                                                            gcGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCa                                                                                                                         <  1:91728/250‑1 (MQ=255)                                                                                                                                     gttcgcgCCCCTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCATAGGG‑AGTa                                                                                                               <  2:406902/248‑1 (MQ=255)                                                                                                                                          cgcCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGG‑AGTAAGa                                                                                                            <  1:472862/249‑1 (MQ=255)                                                                                                                                          cgcCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGAACTGATCATTAAGCTCACTGTGGTACTCAATTTCATCGGACGATTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATATCTGATGAGGCggggtattacacaggggtaaaga                                                                                                            >  2:472862/1‑230 (MQ=255)                                                                                                                                                 gCCAGTGCCTTGTTGGCAAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACGGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGTG‑AGTAAGAATGATTGGc                                                                                                   <  2:558835/251‑1 (MQ=255)                                                                                                                                                  ccAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGG‑AGTAAGAATGATTGGCt                                                                                                  <  1:65442/251‑1 (MQ=255)                                                                                                                                                     gTGCCTTGTTGGCTAGCCCATCTCGGGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTTCCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGG‑AGTAAGAATGATTGGCtata                                                                                               <  2:205983/251‑1 (MQ=255)                                                                                                                                                       gCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGG‑AGTAAGAATGATTGGCtatat                                                                                              <  1:116778/250‑1 (MQ=255)                                                                                                                                                              tgGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCaatgaa                                                                                                                                                                                          >  1:302328/1‑152 (MQ=255)                                                                                                                                                              tgGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCaatgaa                                                                                                                                                                                          <  2:302328/152‑1 (MQ=255)                                                                                                                                                               ggCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGG‑AGTAAGAATGATTGGCTATATAGGGGAaga                                                                                     <  1:313834/251‑1 (MQ=255)                                                                                                                                                                     cccatctcctttCATATCACTGATATGCGAGAAGTAAATCAGATGCACCACCAAAAGCACCGAGAACATTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATTGGACGTTTGCGTAACATGCGGTATGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCTGGTAATTCAAAGGG‑AGTAAGAATGATTGGCTATTTAGGGGAAGATActct                                                                               <  2:143850/240‑1 (MQ=255)                                                                                                                                                                            cgtgTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCTCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGG‑AGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCa                                                                           <  2:330616/245‑1 (MQ=255)                                                                                                                                                                                    atCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGG‑AGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCACCGGAAACTGc                                                                <  2:170873/251‑1 (MQ=255)                                                                                                                                                                                          gATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGG‑AGTATGAATGGTTTGCTTTATAGGGGAAAGGACTCTGGCAAAGGGAACTGCCA‑GTgg                                                         >  1:469312/1‑250 (MQ=255)                                                                                                                                                                                           aTTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTGAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTa                                                                                                                                                >  1:39859/1‑165 (MQ=255)                                                                                                                                                                                           aTTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTGAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTa                                                                                                                                                <  2:39859/165‑1 (MQ=255)                                                                                                                                                                                               gCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGG‑AGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAACTGCCA‑GTGCTGTa                                                     <  1:247123/250‑1 (MQ=255)                                                                                                                                                                                                                    ccaccaATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGGAATTTAAAAGG‑GGGAAGAATGATTTGCTTTATAGGGGGAGAGACTCTGGCAACGGGAAATGCCC‑GTGCTGTATTAAGATTCCCGGGCTATGCtt                               >  1:6308/1‑251 (MQ=255)                                                                                                                                                                                                                              cACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGG‑AGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAACTGCCA‑GTGCTGTATGAAGATTCCGGGGCTATGCTTATAGCGATaa                     <  1:488505/251‑1 (MQ=255)                                                                                                                                                                                                                                ccGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGATATTCAAAGGG‑GGGAAAAAAGATTGGCTAATTTGGGGGAGAGCATCAGGCAACAGAAAACGCCA‑GAGGTGTCAGTAGATTCCGGGGCTATGCTTATAGCGaatatc                   >  1:336759/1‑246 (MQ=255)                                                                                                                                                                                                                                    gaACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCATTGTGTGACAAATTTCTTAAGCATTATCTCTGCTGAGGCGGGGAAATTAAAGGG‑GGTAAAAAAGGTTTGCTATATAGGGGGAGAGACTCTTGCACCGAAAACTGCCA‑ATGCGGTATTAAGATTCAGGGGGTATGCTTATAGCGATAATTATAc               >  1:536964/1‑251 (MQ=255)                                                                                                                                                                                                                                     aaCAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGGAATTCAAAGGGGAG‑AAGAATGATTGGATAAATAGGGGAAGAGGATCTGGCAAAGGGAAATGGCA‑GTGGTGAATGAAGATTCCGAGGGTATGATTATAgaggtaataataat              >  1:206303/1‑238 (MQ=255)                                                                                                                                                                                                                                      aCAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGG‑AGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAACTGCCA‑GTGCTGTATGAAGATTCCGGGGCTATGCTTATAGCGATAATCATACtg             <  1:50198/251‑1 (MQ=255)                                                                                                                                                                                                                                            ccGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCt                                                                                                                                                         >  2:192132/1‑107 (MQ=255)                                                                                                                                                                                                                                            ccGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCt                                                                                                                                                         <  1:192132/107‑1 (MQ=255)                                                                                                                                                                                                                                                aCTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACACATTTCTTAAGCATTATCTCTGATGAGGCGGGTCAATCATAGGG‑AGTGAGATAGATTTGCTATATAGGGGTAGCAACACTGGCTACGGCCACT‑CCACGTGCTGTATAAAATTTCAGGGGCTATGCTTATATCGATAATCCTAATAGTGAGAGAgg   >  1:456669/1‑251 (MQ=255)                                                                                                                                                                                                                                                 cTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTAATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGGAATTCCAAGGG‑AGTAAAAATAATTGGGTATTAAGGGGAAGAGACATTGGCAACGGAAAATGGCA‑GTGGTGTATGAAGATTCCGGGGCGCAGGTTATAGCACTAATCATACtggttagggaggg  >  1:273569/1‑240 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                         AGACTTTGCATATCGGTAACGACAATAAACAGCAGATCGTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGCGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCACCAATAGCACCGAGAACAGTACCGCACTGACCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGG‑AGTAAGAATGATTGGCTATATAGGGGAAGAGACTCTGGCAACGGAAACTGCCA‑GTGCTGTATGAAGATTCCGGGGCTATGCTTATAGCGATAATCATACTGATGAGAGAGGG  >  NC_000913/4351346‑4351836

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