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

ATAGTCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGG‑‑G‑T‑G‑TTATTCATTGTTATCGCCTTCGTGA‑TGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTT  >  NC_000913/3216228‑3216712                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      ataGTCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTGCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCggg                                                                                                                                                                                                                                                  <  1:432864/251‑1 (MQ=255)       aCGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCggggg‑‑ggtg                                                                                                                                                                                                                                          >  1:467805/1‑248 (MQ=255)         gcgACGGGGGGCGTGGTCCTGCCGGCCAGACGATGCCCTCTGCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACTATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGGCAGTTCACGGAAGCGCAGAGCATTTCGAACGCGCTGCGGGTGGCggggg‑‑ggtgtt                                                                                                                                                                                                                                        <  2:191421/251‑6 (MQ=255)           gACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTAt                                                                                                                                                                                                                                    <  1:10110/251‑1 (MQ=255)             gggTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTc                                                                                                                                                                                                                                  <  1:367253/250‑1 (MQ=255)              ggTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACAGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTc                                                                                                                                                                                                                                  >  1:44287/1‑250 (MQ=255)               gTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCAt                                                                                                                                                                                                                                <  2:391943/251‑1 (MQ=255)                 gggCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTGCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTg                                                                                                                                                                                                                              <  2:156450/251‑1 (MQ=255)                   gCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCATTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGtt                                                                                                                                                                                                                            <  2:341303/251‑1 (MQ=255)                              ccgcgcAGACGACGCCCTAATCCGTTACGGATGGCGTCACAAAGTTAGCGTAAGGGTGAGGAAAGAAGAGTTTGCTGTGATCGACAGAGCCACGCGATGTAGAGCCGGCCAGCGCCGCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAACGCAGATCATTGCGAACGGGCTGCGTGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTct                                                                                                                                                                                                                  <  2:259221/245‑2 (MQ=255)                              ccGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAACGTGAGCGTCCGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTGACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGt                                                                                                                                                                                                                 >  2:60343/1‑251 (MQ=255)                                ggCCAGACGATGCCCTCTTCCGTGACCGTTGGTGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGTGa                                                                                                                                                                                                               <  1:13981/251‑1 (MQ=255)                                   cAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGTGA‑Tgg                                                                                                                                                                                                           >  1:395979/1‑251 (MQ=255)                                    aGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGTGA‑TGGt                                                                                                                                                                                                          >  1:195809/1‑251 (MQ=255)                                      aCGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACCGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGTGC‑TGGTa                                                                                                                                                                                                         >  2:220332/1‑250 (MQ=255)                                      aCGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGTGA‑TGGTaa                                                                                                                                                                                                        <  1:21426/251‑1 (MQ=255)                                        gattcccTCTGCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGAAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGTCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGTGA‑TGGTAATc                                                                                                                                                                                                      <  1:80292/247‑1 (MQ=255)                                                    cGTTACCGTTGGCGGCACAAAGGGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCGCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGGGA‑TTGTAATCAGATATATCtaa                                                                                                                                                                                          <  2:406586/251‑1 (MQ=255)                                                      ttACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGACGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGA‑G‑TTATTCATTGTTATCGCCTTCGTGA‑TGGTAATCAGATATATCtaaat                                                                                                                                                                                        >  2:213231/1‑251 (MQ=255)                                                            ttGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGATTGCTGTGATCGTCAAAGCCACGCGATGTAACACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGTGA‑TGGTAATCAGATATATCTAAATAAAACt                                                                                                                                                                                  >  1:64612/1‑251 (MQ=255)                                                                  gcACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGTGA‑TGGTAATCAGATATATCTAAATAAAACTCGCaaa                                                                                                                                                                            <  1:413103/251‑1 (MQ=255)                                                                      aaaGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCAtt                                                                                                                                                                                                                               <  1:366497/197‑1 (MQ=255)                                                                      aaaGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCAtt                                                                                                                                                                                                                               >  2:366497/1‑197 (MQ=255)                                                                             gCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCTCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGTTA‑TGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGAtt                                                                                                                                                                 <  2:411986/251‑1 (MQ=255)                                                                             gCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGGAACACAGTCAGTTCACGGAGGCGCAGATCATTGCGAACGCGCGGCTGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCGTTCGTGA‑TGGTAATCAGATAGATCTAAATAGAACTCGCAAATGATAATGAtt                                                                                                                                                                 <  2:480551/251‑1 (MQ=255)                                                                                  aGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGTGA‑TTGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTaa                                                                                                                                                            <  1:35978/251‑1 (MQ=255)                                                                                    gTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGTGA‑TGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATc                                                                                                                                                          <  1:146489/251‑1 (MQ=255)                                                                                        aGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGTGA‑TGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATg                                                                                                                                                       >  2:78878/1‑250 (MQ=255)                                                                                                  agTTTGCTGTGATCGTCAAAGCCCCGCGATGTAACACCGTCCAGCGCCTCGCAGCCGAGGACAATGCGCTGAAAGCCGGCGCTGATGCGCTTAACGCGTAACACAGTCAGTTCACGGAAGCGCCGATCATTGCGAACGCGCTGCGGGGGGCGGGGG‑‑G‑G‑GGTTCTTCATTGTTATCTCCTTCGTGC‑TGGTAATCAGCTATATCTCAATAAAACTCGCAAATGATAATGTTTGTTAATCATGCTAAATGCa                                                                                                                                              >  2:312937/1‑249 (MQ=255)                                                                                                       gCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGTGA‑TGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGAtt                                                                                                                                       <  2:361959/251‑1 (MQ=255)                                                                                                          gtgATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTAGTCATTTTTATCGCCTTCGTGA‑TGGTAATCAGATATATCTAAATAAAACTCGCAAAGTATAATGATTGTTAATCATGATAAATGCAAGCGAttt                                                                                                                                      <  2:478953/249‑1 (MQ=255)                                                                                                                       ccACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGTGA‑TGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGt                                                                                                                       <  1:281521/251‑1 (MQ=255)                                                                                                                        cACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGCACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGTGA‑TGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTACTCATGATAAATGCAAGCGATTTGTAGAACTGATATGTc                                                                                                                      >  1:469963/1‑251 (MQ=255)                                                                                                                                  aaaaCCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGTGA‑TGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTg                                                                                                             <  1:175566/250‑1 (MQ=255)                                                                                                                                  aaaaCCGCCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAGCCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGGGA‑TGGTAATCAGATATATCTAAAGAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTg                                                                                                             <  1:312937/250‑1 (MQ=255)                                                                                                                                   aaaCCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGGAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGTGA‑TGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGAGTTGTAGAACTGATATGTCTATAGTCTGAt                                                                                                           <  2:389975/251‑1 (MQ=255)                                                                                                                                   aaaCCGGCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGTGA‑TGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGAt                                                                                                           <  2:152093/251‑1 (MQ=255)                                                                                                                                      ccGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGT‑G‑T‑G‑TTATTCATTGTTATCGCCTTCGTGA‑TGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAg                                                                                                        <  1:202902/251‑1 (MQ=255)                                                                                                                                          ccAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCGGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTG‑T‑G‑TTATTCATTGTTATCGCCTTCGTGA‑TGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACg                                                                                                     <  1:479434/251‑1 (MQ=255)                                                                                                                                           cAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGTGA‑TGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGa                                                                                                    >  2:321629/1‑250 (MQ=255)                                                                                                                                            aGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGTGA‑TGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAAc                                                                                                  <  2:467805/251‑1 (MQ=255)                                                                                                                                             gcgcCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGTGA‑TGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAATCGATTTGGAGAACTGATATGTCTATAGGCTGATAAGACGAAcc                                                                                                 <  2:64612/251‑1 (MQ=255)                                                                                                                                              cgcCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGTGA‑TGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCg                                                                                                <  1:392715/251‑1 (MQ=255)                                                                                                                                                    gccgccGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGTGA‑TGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCtctt                                                                                          >  2:283929/1‑251 (MQ=255)                                                                                                                                                           aGGACAATGCGCTGAACACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGTGA‑TGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGc                                                                                   >  1:407034/1‑251 (MQ=255)                                                                                                                                                                        gAAACCCGGCTCTGATGCGCTAAACGGGTAACACAGTCAGTTCACGGAGGCGCAAATCATTGCGACCGCGCTGGGGGGGGGGGGGG‑‑GGT‑G‑TTATTCATTGTTACCGCCTTCGTGA‑TGGTAATCAGATATATCTAAATAAAACGCGAAAATGATAATGTGTGTTAATCAGAATAAAGGCAAGCGATTTGTAGAACTGATATTTTTATGTTCTAATAAGACGAGCCGCCTCTTCTCAGGCATCATTACTCa                                                                        <  2:125673/249‑1 (MQ=255)                                                                                                                                                                                  gcTGATGCGCTCAACGCGAAACACAGTCAGTTCACGGAAGCGCAGATAATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGTGA‑TGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCACGCGATTTGTAGAACTGAAATGTCCATAGTCTGATAAGACGAACCGCCTCTTATCAGGCATCATTACTCAACCCCggcagcg                                                            >  2:326595/1‑246 (MQ=255)                                                                                                                                                                                                          agtcattTCGCGGAAGCGCAGATCATGGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTTTTCATTGTTATCGCTTTC‑TGAATGGTAATCAGAAATATCTAAAAAAACCTCGCAAATGATAATGATTGTAAAGCATGATAAATGCAACCGTTTTGTAGAACTGAAATGTTTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTAc                                    <  1:371084/251‑1 (MQ=255)                                                                                                                                                                                                              ggtTCACGGAAGCGCAGTTTATTGCGATCGCGCTGGGGGGGGGGGGGG‑‑GGT‑G‑TTATTCTTTGTTATCGCCTTCGTTA‑TGGTAATCAGATATATCTAAAAAAAACGCGCAAATGAAAATGATTGTTAATCAGGATAAATGCAAGCGATTGGTAGAACTGATATGTGTATAGTCGGATAAGACGAACCGCCTCTTCTCAGGCATCAGTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACgtgt                                <  2:407034/250‑1 (MQ=255)                                                                                                                                                                                                              agtTCACGGAAGCGCAGATCATGGCGAACGCGCTGCGGGGGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGGGA‑TGGTAATCAGATATATCTAAAAAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGAAAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACgtgt                                <  1:213231/251‑1 (MQ=255)                                                                                                                                                                                                                   ccgaaaGCTCAGATCATTGCGACCCCGCTGCGGGTGGCGGGGG‑‑G‑TTG‑TTATTCATGGTTAGCGCCTTCGTAA‑TGGTAATCAAATATATCTAAATAAACCCCGCAAATGAAAATGATTGTTAATCATGATAAATCCAAGCGGTTGGTAGACCTAATATTTCTATAGTCGGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGAt                           <  1:419865/247‑1 (MQ=255)                                                                                                                                                                                                                   aCGGAAGCTCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGTGA‑TGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGTGAt                           >  2:380169/1‑251 (MQ=255)                                                                                                                                                                                                                           gcAGATGATTGCGAACGCGCGGCGGGTGGCGGGGG‑‑GGT‑T‑TTATTCATTGTTATCGCCTTCGTGA‑GGTTAATAAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCACGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCAc                   <  2:86126/251‑1 (MQ=255)                                                                                                                                                                                                                              gATCATTGCGAACGCGCTGCGGGTGGGGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGTGA‑TGGTAATCAGATATATCAAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAGGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGcatcat                                                                              <  1:74954/189‑1 (MQ=255)                                                                                                                                                                                                                              gATCATTGCGAACGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGTGA‑TGGTAATCAGATATATCTCAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGcatcat                                                                              >  2:74954/1‑189 (MQ=255)                                                                                                                                                                                                                               aTCATTGCGACCGCGCTGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTGTTATCGCCTTCGTGA‑TGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGttt               <  1:92089/251‑1 (MQ=255)                                                                                                                                                                                                                               aTCATTGCGACCGCGCGGCGGGTGGCGGGGG‑‑GGT‑G‑TTATTCATTTTTATCGCCTTCGTGA‑TGGTAACCAGATATATCTAAAAAAAACTCGCAAATGATAATGATTGTTAATCATGAAAAATGCAAGCGATTTGTAAAACTAATTTGTCTATAGTCGGATAAGACAAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGtt                <  1:254944/250‑1 (MQ=255)                                                                                                                                                                                                                                             gcTGCGGGTGGCGGGGG‑‑GGT‑G‑TTCTTTATTGTTATCGCCTTCGTGA‑TGGTAATCAGATATATATAAATAAAACTCGCAACTGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGCTAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGAGAACGCCTTATCCGGCCTACGGGTGAGGTGAGTCCACGTTTAAAAATCAttttt  >  2:306438/1‑250 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      ATAGTCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGG‑‑G‑T‑G‑TTATTCATTGTTATCGCCTTCGTGA‑TGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTT  >  NC_000913/3216228‑3216712

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