BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	498,116	C→T	P21L (CCT→CTT)	hemH →	ferrochelatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	498,116	0	C	T	99.1%	391.7 / NA	111	P21L (CCT→CTT)	hemH	ferrochelatase
Reads supporting (aligned to +/- strand): ref base C (1/0); major base T (60/49); minor base G (0/1); total (61/50)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 111 total aligned reads displayed.

GGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACCT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCTGGGAATGAGCT  >  NC_000913/497877‑498358

ggAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTa                                                                                                                                                                                                                                          >  1:492545/1‑250 (MQ=255)

ggAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTa                                                                                                                                                                                                                                          >  1:210837/1‑250 (MQ=255)

    cAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACg                                                                                                                                                                                                                                     <  1:59513/251‑1 (MQ=255)

    cAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACg                                                                                                                                                                                                                                     >  1:292898/1‑251 (MQ=255)

        cccGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCCATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCCACAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGAACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTAt                                                                                                                                                                                                                                 >  1:851118/1‑251 (MQ=255)

               aaGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAAc                                                                                                                                                                                                                          <  2:434988/251‑1 (MQ=255)

                     aCTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCCCACTT‑GAAGCGGTAAAACGCTATCTGAAACAAtttt                                                                                                                                                                                                                    >  2:259405/1‑251 (MQ=255)

                     aCTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAAtttt                                                                                                                                                                                                                    <  1:331187/251‑1 (MQ=255)

                     aCTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAAtttt                                                                                                                                                                                                                    <  2:301148/251‑1 (MQ=255)

                        tATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACTCTATCTGAAACAATTTTTa                                                                                                                                                                                                                  <  1:252406/250‑1 (MQ=255)

                        tATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAAttttt                                                                                                                                                                                                                   <  1:210979/249‑1 (MQ=255)

                         aTCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAATCCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAg                                                                                                                                                                                                                <  2:242033/251‑1 (MQ=255)

                             gCAATATCAGGACGGATATTCATTCGGCCTTTTACAAAAACCTCACTCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAACACGCTATCTGAAACAATTTTTAAGCGAc                                                                                                                                                                                                            >  2:397889/1‑251 (MQ=255)

                              cAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATGCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACa                                                                                                                                                                                                           <  1:397889/251‑1 (MQ=255)

                                aTCTCCGGCCGGATATTCATTCGGCCTTTTACACAAACATCCATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTCGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTCAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGGTGCCCCCACACTT‑GAAGCGGTCTAAAGCTATCTGAAACCATTTTTAAGCGAc                                                                                                                                                                                                            >  2:847480/1‑248 (MQ=255)

                                   tcAGGCCGGATATTCATTCGGCCTTTTCCAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGAcgc                                                                                                                                                                                                      >  1:717120/1‑251 (MQ=255)

                                    cAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATTGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGAcgcg                                                                                                                                                                                                     >  2:760276/1‑251 (MQ=255)

                                     aGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGAcgcgg                                                                                                                                                                                                    >  1:774079/1‑250 (MQ=255)

                                                       gCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCCCACGgtt                                                                                                                                                                                  >  2:675290/1‑251 (MQ=255)

                                                              aCAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCTCGGTTGTTAtg                                                                                                                                                                            <  1:753661/250‑1 (MQ=255)

                                                              aCAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTAtg                                                                                                                                                                            <  2:210837/250‑1 (MQ=255)

                                                                aaaaaCCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTAtggtg                                                                                                                                                                         <  2:292898/251‑1 (MQ=255)

                                                                     ccTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCAt                                                                                                                                                                    <  2:492545/251‑1 (MQ=255)

                                                                     ccTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCAt                                                                                                                                                                    <  1:855490/251‑1 (MQ=255)

                                                                       tCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTACTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAAACTGGGTCCGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATtg                                                                                                                                                                  >  1:149919/1‑251 (MQ=255)

                                                                        cAATCTGCGCTATTCCTTGTTCTGATTGGACCTCTCACAGCAATTAGTTCTTCTTCCTCACGTTTCCGCTACAATGATCAACAAGTGGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTGAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATtgc                                                                                                                                                                 <  2:31998/251‑1 (MQ=255)

                                                                         aaTCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATtgct                                                                                                                                                                >  1:9978/1‑251 (MQ=255)

                                                                            cTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTgcg                                                                                                                                                             >  1:538248/1‑251 (MQ=255)

                                                                             tGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGgttgt                                                                                                                                                                                <  1:861917/231‑1 (MQ=255)

                                                                             tGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGgttgt                                                                                                                                                                                >  2:861917/1‑231 (MQ=255)

                                                                               cgcTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGAc                                                                                                                                                                                                            >  2:322423/1‑201 (MQ=255)

                                                                               cgcTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGAc                                                                                                                                                                                                            <  1:322423/201‑1 (MQ=255)

                                                                                   aTTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTg                                                                                                                                                      <  1:271451/251‑1 (MQ=255)

                                                                                      ccTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAAAAATTTTTAAGCGACAGACGCGGGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGAtt                                                                                                                                                   >  2:730025/1‑251 (MQ=255)

                                                                                           ttCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑TAAGCGGTAAAACGCTATCTg                                                                                                                                                                                                                              >  1:535531/1‑171 (MQ=255)

                                                                                               gATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCt                                                                                                                                           <  2:326103/250‑1 (MQ=255)

                                                                                               gATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTg                                                                                                                                          >  1:303098/1‑251 (MQ=255)

                                                                                                aTTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTgc                                                                                                                                         >  1:920346/1‑251 (MQ=255)

                                                                                                aTTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGGGATTTTGCCGCTgc                                                                                                                                         >  2:343753/1‑251 (MQ=255)

                                                                                                  ttGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTTGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTgcg                                                                                                                                        <  1:604973/250‑1 (MQ=255)

                                                                                                   tGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTACCCTGCTGGCAAACCTGGGGACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTTAGCGACAGAAGCGGGGTTGATACCTCTCGGTTGTTTTGGCTGCCACTGCTGCgcggcggggtttttgcggctcgcc                                                                                                                                      >  2:698432/1‑233 (MQ=255)

                                                                                                     aCCTCTCACAGCAATAAGTTCTTCTTCCTCACTTTTCCGCTACCATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCCGACTACAAACGGTATCCTGCTGGCAAACCTGGGTACGCCAGCTGCCCCCACACTT‑GCAGCGGTAAAACGCTATCTGAAACAATTTTTCAGCGACAGACGCGTGGTTGATACCTCCCGGTTGTTATGGGGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTcc                                                                                                                                    >  2:613197/1‑250 (MQ=255)

                                                                                                        tctcACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTcgccg                                                                                                                                 >  1:320113/1‑251 (MQ=255)

                                                                                                         ctcACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGCCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTGGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCccgccgg                                                                                                                                >  2:532410/1‑250 (MQ=255)

                                                                                                         ctcACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATAATCAACAAGTTGAATCGATCAGAGGCGGTAATGCGTCCGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGCTGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTcgccgc                                                                                                                                >  1:427770/1‑251 (MQ=255)

                                                                                                          tcACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCcgcg                                                                                                                               >  1:498939/1‑251 (MQ=255)

                                                                                                               gCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTgg                                                                                                                          >  1:537249/1‑251 (MQ=255)

                                                                                                               gCAATTAGTTCTGCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGGAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATTCCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGAttt                                                                                                                                                  <  1:847079/227‑1 (MQ=255)

                                                                                                                   ttAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCCTCAGACCAAAACCGGTATCCTGCTGGCCAACCTGGGTACGCCCGATGCCCCCA‑ACCTTGAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGCCGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGa                                                                                                                       >  1:139044/1‑250 (MQ=255)

                                                                                                                    tAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATTCCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGCTTTTGCCGCTGCGCTCGCCGCGTGTGGCGa                                                                                                                       >  2:264059/1‑249 (MQ=255)

                                                                                                                    tAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGaa                                                                                                                      <  2:261877/250‑1 (MQ=255)

                                                                                                                      gTTCTTATTCCTCACTTTTCCGCTACAGTTATCAACAAGCTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGGATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTTAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATATCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCt                                                                                                                   >  1:651107/1‑251 (MQ=255)

                                                                                                                         cttcttcCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑TAAGCGGTAAAAAGCTATATGAAACAATTTTTAAGCGACAGAAGCGTGGTTGATACATCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTccgatcggaagagcgtcgagtg                                                                                                                >  2:709632/1‑230 (MQ=255)

                                                                                                                         cttcttcCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGt                                                                                                                 >  1:915584/1‑250 (MQ=255)

                                                                                                                         cttcttcCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGAGTTGGCCGCTGCGCTc                                                                                                                                     <  1:709632/230‑1 (MQ=255)

                                                                                                                              tcCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTAGCCGCGTGTGGCGAAGCTGTATGcc                                                                                                            <  1:532410/250‑1 (MQ=255)

                                                                                                                               ccTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCct                                                                                                           <  1:847480/250‑1 (MQ=255)

                                                                                                                                cTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGcc                                                                                                            <  2:81104/248‑1 (MQ=255)

                                                                                                                                cTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGcc                                                                                                            >  1:81104/1‑248 (MQ=255)

                                                                                                                                 tCACTTTTCCGATACAATTATAACCAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTAACCTGCTGGCAAACCTGGGGACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGGGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCGGCGCGGGGTGATTTTGCCGCTGCGCTAGACGCGGGGGGGGAAGCGGTATGCCtctg                                                                                                        >  2:844508/1‑251 (MQ=255)

                                                                                                                                  cACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCctct                                                                                                         <  1:343753/249‑1 (MQ=255)

                                                                                                                                      tttCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTTGGTACGCCCGATGCCCCCTCACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCtctgtctg                                                                                                    <  1:928742/250‑1 (MQ=255)

                                                                                                                                      tttCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCtctgtctg                                                                                                    <  1:259405/250‑1 (MQ=255)

                                                                                                                                              aCAATTATCACCAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAgg                                                                                           >  1:243694/1‑251 (MQ=255)

                                                                                                                                                aaTTATCACCAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCCAACCTGGGAACGCCCGCTGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGAGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTg                                                                                         >  2:353752/1‑251 (MQ=255)

                                                                                                                                                    aTCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTAGGCCTCTGTCTGGATGGAAGGTGGCt                                                                                      <  2:149919/250‑1 (MQ=255)

                                                                                                                                                          aaGTGGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCt                                                                               <  2:320113/251‑1 (MQ=255)

                                                                                                                                                           aGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGGAGCTGGATGCCCCTGTCTGGATGGAAGGTGGCCCGCCGCTg                                                                              >  2:221117/1‑251 (MQ=255)

                                                                                                                                                              tGAATAGATAAGAGGTGGTACGGCGTAAGATTAAAAACGGTACCCTGCTGACAAAACTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTTAAACGAAATCTGACACAATTTTTAAGCGAACGACGCGTGGTTGATACCTCCCGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGGTTTTGCCGCCGCGCTTGCCGCGTGTGGCGAAGCTGTATGCCCCTTACTGGGTGGATGGTGGCTCGCCGGTGATg                                                                           >  2:861777/1‑251 (MQ=255)

                                                                                                                                                               gAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTCGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCCCGGTTGTTATGGTGGCCATTGCTGCGCGGGGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAAGGGGCTCGCCGCTGATgg                                                                          >  2:764158/1‑251 (MQ=255)

                                                                                                                                                               gAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCCGGCAAACCTGGGGACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGGGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATgg                                                                          >  1:47127/1‑251 (MQ=255)

                                                                                                                                                                aaTCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGt                                                                         <  1:221117/251‑1 (MQ=255)

                                                                                                                                                                    gATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGGTGGGAGGGGGCTCGCCGCCGATGGTTTAc                                                                     >  2:756754/1‑251 (MQ=255)

                                                                                                                                                                    gATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTAc                                                                     >  1:288674/1‑251 (MQ=255)

                                                                                                                                                                    gATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTAc                                                                     >  1:598415/1‑251 (MQ=255)

                                                                                                                                                                       aaGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCAGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAg                                                                   <  1:803513/250‑1 (MQ=255)

                                                                                                                                                                          gggCGGTAATGCGTCAGACTAAAACCGGTATCCTGCGGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAgccgc                                                               <  2:538248/250‑1 (MQ=255)

                                                                                                                                                                          agGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAgccgc                                                               >  1:582266/1‑251 (MQ=255)

                                                                                                                                                                            gCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCa                                                             <  2:537249/251‑1 (MQ=255)

                                                                                                                                                                                 aaTGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCaa                                                         <  1:730025/250‑1 (MQ=255)

                                                                                                                                                                                      gTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCGGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAAcagcag                                                   <  1:353752/251‑1 (MQ=255)

                                                                                                                                                                                        cAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACGT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCGCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGgc                                                 <  1:844508/251‑1 (MQ=255)

                                                                                                                                                                                           aCTAAAACCGGTTTCCTGCTGGCCAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCt                                              <  2:774079/251‑1 (MQ=255)

                                                                                                                                                                                                  ccGGTATCCGGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATTGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACaa                                       <  1:53269/251‑1 (MQ=255)

                                                                                                                                                                                                   cGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAAc                                      <  2:598415/251‑1 (MQ=255)

                                                                                                                                                                                                    ggTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTcgcc                                                                                  >  1:386030/1‑206 (MQ=255)

                                                                                                                                                                                                    ggTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCCGGCGCTGGCACAACg                                     >  2:498482/1‑251 (MQ=255)

                                                                                                                                                                                                    ggTATCCGGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATTGAAGGTGGCTcgcc                                                                                  <  2:386030/206‑1 (MQ=255)

                                                                                                                                                                                                       aTCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACACCCGCCAGCAACAGCAGGCGCTGGCACAACGttt                                  >  2:19910/1‑251 (MQ=255)

                                                                                                                                                                                                           tgctgGCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCg                              >  2:736965/1‑251 (MQ=255)

                                                                                                                                                                                                                gCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACCATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGCTTTTGCCGCTGCGCTCGCCGCGTGTgg                                                                                                                          >  2:426183/1‑154 (MQ=255)

                                                                                                                                                                                                                gCAAACCTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTgg                                                                                                                          <  1:426183/154‑1 (MQ=255)

                                                                                                                                                                                                                     ccTGGGTACGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCg                    <  1:675290/251‑1 (MQ=255)

                                                                                                                                                                                                                             cGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGc                                                                                       >  2:372035/1‑176 (MQ=255)

                                                                                                                                                                                                                             cGCCCGATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGc                                                                                       <  1:372035/176‑1 (MQ=255)

                                                                                                                                                                                                                                  gATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGc                                           <  2:295854/215‑1 (MQ=255)

                                                                                                                                                                                                                                  gATGCCCCCACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGc                                           >  1:295854/1‑215 (MQ=255)

                                                                                                                                                                                                                                       ccccACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCTGGGAATGAGCt  >  2:488991/1‑251 (MQ=255)

                                                                                                                                                                                                                                       ccccACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCTGGGAATGAGCt  <  1:488991/251‑1 (MQ=255)

                                                                                                                                                                                                                                       ccccACACTT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCTGGGAATGAGCt  >  1:137198/1‑251 (MQ=255)

GGAACAGGCCCGCACAAGCAAACTTATCAGCAATCTCAGGCCGGATATTCATTCGGCCTTTTACAAAAACCTCAATCTGCGCTATTCCTTTTTCTGATTTGACCTCTCACAGCAATTAGTTCTTCTTCCTCACTTTTCCGCTACAATTATCAACAAGTTGAATCGATAAGAGGCGGTAATGCGTCAGACTAAAACCGGTATCCTGCTGGCAAACCTGGGTACGCCCGATGCCCCCACACCT‑GAAGCGGTAAAACGCTATCTGAAACAATTTTTAAGCGACAGACGCGTGGTTGATACCTCACGGTTGTTATGGTGGCCATTGCTGCGCGGCGTGATTTTGCCGCTGCGCTCGCCGCGTGTGGCGAAGCTGTATGCCTCTGTCTGGATGGAAGGTGGCTCGCCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCTGGGAATGAGCT  >  NC_000913/497877‑498358

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 13 ≤` `ATCG/ATCG` `< 16 ≤` `ATCG/ATCG` `< 33 ≤` `ATCG/ATCG` `< 39 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`