BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_007779	4,017,888	A→C	12.1%	S1112R (AGT→AGG)	Y75_RS19860 ←	rhs element protein RhsB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007779	4,017,888	0	A	C	12.1%	229.1 / 2.7	131	S1112R (AGT→AGG)	Y75_RS19860	rhs element protein RhsB
Reads supporting (aligned to +/- strand): ref base A (34/72); new base C (8/14); total (42/89)
Fisher's exact test for biased strand distribution p-value = 8.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 133 total aligned reads displayed.

CAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGG‑AACAC‑CACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGCTCCTGTCGTTCTGTAT  >  NC_007779/4017651‑4018115

    cGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCctcact                                                                                                                                                                                                                                      <  2:2106312/236‑1 (MQ=38)

             tgctgcAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAg                                                                                                                                                                                                                           >  1:1441610/1‑238 (MQ=35)

                tgcAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGGCAGCCAg                                                                                                                                                                                                                       >  1:796312/1‑239 (MQ=25)

                tgcAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCCCGATGGTCGCCGTGGTACAGGTGGAGTTTTCGCGCCGGCGGGTACACCGGGGCCCTCT‑GGGTTTTCATCTGCGCCAAGGTCAGGCCGCACGCGGCCAGCCAGCGGCGGCGTTTCCCCCTCCCCCCGTCAccccg                                                                                                                                                                                                                       >  2:1542593/1‑237 (MQ=11)

                     cTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCCCTCACCCCGTCAGCCAGGAtt                                                                                                                                                                                                                   >  1:1937203/1‑238 (MQ=21)

                      tgctGCGGGTTTTCTTCATTCAGCAGGTTGCCCCCTTCGTCATATTCTGCGCGCCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTTCAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCT‑GGTTTTGCATCTGCTCCACGCTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGAttt                                                                                                                                                                                                                  <  2:573405/235‑1 (MQ=255)

                                 ctctTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCTCGCTGTTGCCGCTTCCGGGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGcc                                                                                                                                                                                                      <  1:1131261/239‑1 (MQ=18)

                                   ctTCATTCAGCCGGATGCCCCATTCATCATATTATGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGAAGGTCGCATCGGAAAAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCT‑GGTGTTGCATCTGCTCCCCGGTCCGGCCGCCCCATGCCAGCCCGCGGCGAGTTTCTCCCAGCACCCCGTCCGCCCGGGTTTCGCCTTCCAACCgg                                                                                                                                                                                                    >  2:372392/1‑239 (MQ=11)

                                      tattcaGCAGGTTGCCCCATTCATCCTATGCTGCGCCCCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCCGCAGGCCGCGATTGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCg                                                                                                                                                                                                 <  1:1942871/238‑1 (MQ=12)

                                      cattcaGCAGGTTGCCCCATTCATCATGTTCTGCGCACCCCGCGTTTGCCCCTTCCGTGCTGATGAGCGCCCGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCg                                                                                                                                                                                                 <  2:548563/239‑1 (MQ=12)

                                       attcaGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGa                                                                                                                                                                                                <  1:449455/239‑1 (MQ=32)

                                            gcagGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCCTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGGGCa                                                                                                                                                                                             >  1:684640/1‑237 (MQ=14)

                                             cagGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCT‑GGTTTTGCATCTGCTCCCCGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCCCCCTCACCCGGCCAGCCAGGATTTCACTTTCCAGCCGGTCGCGCATc                                                                                                                                                                                           >  2:374196/1‑238 (MQ=9)

                                             cagGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGGCCAGCATCt                                                                                                                                                                                          >  1:1406845/1‑239 (MQ=14)

                                                gTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCCCCCACCCGGCCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGc                                                                                                                                                                                        >  1:1040688/1‑238 (MQ=14)

                                                gTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCCGc                                                                                                                                                                                        >  1:627750/1‑238 (MQ=14)

                                                  tGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCAcc                                                                                                                                                                                     >  1:50989/1‑239 (MQ=31)

                                                   gCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCCGGTCCATCT‑GGTTTTGCATCTGCTCCACCGTCCGGCCGCCCGATTCCAGCCAGCGGCGGCGTTCCTCCCTCCCCCCGGCGCCCAGGGTTTCACTTTTCACCCTGTCGAGCATCTGcccc                                                                                                                                                                                     >  2:1877098/1‑238 (MQ=11)

                                                     cccATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGc                                                                                                                                                                                  >  1:975964/1‑239 (MQ=32)

                                                      ccATTCATCATATTCTGCGCGCCACGGTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGTTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCa                                                                                                                                                                                 <  2:1914755/239‑1 (MQ=18)

                                                      ccATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCa                                                                                                                                                                                 <  1:2133386/239‑1 (MQ=33)

                                                      ccATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCa                                                                                                                                                                                 <  1:383850/239‑1 (MQ=33)

                                                      ccATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGAGGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCa                                                                                                                                                                                 <  1:1542593/239‑1 (MQ=18)

                                                       ccggcatcatATTCTGGCCACACCGCTTTTGCCCCTTCCGGGCGGATGGGCGCCGGCGGAAGGCCGGGATGGTGGCAGTGGTACAGGTGGATTTGTCGCGCCGGCGTGTACACCGGGTCCATCT‑GGTTTTGCAGCTGCTCCCCGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCa                                                                                                                                                                                 <  2:1610057/234‑1 (MQ=11)

                                                         ttCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACcg                                                                                                                                                                              >  1:1610148/1‑239 (MQ=34)

                                                          tcatcaTATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCACCCAGGATTTCACTTTCCAGCCGGTCGAGCCTCTGCACCAGCACcg                                                                                                                                                                              >  1:907135/1‑238 (MQ=21)

                                                           tctcatATTCTGCGCGCCCCGCTGTTGCCCCTCCCGTGCTGATGAGCGCCGGCGGCAGGCCGCGATGGTCGCGGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACcggc                                                                                                                                                                            <  2:431599/237‑1 (MQ=255)

                                                           catcatATTCTGCGCCCCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCGGCGGCAGGCCGCGATGGTCGCAGTGGTACGGGTGGATTTTTCGCGCCGGCGTGTACTCCGGGTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCGCCAGCACcggc                                                                                                                                                                            <  1:2038835/239‑1 (MQ=255)

                                                           catcatATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACcggc                                                                                                                                                                            <  1:1260179/239‑1 (MQ=35)

                                                             tcatATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACcggc‑gg                                                                                                                                                                         <  1:332972/239‑1 (MQ=37)

                                                             tcatATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACcggc‑gg                                                                                                                                                                         <  1:114085/239‑1 (MQ=37)

                                                             tcatATTCGGCGCACCGCGCCGTTGCCCCTTCCGGTCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCTCCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGc                                                                                                                                                                                              <  2:1798483/219‑1 (MQ=9)

                                                              catATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTTCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCCTCTGCACCAGCACCGGC‑ggg                                                                                                                                                                        >  1:1769762/1‑239 (MQ=37)

                                                              catATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑ggg                                                                                                                                                                        <  1:856833/239‑1 (MQ=37)

                                                              catATGCTGCGCCCCCCGGTGTGGCCCCTTCCGGGCTGATGGGCGCCGGCGGCGGGCCGCGATGGTCGCGGGGGTACAGTTGGATTTTTCGCGCCGGCGGGTACACCGGGTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCTCTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑ggg                                                                                                                                                                        <  2:1024869/239‑1 (MQ=11)

                                                                tatTCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGG‑aa                                                                                                                                                                     <  1:1874888/239‑1 (MQ=37)

                                                                  ttCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGG‑aaa                                                                                                                                                                    >  1:1829076/1‑237 (MQ=37)

                                                                   tCTGCGCCCCACGCTGTGGCCCCTTCCGGGCTGGTGAGGGCCAGCGGCAGGCGGGGATGGTGGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGGGTACACCGGGTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGG‑Aacac                                                                                                                                                                  <  2:1241373/239‑1 (MQ=11)

                                                                     tgctcgccacGCTGGTGCCCCTTCGTTGCTGGTGAGCGCCAGCGGCGGGCCGGTATGGTCGCAGTTGTACGGGTGGAGTTTTCGCGCCGGCGGGTACACCGGGTCCATTT‑GGTTTTCCATATGCTCCACCGTCAGGCCGCGCGAGGCCGGCCAGCGGCGGCTTTCCTCACTCACCCGCTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCGCCAGCACCGGC‑GAG‑AAccc‑cg                                                                                                                                                               <  2:568606/233‑2 (MQ=11)

                                                                     tGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTTCTCCCTCACCCCGTCCGCCCAGGTCTCCCTTTCCCGCCGGTCGAGCCTCCGCACCAGCCACGGC‑ggggcgacc                                                                                                                                                                  >  2:1704748/1‑232 (MQ=11)

                                                                       cgcACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGG‑Gacaaccc                                                                                                                                                               >  1:240657/1‑234 (MQ=37)

                                                                       cgcACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCCGTCAGCCAGGGCTTCACTTTCCAGCCGGTCGCGCATCTGCACCAGCACcggc‑ggtgaaaaccc                                                                                                                                                                >  1:603822/1‑229 (MQ=255)

                                                                       cgcACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGGACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCCGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCa                                                                                                                                                                                    >  1:1798483/1‑219 (MQ=14)

                                                                       cgcACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGAAGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCT‑GGTTTTGCATCTGCTCCACCGGCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCCGCCGGTCGAGCATCCGCACCAACACCGGC‑CGG‑GAcac‑ccc                                                                                                                                                              >  1:1795555/1‑236 (MQ=11)

                                                                        gccccccggtGTTGCCCCTTCCGTGTTGAGGAGCGCCGGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCCCCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGG‑AACAC‑caca                                                                                                                                                             <  2:840997/229‑1 (MQ=255)

                                                                        gcACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGG‑AACAC‑caca                                                                                                                                                             <  1:571002/238‑1 (MQ=37)

                                                                        gcACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGG‑AACAC‑Cacac                                                                                                                                                            <  1:1540680/239‑1 (MQ=37)

                                                                            cacGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑CGG‑GACCC‑Cccacccc                                                                                                                                                         >  1:2004330/1‑235 (MQ=255)

                                                                              cGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGCTGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCCGCGGGGACACCGGGTCCATCT‑TGTTTTGCATCTTCTCCACCGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTCTCCTCACTCCCCCGGCCTGCCCGGGCTTCACTTTCCTGCCGGTCTCGCATCTGCACCACCACCCGC‑GGGTAACACACACGCTGc                                                                                                                                                         >  2:710652/1‑238 (MQ=11)

                                                                                      ccccTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGG‑AACCC‑CCCACTcccccctcctc                                                                                                                                                >  1:1382712/1‑233 (MQ=255)

                                                                                         cTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGG‑AACAC‑CACACTGCCACCGTCTTcgccg                                                                                                                                           <  1:374196/239‑1 (MQ=33)

                                                                                           tCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGG‑AACAC‑C‑CACCGCCCCCGTCTTcccccc                                                                                                                                          >  1:1461552/1‑236 (MQ=255)

                                                                                             cGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGG‑GAAAA‑CACACTGCCCCCCTccttcccccccg                                                                                                                                       >  1:757685/1‑228 (MQ=11)

                                                                                              gTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGG‑AACAC‑CACACTGCCACCGTCTTCGCCGCCGGa                                                                                                                                      <  1:615646/239‑1 (MQ=33)

                                                                                              gTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGG‑AACAC‑CACACTGCCACCGTCTTCGCCGCCGGa                                                                                                                                      <  1:644244/239‑1 (MQ=33)

                                                                                              gTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGG‑AACAC‑CACACTGCCACCGTCTTCGCCGCCGGa                                                                                                                                      <  1:794691/239‑1 (MQ=33)

                                                                                               tGCTGATGCCCGCCAGCTGCAGGCTGGGATGTGTGCAGTGTTCCAGGTGGATTTTTCGCGCCAGCGGGTACTCCGGTTCCATCT‑GGTTTTGCATCTGCTCCACGGTCGGGCCGCACGATGCCGGCCAGCGGCGGCTTTCCTCACTCACCCTGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGG‑AACAC‑CACACTGCCACCGTCTTCGCCGCCGGa                                                                                                                                      <  2:1769510/238‑1 (MQ=11)

                                                                                                       aGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGG‑AACAC‑CACACTGCCACCGTCTTCGCCGCCGGACTGCTGaa                                                                                                                              <  2:515192/238‑1 (MQ=31)

                                                                                                          gcCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGG‑AACAC‑CACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGc                                                                                                                          <  1:1469263/239‑1 (MQ=31)

                                                                                                            cAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTCCACCGGGTCCATCT‑GGTGTTGCATCTGCTCCTCGGTCAGGCCGCTCGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGG‑AACAC‑CACACTGCCTCCGTCTTCGCCGCCGGACTGCTGAAGGGCa                                                                                                                         <  1:1619926/238‑1 (MQ=9)

                                                                                                                    ggCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCCCCGGC‑GGG‑AAAAC‑CACACCCCCACCCCCTTCCCCCCCGGCCTGCCGAAgggggtccccccg                                                                                                                 >  1:2131013/1‑228 (MQ=11)

                                                                                                                        gcgATGGTCGCTGTGGTACAGGTGGATTTTTCGTCCCGGCGTGTACCGCGGGTCCCTCT‑GGTTTTGCATCTTCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCCCCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGG‑AACCC‑CACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTg                                                                                                             <  2:172314/238‑1 (MQ=11)

                                                                                                                                cGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGG‑GACAC‑CACACCCCCACCCCCTTCTCCCCCCGCACGCTGGAGGGCACCCGCCCCGCCGCCgcgcgcc                                                                                                    >  1:2025452/1‑236 (MQ=11)

                                                                                                                                  cAGTGGTACCGGGGTATTTTTCGCGCGGGCGTGTCCCCCGGGCCCAGCT‑GGTTTTCCATCTGCCCCCCGGTCGGGCCGCCCGCTGCCAGCCAGCGGCGGCTTTCCTCACTCCCCCGGTCAGCCAGGATTTCCCTTCCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGG‑ACCAC‑CACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCg                                                                                                   <  2:1975960/238‑1 (MQ=11)

                                                                                                                                   cgTGGGACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGG‑AACAC‑CACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGtt                                                                                                 <  2:1136087/238‑1 (MQ=14)

                                                                                                                                       gtgccggTTGATTTTCAGCGCGGACGGGTAAACCGTGTCCATCT‑GGTTTTCCACCTGCTCCACGGTCGGGCCGGACGATGCCAGCCAGCGGCGCCTTTCCCCACCCCCCCGGTCATCCAGGGTTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGG‑AACCC‑CACACTGCCACCGTCTTCGCCGCCGGACTGCTGAATGGCATCCGCCAGGCTGCGGCGCTGCGTTTTc                                                                                              <  2:2092207/233‑1 (MQ=11)

                                                                                                                                       gTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCT‑GGTTTTGCATCTGCTCCTCGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCTCTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGG‑AACAC‑CACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGc                                                                                                                          <  1:2173485/210‑1 (MQ=14)

                                                                                                                                             gTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGG‑AACAC‑CACCCCGCCACCGTCCTCGCCCCCGGACCGCGCAAGGGCACCCCCCACGCTGCGGCGCCGCGCTTTCGCCAGc                                                                                        >  1:1788972/1‑238 (MQ=11)

                                                                                                                                               ggATTTTTCGCGCCGGCGTGTACACCGGGTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGG‑AACAC‑CCCACCGCCCCCGTCTTCCCCCCCGGACTGCTGCAGGGCATCCCCCCCGCTGCGGCCCCCCGGTTTCGCCAGccc                                                                                      >  1:155713/1‑236 (MQ=11)

                                                                                                                                                      tCGCGCCGGCGTGTACCCCGGGTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGG‑AACAC‑CACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTgg                                                                              <  1:304894/239‑1 (MQ=14)

                                                                                                                                                      ccggcccgccgTGTACCCCGGGTCC‑TCTGGGTTTTGCTTCTCCCCCCCGGTCTGGCCGCCCGAGGCCAGCCAGCGGCGGCTTTCCCCCCTCACCCGGTCAGCCAGGATTTCCCTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGG‑AACAC‑CACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTgg                                                                              <  2:1337337/234‑1 (MQ=11)

                                                                                                                                                       gcggcCGGGGTGTCCCCGGGGTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCCCTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGG‑AACAC‑CACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGc                                                                             <  2:2004330/236‑1 (MQ=9)

                                                                                                                                                       cgcgCCGGCGTGTACACCGGGTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGG‑AACAC‑CACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGc                                                                             <  1:1230710/239‑1 (MQ=31)

                                                                                                                                                            cGGCGTGTACCCCGGGTCCATCT‑GGTTTTGCTTCTGCCCCCCGGCCGGGCGGCACGATGCCAGCCAGCGGCGGCTTTCCTCCCCCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGG‑AACAC‑CACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGttt                                                                        <  2:1183232/239‑1 (MQ=11)

                                                                                                                                                            cGGCGTGTACACCGGGTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGG‑AACAC‑CACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGttt                                                                        <  1:1324254/239‑1 (MQ=32)

                                                                                                                                                                gtgtACACCGGGTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGG‑AACAC‑CACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGa                                                                     <  1:327444/238‑1 (MQ=33)

                                                                                                                                                                  gtccacCGGGTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGG‑AACAC‑CACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGAcc                                                                   <  1:870193/235‑1 (MQ=21)

                                                                                                                                                                  gtACACCGGGTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCGCTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGG‑AACAC‑CACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTGCGAcc                                                                   <  1:1402650/238‑1 (MQ=18)

                                                                                                                                                                      acCGGGTCCATCT‑GGTTTTGCATCCCCTCCCCGGCCAGGCCGCCCGATGCCACCCAGCGGCGGCTTTCCTCCCTCACCCGGTCAGCCAGGATGTCCCTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGG‑AACAC‑CACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCtgat                                                              <  2:684640/239‑1 (MQ=11)

                                                                                                                                                                         ggttCCTTCT‑GGTTTTGCGTTGGCCCCCCGGCCAGGCCGCACGATGCCGGCCCGCGGCGGCTTTCCTCCCTCCCCCGGTCGGCCAGGATTTCCCTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGG‑AACAC‑CACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAg                                                           <  2:1377257/236‑1 (MQ=11)

                                                                                                                                                                         gggTCCATCT‑GGTTTTGCATCTGCTCCCCGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGG‑AACAC‑CACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAg                                                           <  1:123182/239‑1 (MQ=25)

                                                                                                                                                                         gggTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGG‑AACAC‑CACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAg                                                           <  1:1234003/239‑1 (MQ=35)

                                                                                                                                                                                         tGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGG‑AACAC‑CACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTccc                                            <  1:236955/239‑1 (MQ=35)

                                                                                                                                                                                            aTCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGG‑AACAC‑CACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCgg                                         <  1:1650327/239‑1 (MQ=35)

                                                                                                                                                                                                    cACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGG‑AACAC‑CACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATaaa                                 <  1:996360/239‑1 (MQ=35)

                                                                                                                                                                                                    cACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGG‑AACAC‑CACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATaaa                                 <  1:510828/239‑1 (MQ=35)

                                                                                                                                                                                                         tCAGGCCGCCCGATGCCAGCCAGCGGCGGCTTTCCTCCCTCCCCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGG‑AACAC‑CACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGt                             <  1:1126035/238‑1 (MQ=18)

                                                                                                                                                                                                         tCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCTCTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGG‑AACAC‑CACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGt                             <  1:710652/238‑1 (MQ=18)

                                                                                                                                                                                                          cAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGG‑AACAC‑CACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCt                           <  1:296166/239‑1 (MQ=33)

                                                                                                                                                                                                              ccGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCCCCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGG‑AACAC‑CACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTg                          <  1:382743/236‑1 (MQ=18)

                                                                                                                                                                                                                   cgaacccggccagcGGCGGCTTTACTCCCTCGCCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGG‑AACGC‑CACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGc                  <  2:1535771/231‑1 (MQ=255)

                                                                                                                                                                                                                        cccgccagcGGCGGCTTCCCTCCCTCACCCGGTCAGCCAGGGTTTCCCTTTCCGGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGG‑AACAC‑CACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGCTCCTg             <  2:302303/236‑1 (MQ=255)

                                                                                                                                                                                                                         cagcccgcGGCGGCTTTCCTCCCTCACCCGGTCGGCCGGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCGGCACCGGC‑GGG‑AACAC‑CACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGCTCCTg             <  2:1382712/238‑1 (MQ=255)

                                                                                                                                                                                                                         cagccagcGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGG‑AACAC‑CACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGCTCCTg             <  1:1318202/238‑1 (MQ=37)

                                                                                                                                                                                                                         cagccagcGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGG‑AACAC‑CACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGCTCCTg             <  1:1749941/238‑1 (MQ=37)

                                                                                                                                                                                                                          ggccagcGGCGGCTGTCCTCACCCCCCCGGGCGGCCAGACTTCCCCCTCCCGGCCGGTCGCGCGTCTGCACCCGCACCGGC‑GGG‑AACAC‑CACACTGCCCCCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTGTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGGGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGCTCCTGTc           <  2:19538/238‑1 (MQ=11)

                                                                                                                                                                                                                           ccccgcggcgGCTTTCCTCCCCCACCCGGTCGGCCAGGATTCCCCTTCCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGG‑AACAC‑CACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGCTCCTGTCg          <  2:50989/235‑1 (MQ=255)

                                                                                                                                                                                                                             cagcGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGG‑AACAC‑CACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGCTCCTGTc           <  1:527324/236‑1 (MQ=37)

                                                                                                                                                                                                                                   cgggTTACCTCCCTCCCCCGGCCAGCCAGGAGTTCCCTTCCCAGCCGGTCGAGCGTCTGCACCGGCACCGGC‑GGG‑AACAC‑CACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGCTCCTGTCGTTCTGTa   <  2:2131013/238‑1 (MQ=11)

                                                                                                                                                                                                                                    ggcTTTCCTCACTCCCCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGG‑AACAC‑CACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGCTCCTGTCGTTCTGTAt  <  1:805148/238‑1 (MQ=37)

CAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCT‑GGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGC‑GGG‑AACAC‑CACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGCTCCTGTCGTTCTGTAT  >  NC_007779/4017651‑4018115

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