Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	172,569	C→T	R370C (CGT→TGT)	fhuB →	iron(III) hydroxamate ABC transporter membrane subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	172,569	0	C	T	100.0%	238.7 / NA	70	R370C (CGT→TGT)	fhuB	iron(III) hydroxamate ABC transporter membrane subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (33/37);  total (33/37)

TTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGATGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGG‑TGGCGGGCTGTATTATTCAGCGACTGA‑CCGGAAACCCGATGGCAAGCCCGG‑AAGTGCTGGGGATTAGCTCCGGCGCGGCGTTTGGCGTGGTGTTGATGCTGTTTCTGGTGCCGGGTAATGCCTTTGGCTGGC  >  NC_000913/172330‑172809                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    tttCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGATGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTTGTGATGCGCAc                                                                                                                                                                                                                                          >  1:340113/1‑251 (MQ=255)      gATCAAATCATCCTCTGGCTGACTCGCGTGTGGATGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTTGTGATGCGCACGGCTg                                                                                                                                                                                                                                     <  2:162730/251‑1 (MQ=255)        tCAAATCATCCTCTGGCTGACTCGCGTGTGGATGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTTGTGATGCGCACGGCTgg                                                                                                                                                                                                                                    >  1:219008/1‑250 (MQ=255)              catcCTCTGGCTGACTCGCGTTTGATTGGAAGTGTCCCACGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGGTGCCGCGTTTACGCCGCATTAGCGCGCCGGATATGAAGGGCAACGATCATGTCGCGGCTGAACGCCAACATGTGCTGGCGTTGGCCCTCGCGGGCGGCGTGCTGCTGTTGTTGGCTGTGGTGGTGGCGCTATCGTTTGGTTGTGATGGGCACGGCTGGACGTggg                                                                                                                                                                                                                             >  2:161455/1‑251 (MQ=255)              catcCTCTGGCTGACGCGCGTGTGGATGGAAGGGGCCTCCGGTTCGGTCACTGCGTTGAGCGGTGCGCAGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAGCGCCAACTTGTGCTGGGGTTTGCTCTCGCGGGCGGCGTGCGGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGGTTGGTTGTGATGCGCAGGGCGGGACGTggg                                                                                                                                                                                                                             <  2:157107/251‑1 (MQ=255)                           aCTCGCGTGTGGATGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTTGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGtt                                                                                                                                                                                                                >  1:75180/1‑251 (MQ=255)                                ggtgtGGATGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCAGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTTGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCg                                                                                                                                                                                                           <  2:261909/250‑1 (MQ=255)                                cgTGTGGATGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTTGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCg                                                                                                                                                                                                           >  2:42422/1‑251 (MQ=255)                                       atggaAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTTGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCGCGAGGAttt                                                                                                                                                                                                    <  1:110328/251‑1 (MQ=255)                                        gggaAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGGGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGAGATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTTGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTa                                                                                                                                                                                                   <  2:333872/250‑1 (MQ=255)                                                tCCACCGGTTCGGTCACTGCGTTGATCGGAGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGGGTCGCGGCTGAACGCCAACATGTGCTGGCGCTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGGGGTGGTGGCGCTGGCGTTTGGTTGTGATGCGCACGGCTGGAGGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAAGGCCCTg                                                                                                                                                                                           >  2:162875/1‑251 (MQ=255)                                                      ggTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGTCGGATATGAAGGGCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGTTTGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTg                                                                                                                                                                                     <  1:163246/251‑1 (MQ=255)                                                      ggTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTTGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTg                                                                                                                                                                                     <  1:45864/251‑1 (MQ=255)                                                          cGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTTGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGcc                                                                                                                                                                                  <  2:90379/250‑1 (MQ=255)                                                          cGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTTGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGcc                                                                                                                                                                                  <  1:41944/250‑1 (MQ=255)                                                            gTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTTGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCcgcg                                                                                                                                                                               >  1:226259/1‑251 (MQ=255)                                                                      tGATCGGAGCGCCGCTGCTACTGCGGCTGTTGCCGCGTTTACGCAGCCTTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGGGGCGCTGTCGTTTGGTTGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCATGGCGCTGGCCGCGAATTATggcg                                                                                                                                                                     >  2:147181/1‑251 (MQ=255)                                                                       gATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTTGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATggcg                                                                                                                                                                     <  1:26704/250‑1 (MQ=255)                                                                           ggTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTTGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATggcggggcg                                                                                                                                                                >  2:168376/1‑247 (MQ=255)                                                                           ggTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTTGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGgcgcc                                                                                                                                                                >  1:326656/1‑250 (MQ=255)                                                                              gcgcCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGATTGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGt                                                                                                                                                              <  2:75180/250‑1 (MQ=255)                                                                                gcggctgctGCGGTGGCTGTTGACGCTTTTACGCGTCATTAGCGCGCCGGTGATGAGGGTCAACGTTGGTGTCGCGGCTGAACGCCAACATGTGGTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTTGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCGCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCTCTGTTTg                                                                                                                                                           <  2:37975/248‑1 (MQ=255)                                                                                gccgcTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTTGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTg                                                                                                                                                           <  2:219008/251‑1 (MQ=255)                                                                                 ccgcTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTTGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGc                                                                                                                                                          >  1:369622/1‑251 (MQ=255)                                                                                  cgcTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCCACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTTGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGGTTTAATGCCCTGGCGCTGGCCGCGAATTATggcggggcggttgtcg                                                                                                                                                         >  2:14636/1‑240 (MQ=255)                                                                                  cgcTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTTGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGc                                                                                                                                                          <  1:235867/250‑1 (MQ=255)                                                                                  cgcTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTTGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGGTTGCg                                                                                                                                                         >  2:235867/1‑251 (MQ=255)                                                                                    atgctACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGCAGGTCAACGATCGTGTCGAGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTTGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGGTTTAATGCCCTGGCGCTGGCCGCGAATTATggcgggcgtttttgtggg                                                                                                                                                       >  2:192831/2‑238 (MQ=255)                                                                                      gctACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCATCATGTGCTGGCGTGTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTTGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCggg                                                                                                                                                       <  2:22818/249‑1 (MQ=255)                                                                                        tACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTTGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGGCGCGAATTATGGCGGCGCTGTTTGCGGGCGTc                                                                                                                                                   >  1:166190/1‑251 (MQ=255)                                                                                         aCTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTTGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCCCTGTTTGGGGGGGTcc                                                                                                                                                  >  2:178513/1‑250 (MQ=255)                                                                                              ggttGTTGCCGCGTTTACGCGGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTTAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTTGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCtg                                                                                                                                             <  2:90111/248‑1 (MQ=255)                                                                                                  gTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTTGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCtggcg                                                                                                                                          <  1:339044/250‑1 (MQ=255)                                                                                                       cgcgTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGAGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTTGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCCCTGTTTTCGGGCGTCATGCtggcgg‑gggcg                                                                                                                                   >  1:2502/1‑251 (MQ=255)                                                                                                          gTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTTGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTc                                                                                                                                                                                                            >  2:125243/1‑176 (MQ=255)                                                                                                          gTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTTGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTc                                                                                                                                                                                                            <  1:125243/176‑1 (MQ=255)                                                                                                              aCGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTTGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGG‑TGGCGGGCTGta                                                                                                                            <  1:168994/251‑1 (MQ=255)                                                                                                                 cagcaTTAGCGCGCCGGCTATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTTGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCTCGGATTATTGCGGCGCGGTGTGCGGGCGTCCTTCTGGGGGATGGCgggtggagtt                                                                                                                          >  2:262756/1‑244 (MQ=255)                                                                                                                    caTTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTTGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGGGGCCTGCtggcgggggggggggggggtattt                                                                                                                       >  2:121334/1‑233 (MQ=255)                                                                                                                      ttAGCGCGCCGGCTATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCCACATGTGCTGGAGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTTGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTTTGGGGGGCATGCtggcgggtgagggatcgtgttatttt                                                                                                                     >  2:228687/1‑231 (MQ=255)                                                                                                                        aGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTTGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGG‑TGGCGGGCTGTATTATTCAGc                                                                                                                   <  1:155232/250‑1 (MQ=255)                                                                                                                        aGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTTGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGG‑GGGGGGGcggtcttttttagcg                                                                                                                  >  1:149228/1‑237 (MQ=255)                                                                                                                              ccGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTTGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGG‑TGGCGGGGCGTATTTTTTAGCGGCTGA‑c                                                                                                           >  1:139594/1‑251 (MQ=255)                                                                                                                                      gAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTTGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGG‑TGGCGGGCTGTATTATTCAGCGACTGA‑CCGGAAAcc                                                                                                   <  1:288082/251‑1 (MQ=255)                                                                                                                                       aaGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTTGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGG‑TGGCGGGCTGTATTATTCAGCGACTGA‑CCGGAAAccc                                                                                                  <  1:216860/251‑1 (MQ=255)                                                                                                                                           tCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTTGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGG‑GGGGGGGGTGGTTTATTTAGCGGCTTA‑ACGGGAACCCGATg                                                                                              >  1:95430/1‑251 (MQ=255)                                                                                                                                                      gtCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTTGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGG‑TGGCGGGCTGTATTATTCAGCGACTGA‑CCGGAAACCCGATGGCAAGCCCgg                                                                                    <  1:178513/250‑1 (MQ=255)                                                                                                                                                         gcgGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTTGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGG‑TGGCGGGCTGTATTATTCAGCGACTGA‑CCGGAAACCCGATGGCAAGCCCGG‑AAGt                                                                               <  1:112768/251‑1 (MQ=255)                                                                                                                                                            gCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTTGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGGCATGCTGGCGG‑GGGCGGGGTGTATTTTTTAGCGACCGA‑ACGGAAAcaccggggcaaaccccggaggggg                                                                             >  2:179022/1‑228 (MQ=255)                                                                                                                                                             cTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTTGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGG‑TGGCGGGCTGTATTATTCAGCGACTGA‑CCGGAAACCCGAt                                                                                               <  1:125733/232‑1 (MQ=255)                                                                                                                                                             cTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTTGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGG‑GGGCGGGGTGTTTTTTTCTGGGACTGG‑CCGGAGAACCGAtagatacggaacggggcggg                                                                            >  2:125733/1‑232 (MQ=255)                                                                                                                                                             cTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTTGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGGCATGCTGGCGG‑GGGGGGGCTGGTTTATTTAGCGGCTGA‑CCGGGAACCCGAAGGCAAGCCCGG‑aggggctg                                                                           >  1:310949/1‑244 (MQ=255)                                                                                                                                                                         aTGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTTGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGG‑TGGCGGGCTGTATTATTCAGCGACTGA‑CCGGAAACCCGATgg                                                                                             <  2:243174/222‑1 (MQ=255)                                                                                                                                                                              cTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTTGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGGGTCATGCTGGCGG‑TGGGGGGCTGTATTATTCAGGGACTGA‑CCGGAAAACCGAGGGGCAGACCGG‑AAGGGGTGGGGGTTAGCTCcggggg                                                          >  1:73662/1‑248 (MQ=255)                                                                                                                                                                                 gCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGGTGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGG‑TGGCGGGCTGTATTATTCAGCGACTGG‑CCGGAAACCCGATGGCAAGCCCGG‑AAGGGCTGGGGATTAGCTCcggcgcgga                                                       <  2:19629/251‑2 (MQ=255)                                                                                                                                                                                     ttGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTTGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGG‑GGGCGGGCTGGATTATTCAGCGACTGA‑CCGGAAACCCGATGGCAAACCCCG‑AAGTGCTGGGGGTTAGCTCcggcgcggcggtt                                                   >  1:248663/1‑248 (MQ=255)                                                                                                                                                                                      tGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTTGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGG‑TGGCGGGCTGTATTATTCAGCGACTGA‑CCGGAAACCCGATGGCAAGCCCGG‑AAGTGCTGGGGATTAGCTCCGGCGCGGCGttt                                                   <  2:139594/250‑1 (MQ=255)                                                                                                                                                                                       gCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTTGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGG‑TGGCGGGCTGTATTATTCAGCGACTGA‑CCGGAAACCCGATGGCAAGCCCGG‑AAGTGCTGGGGATTAGCTCCGGCGCGGCGTTTg                                                  <  1:307487/250‑1 (MQ=255)                                                                                                                                                                                        cccTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTTGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGGTTGCGGGCGTCCAGCTGGGGG‑TGGTGGGGTGTATTATTCCGCGCCTGGGCGGGAAACCC‑ATGGCAAACCCGGGAAG‑GCTGGGGGTTAGCTCccgcgcgggggtttgg                                                >  2:282817/1‑243 (MQ=255)                                                                                                                                                                                            cgcgGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTTGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGG‑TGGCGGGCTGTATTATTCAGCGACTGA‑CCGGAAACCCGATGGCAAGCCCGG‑AAGTGCTGGGGATTAGCTCCGGCGCGGCGTTTGGCgt                                              <  1:118527/249‑1 (MQ=255)                                                                                                                                                                                                        gctgctGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTTGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGG‑TGGCGGGCTGTATTATTCAGCGACTGA‑CCGGAAACCCGATGGCAAGCCCGG‑AAGTGCTGGGGATTAGCTCCGGCGCGGCGTTTGGCGTGGTGTGGATGCTGt                                <  1:265741/251‑1 (MQ=255)                                                                                                                                                                                                              gTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTTGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGG‑TGGCGGGCTGTATTATTCAGCGACTGA‑CCGGAAACCCGATGGCAAGCCCGG‑AAGTGCTGGGGATTAGCTCCGGCGCGGCGTTTGGCGTGGTGTTGATGCTGTTTCt                            <  2:248663/249‑1 (MQ=255)                                                                                                                                                                                                              gTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTTGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGG‑TGGCGGGCTGTATTATTCAGCGACTGA‑CCGGAAACCCGATGGCAAGCCCGG‑AAGTGCTGGGGATTAGCTCCGGCGCGGCGTTTGGCGTGGTGTTGATGCTGTTTCt                            <  1:93672/249‑1 (MQ=255)                                                                                                                                                                                                              gTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTTGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGG‑TGGCGGGCTGTATTATTCAGCGACTGA‑CCGGAAACCCGATGGCAAGCCCGG‑AAGTGCTGGGGATTAGCTCCGGCGCGGCGTTTGGCGTGGTGTTGATGCTGTTTCTg                           <  2:183961/250‑1 (MQ=255)                                                                                                                                                                                                                    ggCTGTGGTGGTGGCGCTGTCGTTTGGTTGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGG‑TGGCGGGCTGTATTATTCAGCGACTGA‑CCGGAAACCCGATGGCAAGCCCGG‑AAGTGCTGGGGATTAGCTCCGGCGCGGCGTTTGGCGTGGTGTGGATGCTGTTTCTGGTGCCgg                    <  1:272722/251‑1 (MQ=255)                                                                                                                                                                                                                              gtggCGCTGTCGTTTGGTTGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGG‑TGGCGGGCTGTATTATTCAGCGACTGA‑CCGGAAACCCGATGGCAAGCCCGG‑AAGTGCTGGGGATTAGCTCCGGCGCGGCGTTTGGCGTGGTGTTGATGCTGTTTCTGGTGCCGGGTAATGCCtt          <  1:168376/251‑1 (MQ=255)                                                                                                                                                                                                                                ggCGCTGTCGTTTGGTTGTGATGCGCACGGCTGGACGTGGGCGAGAGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGG‑TGGCGGGCTGTATTATTCAGCGACTGA‑CCGGAAACCCGATGGCAAGCCCGG‑AAGTGCTGGGGATTAGCTCCGGCGCGGCGTTTGGCGTGGTGTTGATGCTGTTTCTGGTGCCGGGTAATGccgtt         <  2:91594/250‑4 (MQ=255)                                                                                                                                                                                                                                    cTGTCGTTTGGTTGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGG‑TGGCGGGGTGTATTATTTAGCGACTGA‑CCGGAAAACCGATGGCAAACCCGG‑GAGTGCTGGGGGTTAGCCCCGGCGCGGCGGTTGGCGTGGTGGTGGTGCTGTTTTTGGTGCCGGGTGATGCCTTttgctg    >  1:258315/1‑251 (MQ=255)                                                                                                                                                                                                                                    cTGTCGTTTGGTTGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGG‑GGGCGGGGTGGATTTTTTCGCGACTGA‑CCGGGAAACCGATGGGAAACCCGG‑GAGGGCTGGGGACTAGCCCCCGCGCGGCGCTTGGCGTGGGGTTGTTGCTGTTTCTCGTGcccggggacgcccttggccg    >  1:2985/1‑233 (MQ=255)                                                                                                                                                                                                                                      gTCGTTTGGTTGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGG‑GGGCCGGCTGCATTATTCAGCGGCTGA‑CCGGAAACCCCCTGGCAAACCCCG‑ACGTGCTGGGGATTTGCTCCGGCGCGGCGGTTGGGGGGGGGTTGATGGTTTTTCTTGTGCCGGGTAATGCCTTtgggcggc  >  1:369595/1‑246 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    TTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGATGGAAGTGTCCACCGGTTCGGTCACTGCGTTGATCGGTGCGCCGCTGCTACTGTGGCTGTTGCCGCGTTTACGCAGCATTAGCGCGCCGGATATGAAGGTCAACGATCGTGTCGCGGCTGAACGCCAACATGTGCTGGCGTTTGCCCTCGCGGGCGGCGTGCTGCTGTTGATGGCTGTGGTGGTGGCGCTGTCGTTTGGTCGTGATGCGCACGGCTGGACGTGGGCGAGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGG‑TGGCGGGCTGTATTATTCAGCGACTGA‑CCGGAAACCCGATGGCAAGCCCGG‑AAGTGCTGGGGATTAGCTCCGGCGCGGCGTTTGGCGTGGTGTTGATGCTGTTTCTGGTGCCGGGTAATGCCTTTGGCTGGC  >  NC_000913/172330‑172809

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