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%	134.1 / NA	48	Y7H (TAC→CAC)	nfeF	NADPH‑dependent ferric‑chelate reductase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (14/34);  total (14/34)

CACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCG‑TGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACG  >  NC_000913/3216233‑3216696                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                 cACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACACAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCATGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGCACGCGCTGCGGGTGGCGGGGGgtg                                                                                                                                                                                                                        >  2:133827/1‑251 (MQ=255)  ccgcgACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGgtgt                                                                                                                                                                                                                       <  2:184071/250‑1 (MQ=255)             ggCGTGGTCCTTCCGGCCAGACGATGCCCTATTCCGTTACCGTTGGCGGCACAAAGTGAGAGTCAGGTTGCGGACCGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGCTGAAAACCCGTCCCGCGCATCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGtt                                                                                                                                                                                                            >  1:56357/1‑251 (MQ=255)                gTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATc                                                                                                                                                                                                         <  1:10166/251‑1 (MQ=255)                gTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATc                                                                                                                                                                                                         <  2:323012/251‑1 (MQ=255)                gTGGTCCGTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATc                                                                                                                                                                                                         <  2:328277/251‑1 (MQ=255)                  ggTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGATGGCGGCACAAAGTGAGCGTCAGGTTGAGGAACGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGCGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGc                                                                                                                                                                                                       >  1:208422/1‑251 (MQ=255)                     ccTTCCGGCCAGACGATTCCCTCTGCCGTTACCGTTGGCGTCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGAGGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCt                                                                                                                                                                                                     <  1:265548/250‑1 (MQ=255)                     ccTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGGGTTATTCATTGTTATCGCCtt                                                                                                                                                                                                    >  1:126576/1‑251 (MQ=255)                         ccGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCG‑t                                                                                                                                                                                                <  1:297273/250‑1 (MQ=255)                         ccGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTGATCGCCTTCG‑t                                                                                                                                                                                                <  1:128562/250‑1 (MQ=255)                              cAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCTGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCG‑TGATGGt                                                                                                                                                                                          <  1:405383/251‑1 (MQ=255)                                    aggCCCTCGTCCGGTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTGGCTGTGATCGTCAAATCCACTCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAGACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTAGTCATTGGTATCGCCTTCG‑TTATTGTAATCAg                                                                                                                                                                                    <  2:155269/249‑1 (MQ=255)                                             gccGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGGTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCG‑TGATGGTAATCAGATATATc                                                                                                                                                                             <  1:126675/248‑1 (MQ=255)                                                gTTACCGTTGGCGGCACCCAGCGAGCGTCAGGTTGAGGCAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTACAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCCTTGCGAACGCGCTGCGGGTGGCGGGGGGGGTTATTCATTGTTATCGCCTTCG‑TGGTGGTACTCCGCTATATCtcaat                                                                                                                                                                        >  2:400766/1‑251 (MQ=255)                                                   aCCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCG‑TGATGGTAATCAGATATATCtaaataa                                                                                                                                                                      <  2:271277/250‑1 (MQ=255)                                                            ggcACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCG‑TGATGGTAATCAGATATATCTAAATAAAACTCGCaa                                                                                                                                                             >  2:101657/1‑250 (MQ=255)                                                             gcACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGATGTGATAGACAAAGGAACGCGAGGGAAAACCGTCCAGCGGCTCGACGCCGAGGACAATGCGCTGAAAAACGGCGCTGAAGAGGTCAAAGCGAAACACAGTCAGTTCACGGAAGCGAAGATCATTGCGAACGCGCTGCGGGTGGGGGGGGGGGGTATTCAATGGTATCGCCTTCG‑TGATGGGAATCAGAGATAACTAAAAAAAACTCGCaaa                                                                                                                                                            >  2:407788/1‑250 (MQ=255)                                                              cacaAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCATGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCG‑TGATGGTAATCAGATATATCTAAATAAAACTCGCAaatg                                                                                                                                                          <  1:133827/251‑1 (MQ=255)                                                               acaAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCG‑TGATGGTAATCAGATATATCTAAATAAAACTCGCAaatg                                                                                                                                                          >  1:127977/1‑250 (MQ=255)                                                               acaAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGGTCACGGAAGCGCAGATCATTGCGAACGCGCTTCGGGTGGCGGGGGGTGTTATTCATGGTTAGCGCCTTCG‑TGATGGTATTCAGATATATCTAAATAAAACTCGCAaatga                                                                                                                                                         <  1:436794/251‑1 (MQ=255)                                                                caAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTAGCGCCTTCG‑TGATGGTAATCAGAGATATCTAAATAAAACTCGCAaatgat                                                                                                                                                        <  2:405950/251‑1 (MQ=255)                                                                    gTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCG‑TGATGGTAATCAGATATATCTAAATAAAACTCGCAaatgataatg                                                                                                                                                    <  1:400766/251‑1 (MQ=255)                                                                         cGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCG‑TGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGt                                                                                                                                               >  1:322227/1‑251 (MQ=255)                                                                               gTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCG‑TGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCa                                                                                                                                         <  2:289030/251‑1 (MQ=255)                                                                                           agagTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCG‑TGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCaa                                                                                                                             <  2:127977/251‑1 (MQ=255)                                                                                                  gCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCG‑TGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGAttt                                                                                                                      <  1:390910/251‑1 (MQ=255)                                                                                                       gATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCG‑TGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGt                                                                                                                    >  1:273594/1‑248 (MQ=255)                                                                                                           gTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGTACAATGCGCTGAAAACCGGCGCTGAGGCGCTCAACGCTTAACACAGTCAGTTCACGGAAGCTCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCG‑TGATGGTAATCAGATATATCTAAATAAAACTCGCAAATTATAATTATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTg                                                                                                             <  1:239558/251‑1 (MQ=255)                                                                                                           gTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAGGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCG‑TGATGGTAATCAGATATATCTAAAGAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTg                                                                                                             <  2:280956/251‑1 (MQ=255)                                                                                                                             aaaaCCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCG‑TGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGAt                                                                                           <  1:365046/251‑1 (MQ=255)                                                                                                                               aaCCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCG‑TGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGa                                                                                            <  2:83676/248‑1 (MQ=255)                                                                                                                               aaCCGTCCAGCGCCTCGCCGACGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCG‑TGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGa                                                                                            >  1:83676/1‑248 (MQ=255)                                                                                                                                           ccTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCG‑TGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCct                                                                             <  1:396824/251‑1 (MQ=255)                                                                                                                                           ccTCCCCGGCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGGTATCGCCTTCG‑TGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCct                                                                             <  2:273594/251‑1 (MQ=255)                                                                                                                                              cgccgcCGAGGACAATGCGCTGAAAACCGGCGGTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCG‑TGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACGGATATGTCTATAGTCTGATAAGACGAACCGCCtctt                                                                          <  2:322227/251‑1 (MQ=255)                                                                                                                                                 cgcggAGGACAGTGCGCTGAAAACCGGCTCTGAGGCGCTCAACGCGTACCAACGTCATTTCACGGAAGCGCAGATCATTGGGACCGCGCGGGGGGTGGGGGGGGGTGTTATTCATTGTTATCGCCTTTG‑TGAGGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGAAAATTGCAAGCGATTTGTAGACCTGATATGTATATAGTTGGATAAGACGAACCGCCTCTtctc                                                                       <  2:398697/247‑1 (MQ=255)                                                                                                                                                 cgccgAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCG‑TGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTtctc                                                                       <  1:16843/251‑1 (MQ=255)                                                                                                                                                              gcgcgaaaagccgccgcTGATGCCCTCAACGCGTACCACGGACAGTTCCCGGAATCGCAGATCTCTGCGACCGCGCTGCGGGGGGCGGGGGGTTTTATTCTTTGTTATCGCCTTCG‑TGCTGGTAATCAAATATATCGAAAAAAGACTCTCAAATTATAATGATTGTCAACCTTGACAAAAGCAGGCTCTTCGTAGAACTGATATGTCGATTCTCTGATAGGACGAACCGCCTCTTCTCAGGCATCATTACt                                                          <  1:300557/237‑1 (MQ=255)                                                                                                                                                                 cTGAATACCGGCGCTGATGCGCTCAACGCGTAACACCGTCAGTTCACGGAAGCGCAGCTCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCG‑TGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTACTCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCa                                                        >  2:273970/1‑250 (MQ=255)                                                                                                                                                                          ggCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCG‑TGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATg                                              >  1:205497/1‑251 (MQ=255)                                                                                                                                                                           gcccTGATGCCCTCAACGCGTAACACATTCATTTCCCGGAAGCTCAGATCATTGCGAACGCGCTGCGGGGGGCGGGGGGTTTTATCCATTGTTATCTCCTTCG‑TAATGGTAATCAGATATATCGAAATAAACCTCTCAAATTATAATGATTTTTAATCATGATAAATGTAAGGTTTTTGTATAAGTGATATGTCTATATTCAGATAATACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATgc                                             <  1:407788/248‑1 (MQ=255)                                                                                                                                                                                        aaccgcgCAACACAGTCAGTTCACGGAAGCGCAGATCATTGCAAACGCGCTGGGGGTGGCGGGGGGTGTTATTCATTTTTATCGCCTTCG‑TGATGGTAATCAGATATATCAAAATAAAACTCGCAAATGATAATTATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACGGATATGTCTAAAGTCTGAAAAGGCGAGCCGCCTCTTCGCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGc                                  <  2:437592/246‑1 (MQ=255)                                                                                                                                                                                          aCGCGTAACTCCGTAAGTTACCGGAGGCGCAGAGCATTGAGACCGCGCGACGGGGGGGGGGGTGTGTTATTCATTGTTATCGCCTCCG‑GGATGGTAATCACATAAATCTAAAAAAACCGCGGAAATGATAAAGATTGTTAAGCATGACAAAAGCAAGCGATTTGAATACCGGATATCTCTAAAGTCTGATAAGCCGACCCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCATCAACGCCt                                <  2:320595/249‑1 (MQ=255)                                                                                                                                                                                                  cacaGTCAGTTCGGGGAAGCGCAGGTCATGTCCAACGCGATGCGGGGGGCGGGGGGTGTTATTCATTGTTATCCGCTTCGGTG‑TGGTAATCAGATATATCAAAATAAAACTCGAAAATGATAATGATTGTTAATCATGATAAATGCAAGCGTTTGGGAGAACTGAGATGTTTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGGGAACGCCTTATCCGGCCt                      <  2:219752/251‑1 (MQ=255)                                                                                                                                                                                                                      tcggatttttGCGAACGCGTGGCGGGGGGGGGGGGGTTTTATTCATTTTTATCGCCTTCG‑TGATGGAAATCAGATATATCTAAAAAAACCTCGCAAATAAAAATGATTGTAAATCATGATAAATGCAAGCGATTTGTAGAACTGATATTTCAAGAGTCTGATAAGCCGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCa    <  2:380812/241‑1 (MQ=255)                                                                                                                                                                                                                      gcAGATCATTGCGACCGCGCGGCGGGTGGGGGGGGGTTTTATTCATTTTTATCGCCTTCG‑TGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTTTTAATCATGAAAAATGCAAGCGATTTGTAGAACTGATATGTTTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACg  <  1:101657/251‑1 (MQ=255)                                                                                                                                                                                                                              ttGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCG‑TGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCa    >  1:380812/1‑241 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                 CACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCG‑TGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACG  >  NC_000913/3216233‑3216696

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑