Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	784,677	G→A	A49A (GCC→GCT)	zitB ←	Zn(2(+))/Cd(2(+))/Ni(2(+))/Cu(2(+)) exporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	784,677	0	G	A	100.0%	255.3 / NA	73	A49A (GCC→GCT)	zitB	Zn(2(+))/Cd(2(+))/Ni(2(+))/Cu(2(+)) exporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (32/41);  total (32/41)

ACCCAGCCACGGCAATTGCCATCATCATGCCGCCCTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCGGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATACTGCTTTTAGAAGTGAAAAATTAAAAGGGAGAGACTCCGCTCTCCC  >  NC_000913/784430‑784915                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                               aCCCAGCCACGGCAATTGCCATCATCATGCCGCCCTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGAGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGcc                                                                                                                                                                                                                                             >  1:664415/1‑251 (MQ=255)    cAGCCACGGCAATTTCCATCATCATGCCGCCCTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATcagccagc                                                                                                                                                                                                                                          <  2:541821/251‑1 (MQ=255)     aGCCACGGCCATTGCCATCATCATGCCGCCCTCGACCGGAAGCGCCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAACATGGTAATTAACACGAGCGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGAAGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGAACCGCATcagccagc                                                                                                                                                                                                                                          >  2:711572/1‑250 (MQ=255)     aGCCACGGCAATTGCCATCATCATGCCGCCCTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCa                                                                                                                                                                                                                                         <  1:164101/251‑1 (MQ=255)      gCCACGGCAATTGCCATCATCATGCCGCCCTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAATCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCaa                                                                                                                                                                                                                                        <  2:329644/251‑1 (MQ=255)         aCGGCAATTGCCATCATCATGCCGCCCTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGc                                                                                                                                                                                                                                     >  1:17668/1‑251 (MQ=255)               aTTGCCCTCATCATCACGCCCTCGACCGGACGCGGCGTGCGGAAGCGTTCTACCGCCTACCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCAAAAACGCCGCGCGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGGAGGGCGAAAAGCAAGGCGGCAGTAGCGGTCAACATATGACCCGCATCAGCCAGCAATGCCag                                                                                                                                                                                                                                  >  2:464790/1‑248 (MQ=255)               aTTGCCATCATCATGCCGCCCTCGACCGGACGCGGCGTGCGGCAGCGTTCTATCGCCTCCCAGAAACTCAAAATGGTAATTACCACGAGGGCTATGGCCTTCACAAACGCCGCGAGCGTGGTAAGTCTTAGCCAACCGAAAGTGGGGCGCATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACCTATGACCCGCATCAGCAAGCCATGCCag                                                                                                                                                                                                                                  >  1:296593/1‑248 (MQ=255)                      tcatcatGCCGCCCTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCaga                                                                                                                                                                                                                         >  1:97168/1‑250 (MQ=255)                         tcatGCCGCCCTCGACCGGACGCGGCGTGCGGCAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCagaa                                                                                                                                                                                                                        >  1:281477/1‑248 (MQ=255)                             gccgccCTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAac                                                                                                                                                                                                                 <  2:297051/251‑1 (MQ=255)                                 cccTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGATGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAaccacc                                                                                                                                                                                                             <  1:330555/251‑1 (MQ=255)                                  ccTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCCGAAAGAAAaccacc                                                                                                                                                                                                             >  2:610577/1‑250 (MQ=255)                                     cGACCGGACGCGGCGTGCGGAAGCGTTCTCGCGCCGCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCTTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAAGGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACg                                                                                                                                                                                                         <  1:22923/251‑1 (MQ=255)                                      gACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAACGAAAACCAACAACg                                                                                                                                                                                                         >  2:410273/1‑250 (MQ=255)                                      gACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACg                                                                                                                                                                                                         >  1:583942/1‑250 (MQ=255)                                         cGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAACCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCTAGCGTGGGCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACtt                                                                                                                                                                                                     <  1:423262/251‑1 (MQ=255)                                         cGGACGCGGCGTGCGGAAGCGTTCTAGCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACtt                                                                                                                                                                                                     <  1:218145/251‑1 (MQ=255)                                          ggACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCCTTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGGCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGCAAGAAAAccaccaccgcctt                                                                                                                                                                                                     >  2:325761/1‑243 (MQ=255)                                             cgcgGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTAc                                                                                                                                                                                                 <  2:343903/251‑1 (MQ=255)                                                 gcgTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGc                                                                                                                                                                                             <  1:19432/251‑1 (MQ=255)                                                         aaGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAAccc                                                                                                                                                                                     <  1:315491/251‑1 (MQ=255)                                                           gCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCgg                                                                                                                                                                                   <  2:481892/251‑1 (MQ=255)                                                             gTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGc                                                                                                                                                                                  >  1:398012/1‑250 (MQ=255)                                                             gTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCa                                                                                                                                                                                 <  1:446240/251‑1 (MQ=255)                                                             gTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCa                                                                                                                                                                                 <  1:569031/251‑1 (MQ=255)                                                                   ccGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACg                                                                                                                                                                           <  2:139643/250‑1 (MQ=255)                                                                    cGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACg                                                                                                                                                                           <  1:216234/250‑1 (MQ=255)                                                                    cGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGc                                                                                                                                                                          <  2:635380/251‑1 (MQ=255)                                                                    cGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGc                                                                                                                                                                          <  1:300706/251‑1 (MQ=255)                                                                         cccAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGaa                                                                                                                                                                      >  1:394278/1‑250 (MQ=255)                                                                                   cAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATacaaca                                                                                                                                                           <  1:385598/251‑1 (MQ=255)                                                                                              ttACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGACAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGc                                                                                                                                                 >  2:613262/1‑250 (MQ=255)                                                                                                  cACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCattat                                                                                                                                            <  2:414352/251‑1 (MQ=255)                                                                                                    cGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCattatt                                                                                                                                           >  1:637064/1‑250 (MQ=255)                                                                                                     gaggtcttaggCATTCGCAAACGCCGCGGTCGTGGTCAGTCGTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGGTGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGTCCAACATATGACCCGCATCAGCCAGCAATGCCCGAGAACCAGAAAGAAAAACACCATCGTCATCTGCCAGCACAATTCCAGCAGTCATGCCGAAAGCATACAACAAGCGCCGAGCattatgat                                                                                                                                         <  2:542738/242‑4 (MQ=255)                                                                                                      aGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTAtc                                                                                                                                        >  1:291331/1‑251 (MQ=255)                                                                                                        ggCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTAtctt                                                                                                                                      >  1:484919/1‑251 (MQ=255)                                                                                                                 aTTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCAGgt                                                                                                                             <  1:689386/251‑1 (MQ=255)                                                                                                                     acaAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCAGGTGTGa                                                                                                                         <  2:729150/251‑1 (MQ=255)                                                                                                                          cgccgcGAGCGTGGTCAGTCTTAGCCAACCGACAGTGTGGCGAATGGCGGGAGGCCGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGAAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAAACTTATTATCATCTGGCAGGtgtgcagccg                                                                                                                    >  2:294051/1‑245 (MQ=255)                                                                                                                           gccgcGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCAGGTGTGAGGACgt                                                                                                                   <  1:711572/251‑1 (MQ=255)                                                                                                                           gccgcGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCAGGTGTGAGGACgt                                                                                                                   >  1:348233/1‑251 (MQ=255)                                                                                                                                               tAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTgcg                                                                                                  <  1:103559/248‑1 (MQ=255)                                                                                                                                               tAGCCAACCGAAAGTGTGGCGAATGGGGTGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGcc                                                                                                <  2:680915/250‑1 (MQ=255)                                                                                                                                                    aaCCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGa                                                                                            >  1:332916/1‑249 (MQ=255)                                                                                                                                                      ccGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCAGGTGTGAGGACGTGTGTGAGTGTGGGTGCGCCATGACg                                                                                          >  2:282903/1‑249 (MQ=255)                                                                                                                                                         aaaGGGTGGCGAATGGTGGGAGGACGACGTGAAAATGGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCa                                                                                     <  2:309436/251‑1 (MQ=255)                                                                                                                                                                ggCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCAGCATCAGCCAGCAATGCCAGAGAACCCGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGAATTATTATCT‑CTGGCAGGGGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATCCCTTa                                                                              >  2:451829/1‑250 (MQ=255)                                                                                                                                                                ggCGAATGGTGGGAGGACGACGTGAAAAGTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTAGCGGGCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCTCCAGCGACTTCTACCAGCATAACCCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTAACTTCTGGCAGGTGTGTGGGCGTGTGTGAGTGTGAGTTCGCCATGACGAGGCAACCCTTa                                                                              <  2:664415/251‑1 (MQ=255)                                                                                                                                                                   gAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACGTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCAGGTGTGAGGACTTGTGTGAGTGTCAGTGCGCGATGACGATTCATCCCTTAatt                                                                           <  2:373086/251‑1 (MQ=255)                                                                                                                                                                    aaTGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATCCCTTAatt                                                                           >  1:162521/1‑250 (MQ=255)                                                                                                                                                                     aTGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCT‑CTGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATCCCTTAattatt                                                                        <  1:451829/251‑1 (MQ=255)                                                                                                                                                                           ggaggaCGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATCCCTTAATTATTCTCAc                                                                   >  1:638156/1‑251 (MQ=255)                                                                                                                                                                            gaggaCGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATCCCTTAATTATGCTCACt                                                                  <  1:555932/251‑1 (MQ=255)                                                                                                                                                                             aggaCGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCatta                                                                                                                                             <  2:408969/175‑1 (MQ=255)                                                                                                                                                                             aggaCGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCatta                                                                                                                                             >  1:408969/1‑175 (MQ=255)                                                                                                                                                                                acgactTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTAGCGGTCAACATATGACCCGCCTCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGACATCCCTTAATTATTCTCACTTAAc                                                              <  2:436781/251‑1 (MQ=255)                                                                                                                                                                                          aaTTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATc                                                                                   >  1:626603/1‑220 (MQ=255)                                                                                                                                                                                          aaTTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATc                                                                                   <  2:626603/220‑1 (MQ=255)                                                                                                                                                                                           aTTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGAcatc                                                   <  2:17668/251‑1 (MQ=255)                                                                                                                                                                                              gCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATa                                                <  2:414660/251‑1 (MQ=255)                                                                                                                                                                                                aCGGCGAGTAGGGCGAAAAGCAGGGCGTCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATACt                                              <  1:656270/251‑1 (MQ=255)                                                                                                                                                                                                aCGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATACt                                              >  1:453719/1‑251 (MQ=255)                                                                                                                                                                                                  ggCGAGTCGGGCGAAAAGCACGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCCGAGAACCAGAAAGAAAACCACCCACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGCCGAGCATTATTATCTTCTGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATCCCTTAATTCTTCTCACTTAACTTTATGACATCATACTGc                                            >  2:83561/1‑251 (MQ=255)                                                                                                                                                                                                    cGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCGCCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATCCCTTAATTATTCTCACTTACCTTTATGACATCATACTGCtt                                          <  2:557899/251‑1 (MQ=255)                                                                                                                                                                                                          gggCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCAGGTGTGAGGACGTGTGTGAGGGTGAGTGCGCCATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATACTGCTTTTAGaa                                    <  1:325761/251‑1 (MQ=255)                                                                                                                                                                                                          gggCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACACGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATACTGCTTTTAGa                                     >  1:692871/1‑250 (MQ=255)                                                                                                                                                                                                               aaaaGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATACTGCTTTTAGAAGTGa                                >  1:479397/1‑250 (MQ=255)                                                                                                                                                                                                                      gggCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATACTGCTTTTAGAAGTGAAAAATTa                         >  2:292149/1‑250 (MQ=255)                                                                                                                                                                                                                            cAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATACTGCTTTTAGAAGTGAAAAATTAAAAggg                   >  2:274238/1‑250 (MQ=255)                                                                                                                                                                                                                                      cAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATACTGCTTTTAGAAGTGAAAAATTAAAAGGGAGAGACTCCg         <  1:553367/250‑1 (MQ=255)                                                                                                                                                                                                                                       aaCATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCATGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATACTGCTTTTAGAAGTGAAAAATTAAAAGGGAGAGACTCCGct       <  2:620432/251‑1 (MQ=255)                                                                                                                                                                                                                                            atGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATACTGCTTTTAGAAGTGAAAAATTAAAAGGGAGAGACTCCGCTCTccc  <  1:274238/251‑1 (MQ=255)                                                                                                                                                                                                                                            atGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATACTGCTTTTAGAAGTGAAAAATTAAAAGGGAGAGACTCCGCTCTccc  >  2:479577/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                               ACCCAGCCACGGCAATTGCCATCATCATGCCGCCCTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCGGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATACTGCTTTTAGAAGTGAAAAATTAAAAGGGAGAGACTCCGCTCTCCC  >  NC_000913/784430‑784915

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