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%	171.0 / NA	51	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 (22/29);  total (22/29)

GCCACGGCAATTGCCATCATCATGCCGCCCTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCGGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATACTGCTTTTAGAAGTGAAAAATTAAAAGGGAGAGACTCCGCTCTCCCATTATTGG  >  NC_000913/784435‑784923                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                       gCCACGGCAATTGCCATCATCATGCCGCCCTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCATACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGGGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGTGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCaa                                                                                                                                                                                                                                                <  2:533798/251‑1 (MQ=255)   cACGGCAATTGCCATCATCATGCCGCCCTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATg                                                                                                                                                                                                                                              >  1:50016/1‑251 (MQ=255)          aTTGCCATCATCATGCCGCCCTCGACCGGCCGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCagcg                                                                                                                                                                                                                                        >  2:524701/1‑248 (MQ=255)                     catGCCGCCCTCGACCGGACGAGGCGTGCGGATGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCCTTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATTTGACCCGCATCAGCCAGCCATGCCagagcacacgaacgg                                                                                                                                                                                                                             >  2:303268/1‑239 (MQ=255)                            cccTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAAGTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAaccacc                                                                                                                                                                                                                     <  2:558754/251‑1 (MQ=255)                             ccTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAAccacca                                                                                                                                                                                                                    <  1:533752/251‑1 (MQ=255)                                               tGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCAt                                                                                                                                                                                                   >  1:182595/1‑250 (MQ=255)                                                  ggAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAAc                                                                                                                                                                                               >  1:10145/1‑251 (MQ=255)                                                             aTCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACTCCTCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAAGTGCACGGCGAGTAGGGCGAAAAGCAGGTCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCa                                                                                                                                                                                     <  1:243889/250‑1 (MQ=255)                                                               cGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCTCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACg                                                                                                                                                                                   <  1:320682/250‑1 (MQ=255)                                                                     ccAGACAAACAAAATGGTAATTACCACGAGGGCGATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGGGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTAGCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAg                                                                                                                                                                            <  1:60959/251‑1 (MQ=255)                                                                              cAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATacaac                                                                                                                                                                    <  2:50016/250‑1 (MQ=255)                                                                                           aCCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCat                                                                                                                                                       <  1:524701/250‑1 (MQ=255)                                                                                             cACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAACACCACAAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCcttat                                                                                                                                                    >  2:478943/1‑251 (MQ=255)                                                                                                gAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGAGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGACAGTCACGCCGAAAGAATACAACAAGCGCCGAGCATtattat                                                                                                                                                 >  2:308560/1‑251 (MQ=255)                                                                                                gAGGGCTATAGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGCAAGCATACAACAAGCGACGAGCATtattat                                                                                                                                                 >  1:452905/1‑251 (MQ=255)                                                                                                 aGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTAtc                                                                                                                                                <  2:272560/251‑1 (MQ=255)                                                                                                         ggCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCa                                                                                                                                        <  1:112692/251‑1 (MQ=255)                                                                                                                            aGCGTGGTCAGTCTTAGCCAACCGAAAGAGTGGCGAATGGGGGGAGGACGACGTGAAAATTGCAAGGCGCGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGCAAGCATACAACAAGCGACGAGCATTATTATCTTCTAGCAGGTGTGAGGACGtgtgtg                                                                                                                      >  2:290744/1‑250 (MQ=255)                                                                                                                              cGTGGTCAGTCTTAGACAACCGAAAGTGTGGCGCATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCAGGTGTGAGGACGTgtgtgagt                                                                                                                   >  2:343658/1‑251 (MQ=255)                                                                                                                               gTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGCGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCAGGTGTGAGGACGTgtgtgagtg                                                                                                                  >  2:520211/1‑251 (MQ=255)                                                                                                                                      gTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCTGTTAACATATGACCCGCATCAGTCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACGGCTACCATTTTAAACCCGGCAGTCACGCCGATAGCATACTACAAGTGACTAGGATTATTATCTTCTGGCAGGTGTGAGGACGTGTGTGAGTGTTAGTgc                                                                                                           <  1:87718/251‑1 (MQ=255)                                                                                                                                              cAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGAc                                                                                                   <  1:254931/251‑1 (MQ=255)                                                                                                                                              cAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGAc                                                                                                   >  1:142488/1‑251 (MQ=255)                                                                                                                                               aaCCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCCGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACg                                                                                                  >  1:299112/1‑251 (MQ=255)                                                                                                                                               aaCCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACg                                                                                                  <  2:444208/251‑1 (MQ=255)                                                                                                                                                           ggCGAATGGGGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCATCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATcccgg                                                                                       <  2:331362/250‑3 (MQ=255)                                                                                                                                                            gCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCAGGTGTGAGGACGTGTGGGAGTGTGAGTGCGCCAGGACGAGTCATCCCTTaa                                                                                     <  2:109235/251‑1 (MQ=255)                                                                                                                                                             cGAATGGTGGGAGGACGACGTGAAAAGTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCAGGTGTGAGGACGGGTGTGAGTGTGAGGGCGCCATGACGAGTCAGCCCTTaag                                                                                    <  2:350342/251‑2 (MQ=255)                                                                                                                                                              gAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCCGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCCTCCCTTAatt                                                                                   >  2:262792/1‑251 (MQ=255)                                                                                                                                                                   gtggGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATCCCTTAATTATtct                                                                              <  2:142488/251‑1 (MQ=255)                                                                                                                                                                      ggaggaAGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGCGAACCAGAAAGAAAACCCCCAACGACTACTACCAGCATAAACCCGGCAGTAACGCCGAAAGAATACAACAAGCGACGAGCATTATTATCTTCTGGCAGGTGTGAGGACGTGTGTGAGTGAGCGTGCGCCATGACGAGCCATCCCTTAATTATTCTCAc                                                                           >  2:376711/1‑251 (MQ=255)                                                                                                                                                                                   aaaaTTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCAGCCCTTAATTATTCTCACTTAACTTTATGAc                                                              <  1:478943/251‑1 (MQ=255)                                                                                                                                                                                               cGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCAGACAACAAGCGACGAGCATTATGAGCTTCTGGCAGGTGTGAGGACGTGGGTGAGTGTGAGTGCGCCATGACGAGTCATCCCGTAAGTAGGCTCACTTAACTTTATGACAGCATACTGCtt                                                  <  2:80851/251‑1 (MQ=255)                                                                                                                                                                                                  gTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGGCATCCCTTAATTATTCTCACTTAACTTTATGACATCATACTGCTTTTa                                               <  2:299112/251‑1 (MQ=255)                                                                                                                                                                                                       gCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTATCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTAGCTTCTGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATACTGCTTTTAGAAGt                                          <  1:521250/251‑1 (MQ=255)                                                                                                                                                                                                                 gggCGGCAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATACTGCTTTTAGAAGTGAAAAATTaa                                <  1:280202/251‑1 (MQ=255)                                                                                                                                                                                                                  ggcggcAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATACTGCTTTTAGAAGTGAAAAATTaaa                               <  1:520211/251‑1 (MQ=255)                                                                                                                                                                                                                     ggcAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATACTGCTTTTAGAAGTGAAAAATTAAAAgg                            <  1:25310/251‑1 (MQ=255)                                                                                                                                                                                                                       cAGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATACTGCTTTTAGAAGTGAAAAATTAAAAGGga                          <  2:10145/251‑1 (MQ=255)                                                                                                                                                                                                                        aGTATCGGTCAACATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATACTGCTTTTAGAAGTGAAAAATTAAAAggg                           >  2:169740/1‑249 (MQ=255)                                                                                                                                                                                                                                   aCATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCCCCAACGCCTTCTACCAGCATAAACCCGGCAGTCACGCCGCACGCATACAACAAGCGACGAGCATTATTATCTTCTGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATCCCTAAATTCTTCTCCCTTACCTTTATGCCATCATACTGCTTTTAGAAGTGAAAACTTAAAAGGGAGAGACTCCGctc              >  2:551103/1‑251 (MQ=255)                                                                                                                                                                                                                                   aCATATGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCAGGTGGGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATCCCTTATTTATTCTCACTTAACTTTATGACATCATACTGCTTTTAGAAGTGAAAAATTAAAAGGGAGAGACTCCGctc              >  1:130821/1‑251 (MQ=255)                                                                                                                                                                                                                                      tatGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATACTGCTTTTAGAAGTGAAAAATTAAAAGGGAGAGACTCCGCTCTcc           >  1:316854/1‑251 (MQ=255)                                                                                                                                                                                                                                       atGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCTGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATACTGCTTTTAGAAGTGAAAAATTAAAAGGGAGAGACTCCGCTCTccc          <  1:553594/251‑1 (MQ=255)                                                                                                                                                                                                                                       atGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATACTGCTTTTAGAAGTGAAAAATTAAAAGGGAGAGACTCCGCTCTccc          >  1:454470/1‑251 (MQ=255)                                                                                                                                                                                                                                       atGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATACTGCTTTTAGAAGTGAAAAATTAAAAGGGAGAGACTCCGCTCTcc           <  2:182595/250‑1 (MQ=255)                                                                                                                                                                                                                                        tGACCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAAAGACTTCTACCAGAATAAACACGGCCGTCAACACGAAAGCATACAACAAGCGACGAGAATTATTATCTTCTGGCAGGTGAGAGGACGTGTGTGAGTGTGAGTGCGCAATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATACTGATTTTAGAAGTGAAAAATTAAAAGGGAGAGACTCCGCTCTccc          >  1:543614/1‑250 (MQ=255)                                                                                                                                                                                                                                          aCCCGCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATACTGCTTTTAGAAGTGAAAAATTAAAAGGGAGAGACTCCGCTCTCCCatt       >  2:156444/1‑251 (MQ=255)                                                                                                                                                                                                                                              gCATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATACTGCTTTTAGAAGTGAAAAATTAAAAGGGAGAGACTCCGCTCTCCCATTATTg   <  1:290744/251‑1 (MQ=255)                                                                                                                                                                                                                                               cATCAGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACCTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACTAGCATTATGATCTTCTGTCAGGTGCGAGGACGTGTGTGAGTGTTAGTGCGCCATGACGAGTCATCCCTTAATTATACTCACTTAACTTTATGACATCATACTGCTTTTAGAAGTGAAAAATTAAAAGTCAGAGACACCGCTCTCCCATTATTgg  <  2:543614/251‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                       GCCACGGCAATTGCCATCATCATGCCGCCCTCGACCGGACGCGGCGTGCGGAAGCGTTCTATCGCCTCCCAGACAATCAAAATGGTAATTACCACGAGGGCTATGGCATTCACAAACGCCGCGAGCGTGGTCAGTCTTAGCCAACCGAAAGTGTGGCGAATGGTGGGAGGACGACGTGAAAATTGCACGGCGAGTAGGGCGAAAAGCAGGGCGGCAGTATCGGTCAACATATGACCCGCATCGGCCAGCAATGCCAGAGAACCAGAAAGAAAACCACCAACGACTTCTACCAGCATAAACCCGGCAGTCACGCCGAAAGCATACAACAAGCGACGAGCATTATTATCTTCTGGCAGGTGTGAGGACGTGTGTGAGTGTGAGTGCGCCATGACGAGTCATCCCTTAATTATTCTCACTTAACTTTATGACATCATACTGCTTTTAGAAGTGAAAAATTAAAAGGGAGAGACTCCGCTCTCCCATTATTGG  >  NC_000913/784435‑784923

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: ‑