Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,297,469	C→T	G633D (GGC→GAC)	fdhF ←	formate dehydrogenase H

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,297,469	0	C	T	100.0%	227.9 / NA	73	G633D (GGC→GAC)	fdhF	formate dehydrogenase H
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (43/30);  total (43/30)

CGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACT‑‑GGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCAC‑GGGTGATAATTTTGCCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACC  >  NC_000913/4297222‑4297710                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                  cgccatggCCGCTTCGGGAAGGCGAGGTTGCAACTGGTTGTACTCGTCAATCACGTACTGCTCGGCGGGGCTCGGAGCGGCGATCTGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTGTCGGTAACCAGCTCGTTACAGGCACCAATCCTCCACT‑‑GGTAGGTCAGGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCTCGC‑GGGTGAGAATTTTGTCtt                                                                                                                                                                                                                                                <  2:527450/244‑1 (MQ=255)  gCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACT‑‑GGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGCTCGCTGACCTGCG‑ACAGGGTGATAATTTTGTCttt                                                                                                                                                                                                                                               >  2:393416/1‑251 (MQ=255)   ccAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACT‑‑GGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCAC‑GGGTGATAATTTTGTCttt                                                                                                                                                                                                                                               >  1:526947/1‑250 (MQ=255)     aGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACT‑‑GGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCAC‑GGGTGATAATTTTGTCTTTAcg                                                                                                                                                                                                                                            <  2:99673/251‑1 (MQ=255)          cGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGCTCGGCTAGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCACGTTACAGGCACCAATCCACCACT‑‑GGTAGGTCCTGTAAATCGCCCCTTTGTTCGGACGATCGCTGCCCTGCGCAC‑GGGTGATAATTTTGTCTTTAcgcgc                                                                                                                                                                                                                                         >  2:544982/1‑248 (MQ=255)          cGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTAAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACT‑‑GGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCAC‑GGGCGATAATTTTGTCTTTATGCGAGt                                                                                                                                                                                                                                       >  2:64527/1‑251 (MQ=255)           gCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACT‑‑GGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCAC‑GGGTGATAATTTTGTCTTTACGCGAGTg                                                                                                                                                                                                                                      <  1:129098/251‑1 (MQ=255)              tCGCGCAGGCGAGTTTTCAACTTGTTGTACCCGTCAATCACGTACGGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACT‑‑GGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCAC‑TGGTGATAATTTTGTCTTTACGCGAGTGCa                                                                                                                                                                                                                                    <  2:482496/250‑1 (MQ=255)                   cAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGACGGCGCACTGATCGGCGATCGGCTCGACGCAAACGACGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGCCTTACGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACT‑‑GGTAGGTAATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCCC‑GGGTGATAATTTTGTCTTTACGCGAGTGCaacaac                                                                                                                                                                                                                               >  1:81973/1‑245 (MQ=255)                      gCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACT‑‑GGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCAC‑GGGTGATAATTTTGTCTTTACGCGAGTGCACCCaaaac                                                                                                                                                                                                                            >  1:547690/1‑248 (MQ=255)                          gTTTTCAACTTGTTGTACTCGGCAATCACGTACTGCGCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCATTATTTGTACTCCGGCGTTTTCGTAAGCGGGCTTAAGGTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACT‑‑TGTAGGTCATGTAAATCGCCCCTTTTTTCGGACGATCGCTGACCTGCGCAC‑GGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGc                                                                                                                                                                                                                       <  2:336353/251‑1 (MQ=255)                            tttCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACT‑‑GGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCAC‑GGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGcc                                                                                                                                                                                                                      <  1:66238/250‑1 (MQ=255)                             ttCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACT‑‑GGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCAC‑GGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTc                                                                                                                                                                                                                    >  1:483318/1‑251 (MQ=255)                             ttCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTGCGGTAACCAGCTCGTTACAGGCACCAATCCACCACT‑‑GGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCAC‑GGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCt                                                                                                                                                                                                                     <  2:425612/250‑1 (MQ=255)                                aaCTTTTTGTACTAGTCAATCACGTACTGATCGGCGGCGCGCTGCTCGGCGCTCGGCACGGAGAGAACGGCCCAGAATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCCGCTCGTTACCGGCACCAATCCACAACT‑‑GGTACGTCATGTCAATCACCCCTTTGTTCGGACGATAGCTGACCTGCGCAC‑GGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATc                                                                                                                                                                                                                 >  2:575927/1‑251 (MQ=255)                                 aCTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACT‑‑GGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCAC‑GGGCGATAAGTTTGTCTTTACACGAGTGCACCCAAACCAATGCCTCATCt                                                                                                                                                                                                                >  1:318410/1‑251 (MQ=255)                                    tgttgtACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACT‑‑GGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCAC‑GGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTc                                                                                                                                                                                                              >  1:537706/1‑250 (MQ=255)                                              gTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACT‑‑GGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCAC‑GGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGAc                                                                                                                                                                                                   >  1:458941/1‑251 (MQ=255)                                                cAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACT‑‑GGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCAC‑GGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGt                                                                                                                                                                                                 >  1:164676/1‑251 (MQ=255)                                                 aaTCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACT‑‑GGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCAC‑GGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGtt                                                                                                                                                                                                >  2:19677/1‑251 (MQ=255)                                                    cACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACT‑‑GGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCAC‑GGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTg                                                                                                                                                                                              <  1:428001/250‑1 (MQ=255)                                                    cACGTACTGCTCGGCGGCGCGCAGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACT‑‑GGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCAC‑GGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTgg                                                                                                                                                                                             >  1:398797/1‑251 (MQ=255)                                                      cggaCTGCTCGGCGGCGCGCTGAGCGGCGATCGGCTCGACGCGAACGGCGCAGGATTTTTACTCCGGCGTTTTAGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCAACACT‑‑GGTAGGTCATGTAAATCGCACCTTTGTTCGGACGATCGCTGACCTGAGCCC‑GGGTGATAATTTTGTCTTTACGCGAGTGCACCCATACCAATGCCACATCTTCAATACCCCGACGTTTGGct                                                                                                                                                                                           >  2:280592/4‑250 (MQ=255)                                                            gCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACT‑‑GGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCAC‑GGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCtt                                                                                                                                                                                       <  1:146668/249‑1 (MQ=255)                                                            gCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACT‑‑GGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCAC‑GGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTc                                                                                                                                                                                      <  1:88558/250‑1 (MQ=255)                                                                      gcgcTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGCACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACT‑‑GGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCAC‑GGGGGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGAttt                                                                                                                                                                           >  2:39108/1‑251 (MQ=255)                                                                          tGATCTGCGATCGGCTGTCCGCGAACGGCCCAGCATTTGTACTCCGGAGTTTTCGTCATCGGGCTTACGTTTTCAATAACCAGCTCGTTACAGGCACCAATCCACCACT‑‑GGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCCC‑GGGTGATAATTTTGTCTTTACGCGCGTGCCCCCCACCCACTGCATCATCTTCAATACACAGCCGTTTTGCCTCTTCGGTATTGTTTTGTGc                                                                                                                                                                       >  2:289622/1‑251 (MQ=255)                                                                           gatcggcgatcggcTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACT‑‑GGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCAC‑GGGTGATAATTTTGTCTTTACGCGAGTGCACCCAATCCAATGCCTCATCTTCAACACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCg                                                                                                                                                                      <  2:89335/251‑1 (MQ=255)                                                                            atcggcgatcgccTCGACGAGAACGGCGCAGTATTTGTACTCCGGCCTTATAGCAATCCGGCTTAAGTTTACGGAAACCAGCTCGATACAGGCCCCAA‑‑CACCAAAAAGGTAGGACATGTAAATCGCCACTTTGTTCGGACGAACGCTGACCTGAACAC‑GGGAGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTACATACCAAGACATTTCGCGTCTTCGGAATTGATTTGTGCGt                                                                                                                                                                     >  2:387322/1‑251 (MQ=255)                                                                                gcgatcggcTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACT‑‑GGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCAC‑GGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGcc                                                                                                                                                                 <  1:42017/251‑1 (MQ=255)                                                                                 cgatcggcTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACT‑‑GGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCAC‑GGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCa                                                                                                                                                                <  1:278556/251‑1 (MQ=255)                                                                                          cGACGCGAACGGCGCAGTCTTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACAAGCTCGTTACAGGCACCAATCCACCACT‑‑GGTAGGTCATGTAAATCGCCCCTTTGTTCGGCCGATCGCTGACCTGCGCAC‑GGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTtcatc                                                                                                                                                        >  2:580913/1‑250 (MQ=255)                                                                                                gAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACT‑‑GGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCAC‑GGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATcagccagc                                                                                                                                                 >  1:425580/1‑251 (MQ=255)                                                                                                   cGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACT‑‑GGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCAC‑GGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGcgcc                                                                                                                                              >  1:108163/1‑251 (MQ=255)                                                                                                        cAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCCCT‑‑GTTAGGTCGTGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCAC‑GGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATACCTCATCTTCAATACCCAGACGCTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAACGCCGCCat                                                                                                                                         >  1:210920/1‑250 (MQ=255)                                                                                                             tttGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACT‑‑GGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCAC‑GGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCg                                                                                                                                    >  1:554410/1‑251 (MQ=255)                                                                                                                                 tAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACT‑‑GGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCAC‑GGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGa                                                                                                                 >  1:475533/1‑250 (MQ=255)                                                                                                                                  aaTCGCGCTTAAATTTACGGTAACCAGCTCGTTCCAGACACCAATCCAAAACT‑‑GGTCGGCCCTGTAACACGACCCTTTGCTCGGACGCTCGATGACTTGCGCAC‑GGGTCCTAATTTTGTCTCTACGCGAGTGCACCCAAACAAATGCAGAATCTTCAATAACCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCAGCCAGTGCCGCATAGCTACCGGTCATCGAAc                                                                                                               >  2:572955/1‑251 (MQ=255)                                                                                                                                   aTCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACT‑‑GGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCAC‑GGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACg                                                                                                              >  1:352040/1‑251 (MQ=255)                                                                                                                                         cTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCAAT‑‑GGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCAC‑GGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGCATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCTACGCCAGTGCCGCACAGTTACCGATCATCGAACGGCa                                                                                                           >  2:189170/1‑248 (MQ=255)                                                                                                                                          ttAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACT‑‑GGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCAC‑GGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGCACGGCAag                                                                                                         >  2:198876/1‑249 (MQ=255)                                                                                                                                           tAAGTTTTGGGTAACCAGCTCGTTACCGGCACCAATCCACCACT‑‑GGTAGGTCATGTAAATCCCCCCTTTGTTCGGCAGATCGCTGACCTGGGCAC‑GGGTGATAATTTTGTCTTGACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATCCCCAGACGTTTGGCGTCTTCGGTATGGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAag                                                                                                         <  2:514164/248‑1 (MQ=255)                                                                                                                                           tAAGTTTTCGGTAACCAGCTCGTTCCAGGCACCACGCCACCACT‑‑GGTAGGTCATGTAAATCGCCCCTTTGTCCGGACGATCGCTGACCTGCGCAC‑GGGTGATAATGTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCGTCAAGACCCAGACGTTGGGCGTCTTCGATATTGATTTGGGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGGCATCGAACGGCAag                                                                                                         <  2:392912/248‑1 (MQ=255)                                                                                                                                                    ggTAACCAGCGCGTTACAGGCACCAATCCACCACT‑‑GGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCAC‑GGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCaa                                                                                             <  1:25942/251‑1 (MQ=255)                                                                                                                                                           aGCTCGTTACAGGCACCAATCCACCACT‑‑GGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCCC‑GGTTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTACGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACg                                                                                      <  1:4330/251‑1 (MQ=255)                                                                                                                                                               cgtaacagtaccaaaTCCCCCACT‑‑GTTAGGTCATGTAAATCGCCCCTTTGTTCGGACTATCGCGTTCCTGCGCAC‑GGGTTATAATTTTGTCTTTACGCGAGTGCCCCCAACCCAAGGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTGGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCAcc                                                                                  <  2:92379/238‑1 (MQ=255)                                                                                                                                                                     aGGCACCACACCACCACG‑‑GGGAGGTCATGGAAATCGCCCTTTTTTTCGGACGTGCGCTGACCTGCGCAC‑GGTTGATAATTTTTTTTTGACGCGAGTGCACCCAAACCAATGCCTCATCTTAAATCCCCAGACGTGTGGCGTCTCCGGTATTGAGTTGTGCGTAGCCAGGTTCATCAGCCAGGGTCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGAc                                                                            <  2:526947/251‑1 (MQ=255)                                                                                                                                                                       gCACCAATCCCCCACT‑‑GGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCAC‑GGGTTATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCAGCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAg                                                                          <  1:279957/251‑1 (MQ=255)                                                                                                                                                                        ccccAATCCACCACT‑‑GGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCAC‑GGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGt                                                                         <  2:108163/249‑1 (MQ=255)                                                                                                                                                                               ccaccaCT‑‑GGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCAC‑GGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCg                                                                  >  1:454233/1‑251 (MQ=255)                                                                                                                                                                                 accacT‑‑GGTCGGTCCTGTACATAGCCACTTTGTTCGGACGATCGCTGACCTGCGCAC‑GGGTGATAATTTTGTCTTTACGCGCGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTCCCGGTCATCGAACGGCAAGAGTCGTGACCAACTTCACGCACCGTTGACAGTACCATCggg                                                                >  2:396124/1‑251 (MQ=255)                                                                                                                                                                                  cctcT‑‑GGTAGGTCATTTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCAC‑GGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGt                                                               <  2:494701/248‑1 (MQ=255)                                                                                                                                                                                   cacT‑‑GGTAGGGCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCAC‑GGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTGCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTa                                                              <  1:198876/251‑1 (MQ=255)                                                                                                                                                                                    acT‑‑GGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCAC‑GGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTGCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTAc                                                             <  1:393416/251‑1 (MQ=255)                                                                                                                                                                                            aGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCAC‑GGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACtcgtcg                                                       >  1:478065/1‑251 (MQ=255)                                                                                                                                                                                                    tAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCAC‑GGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGttt                                               >  1:250247/1‑251 (MQ=255)                                                                                                                                                                                                         cGCCCCTTTGTTCGGACGATCGCTGACCTGCGCAC‑GGGAGATAATTTTGTCTTTACGCGAGTGCACCCACACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGAATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCg                                           >  1:517471/1‑250 (MQ=255)                                                                                                                                                                                                              cTTTGTTCGGACGATCGCTGACCTGCGCAC‑GGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGgc                                     >  1:472449/1‑251 (MQ=255)                                                                                                                                                                                                              cTTTGTTCGGACGATCGCTGACCTGCGCAC‑GGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGgc                                     >  2:239930/1‑251 (MQ=255)                                                                                                                                                                                                               tttGTTCGGACGATCGCTGACATACGCAC‑GGGTGATAATTTTGTCTTTACGTGAGTGCACCCAAACCAATGCCTCATCTTCAATACACAGACGCTTGTCGTCTTCGGTATTGATGTGTGCGTAGCCAGGTTCATCAGGCAGAGCCGCCCGTGCCGCACAGTTAGCGGTCCTCGCAAGGAAAGAGTAGTGACGACCTTCACGCACCTATGAAAGTCCCATCGGGTACTCGTCGGAGCGTTTGTGCATTGgcg                                    >  2:190038/1‑251 (MQ=255)                                                                                                                                                                                                                     cGGACGATCGCTGACCTGCGCAC‑GGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTAccc                              >  1:8782/1‑251 (MQ=255)                                                                                                                                                                                                                      ggACGATCGCTGACCTGCGCAC‑GGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTAccc                              >  1:324383/1‑250 (MQ=255)                                                                                                                                                                                                                       gACGATCGCTGACCTGCGCAC‑GGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAg                            <  1:508886/251‑1 (MQ=255)                                                                                                                                                                                                                                 gACCTGCGCAC‑GGGTGATAATTTTGTCTTTACGCGAGTGCACCCAACCCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGCTTGGCGCTACCCCGTCGCAGGTg                   >  2:384475/1‑250 (MQ=255)                                                                                                                                                                                                                                   ccTGCGCAC‑GGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTgaag                <  1:174346/251‑1 (MQ=255)                                                                                                                                                                                                                                   ccTGCGCAC‑GGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTga                  >  1:117962/1‑249 (MQ=255)                                                                                                                                                                                                                                   ccTGCGCAC‑GGGGGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTga                  >  2:408730/1‑249 (MQ=255)                                                                                                                                                                                                                                      gcgcAC‑GGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGTCGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTgaagaa              <  1:183843/250‑1 (MQ=255)                                                                                                                                                                                                                                       cgcAC‑GGGTGATAATTTTGTCTTTACGCGGGTGCCCCCAAACCAATTCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCACCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGAtg                                        <  2:495820/223‑2 (MQ=255)                                                                                                                                                                                                                                       cgcAC‑GGGTGATAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGAtt                                        >  1:495820/1‑223 (MQ=255)                                                                                                                                                                                                                                           aggggTGATAATTTTTTCTTTACGCGAGTGCACCCAAACCAATGCCTCACCTTCAATCCCCAGACGTTTGGCGCCTTCGGTATTGATTTGTGCGTAGCCAGGTTCACCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTgcg         <  2:571086/249‑1 (MQ=255)                                                                                                                                                                                                                                             ggttGATATTTTTGTCTTAACGCGAGTGCACCCAAGCCAATGCTTCATCTCCAACACCCAGAGGTTTGGAGTCTGCGGTTTTGATTTGTGCGTAGCCAGGTTCATCGGCCAGGTCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACTGTACCAGCGGGTACTCGTCGGTGAGTTTGTCGATTGGAGCTACCCAGTGGCAGGTGAAGAACTGCGcc       <  2:35505/248‑1 (MQ=255)                                                                                                                                                                                                                                                  aTAATTTTGTCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGAcc  >  1:40263/1‑251 (MQ=255)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                  CGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGGCGTTTTCGTAATCGGGCTTAAGTTTTCGGTAACCAGCTCGTTACAGGCACCAATCCACCACT‑‑GGTAGGTCATGTAAATCGCCCCTTTGTTCGGACGATCGCTGACCTGCGCAC‑GGGTGATAATTTTGCCTTTACGCGAGTGCACCCAAACCAATGCCTCATCTTCAATACCCAGACGTTTGGCGTCTTCGGTATTGATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACC  >  NC_000913/4297222‑4297710

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