BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_007779	3,874,921	A→C	13.3%	S1104R (AGT→AGG)	Y75_RS19225 ←	RHS element protein RhsA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007779	3,874,921	0	A	C	13.3%	237.0 / 2.8	125	S1104R (AGT→AGG)	Y75_RS19225	RHS element protein RhsA
Reads supporting (aligned to +/- strand): ref base A (24/80); new base C (4/17); total (28/97)
Fisher's exact test for biased strand distribution p-value = 7.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 126 total aligned reads displayed.

GATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCGTATTCTGCGCACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAG  >  NC_007779/3874683‑3875158

gATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCGTATTCTGCGCACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCa                                                                                                                                                                                                                                               <  1:1856633/239‑1 (MQ=37)

           tctTCATTCAGCAGGTTGCCCCATTCATCGTATTCTGCGCACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCCGTGGTACAGGTGGCTTTTTCGCGCCGGCGAGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGAACGATGCCAGCCAGCGGCGGCTTTCCTCACACACCCGGTCCGCCAGGACTTGACTTTCCAGCCg                                                                                                                                                                                                                                    >  2:1063888/1‑239 (MQ=11)

                attcaGCAGGTTGCCCCATTCATCGTATTCTGCGCGCCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGa                                                                                                                                                                                                                               <  1:1009562/239‑1 (MQ=37)

                    agcagGTTGCCCCATTCATCGTATTCTGCGCACCACTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCCGGATTTCCCTTTCCAGCCGGTCggtca                                                                                                                                                                                                                            >  2:656898/1‑234 (MQ=255)

                     gcgggTTGCCCCATTCAGCGTATTCTGCGCCCCATTCTTTTGTCCCTTCCTGGCTGATAAGGGCCAGCGCCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACCCCGGGTCCATCTGGTTTTGCATCTGCTCCCCGGTCAGGCCGCACGAGGCCAGCCAGCGGCGGCTTTCCTCACTCCCCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATc                                                                                                                                                                                                                          <  2:312817/236‑1 (MQ=11)

                     gcagGTTGCCCCATTCATCGTATTCTGCGCACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCTTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCa                                                                                                                                                                                                                            >  1:1092451/1‑237 (MQ=37)

                      cagGTTGCCCCATTCATCGTATTCTGCGCACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCCGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGGTTTCACTTTCCAGCCGGTCGAGCATcc                                                                                                                                                                                                                         >  1:2237831/1‑238 (MQ=37)

                           tGCCCCATTCATCGTATTCTGCGCACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACCCACCCGGTCCGCCAGGATTTCACTTTCCAGCCGGGCGAGCATCTGCAcc                                                                                                                                                                                                                    >  1:1487173/1‑239 (MQ=255)

                              cccATTCATCGTATTCTGCGCACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGA‑TTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCCCCCGGTCAGCCCGGAGTTCACCTTCCCGCCCGTCGGGCCACCGCACCCGc                                                                                                                                                                                                                 >  2:526345/1‑238 (MQ=11)

                               ccATTCATCGTATTCTGCGCACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCa                                                                                                                                                                                                                <  1:1742028/239‑1 (MQ=37)

                                 aTTCATCGTATTCTGCGCACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCAcc                                                                                                                                                                                                              <  1:1779773/239‑1 (MQ=37)

                                    cATCGTATTCTGCGCACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACcggc                                                                                                                                                                                                           <  1:562072/239‑1 (MQ=37)

                                    cATCGTATTCTGCGCACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACcggc                                                                                                                                                                                                           <  1:2113120/239‑1 (MQ=37)

                                      tCGTATTCTGCGCACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACcggc                                                                                                                                                                                                           <  1:1889817/237‑1 (MQ=37)

                                          aTTCTGCGCACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGaa                                                                                                                                                                                                      <  1:2084861/238‑1 (MQ=37)

                                          aTTCTGCGCACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGaa                                                                                                                                                                                                      <  1:162543/238‑1 (MQ=37)

                                                cgcACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGGTTTCCTCACTCACCCGGTCAGCCAGGGTTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGGAcacccc                                                                                                                                                                                                >  1:1645873/1‑236 (MQ=255)

                                                cgcACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTTCTCACCCCCCCCGTCCGCCAGGGTTTTACTTTTCCGGCCGTCGAGCAACTGCACCTGCACCGGCGGGAAaacccc                                                                                                                                                                                                >  2:2084861/1‑236 (MQ=11)

                                                cgcACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCCCCCGGTCAGCCAGGGATTCCCTTTCCAGCCGGGTCAGCCTCCGCACCCGCCACCGCggggggacacc                                                                                                                                                                                                >  2:142179/1‑230 (MQ=11)

                                                cgcACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGGTTTCACTTTCCAGCCGGTCGAGCAACCGCCCCAGCACCGGCGGGGAaaccac                                                                                                                                                                                                >  2:2266653/1‑238 (MQ=255)

                                                cgcACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCCGCCCCCGCACCCGCGGGagacacac                                                                                                                                                                                                >  2:1179502/1‑231 (MQ=255)

                                                cgcACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGGCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTTCTCACTCACCCGGTCAGCCCGGATTTCACCTTCCCGCCGGTCGGGCATCCTCCCCCGCCACCGCGGGGAcaccca                                                                                                                                                                                                >  2:1419698/1‑236 (MQ=11)

                                                cgcACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGGTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCCCCCTCACCCGGGCAGCCCGGATTTCACTTTCCCGCCCGGCGAGCCTCTCCCCCAGCCCCCGCGGGGAaacccc                                                                                                                                                                                                >  2:1889817/1‑236 (MQ=11)

                                                cgcACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGCTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCCCCCTCACCCGGTCCGCCAGGATTTCACCTTCCAGCCGGCCGAGCATCCGCACCAGCCCcggcggtgaccacac                                                                                                                                                                                                >  1:231959/1‑229 (MQ=11)

                                                cgcACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGGGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCCGCCAGGATTTCACCTTCCAGCCGGTCGGACATCCTCCCCAGCAACCGCGGGGaccacca                                                                                                                                                                                                >  2:2177256/1‑233 (MQ=11)

                                                 gcACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACcaca                                                                                                                                                                                               <  1:853685/238‑1 (MQ=37)

                                                   aCCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCCGCCGGTCGAGCATCTGCACCAGCACCCGCGGGGACACCccacc                                                                                                                                                                                             >  1:1271583/1‑237 (MQ=255)

                                                     cATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCAc                                                                                                                                                                                                                     <  1:331898/212‑1 (MQ=38)

                                                     cATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGaa                                                                                                                                                                                                      <  1:515112/227‑1 (MQ=37)

                                                      aTTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCa                                                                                                                                                                                         <  1:705402/239‑1 (MQ=37)

                                                        tCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCAcc                                                                                                                                                                                       <  1:2177256/239‑1 (MQ=35)

                                                          tgttgtCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTCCACCCGCACCGGCGGGGACACCCAACCGCCAccc                                                                                                                                                                                      >  1:1967175/1‑237 (MQ=255)

                                                            ttgtCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAg                                                                                                                                                                                                                                       <  1:401597/187‑1 (MQ=37)

                                                               tCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAAAACCACACCGCcccccgctgc                                                                                                                                                                                 >  1:554448/1‑232 (MQ=255)

                                                                 ccttccttGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCCCCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACAACACACTCCCCCCctcctccc                                                                                                                                                                               >  1:283896/1‑236 (MQ=255)

                                                                 ccttccttGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGGACACCCCACTGCcccccgcctcgc                                                                                                                                                                               >  1:138087/1‑230 (MQ=255)

                                                                   gtccttGCTGATAAGGGCCCGGGGCAGGCCGCGATGGTCGCGGTGGTGCAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTcgccgc                                                                                                                                                                            <  2:1767016/238‑1 (MQ=16)

                                                                     ccttGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCCCCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGccgccg                                                                                                                                                                          <  1:135388/239‑1 (MQ=18)

                                                                      cttGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCCGGGTTTCCCTTTCCAGCCGGTCgcgc                                                                                                                                                                                                                             >  2:401597/1‑187 (MQ=255)

                                                                         gCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAg                                                                                                                                                                                                                                       <  1:1635717/174‑1 (MQ=33)

                                                                         gCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAg                                                                                                                                                                                                                                       >  2:1635717/1‑174 (MQ=33)

                                                                              tAAGGGCCAGCGGCAGGCCGCGAGGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGActgctg                                                                                                                                                                  <  1:2167052/238‑1 (MQ=17)

                                                                                    ccggCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACGCCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAgcgc                                                                                                                                                            <  1:855296/235‑1 (MQ=14)

                                                                                    ccAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAgcgc                                                                                                                                                            <  1:1395299/238‑1 (MQ=31)

                                                                                         ggccggccGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTTCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCg                                                                                                                                                       <  1:1677461/238‑1 (MQ=14)

                                                                                                      ggTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCCCCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGTCTGCTGAAGCGCATCCGCCAGGCTgcggcg                                                                                                                                          <  2:1017183/238‑1 (MQ=14)

                                                                                                      ggTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGTAACGCCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGgcgc                                                                                                                                         <  1:242423/239‑1 (MQ=14)

                                                                                                          gCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGGCCACCACACCTCCACCGCCCTTGCCCCCGGACTCCTGGAGCGCACCCGCCCCGCCGCGGCGCcgcg                                                                                                                                     >  1:1225876/1‑239 (MQ=11)

                                                                                                                       tGGATTTTTCGTGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAg                                                                                                                           <  1:1849096/236‑1 (MQ=14)

                                                                                                                                cgcgcgGGCGTGTACCCCGGGTCCATCTGGTTTTGCATCTGCTCCCCGGTCGGGCCGCACGATGCCGGCCAGCGCCGGCTTTCCTCCCTCACCCGGTCAGCCAGGATTTCCCTTCCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGc                                                                                                               <  2:1487173/239‑1 (MQ=11)

                                                                                                                                cgcgcgGGCGGGTACACCGGGTCCATCTGGTTTTGCCTCTGCTCCCCGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGc                                                                                                               <  2:1645873/239‑1 (MQ=9)

                                                                                                                                cccgccgcCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGc                                                                                                               <  1:887822/237‑1 (MQ=14)

                                                                                                                                   gcCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCCCGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCCCCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGt                                                                                                            <  1:2021035/239‑1 (MQ=17)

                                                                                                                                     cggggggtACACCGGGTCCATCTGGTTTTGCCTCTGCTCCACGGCCGGGCCGCAGGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCCCTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCGCACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGttt                                                                                                          <  2:231959/233‑1 (MQ=11)

                                                                                                                                         gtgtACACCGGGTCCGTCTGGTGTTGCATATGCTCCACGGTCAGGCCGCGCGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCaa                                                                                                       <  1:1557536/238‑1 (MQ=16)

                                                                                                                                         gtggACCCCGGGTCCATCTGGTTTTGCATCGGCTCCCCGGTCGGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCaa                                                                                                       <  2:1681906/238‑1 (MQ=16)

                                                                                                                                         gggagcacCGGGTCCATCTGGTTTTGCACCTGCGCCCCGGTCGGGCCGCCCGCTGCCAGCCAGCGGCGGCTTTCCCCCCTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCaa                                                                                                       <  2:2034975/233‑1 (MQ=11)

                                                                                                                                           gtACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAAct                                                                                                     <  1:693440/238‑1 (MQ=34)

                                                                                                                                             acacCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCtg                                                                                                  <  1:661662/239‑1 (MQ=35)

                                                                                                                                               acCGGGTCCATCTGGTTTTGCATCTGCTCCACGGCCAGGCCGCACGATGCCAGCCCGCGGCGGCTTTCCTCACTCCCCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGttt                                                                                                          <  1:1479099/229‑1 (MQ=17)

                                                                                                                                                  gggTCCATCTGGTTTTGCATCTGCTCCCCGGTCAGGCCGCCCGATGCCAGCCAGCGGCGGCTTTCCTCCCTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAg                                                                                             <  1:1971681/239‑1 (MQ=25)

                                                                                                                                                  gggTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAg                                                                                             <  1:1179502/239‑1 (MQ=35)

                                                                                                                                                     tCCATCTGGTTTTGCATCTGCCCCACGGTCAGGCCGCACGATGCCAGCCCGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCt                                                                                                                                        <  1:929446/193‑1 (MQ=14)

                                                                                                                                                          cTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGa                                                                                     <  1:1703180/239‑1 (MQ=35)

                                                                                                                                                             gtgttGCATCTGCCCCACGGTAAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAg                                                                                   <  2:1342787/235‑1 (MQ=25)

                                                                                                                                                                  gCATCTCCCCCCCGGCCGGGCCGAACGATGCCAGCCAGCGGCGGCGTTCCTCACTCACCCGGTCAGCCAGGATTTCCCTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTcccc                                                                             <  2:88326/239‑1 (MQ=11)

                                                                                                                                                                    aTCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCtgat                                                                                                <  1:670069/218‑1 (MQ=35)

                                                                                                                                                                    aTCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCgg                                                                           <  1:88531/239‑1 (MQ=35)

                                                                                                                                                                    aTCTGCCCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGCTTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACCCTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCgg                                                                           <  1:569484/239‑1 (MQ=25)

                                                                                                                                                                            cACGGCCGGGCCCCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATaaa                                                                   <  2:279392/239‑1 (MQ=21)

                                                                                                                                                                             aCGGTCAGGCCGCACGATGCAAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCCCTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCGCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAAt                                                                  <  1:2175619/239‑1 (MQ=25)

                                                                                                                                                                                 tCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGt                                                               <  1:1672276/238‑1 (MQ=33)

                                                                                                                                                                                 tCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGt                                                               <  1:1419698/238‑1 (MQ=33)

                                                                                                                                                                                 cctggCCTCCCGATGCCAGCCAGCGGCGCCTTCCCTCCCTCACCCGGTCGGCCGGGATTCCACTTCCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGt                                                               <  2:86571/235‑1 (MQ=11)

                                                                                                                                                                                  cAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTc                                                              <  1:1375930/238‑1 (MQ=33)

                                                                                                                                                                                              tGCCAGCCAGCGGCGGCTTTCCTCACTCCCCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGCCTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTTATGAGTGGCGTGAAGCTCCCCgg                                                                           <  1:976508/213‑1 (MQ=21)

                                                                                                                                                                                              tGCCAGCCAGCGGCGGCTTTCCTCACTCCCCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGtt                                                  <  1:1629730/238‑1 (MQ=18)

                                                                                                                                                                                              tGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCCCTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCaa                                                                                                       <  1:752889/185‑1 (MQ=18)

                                                                                                                                                                                                 cagccagcGGCGGCTTTCCTCCCTCCCCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACGCCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTg                                               <  1:1409527/238‑1 (MQ=18)

                                                                                                                                                                                                 cagccagcGGCGGCTTTCCTCACTCACCCGGTGAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTg                                               <  1:1736093/238‑1 (MQ=18)

                                                                                                                                                                                                     cggcggcgGCTTTCCTCCCCCACCCGGTCAGCCAGGATTTCCCTTCCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCACCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGtt                                          <  2:554448/237‑1 (MQ=255)

                                                                                                                                                                                                     cggcgccgGCTTTCCTCCCTCACCCGGTCAGCCAGGATTTCCCTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACGGCCACCGTCTTCGCCGCCGGACGGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGTCTGATAAATCGTCTGGATGCGGGTTCTGTCGtt                                          <  1:2051367/237‑1 (MQ=255)

                                                                                                                                                                                                     cagcGGCTGCTTCCCCCCCACACCCGGTCAGCCAGCATTTACCCTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGATCACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGtt                                          <  2:154530/239‑1 (MQ=11)

                                                                                                                                                                                                     cagcGGCGGCTTTCCTCCCTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACTCCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTc                                             <  1:2031686/236‑1 (MQ=18)

                                                                                                                                                                                                     cagcGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCCCTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGtt                                          <  2:1894757/239‑1 (MQ=18)

                                                                                                                                                                                                     cagcGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCCCCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGtt                                          <  1:1159997/239‑1 (MQ=18)

                                                                                                                                                                                                     cagcGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTc                                             <  1:859552/236‑1 (MQ=33)

                                                                                                                                                                                                     cagcGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGtt                                          <  1:560413/239‑1 (MQ=33)

                                                                                                                                                                                                         ggcggcTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCt                                        <  1:1044105/237‑1 (MQ=33)

                                                                                                                                                                                                           cggcTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTa                                     <  1:874862/238‑1 (MQ=33)

                                                                                                                                                                                                                          acCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCaa                                                                                                       <  1:654722/157‑1 (MQ=33)

                                                                                                                                                                                                                          acCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCCTCTTCGCCCCCGGACCGCTGAAGCGCCTCCGCCAGGCCGCGGCGCTGCGTTTTCGCCAGCTCACCCGTGGCGGGTTCCACTCCGCTGGGTGGCGTGGATCTCCCCCGCGGGTAAATCGTCCTTGTGCGGGGTTTGTCTGTCTGTCTTTTgggccaccagg                      >  2:703167/1‑229 (MQ=11)

                                                                                                                                                                                                                           cccGGTCAGCCAGGATTCCCCTGTCCAGCCGGTCGAGCATCTGCTCCAGCACCGGCGGGAACACCCCACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCg                    <  1:703167/239‑1 (MQ=255)

                                                                                                                                                                                                                           cccGGTCAGCCAGGATGTCACTTCCCAGCCGGTCGAGCATCTGCACCAGCCCCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCg                    <  1:2105459/239‑1 (MQ=255)

                                                                                                                                                                                                                            ccGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCCCCCCCCGACCGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCCGGGGCGGGTTTCACGCTGGTGAGGGGCGTGAAGCTCCCCGGCTGATAAAAAGTCTTGTTGGGGGGTTTGTCGTTCTGTATTTTGGGCTGcccggcg                    >  2:2105459/1‑233 (MQ=11)

                                                                                                                                                                                                                                tCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCCCTCTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCTTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGt                <  1:1063888/238‑1 (MQ=255)

                                                                                                                                                                                                                                     cAGGATTTCCCTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAg           <  1:526345/238‑1 (MQ=37)

                                                                                                                                                                                                                                     cAGGATTTCCCTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGc          <  1:191269/239‑1 (MQ=37)

                                                                                                                                                                                                                                     cAGGATTTCACTTCCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGc          <  1:656898/239‑1 (MQ=37)

                                                                                                                                                                                                                                             cACTTCCCAGCCGGTCGAGCATCCGCACCGGCCCCGGCGGGACCCCCCCCCGCCCCCCGTCTCCCCCGCCGGCCTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTCCAACTCTGATGGGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAg  <  2:1786102/239‑1 (MQ=11)

GATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCGTATTCTGCGCACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAG  >  NC_007779/3874683‑3875158

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