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	98.6%	187.2 / ‑4.7	70	coding (19/765 nt)	nfeF	NADPH‑dependent ferric‑chelate reductase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (22/47);  minor base . (0/1);  total (22/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

GTATAGTCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACG‑GAAGCGCAGATCATTGCGAACGCGCTGCGGGTAG‑CGG‑GGGG‑TGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGA‑CGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTAC  >  NC_000913/3216226‑3216714                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                      gTATAGTCCCGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACG‑GAAGCGCAGATCATTGCGAACGCGCTGCGGGTGG‑Cg                                                                                                                                                                                                                                                   <  2:192736/251‑1 (MQ=255)     aGTCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACG‑GAAGCGCAGATCATTGCGAACGCGCTGCGGGTGG‑Cgg‑gg                                                                                                                                                                                                                                               <  1:148046/250‑1 (MQ=255)          cgcgACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACG‑GAAGCGCAGATCATTGCGAACGCGCTGCGGGTGG‑CGG‑GGGG‑TGtt                                                                                                                                                                                                                                        <  1:168058/251‑1 (MQ=255)           gcgACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACG‑GAAGCGCAGATCATTGCGAACGCGCTGCGGGTGG‑CGG‑GGGG‑TGTTa                                                                                                                                                                                                                                       <  1:406456/251‑1 (MQ=255)             gACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACG‑GAAGCGCAGATCATTGCGAACGCGCTGCGGGTGG‑Cgg‑g                                                                                                                                                                                                                                                >  1:47693/1‑241 (MQ=255)             gACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACG‑GAAGCGCAGATCATTGCGAACGCGCTGCGGGTGG‑Cgg‑g                                                                                                                                                                                                                                                <  2:47693/241‑1 (MQ=255)              gcGGTGTGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGGTGCCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAACCCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTGCACG‑GAAGCGCAGAGCATTGCGAACGCGCTGCGGTTGG‑CGG‑GGGG‑TGTTATTc                                                                                                                                                                                                                                    <  1:460536/250‑1 (MQ=255)                ggTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACG‑GAAGCGCAGATCATTGCGAACGCGCTGCGGGTGG‑CGG‑GGGG‑TGTTATTCAt                                                                                                                                                                                                                                  <  1:431363/251‑1 (MQ=255)                 gTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTGCCGTTACCGTTGGCTGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACG‑GAAGCGCAGATCATTGCGAACGCGCTGAGTGTGG‑CGG‑GGGG‑TGTTATTCAtt                                                                                                                                                                                                                                 <  1:496875/251‑1 (MQ=255)                    ggCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTGTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACG‑GAAGCGCAGATGAGTGCGAACGCGCTGCGGGTGG‑CGG‑GGGG‑TGTTATTCATTGtt                                                                                                                                                                                                                              <  1:388099/251‑1 (MQ=255)                      cGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACG‑GAAGCGCAGATCATTGCGAACGCGCTGCGGGTGG‑CGG‑GGGG‑TGTTATTCATTGTTAt                                                                                                                                                                                                                            <  1:99761/251‑1 (MQ=255)                            ccTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACG‑GAAGCGCAGATCATTGCGAACGCGCTGCGGGTGG‑CGG‑GGGG‑TGTTATTCATTGTTATCGCCtt                                                                                                                                                                                                                      <  2:433414/251‑1 (MQ=255)                                ccGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACG‑GAAGCGCAGATCATTGCGAACGCGCTGCGGGTGG‑CGG‑GGGG‑TGTTATTCATTGTTATCGCCTTCGt                                                                                                                                                                                                                   <  1:151147/250‑1 (MQ=255)                                   gCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACG‑GAAGCGCAGATCATTGCGAACGCGCTGCGTGTGG‑CGG‑GGGG‑TGTTATTCATTGTTATCGCCTTCGTGATg                                                                                                                                                                                                               <  1:459833/251‑1 (MQ=255)                                        aCGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACG‑GAAGCGCAGATCATTGCGAACGCGCTGCGGGTGG‑CGG‑GGGG‑TGTTATTCATTGTTATCGCCTTCGTGATGGTAAt                                                                                                                                                                                                          <  1:245937/251‑1 (MQ=255)                                        aCGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACG‑GAAGCGCAGATCATTGCGAACGCGCTGCGGGTGG‑CGG‑GGGG‑GGGTATTCATTGTTATCGCCTTCGTGATGGTAAt                                                                                                                                                                                                          >  2:7702/1‑251 (MQ=255)                                         cGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAATAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGTGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACG‑GAAGCGCAGATCATTGCGAACGCGCTGCGGGTGG‑CGG‑GGGG‑TGTTATTCATTGTTATCGCCTTCGTGATGGTAATc                                                                                                                                                                                                         <  2:49363/251‑1 (MQ=255)                                             gCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACG‑GAAGCGCAGATCATTGCGAACGCGCTGCGGGTGG‑CGG‑GGGG‑TGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGat                                                                                                                                                                                                     <  1:412412/251‑1 (MQ=255)                                              cccTCTTCCGTTACCGTTGTCGGCACAAAGCGAGCGTCATGTTGAGGAAAGAATAGTTTGCTGTGATCGTCAAATCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGATAATGCTCTGAAAGCCGGCTCTGATGCGCTCAACGCGTAACGCAGTCAGTTCATG‑GAAGCGCAGATCATTGCGAACGCGCTGCGGGTGG‑CGG‑GGGG‑TGTTATTCATTGTTAGCGCCTACGTGATGGTAATCAGata                                                                                                                                                                                                    <  1:216178/251‑1 (MQ=255)                                              cccTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACG‑GAAGCGCAGATCATTGCGAACGCGCTGCGGGTGG‑CGG‑GGGG‑TGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGata                                                                                                                                                                                                    <  1:7702/251‑1 (MQ=255)                                                    tCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACG‑GAAGCGCAGATCATTGCGAACGCGCTGCGGGTGG‑CGT‑GGGG‑TGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCt                                                                                                                                                                                               <  1:221486/250‑1 (MQ=255)                                                     ccGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACG‑GAAGCGCAGATCATTGCGAACGCGCTGCGGGTGG‑CGG‑GGGG‑TGTTATTCAtt                                                                                                                                                                                                                                 <  2:74675/215‑1 (MQ=255)                                                     ccGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACG‑GAAGCGCAGATCATTGCGAACGCGCTGCGGGTGG‑CGG‑GGGG‑TGTTATTCAtt                                                                                                                                                                                                                                 >  1:74675/1‑215 (MQ=255)                                                          aCCGTTGGCGGCACAACGTGAGCGTCAGGTTGAGGAAAGAAGCGTTTGCTGTGATCGTCAAAGCCACGCGCTGTCAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACG‑GAAGCGCAGATCATTGCGAACGCGCTGCGGGTGG‑CGG‑GGGG‑TGTTATTCCTTGTTATCGCCTTCGTGTTGGTAATCAGATTTATCtaaataa                                                                                                                                                                                         >  2:61584/1‑250 (MQ=255)                                                                    gcACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACG‑GAAGCGCAGATCATTGCGAACGCGCTGCGGGTGT‑CGG‑GGGG‑TGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAAGAAAACTCGCAaat                                                                                                                                                                              <  2:572196/251‑1 (MQ=255)                                                                           gTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACG‑GAAGCGCAGATCATTGCGAACGCGCTGCGGGTGG‑CGG‑GGGG‑TGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAaatgataatg                                                                                                                                                                       <  2:351880/251‑1 (MQ=255)                                                                                 gTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACG‑GAAGCGCAGATCATTGCGAACGCGCTGCGGGTGG‑CGG‑GGGG‑TGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGtt                                                                                                                                                                 >  1:359491/1‑251 (MQ=255)                                                                                    aggttgccgaaagaaaAGTTTGCTGTGCACGTCAAAGCCACGCGCTGTAACACCGTCCAGCGCGTCGTCGCCGAGGACAATGCGATGAAAACAGACGCTGCTGCGCCCAAAGCGTCACACAGTCAGTTCACG‑GAAGCGCAGATCATTGCGAACGCGCTGCGGGTGG‑CGG‑GGGG‑GGCTATTCATTGTTATCGCCTTCGTGCTGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAAt                                                                                                                                                              >  2:222967/9‑251 (MQ=255)                                                                                         gAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACG‑GAAGCGCAGATCATTGCGAACGCGCTGCGGGTGG‑CGG‑GGGG‑CGGTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATgt                                                                                                                                                         >  2:138524/1‑250 (MQ=255)                                                                                          aGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACG‑GAAGCGCAGATCATTGCGAACGCGCTGCGGGTGG‑CGG‑GGGG‑TGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGAt                                                                                                                                                        <  2:418190/251‑1 (MQ=255)                                                                                            gAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGTCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACG‑GAAGCGCAGATCATTGCGAACGCGCTGCGGGTGG‑CGG‑GGGG‑TGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATaa                                                                                                                                                      <  2:357199/251‑1 (MQ=255)                                                                                               agaagaGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAACCCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACG‑GAAGCGCAGATCATTGCGAACGCGCTGCGGGTGG‑CGG‑GGGG‑TGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAACTAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATg                                                                                                                                                   >  1:355923/1‑251 (MQ=255)                                                                                                        tGCTGTGATCATCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACG‑GAAGCGCAGATCATTGCGAACGCGCTGCGGGTGG‑CGG‑GGGG‑TGTTATTCATTGTTATCTCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCa                                                                                                                                                                                 <  1:91497/212‑1 (MQ=255)                                                                                                        tGCTGTGATCATCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACG‑GAAGCGCAGATCATTGCGAACGCGCTGCGGGTGG‑CGG‑GGGG‑TGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATATAAATAAAACTCGCa                                                                                                                                                                                 >  2:91497/1‑212 (MQ=255)                                                                                                            gtgATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACG‑GAAGCGCAGATCATTGCGAACGCGCTGCGGGTGG‑CGG‑GGGG‑TGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGt                                                                                                                                       >  1:18228/1‑250 (MQ=255)                                                                                                               aTCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGGGGTCAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTACCACAGTCAGTTCACG‑GCAGCGCAGGTCATTGAGAAAGCGATGCGGGTGG‑CGG‑GGGG‑GGTTATTCATTGTTAACGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAaatgataatgat                                                                                                                                                                     >  1:529282/1‑217 (MQ=255)                                                                                                               aTCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGGCAGTGCACG‑GTAGCGCAGATCATTGCGAACGCGCTGCGGTTGG‑CGGTGTGT‑TGTT‑GTCAGTGTTATCGCCATCGGTATTGTAATCAGATACATCTAAATAAATCTCGCAaatgataatgat                                                                                                                                                                     <  2:529282/217‑1 (MQ=255)                                                                                                                  gTCAAAGACACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACG‑GAAGCGCAGATCATTGCGAACGCGCTGCGGGTGG‑CGG‑GGGG‑TGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTg                                                                                                                                >  1:469826/1‑251 (MQ=255)                                                                                                                    cAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACG‑GAAGCGCAGATCATTGCGAACGCGCTGCGGGTGG‑CGG‑GGGG‑TGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGat                                                                                                                              <  2:190569/251‑1 (MQ=255)                                                                                                                    cAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACG‑GAAGCGCAGATCATTGCGAACGCGCTGCGGGTGG‑CGG‑GGGG‑TGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAAGGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGa                                                                                                                               <  2:562920/250‑1 (MQ=255)                                                                                                                    aaaaGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACG‑GAAGCGCAGATCATTGCGAACGCGCTGCGGGTGG‑CGG‑GGGG‑TGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGat                                                                                                                              <  2:445928/250‑1 (MQ=255)                                                                                                                      aaGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACG‑GAAGCGCAGATCATTGCGAACGCGCTGCGGGTGG‑CGG‑GGGG‑TGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGatat                                                                                                                            >  1:284891/1‑251 (MQ=255)                                                                                                                       aGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGATGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACG‑GAAGCGCAGATCATTGCGAACGCGCTGCGTGTGG‑CTG‑GGTG‑TGGTAATCATTGTTATCGCCTTCGTGATGGTAATCAGATATAACTAGATATAACGCGCTAATGATAATGATTGGGAAGCAGAATAAGTTCAAGCGATTTGTAGAGCTGATATg                                                                                                                           <  2:217254/251‑1 (MQ=255)                                                                                                                              cgATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACG‑GAAGCGCAGATCATTGCGAACGCGCTGCGGGTGG‑CGG‑GGGG‑TGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAg                                                                                                                    <  1:225054/251‑1 (MQ=255)                                                                                                                                    aaaacctcccAGCCCCTCGCCGCCGGGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAAGACAGTCAGTTCACG‑GAAGCGCAGATCATTGCGAACGCGCTGCGGGTGG‑CGG‑GGGG‑TGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATTATTGATAATGATGATAAATGCAAGCGATTTGTAGAACAGAGATGTCTATAGTCTGAt                                                                                                              <  2:338648/243‑1 (MQ=255)                                                                                                                                      aaCCGTCCAGCGCCTCGCCGCCGAGTACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACG‑GAAGCGCAGATCATTGCGAACGCGCTGCGGGTGG‑CGG‑GGGG‑TGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACGGATATGTCTATAGTCTGATaa                                                                                                            <  2:107002/251‑1 (MQ=255)                                                                                                                                         cGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACG‑GAAGCGCAGATCATTGCGAACGCGCTGCGGGTGG‑CGG‑GGGG‑TGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGA‑c                                                                                                        <  1:326194/251‑1 (MQ=255)                                                                                                                                          gTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACG‑GAAGCGCAGATCATTGCGAACGCGCTGCGGGTGG‑CGG‑GGGG‑TGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGA‑Cg                                                                                                       >  1:171197/1‑251 (MQ=255)                                                                                                                                             cAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACG‑GAAGCGCAGATCATTGCGAACGCGCTGCGGGTGG‑CGG‑GGGG‑TGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACGCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGA‑CGAAc                                                                                                    <  1:571468/251‑1 (MQ=255)                                                                                                                                              aGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCCCAACGCGTAACACAGTCAGTTCACG‑GAAGCGCAGATCATTGCGAACGCGCTGCGGGTGG‑CGG‑GGGG‑TGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGA‑CGAAcc                                                                                                   >  2:491809/1‑251 (MQ=255)                                                                                                                                                 gcCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACG‑GAAGCGCAGATCATTGCGAACGCGCTGCGGGTGG‑CGG‑GGGG‑TGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAaatgataat                                                                                                                                                                        >  1:41099/1‑180 (MQ=255)                                                                                                                                                 gcCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACG‑GAAGCGCAGATCATTGCGAACGCGCTGCGGGTGG‑CGG‑GGGG‑TGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAaatgataat                                                                                                                                                                        <  2:41099/180‑1 (MQ=255)                                                                                                                                                             gggccAATGCGCTGAAACCCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACG‑GAAGCGCAGATCATTGCGAACGCGCTGCGGGGGG‑CGG‑GGTG‑TGTTATTCATTGTTATCGCTTTCGTGAGGGTAATCAGATATATCAAAAGAAAACTCGCAAATGAGAATGATTGTTAATCATGATAAATGCAAGCGATTTGTATATCAGATGGGTCTATAGGCTGATAAGG‑CGAACCGCCTCTTCTCAGGca                                                                                    <  1:61584/247‑1 (MQ=255)                                                                                                                                                                 cAATGCGCTAAAACCCGGCGCGGATGCGCTCAACGGGTAACACAGCCAGTTC‑CGTGAACCGCAAATCATGGCGAACGCG‑TGGGGGTGGGCGG‑GGGT‑TGTTATTTAGTGTTATCGCCTCCGTGATGGTAATCAGATATATCTAAAAAAAACTCGAAAAGGATGATGATGGTGAATCATGAAAAATGCACGCTGTTTGTAGACCTAATATGTTTAGTGTTTAAAAAGA‑CGAACCGCCTCTTCTCAGGcatcat                                                                                <  2:469826/251‑1 (MQ=255)                                                                                                                                                                  aaTGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACCCAGTCAGTTCACG‑GAAGCGCAGATCATTGCGAACGCGCTGCGGGTGG‑CGG‑GGGG‑TGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGa                                                                                                          >  1:565931/1‑225 (MQ=255)                                                                                                                                                                  aaTGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACG‑GAAGCGCAGATCATTGCGAACGCGCTGCGGGTGG‑CGG‑GGGG‑TGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATGTGTAGAACTGATATGTCTATAGTCTGATAAGa                                                                                                          <  2:565931/225‑1 (MQ=255)                                                                                                                                                                     gcgcTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACG‑GAAGCGCAGATCATTGCGAACGCGCTGCGGGTGG‑CGG‑GGGG‑TGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAAACATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGA‑CGAACCGCCTCTTCTCAGGCATCATTACt                                                                            >  1:412049/1‑251 (MQ=255)                                                                                                                                                                          gAAAACCGGCGCTGATGCGCCCAACGCGTAACACAGTCAGTTCACG‑GAAGCGAAGATCATTGCGAGCGCGCTGCGGGTGG‑CGG‑GGGG‑TGTTATTCATTGTTATCGCCTGCGTGATGGTAATCAGATATATCTAAATAAACCCCGCAAATTATAAGGATTGTTAATCATGAAAAATGCAAGCGATTTGTAGAACGAATATTTTGATAGTCTGATAAGA‑CGAACCGCCTCTTCTCAGGCATCATTACTCAAc                                                                        <  2:478241/250‑1 (MQ=255)                                                                                                                                                                              cccGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACG‑GAAGCGCAGATCATTGCGACCGCGCTGCGGGTGG‑CGG‑GGGG‑TTTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGTTTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGA‑CGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCgg                                                                   <  2:171197/250‑1 (MQ=255)                                                                                                                                                                                    gcTGATGCGCTCAACGCGTAACACAGTCAGTTCACG‑GAAGCGCAGATCATTGCGAACGCGCTGCGGGTGG‑CGG‑GGGG‑TGTTATTCATTGTTATCGCATTCGTGATGGTATTCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGAAGACCTGTTATGTTTATATTCTGATAAGA‑CGAACCGCCTCTTCTCAGGCATCATTACTCAACGCGGGATgcgg                                                             >  1:40722/1‑251 (MQ=255)                                                                                                                                                                                                aacccatAACACAGTCAGTCCACG‑GAGGCGCAGATCATTGAGAACGCGCGGGGGGTGG‑GGG‑GGGG‑TTTTATTCATTGTTATCCACTTCGTGAGGGTAACCAGATATATTTAAATAAAACTCGAAAATGAAAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATAGGTCTATAGTCTAATAAGA‑CGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCtt                                                 <  1:153558/245‑1 (MQ=255)                                                                                                                                                                                                  cgcgTACCACAGTCAGTTCACG‑GAAGCGCAGATCATTGCGACCGCGCTGGGGGT‑G‑CGG‑GGGGGTGTTATTCATTGTTATCGCTTTCGGGATGGTAATCAAATATATCTAAAAAAAACTCGCAAATGATAATGTTTGTTAATCATGAAAAATGCAAGCGATTTTAAAAACTGATATGTCTATAGTCTGA‑AAGAACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTAt                                               <  2:294410/251‑1 (MQ=255)                                                                                                                                                                                                        acacAGTCAGTTCACG‑GAAGCGCAGATCATTGCGAACGCGCTGCGGGTGG‑CGG‑GGGG‑TGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGA‑CGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                         >  2:443107/1‑251 (MQ=255)                                                                                                                                                                                                            agtcagtTCACG‑GAAGCGCAGATCATTGCGAACGCGCTGCGGGTGG‑CGG‑GGGG‑TGTTATTCATTGTTATCGCATTCGTGATGGTAAACAGACATATCTAAATAAAACTCGCAAATGATAATGATTGCTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGA‑CGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCtccg                                     >  1:360182/1‑248 (MQ=255)                                                                                                                                                                                                                ggtTCACG‑GAAGCGCAGATCATGGCGAACGCGCTGGGGGTGG‑CGG‑GGGG‑TGTTATTCTTTGTTATCGCCTGCTTGATGGTAATAAGATATATCAAAAAAAAACTCGCAAATGATAATGATTTTTAAACATGATAAATGAAAGCGATTTGTAGAACTAATATGTCTATAGTCTGAAAAGA‑CGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGtgtg                                 <  2:412049/250‑1 (MQ=255)                                                                                                                                                                                                                         gaagcacaattaaTTGCGAACGGGCTGGGGGTGG‑CGG‑GGGG‑TTTTATTCATTTTTGTCGCTTTGGTGTGGGTAATCAAATATATCAAAAAAAACCTCGCAAATAATAATGATTGTTAATCAGGAGAAATGCAAGCGATTTGTAGAACTGATATTTCTAGAGTTTGATAAGC‑CAAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGt                         <  2:195980/239‑1 (MQ=255)                                                                                                                                                                                                                                       ggcgaaggcgtggCGGGTGG‑CGG‑GGGG‑GTTTATTCAGTGTTAGCGCCT‑‑G‑GATGGTAATCAGATATATCTAAAAAAACCTCGCAAATGAAAATGATTTTTAATCATGAAAAAGGCAGGCGATTTGTAGAACGGATATGTCTATAGTCGTATAAGA‑CGAGCCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTAGCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCAt        <  2:284593/239‑1 (MQ=255)                                                                                                                                                                                                                                        gCGAACGCGCGGCGGGTGG‑CGG‑GGGG‑TGTTTTTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAAAAAAACTCGAAAATGATAAGGTTTTTAAACCGTGATAAATGCAAGCGATTTGTAGAACTGAAATGTCTATAGTCGGATAAGA‑CGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAAt           <  1:188950/250‑1 (MQ=255)                                                                                                                                                                                                                                           aacccgcTGCGGGTGG‑GGG‑GTGG‑TGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGTTATATCTAAATAAAACTCGCAAATAATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGA‑CGATCCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCGTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCAtt       <  1:93150/247‑1 (MQ=255)                                                                                                                                                                                                                                                gcTGCGGGTGG‑CGG‑GGGG‑TGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGA‑CGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGCTGAGTCCACGTTTAAAAATCATTTTTac  >  2:151731/1‑251 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                      GTATAGTCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACG‑GAAGCGCAGATCATTGCGAACGCGCTGCGGGTAG‑CGG‑GGGG‑TGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGA‑CGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTAC  >  NC_000913/3216226‑3216714

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