Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,216,473	A→G	Y7H (TAC→CAC)	nfeF ←	NADPH‑dependent ferric‑chelate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,216,473	0	A	G	100.0%	238.8 / NA	65	Y7H (TAC→CAC)	nfeF	NADPH‑dependent ferric‑chelate reductase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (25/40);  total (25/40)

ATAGTCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACC‑GGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACAC  >  NC_000913/3216228‑3216716                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    atatgcACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACC‑GGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCggg                                                                                                                                                                                                                                                <  2:297628/246‑1 (MQ=255)      cACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACC‑GGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGgtg                                                                                                                                                                                                                                           <  1:541872/251‑1 (MQ=255)      cACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACC‑GGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGgt                                                                                                                                                                                                                                            <  1:417984/250‑1 (MQ=255)      cACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACC‑GGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCCGCGGGTGGCGGGGGgtg                                                                                                                                                                                                                                           >  2:503886/1‑251 (MQ=255)              ggTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACC‑GGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCAt                                                                                                                                                                                                                                   <  1:156781/251‑1 (MQ=255)                          ccTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACC‑GGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCtt                                                                                                                                                                                                                       <  1:414440/251‑1 (MQ=255)                             tCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACC‑GGCGCTGATGCGCTCAACGCGTAACACAGGCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCg                                                                                                                                                                                                                     <  2:97162/250‑1 (MQ=255)                                   cAGACCATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGTTCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACC‑GGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCTGGGGGTGTTATTCATTTTTATCGCATTCTTGATGGt                                                                                                                                                                                                              <  2:140657/251‑1 (MQ=255)                                       cGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGGGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCGCTCCGCCGAGGACAATGCGCTGAAAACC‑GGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAACCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGCGGGTGTTATTCATTGTTATCGCCTTTGTGATGGTAATc                                                                                                                                                                                                          <  1:383141/251‑1 (MQ=255)                                          tGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACC‑GGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATc                                                                                                                                                                                                          >  1:369927/1‑248 (MQ=255)                                           gCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACC‑GGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGat                                                                                                                                                                                                      <  2:107192/251‑1 (MQ=255)                                                  tCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACC‑GGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATc                                                                                                                                                                                                 <  1:500754/249‑1 (MQ=255)                                                     gTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACC‑GGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCtaaat                                                                                                                                                                                            <  1:276429/251‑1 (MQ=255)                                                     gTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACC‑GGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCtaaat                                                                                                                                                                                            >  2:511545/1‑251 (MQ=255)                                                             tGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACC‑GGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTc                                                                                                                                                                                     <  2:503314/250‑1 (MQ=255)                                                             tGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACC‑GGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTc                                                                                                                                                                                     <  2:551904/250‑1 (MQ=255)                                                              ggcggcACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACC‑GGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGc                                                                                                                                                                                   <  1:147669/251‑1 (MQ=255)                                                                     caAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACC‑GGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAaatgat                                                                                                                                                                            <  2:490659/251‑1 (MQ=255)                                                                        aGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACC‑GGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAaatgataat                                                                                                                                                                         <  2:436058/251‑1 (MQ=255)                                                                              cGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACC‑GGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGt                                                                                                                                                                   >  1:24399/1‑251 (MQ=255)                                                                                   ggTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACC‑GGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATc                                                                                                                                                              <  2:496576/251‑1 (MQ=255)                                                                                       gAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACC‑GGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGa                                                                                                                                                          <  1:511545/251‑1 (MQ=255)                                                                                       gAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACC‑GGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGa                                                                                                                                                          <  1:279082/251‑1 (MQ=255)                                                                                        aGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACC‑GGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGAt                                                                                                                                                         <  2:328534/251‑1 (MQ=255)                                                                                         ggAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACC‑GGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGGGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGAt                                                                                                                                                         >  2:205537/1‑250 (MQ=255)                                                                                           aaaGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACC‑GGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGGAAATGATAATGATTGTTAATCATGATaaa                                                                                                                                                      <  1:448584/251‑1 (MQ=255)                                                                                           aaaGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACC‑GGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAA‑CTCGCAAATGATAATGATTGTTAATCATGATAAAt                                                                                                                                                     >  1:225968/1‑251 (MQ=255)                                                                                                   gTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACC‑GGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCg                                                                                                                                              >  2:235516/1‑251 (MQ=255)                                                                                                     ttGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACC‑GGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGAt                                                                                                                                            >  1:270737/1‑251 (MQ=255)                                                                                                       gCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACC‑GGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGAttt                                                                                                                                          >  2:158324/1‑251 (MQ=255)                                                                                                          gtgATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACC‑GGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGt                                                                                                                                        >  1:560720/1‑250 (MQ=255)                                                                                                          gtgATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACC‑GGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTa                                                                                                                                       <  2:217465/251‑1 (MQ=255)                                                                                                            gATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACC‑GGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAg                                                                                                                                      <  2:461032/250‑1 (MQ=255)                                                                                                               cGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACC‑GGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACt                                                                                                                                  >  2:297670/1‑251 (MQ=255)                                                                                                                     aGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACC‑GGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATg                                                                                                                            >  1:230634/1‑251 (MQ=255)                                                                                                                                  naaaCCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACC‑GGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGAt                                                                                                               >  2:250840/2‑251 (MQ=255)                                                                                                                                  aaaaCCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACC‑GGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGAt                                                                                                               <  1:158324/251‑1 (MQ=255)                                                                                                                                   aaaCCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACC‑GGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATa                                                                                                              <  2:330376/251‑1 (MQ=255)                                                                                                                                      ccGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACC‑GGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATa                                                                                                              >  1:521086/1‑248 (MQ=255)                                                                                                                                       cGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACC‑GGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGAc                                                                                                          <  1:250840/251‑1 (MQ=255)                                                                                                                                          ccAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACC‑GGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGa                                                                                                        >  2:34647/1‑250 (MQ=255)                                                                                                                                             gcgcCTCGCCGCCGAGGACAATGCGCTGAAAACC‑GGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCg                                                                                                    <  1:511634/251‑1 (MQ=255)                                                                                                                                                 cTCGCCGCCGAGGACAATGCGCTGAAAACC‑GGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCctc                                                                                                >  1:520503/1‑251 (MQ=255)                                                                                                                                                     ccgccgAGGACAATGCGCTGAAAACC‑GGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCtcttct                                                                                            >  1:225110/1‑251 (MQ=255)                                                                                                                                                               cAATGCGCTGAAAACC‑GGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGcatcat                                                                                  >  1:94470/1‑251 (MQ=255)                                                                                                                                                                aaTGCGCTGAAAACC‑GGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCAtt                                                                                 >  1:514809/1‑251 (MQ=255)                                                                                                                                                                 aTGCGCTGAAAACC‑GGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCAtt                                                                                 >  2:128847/1‑250 (MQ=255)                                                                                                                                                                     gcTGAAAACC‑GGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTc                                                                             >  1:534458/1‑250 (MQ=255)                                                                                                                                                                      gcgAAAACCCGGCGCTGATGCGCTCAAGGCGTAACACATTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTTTTATTCATTGTTATCGCCTTCGTGTTGGTAATCAGATATATCTAAATAAAACTCGCAAAGGATATTGATTGTTAATCATGATAAATGCAAGCGATTTTTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCa                                                                            <  1:110380/249‑1 (MQ=255)                                                                                                                                                                        gAAAACC‑GGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTTTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCAAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAAc                                                                          <  1:99914/250‑1 (MQ=255)                                                                                                                                                                          aaaCC‑GGCGCTGATGCGCTAAACGCGTAACACAGTCAGTTCACGGAAGAGAAGATCATGGCGAACCCGCGGGGGGTGGCGGGGGTTGTTATTCATTGTTATCGCTTTGGTAATGGTAATCAGATATATCTAAAAAAAACTCGCAAATAATAATGATTGTAAACCATGATAAATGCAACCGATTTGTAGATTTAATATGTTTATGGTCTAATAGGGCTAACCGCCTCTTCTCAGGCATCATTACTCAACGcc                                                                       <  1:238858/251‑1 (MQ=255)                                                                                                                                                                             tctggCGCGGATGCGCTCAACGCGTAACTCAGCCAGTTCACGGATGCGCAGATCATTGCAACCGCGCTGCGGGGGGTGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAAGAAAACACGCAAATAATAATGATTTTAAATCATGATAAATTCAAGCGATTTGTAGACCTGATATGTCTATGGTCGAATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCg                                                                      <  1:50790/247‑1 (MQ=255)                                                                                                                                                                                ggCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCg                                                                      >  2:50790/1‑247 (MQ=255)                                                                                                                                                                                         gcgcTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGa                                                          >  1:336053/1‑250 (MQ=255)                                                                                                                                                                                               aaCGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGGGGGTGGGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCtt                                                   <  2:369927/251‑1 (MQ=255)                                                                                                                                                                                                    gTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCg                                              >  1:560813/1‑251 (MQ=255)                                                                                                                                                                                                       acacAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGGGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGCTGCGGCGTGAACGCCTTATACGGcc                                           >  2:107549/1‑251 (MQ=255)                                                                                                                                                                                                                 ttaaCGGAAGGGCAGATCTTTGCGACCGCGCTGGGGGTGGCGGGGGGTTTTATTCTTTGTTATCGCCTTCGTGATGTTAATCAGATATACCAAAATAAAACTCGAAAATGAAAATGATTGTTAATAATGATAAATGCAGGAGATTTGTAGAACTGATATGTCTGTAGTCAGATAGGACGAACCGCCTCTTCTCAGGCATCATGACTCAACGCCGGATGCGGCGTGAACGCCTTATCCAGCTTACGTGTga                                  <  2:485490/247‑1 (MQ=255)                                                                                                                                                                                                                       gAAGCGCAGATCTTTGCGACCGCCCTGCGTGTGGCGGGTGGTTTTTTTCATTGTTTTCGCTTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATATTGATTTTTAATCATGATAAATGCAAGCGTTTTGTAGAACTGATATTTCTATAGTCTGATAAGCCGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGt                           <  1:95685/251‑1 (MQ=255)                                                                                                                                                                                                                          gcgcAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTCCGTGATGGTAATCAGATATATCTAAATAAAACTCGAAAATGAAAATGATTGTTAATCAGGATAAATGCAAGCGATTTGTAGAACTGATATGTCAATAGTCTGAGAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCa                        <  2:270737/251‑1 (MQ=255)                                                                                                                                                                                                                               gatatttGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGATAATCAGATATATCTAAAAAAAACTCGAAAATGATAATGATTGTTAATCATGAGAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTTTGATAAGACAAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGttt                   <  1:205537/246‑1 (MQ=255)                                                                                                                                                                                                                               gATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGttt                   <  1:503886/251‑1 (MQ=255)                                                                                                                                                                                                                                            gctcTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAAAACTGATATTTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCAttttt      <  2:520503/248‑1 (MQ=255)                                                                                                                                                                                                                                            gcgcTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGAAAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCAttttt      <  1:362714/251‑1 (MQ=255)                                                                                                                                                                                                                                                 gCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCTTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAACCATGATAAATGCAAGCGATTTGTAGACCTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTacac  <  1:304524/250‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    ATAGTCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACC‑GGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACAC  >  NC_000913/3216228‑3216716

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