BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_007779	4,017,912	A→C	13.8%	S1104R (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,912	0	A	C	13.8%	195.1 / 5.3	110	S1104R (AGT→AGG)	Y75_RS19860	rhs element protein RhsB
Reads supporting (aligned to +/- strand): ref base A (28/64); new base C (3/14); total (32/78)
Fisher's exact test for biased strand distribution p-value = 3.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 111 total aligned reads displayed.

GGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCA‑CCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGCTCCTGTCGTTCTGTATCGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAG  >  NC_007779/4017679‑4018149

gggTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCGCCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACttt                                                                                                                                                                                                                                             <  1:1419005/238‑1 (MQ=18)

   ttCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCCCACCCCCCCCGGCAGCCCGGGTTTCACTTTTCa                                                                                                                                                                                                                                          >  2:1997453/1‑238 (MQ=11)

      tctTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGGCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGGACACCGGGTCCATCTGGGTTTGCATCTGCTCCACGGGCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTGACTTTCCAGCCg                                                                                                                                                                                                                                      >  1:1483280/1‑239 (MQ=255)

           attcaGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGa                                                                                                                                                                                                                                 <  2:1488446/239‑1 (MQ=32)

           attcaGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGa                                                                                                                                                                                                                                 <  2:422381/239‑1 (MQ=32)

            ttcaGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGCCAGCCAGGATTTCACTTTCCCGCCCGTCGAg                                                                                                                                                                                                                                >  1:569514/1‑239 (MQ=14)

             tcaGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCCCCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGc                                                                                                                                                                                                                               >  1:146764/1‑239 (MQ=14)

              cagcagGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACCCACCCGGTCCGCCCGGAGTTCACCTTCCAGCCGGTCgcgc                                                                                                                                                                                                                               >  2:231100/1‑238 (MQ=9)

                gcagGTTGCCCCATTCATCATATTCTGCGCTCCTCGCTTTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGTTGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATTCCAGCCAGCGGCTGCTTTCCTCACTCCCCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATc                                                                                                                                                                                                                            <  1:2073593/239‑1 (MQ=255)

                gcagGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCCGTGGTACAGGTGGATTTTTCGCGCCGGCGGGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCCCGATGCCAGCCAGCGGCGGGCTTCCTCCCTCCCCCGGTCAGCCCGGCGTTCACCTTTCCGCCGGTCGAGccc                                                                                                                                                                                                                             >  2:1440994/1‑236 (MQ=11)

                    gTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCCTCTGGTTTTGCATCTTCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCCTTCCCCCCACACCCCGTCTGCCCGGGTTTCCCTTTCCCGCCCGTCGCGCATcctc                                                                                                                                                                                                                         >  2:1276585/1‑235 (MQ=11)

                    gTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCCCAGTTGCACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGACCGCCAGCGGCGGCCTTTCTCCCCCCCCCGGGCACCCACGATTTCACTATTCAGGTGGTCGAGccttctc                                                                                                                                                                                                                         >  2:271063/1‑232 (MQ=11)

                       gCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGGCTTCCCCCCCCCCCCCGCCAGCCAGGGTTTTACTTTTCAGCCCGTCGCGCATCTGCGcc                                                                                                                                                                                                                      >  2:421936/1‑238 (MQ=11)

                          ccgttCATCATATTCGGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCGGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACCCCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCTCTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCa                                                                                                                                                                                                                  <  1:126337/236‑1 (MQ=16)

                              tcatcaTATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCCGCGTGTACCCCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCCCGACGCCCGCCAGCGGCGGCTTTCCTCACTCACCCCGCCAACCAGGGGTTCCCTTTTCACCCGGTCCAGCATCCCCACCAGGA‑ccc                                                                                                                                                                                                              >  2:1243257/1‑237 (MQ=11)

                               catcatTTTCTGCGCCCCACGCTGTTGCCCCTTCCGTGCTGATGGGCGCCGGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCA‑Ccggc                                                                                                                                                                                                            <  2:403890/239‑1 (MQ=255)

                               catcatATTCTGCGCGCCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCA‑Ccggc                                                                                                                                                                                                            <  1:977436/239‑1 (MQ=25)

                               catcatATTCTGCCCACCCCGCTGTTGCCCCTTCCGTGCGGATGAGCGCCAGCGGCGGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCCCTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCA‑Ccggc                                                                                                                                                                                                            <  2:1009630/239‑1 (MQ=255)

                                 tcatTTTCTGCGCACCGCGCTGTTGCCCCTTCCGGGCGGAGGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTGTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCA‑Ccggc‑gg                                                                                                                                                                                                         <  2:1697595/239‑1 (MQ=255)

                                 tcatATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCA‑Ccggc‑gg                                                                                                                                                                                                         <  2:1456125/239‑1 (MQ=37)

                                    tttgcTGCGCCCCGCGCTGTTCCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCGGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCA‑CCGGC‑GGGaa                                                                                                                                                                                                      <  2:63512/235‑1 (MQ=255)

                                    tatTCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCA‑CCGGC‑GGGaa                                                                                                                                                                                                      <  1:267713/239‑1 (MQ=37)

                                       tCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCA‑CCGGC‑GGGAacac                                                                                                                                                                                                   <  1:1904223/239‑1 (MQ=37)

                                       tCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCA‑CCGGC‑GGGAaca                                                                                                                                                                                                    >  1:1216142/1‑238 (MQ=37)

                                       tCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCCTTCCTCACTCCCCCGGTCAGCCAGGATTTCACTTTCCAGCCCGTCCAGCCTCTCCCCCAGCA‑CCGCC‑ggggaaa                                                                                                                                                                                                    >  2:664659/1‑234 (MQ=11)

                                         tGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCA‑CCGGC‑GGGAAcacca                                                                                                                                                                                                 <  1:52122/239‑1 (MQ=37)

                                         ggcgcACCCCGCTGTTGCCCCTTCCGTGCGGAGGAGCGCCCGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACGCCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCA‑CCGGC‑GGGAAcacca                                                                                                                                                                                                 <  2:1422357/238‑1 (MQ=255)

                                         gccgcACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCA‑CCGGC‑GGGAccacca                                                                                                                                                                                                 <  1:1997453/237‑1 (MQ=37)

                                          gcgcACCACGCTGTTGCCCCTTCCGTGCTGATGAGTGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCCCCCCCACCCGGCCAGCCCGGATTTCCCCTTTCAGCCGGTCGAGCATCCGCACCAGCA‑CCGGC‑gggggaacca                                                                                                                                                                                                 >  2:654028/1‑231 (MQ=11)

                                          gcgcACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCCCTCACCCGGTCCGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCCCCAGCA‑CCGGC‑ggggaacccc                                                                                                                                                                                                 >  1:934925/1‑231 (MQ=255)

                                           cgcACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCCATCCCCCGGGCAGCCAGGGTTTTACTTTCCCACCGGGCGGGCCTCTGCACCCGCACCCGGG‑GGGAAcaccc                                                                                                                                                                                                 >  2:52122/1‑237 (MQ=11)

                                           cgcACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCCGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCA‑CCGGC‑GGGAAcaccac                                                                                                                                                                                                >  1:1975327/1‑238 (MQ=37)

                                           cgcACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATc                                                                                                                                                                                                                            <  2:1934932/212‑1 (MQ=31)

                                           cgcACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCA‑CCGGC‑GGGAAccccac                                                                                                                                                                                                >  1:1408665/1‑238 (MQ=37)

                                           cgcACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGGTTTGCATCTGCTCCCCGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCCCCCCCCCCGGCAGCCCGGATTTCACTTTCCAGCCCGTCGAGCATCTGCACCAGCC‑CCGGC‑gggcaaccccc                                                                                                                                                                                                >  2:2125517/1‑230 (MQ=11)

                                           cgcACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGGGTACACCGGGTCCATCTGGTTTTGCATCTGCGCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCC‑ACGGC‑GGGGAcaccac                                                                                                                                                                                                >  1:707856/1‑238 (MQ=255)

                                                cccGCTGTTGCCCCTTCCGTGCTGAGGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCCCCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGATCATCTGCACCAGCA‑CCGGC‑GGGAACACCACACTGcc                                                                                                                                                                                          <  1:2087093/237‑1 (MQ=18)

                                                cacGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCa                                                                                                                                                                                                                  <  1:1973116/217‑1 (MQ=33)

                                                cacGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCA‑cc                                                                                                                                                                                                               <  1:759557/219‑1 (MQ=34)

                                                 acGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCA‑CCGGC‑GGGAACACCACACTGCCa                                                                                                                                                                                         <  2:773189/239‑1 (MQ=35)

                                                  cGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGCTGGTCGCAGGGGTAAAGGTGGCTTTTTCGCGCCGGCGTGTACACCGGGTCCCTCTGGTTTTGCATCTGCCCCACGGTCAGGCCGCACGACGCCAGCCCGCGGCGGCCTTCCTCCCTCACCCCGCCAGCCACGATTTCACCTTTCCGCCGGGCGCGCATCTGCCCCCCGC‑CCGGCGGGGAAcaccaacccccc                                                                                                                                                                                          >  2:941482/1‑231 (MQ=11)

                                                  cGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTTCTCACTCACCCGGTCAGCCAGGGTTTTACTTTCCAGCCGGTCGAGCATCTGCCCCAGCA‑CCGGC‑GGGGACAACccacccccc                                                                                                                                                                                         >  1:1891447/1‑233 (MQ=11)

                                                             cgtCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCA‑CCGGC‑GGGAACACCACACTGCCACCGTCTTcgccg                                                                                                                                                                             <  1:1120735/237‑1 (MQ=18)

                                                               tCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGGCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGGTTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCA‑CCGGC‑GGGGACACCCCACCCCCACCGCCTTcgcccc                                                                                                                                                                            >  1:169481/1‑236 (MQ=255)

                                                                ccGGGCTGATGGGCGCCAGGGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTCCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGTTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCA‑CCGGC‑GGGAACACCACACTGCCACCGTCTTCGccgccg                                                                                                                                                                          <  2:1211141/239‑1 (MQ=16)

                                                                 cgggCTGAGGAGCGCCCGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCA‑CCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCgg                                                                                                                                                                         <  2:1992424/236‑1 (MQ=17)

                                                                  gTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCA‑CCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGa                                                                                                                                                                        <  1:353778/239‑1 (MQ=33)

                                                                  gTGCTGATGAGCGCCAGCGGCAGGCCGCGAGGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCA‑CCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGa                                                                                                                                                                        <  2:788917/239‑1 (MQ=18)

                                                                      tgatgaGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCCCCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCA‑CCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGActg                                                                                                                                                                     <  1:1776075/238‑1 (MQ=17)

                                                                        atgaGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCCCTTTCCAGCCGGTCGAGCATCTGCACCAGCA‑CCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGa                                                                                                                                                                        <  1:786499/233‑1 (MQ=17)

                                                                                     gcggggcggggtggtggCAGTGGTACGGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTGGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGGGCATCTGCACCAGCA‑CCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGc                                                                                                                                                      <  2:12893/227‑1 (MQ=255)

                                                                                           cgcgATGGTCGCAGTGGTACAGGTGTATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCCCGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCCCTTTCCAGCCGGTCGAGCATCTGCTCCAGCA‑CCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTg                                                                                                                                               <  1:2223150/239‑1 (MQ=9)

                                                                                              gATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCA‑CCGGC‑GGGGACACCACAACGCCCCCCTCCTCGCCCCCGGACAGCTGGAGGGGATCCGCCCGGCCgcgg                                                                                                                                            >  1:560395/1‑239 (MQ=11)

                                                                                                            tACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCA‑CCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCt                                                                                                                                        <  1:1983001/229‑1 (MQ=31)

                                                                                                             aCAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCA‑CCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGc                                                                                                                              <  1:1737372/238‑1 (MQ=31)

                                                                                                                      tttttCGCGCCGGCTTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCA‑CCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCg                                                                                                                    <  1:1848567/239‑1 (MQ=14)

                                                                                                                          tCGCGCCGGCGTGTCCACCGGGTCCATCTGGTTTTGCCTCTGCTCCACGGTCAGGCCGCCCGATGCCAGCCAGCGGCGGCTTTCCTCACTCCCCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCA‑CCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTgg                                                                                                                <  1:953516/239‑1 (MQ=9)

                                                                                                                          tCGCGCCGGCGTGTACCCCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCA‑CCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTgg                                                                                                                <  1:355989/239‑1 (MQ=14)

                                                                                                                          tCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCA‑CCGGC‑GGGAACTCCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTgg                                                                                                                <  1:1956900/239‑1 (MQ=14)

                                                                                                                           cgccCCGGGGGGTACACCGGGTCCATCTGGTTTTGCATCTCCTCCCCGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCCCCCGGTCAGCCAGGATTTCCCTTTCCAGCCGGTCGAGCATCTGCACCAGCA‑CCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGc                                                                                                               <  2:569514/239‑1 (MQ=255)

                                                                                                                                cggcggGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCCCCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCA‑CCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGttt                                                                                                          <  2:1975327/239‑1 (MQ=17)

                                                                                                                                cGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCCCGATGCCAGCCAGCGGCGGCTTTCCTCCCTCACCCGGTCAGCCAGGATTTCCCTTTCCAGCCGGTCGAGCATCTGCACCAGCA‑CCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGttt                                                                                                          <  1:231100/239‑1 (MQ=17)

                                                                                                                                cGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCA‑CCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGttt                                                                                                          <  2:1408665/239‑1 (MQ=32)

                                                                                                                                   cGTGTACCCCGGGTCCTTCTGGTTTTGCATCTGCTCCCCGGTCAGGCCGCCGGAGGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCA‑CCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCg                                                                                                        <  2:1216142/238‑1 (MQ=255)

                                                                                                                                    gtgtACCCCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCA‑CCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGa                                                                                                       <  1:894042/238‑1 (MQ=18)

                                                                                                                                     tgtACACCGGGTCCATCTTGTGTTGCATCTGCTCCCCGGTCAGTCCGCACGAGGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCA‑CCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGAcc                                                                                                     <  1:1901320/239‑1 (MQ=255)

                                                                                                                                      gtACACCGGGTCCATCTGGTTTTGCCTCTGCTCCACGGTCAGGCCGCACGATGCCATCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCA‑CCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCAGCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGAcc                                                                                                     <  1:1907265/238‑1 (MQ=21)

                                                                                                                                            cGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCA‑CCGGC‑GGGGACACCCCACCGCCCCCCCCCTCTCCCCCCGCCCGCTGAAAGGGATCCCCCCGGCCGCGGCGCGGCGTTTTTCCCCCCCCACAGGGGGCGGTTTCGcccccctca                                                                                               >  2:1967945/1‑233 (MQ=11)

                                                                                                                                             tgttCCATTTGGTTTTGCATCTGCTCCCCGGTCAGGCCGCACGATGCCCGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTCCCCTTTCCAGCCGGTCGGGCATCTGCACCAGCA‑CCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAg                                                                                             <  1:986148/236‑1 (MQ=255)

                                                                                                                                             gggTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCCGCCAGCGGCGGCTTTCCTCACTCCCCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCA‑CCGGC‑GGGAACACCCCACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAg                                                                                             <  2:146764/239‑1 (MQ=25)

                                                                                                                                             gggTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCCCTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCA‑CCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAg                                                                                             <  1:279922/239‑1 (MQ=25)

                                                                                                                                             gggTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCA‑CCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAg                                                                                             <  1:487883/239‑1 (MQ=35)

                                                                                                                                                 ccATCTGGTTTTGCATCGGCGCCGCGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACCCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCA‑CCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGc                                                                                         <  1:1243257/239‑1 (MQ=255)

                                                                                                                                                    tCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCA‑CCGGC‑GGGGACACCCCACCGCCCCCCCCCTCCCCCCCCGGCCGCTGGAGGGGATCCGCCAGGCCGCGGCGCCGGGGTTTCCCCACCCCACCCGGGGCAGGTTCGCCCCCGATGAGTGGCggg                                                                                      >  2:91022/1‑237 (MQ=11)

                                                                                                                                                     cTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCCCTTTCCAGCCGGTCGAGCATCTGCACCAGCA‑CCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGa                                                                                     <  1:1253182/239‑1 (MQ=25)

                                                                                                                                                     cTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCA‑CCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGa                                                                                     <  1:654028/239‑1 (MQ=35)

                                                                                                                                                             gCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCA‑CCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTcccc                                                                             <  1:1542919/239‑1 (MQ=35)

                                                                                                                                                              cATCTGCTCCACGGGCAGGCCGCACGATGCCGGCCAGCGGCGGCTTTCCTCACTCCCCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCA‑CCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTcccc                                                                             <  1:941482/238‑1 (MQ=25)

                                                                                                                                                               aTCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCA‑CCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCgg                                                                           <  1:41368/239‑1 (MQ=35)

                                                                                                                                                                       ccgggTCAGGCCGACCGAGGCCGGCCAGCGGTGGCTTTCCTCCCTCACCCGTTCAGCCGGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCA‑CCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATaaa                                                                   <  2:735458/236‑1 (MQ=11)

                                                                                                                                                                             cAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCA‑CCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCt                                                             <  1:1331552/239‑1 (MQ=33)

                                                                                                                                                                                  cGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGCTTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCA‑CCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGAt                                                         <  1:1787100/238‑1 (MQ=18)

                                                                                                                                                                                         tGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCA‑CCGGC‑GGGAACACCACCCTGCCCCCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGCTcc                                                 <  1:2033295/239‑1 (MQ=25)

                                                                                                                                                                                         tGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCA‑CCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTgg                                                                                                                <  1:1427099/176‑1 (MQ=31)

                                                                                                                                                                                           ccggccagcGGCGGCTTTCCTCCCTCACCCGGTCAGCCAGGATTTCCCTTTCCAGCCGGTCGAGCATCTGCACCAGCA‑CCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGCTCCTg                                               <  2:2132193/236‑1 (MQ=255)

                                                                                                                                                                                             cgccatcGGCGTCTTTCCTCGCTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGGGCATCTGCACCAGCA‑CCGGC‑GGGAACACCGCACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGCTCCTGTc                                             <  1:915061/238‑1 (MQ=255)

                                                                                                                                                                                                cggcggcgGCTTTCCTCACTCACCCGGTCAGCCAGGATTCCACTTTCCAGCCGGTCGCGCATCTGCACCAGCA‑CCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGCTCCTGTc                                             <  1:1471619/234‑1 (MQ=37)

                                                                                                                                                                                                cagcGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACAAGCA‑CCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGCTCCTGTc                                             <  1:1440994/236‑1 (MQ=37)

                                                                                                                                                                                                      cggcTTTCCTCCCTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCA‑CCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGCTCCTGTCGTTCTGTa                                     <  1:1401245/238‑1 (MQ=37)

                                                                                                                                                                                                           ttCCTCCCTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCCGCA‑CCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGCTc                                                  <  1:1477539/220‑1 (MQ=25)

                                                                                                                                                                                                             ccTCACTCACCCGGTCGCCCAGGATTTCCCTTTCCAGCCGGTCGAGCATCTGCACCAGCA‑CCGGC‑GGGACCACCACCCTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGCTCCTg                                               <  1:1546449/221‑1 (MQ=255)

                                                                                                                                                                                                                       ccGCTCAGCCAGGATTTCCCTTTCCAGCCGGTCGAGCATCTGCACCAGCG‑CCGGC‑GGGAACACCCCACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGCTCCTGTCGTTCTGTATCGTGGTCAGCCGGTCGc                   <  1:1276585/239‑1 (MQ=255)

                                                                                                                                                                                                                           tcccccggtttttCCCTTCCCCGCCGGCCGACCATCTGCCCCAGCC‑CCGGC‑GGGACCCCCACACTGCCCCCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGCTCCTGTCGTTCTGTATCGTGGTCAGCCGGTCGCCGTc               <  2:506063/229‑1 (MQ=11)

                                                                                                                                                                                                                           tCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCGCCAGCA‑CCGGC‑GGGAACACCCCACTGCCCCCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGCTCCTGTCGTTCTGTATCGTGGTCAGCCGGTCGCCGt                <  1:2073490/238‑1 (MQ=255)

                                                                                                                                                                                                                            cAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCA‑CCGGC‑GGGAACACCCCACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGCTCCTGTCGTTCTGTATCGTGGTCAGCCGGTCGCCGTc               <  1:664659/238‑1 (MQ=38)

                                                                                                                                                                                                                               ccgggATCTCACTGCCCAGCCGGTCGGGCATCTGCACCGGCA‑CCGGC‑GGGAACCCCACCCTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGCTCCTGTCGTTCTGTATCGTGGTCAGCCGGTCGCCGTCCCAg           <  2:169481/236‑1 (MQ=11)

                                                                                                                                                                                                                               ccAGGCTTTCCCTGTCCAGCCGGTCGAGCATCTGCACCAGCT‑CCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGCTCCTGTCGTTCTGTATCGTGGTCAGCCGGTCGCCGTCCCAg           <  1:1102734/239‑1 (MQ=255)

                                                                                                                                                                                                                               ccAGGATTTCACTTTCCAGCCGGTCGAGCATCGGCACCAGCA‑CCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGCTCCTg                                               <  1:785463/203‑1 (MQ=38)

                                                                                                                                                                                                                                     tttCACTTTCCAGCCGGTCGAGCATCTGCACCAGCA‑CCGGC‑GGGAACACCCCACTGCCACCGTCTTCGCCCCCCGACCGCTGAAGGGCATCCCCCACGCTGCCGCGCCGCGTTTTCCCCCGCCCCCCGGGGGCAGGTTCGACCCCGATGAGTGGCGTGGAGCCCCCCCGCTGGTAAAACGTCCGGGTGGGGCCCCCTGCGGTCTGTTTCCTGGGCCGCCGGGCGCCGTCCCCGCCCTa      >  2:1966670/1‑238 (MQ=2)

                                                                                                                                                                                                                                         aCTTTCCAGCCGGTCGAGCATCTGCACCAGCA‑CCGGC‑GGGAACCCCCCACTGCCACCGGCCTCGCCCCCCGGCCGCTGAAGGGCATCCCCCAGGCCGCGGCGCCGCGTTTTTGCCCCGTCACCGGGGGCCGGTTCGACCCCGCTGAGTGGGGTGAACCTCCCCGGCGGATAAATCGCCCTGCTGGGGCTCCCGTCGTTTTGTTTCGTGGGCCGCCGGGCGCCGTCCCAGcccctcccg  >  2:165244/1‑231 (MQ=2)

GGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGTGCTGATGAGCGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCA‑CCGGC‑GGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGGGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCAGTTTCGACCCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGCTCCTGTCGTTCTGTATCGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAG  >  NC_007779/4017679‑4018149

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