BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_007779	3,874,921	A→C	12.0%	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	12.0%	663.8 / 11.5	344	S1104R (AGT→AGG)	Y75_RS19225	RHS element protein RhsA
Reads supporting (aligned to +/- strand): ref base A (84/201); new base C (10/45); total (94/250)
Fisher's exact test for biased strand distribution p-value = 1.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 344 total aligned reads displayed.

TGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCGTATTCTGCGCACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGT‑CAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCA‑CCGGTGGCGGTTTCAACTCTGATGAGTGG‑CGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAG  >  NC_007779/3874685‑3875158

      ttCTCTTCATTCAGCAGGTTGCCCCATTCATCGTATTCTGCGCACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCCCCCGGT‑CAGCCAGGATTTCACTTT‑CCa                                                                                                                                                                                                                                           >  1:148142/1‑238 (MQ=37)

       tttctTCATTCAGCGGGTTGCCCCATTCATCGTATTCTGCGCTCCATTCTGTTGTCCCTTCCTTGCTGATATGGGCCAGCGGCAGGCCGCGGTGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCCCGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCCCCCGTT‑CAGCCAGGATTTCACTTT‑CCAGc                                                                                                                                                                                                                                         <  2:1595790/237‑1 (MQ=255)

        ctattccttcaGCATGTTGCCCCATTCATCGTATTCTGCGCGCCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCCCCCGGT‑CAGCCAGGATTTCACTTT‑CCAGcc                                                                                                                                                                                                                                        <  2:1812234/236‑1 (MQ=255)

         tctTCATTCAGCAGGTTGCCCCATTCATCGTATTCTGCGCACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGT‑CAGCCAGGATTTCACTTT‑CCAg                                                                                                                                                                                                                                          <  1:169412/236‑1 (MQ=37)

           ttcattcaGCAGGTTGCCCCATTCATCGTATTCTGCGCACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCCCCCGGT‑CAGCCAGGGTTTCCCTTT‑CCAGCCgg                                                                                                                                                                                                                                      >  1:1062731/1‑238 (MQ=255)

               ttcaGCAGGTTGCCCCATTCATCGTATTCTGCGCACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGAACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGT‑CAGCCAGGATTTCACTTT‑ACAGCCGGTCGa                                                                                                                                                                                                                                  >  1:544750/1‑238 (MQ=37)

                  agcagGTTGCCCCATTCATCGTATTCTGCGCACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGT‑CAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCa                                                                                                                                                                                                                               >  1:307296/1‑238 (MQ=37)

                  agcagGTTGCCCCATTCATCGTATTCTGCGCACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCCTTTCTCC‑CTCCCCCCGG‑CCGCCAGGGTTTCCCTTT‑CCCGCCGGTCgctca                                                                                                                                                                                                                               >  2:1407857/1‑234 (MQ=11)

                   gcagGTTGCCCCATTCATCGTATTCTGCGCACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACCCCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCCCACCTCACCCGG‑‑CCGCCAGGGTTTCACTTT‑Cccgccgtgctcttat                                                                                                                                                                                                                              >  2:671621/1‑229 (MQ=11)

                   gcagGTTGCCCCATTCATCGTATTCTGCGCACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGT‑CAGCCAGGATTTCACTTT‑CCAGCCGGTCGCGCa                                                                                                                                                                                                                               >  1:1242624/1‑237 (MQ=37)

                    cagGTGGCCCCATTCATCGTTGTGTGCACACCAGTCGTGTGTCCTTGCCTGCCGGATAAGGGCCAGCGGCAGGTCGAGAGGGTCGCAGTGGGACAGGTGGATTTTTCGGGCCGTCGTGTACACCGGGTCCATCTGGTTGTGCATCTTCTCCGCGGTCAGGCGTCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGTT‑CAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCAta                                                                                                                                                                                                                             <  2:2069811/238‑2 (MQ=11)

                       gTTGCCCCATTCATCGTATTCTGCGCACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGG‑CAGCCCGGATTTCACTTT‑CCAGCCGGTCGAGGATCTGc                                                                                                                                                                                                                          >  1:539531/1‑238 (MQ=255)

                       gTTGCCCCATTCATCGTATTCTGCGCACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGGGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGCTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCC‑CTCACCCGGT‑CAGCCCGGCTTTCCCTTT‑CCCGCCGGCCGAGCATCTGc                                                                                                                                                                                                                          >  1:282599/1‑238 (MQ=11)

                               aTTCATCGTATTCTGCGCCCCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCGCCCGGT‑CAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCAcc                                                                                                                                                                                                                 <  1:617063/239‑1 (MQ=37)

                               aTTCATCGTATTCTGCGCACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCC‑CTCACCCGGT‑CAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCAcc                                                                                                                                                                                                                 <  1:688823/239‑1 (MQ=37)

                                  cctCGTGTCCTGCGCCCCATTCGTTTGTCCCTTCCTGGCTGATAAGGGCCAGCGGCGGGCCGCGATGGTCGCAGGGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCGGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGT‑CAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACcggc                                                                                                                                                                                                              <  2:1326119/237‑1 (MQ=11)

                                  cATCGTATTCTGCGCACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGT‑CAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACcggc                                                                                                                                                                                                              <  1:1501273/239‑1 (MQ=37)

                                  cATCGTATTCTGCGCACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGT‑CAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACcggc                                                                                                                                                                                                              <  1:715977/239‑1 (MQ=37)

                                  cATCGTATTCTGCGCACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGT‑CAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACcggc                                                                                                                                                                                                              <  1:1414514/239‑1 (MQ=37)

                                   aTCGTATTCTGGGCACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGT‑CAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACcggc                                                                                                                                                                                                              <  1:1299218/238‑1 (MQ=37)

                                    tctatttCTGCGCCCCATTTTGTTGTCCTTTCCTTGCGGATGAGGGCCTGCGGCAGGCCGCAATGGTGGTAGTGGTTCAGGTGGATTTTTCGCGCCGGCGTTTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCGGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGT‑CAGCCAGGATTTCCCTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACcggcgg                                                                                                                                                                                                            <  2:1717084/234‑1 (MQ=11)

                                     cGTTTTCTGCTCACCTATCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCAGCGATGGTCGCAGTGGTGCAGGTGTATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGT‑CAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCggg                                                                                                                                                                                                           <  2:718843/239‑1 (MQ=255)

                                        tttgtGCGCACCATTCTGTTGGCCCTTCCTGGTTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTCCAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGT‑CAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGaa                                                                                                                                                                                                         <  2:1775017/234‑1 (MQ=255)

                                        tttCTGCCCCCCATTTGGTGGTCCCTTCCTTGCTGATAAGGGCCGGCGGCAGGCCGCGATGGTCGCAGTGGTCCAGGTGGATTTTTCGCGCCGGCGTGTACCCCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGT‑CAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGaa                                                                                                                                                                                                         <  2:1944508/237‑1 (MQ=11)

                                          tCTGCGCACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGT‑CAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAacac                                                                                                                                                                                                      <  1:489902/239‑1 (MQ=37)

                                             gcgcACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGGGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTTCTCA‑CCCACCCGGC‑CAGCCAGGGTTTCACTTT‑CCAGCCGGTTCAGCATCTCCACCAGCCCcgggcggagacacc                                                                                                                                                                                                    >  2:220152/1‑227 (MQ=11)

                                              cgcACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGT‑CAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCCCCGGCGGGGaccacac                                                                                                                                                                                                   >  1:1963785/1‑233 (MQ=255)

                                              cgcACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGT‑CAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAAcaccac                                                                                                                                                                                                   >  1:390195/1‑238 (MQ=37)

                                                   cATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGT‑CAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAAcaccaccccgc                                                                                                                                                                                              >  1:1357387/1‑233 (MQ=37)

                                                     ttCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGT‑CAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGGAaaccaccaccccc                                                                                                                                                                                            >  1:1203616/1‑231 (MQ=255)

                                                           tgtCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCGCCCGGT‑CAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑Ctt                                                                                                                                                                                    <  1:1571288/239‑1 (MQ=25)

                                                           tgtCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCCCCCGGT‑CAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑Ctt                                                                                                                                                                                    <  1:2026055/239‑1 (MQ=25)

                                                             tCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGT‑CAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTCCACCAGCACCGGCGGGCACACCCCACTGCCCCCGT‑Cttt                                                                                                                                                                                   >  1:1093305/1‑237 (MQ=19)

                                                               ccttccttGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGT‑CAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACAACCCACATCCACCgcccttcc                                                                                                                                                                                  >  1:641763/1‑231 (MQ=255)

                                                                 ttccttGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACTCCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGT‑CAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CTTcgccgc                                                                                                                                                                              <  1:1987231/239‑1 (MQ=18)

                                                                  tccttGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGT‑CAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGGACACCACACTGCCCCCCT‑CTTcccccc                                                                                                                                                                              >  1:1376203/1‑236 (MQ=16)

                                                                  tccttGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGT‑CAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGCCGGGGACAACACACCCCCACCGC‑CCTccccgc                                                                                                                                                                              >  1:913820/1‑238 (MQ=11)

                                                                   ccttGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGT‑CAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CTTCGccgccg                                                                                                                                                                            <  1:1554513/239‑1 (MQ=33)

                                                                   ccttGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTGCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGT‑CAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CTTCGccgccg                                                                                                                                                                            <  1:2150997/239‑1 (MQ=18)

                                                                      tGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGT‑CAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGGACACCCCACCGCCCCCCTCCTTCCCCCCCgg                                                                                                                                                                           >  1:183982/1‑238 (MQ=11)

                                                                                  cgaGCGGCAGCCAGTGCGGGTCGCAGTGGTCCAGGTGGATTTTGCGCGCCGGCGTTTACCCCGGGTCCATCTGGTTTTGCATCTGCTCCCCGGTCAGGCCGCCCGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGTT‑CAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACCCCACACTGCCACCGT‑CTTCGCCGCCGGACTGCAGAAgcgc                                                                                                                                                              <  2:807933/236‑1 (MQ=11)

                                                                                  ccAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGT‑CAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CTTCGCCGCCGGACTGCTGAAgcgc                                                                                                                                                              <  1:1823006/238‑1 (MQ=31)

                                                                                   cAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGT‑CAGCCAGGACTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CTTCGCCGCCGGACTGCTGAAGCGCa                                                                                                                                                             <  1:1818090/238‑1 (MQ=14)

                                                                                       ggCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTGGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCCCCCGGT‑CAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CTTCGCCGCCGGACTGCTGAAGCGCATCCg                                                                                                                                                         <  1:1014886/238‑1 (MQ=14)

                                                                                             ccGCGATGGTCGCAGTGGTACAGGTGGA‑TTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGT‑CAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCGCCGT‑CTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAgg                                                                                                                                                    <  1:2219384/236‑1 (MQ=14)

                                                                                                      tcgcgttGGTACAGGGGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGT‑CAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCt                                                                                                                                          <  1:1469170/232‑1 (MQ=14)

                                                                                                             ggTACAGGTGGAGTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGT‑CAGCCAGGATCTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGACCACCACACTGCCACCGT‑CTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGtt                                                                                                                                     <  1:1911120/236‑1 (MQ=14)

                                                                                                               tACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGT‑CAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTCCCACCGT‑CCTCGCCCCCGGGCTCCTGAAGCGCCTCCGCCAGGCCCCGGCGCCGCGTTTTcc                                                                                                                                 >  1:1641542/1‑237 (MQ=11)

                                                                                                               tACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCGCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGT‑CAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCg                                                                                                                                 <  1:1729849/238‑1 (MQ=14)

                                                                                                                aCAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGT‑CAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGc                                                                                                                                <  1:669823/238‑1 (MQ=31)

                                                                                                                    gTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGT‑CAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGGACACCACACCGCCCCCGG‑CCTCTCCCCCGGCGCGCCGGAGCGCCGCCCCCACGCGGCGGCGCCGCGGTTTCCCCAGc                                                                                                                            >  2:418437/1‑238 (MQ=11)

                                                                                                                            tccgggcCGGAGTGTCCACCGGTTCCATCCGGTTTGGCATCTGCTCCACGGTCGGGCCGCAGGAGGCCACCCGGCGGCGGCTTTCCTCC‑CTCCCCCGGT‑CAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCA‑CCGGTg                                                                                                                  <  2:1179364/234‑1 (MQ=11)

                                                                                                                              cgcgCCGGCGTGTCCACCGGGTCCATCGGGTTTTGCCTCTCCTCCGCGGTCAGGCCGCCGTATGCCAGCCAGCGGCGGCTTTCCTCA‑CCCACCCTGT‑CAGCCAGGATTTCCCTTC‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCA‑CCGGTGGc                                                                                                                <  2:1275174/239‑1 (MQ=11)

                                                                                                                              cgcgCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGT‑CAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGGACACCACACTGCCACCGT‑CTTCTCCCCCGGGCTGCTGAAGCGCATCCCCCAGGCCGCGGCGCCGCGCTTTCGCCAGCCCA‑CCGggggg                                                                                                                >  1:757414/1‑238 (MQ=11)

                                                                                                                              cgcgCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGT‑CAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCA‑CCGGTGGc                                                                                                                <  1:2163491/239‑1 (MQ=31)

                                                                                                                                   cgtggtgtACCCCGGGTCCATCTGGTTTTGCATCTGCTCCGCGGTCAGGCCGCACGATGCCAGCCGGCGGCGGCTTTCCTCA‑CTCACCCGGT‑CAGCCAGGATTTCGCTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCA‑CCGGTGGCGGttt                                                                                                           <  2:705580/235‑1 (MQ=12)

                                                                                                                                   cggcggGTCCTCCGGTCCCACCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGAGGCCAGCCAGCGGCGGCTTTCCCCC‑CTCACCCGGT‑CAGCCAGTGTTTCCCTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCAGCGCTGCGTTTTCGCCAGCTCA‑CCGGTGGCGGttt                                                                                                           <  2:540133/239‑1 (MQ=11)

                                                                                                                                   cGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGT‑CAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCA‑CCGGTGGCGGttt                                                                                                           <  1:2063375/239‑1 (MQ=32)

                                                                                                                                       gtgtACCCCGGGTCCATCTGGTTTTGCTCCTGCTCCCCGGTCGGGCCGCACGAGGCCAGCCAGCGGCGGCTTTCCTCA‑CTCCCCCGGT‑CAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCA‑CCGGTGGCGGTTTCaa                                                                                                        <  2:183982/238‑1 (MQ=255)

                                                                                                                                       gtgtACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGT‑CAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCA‑CCGGTGGCGGTTTCaa                                                                                                        <  2:1577933/238‑1 (MQ=33)

                                                                                                                                       gtgtACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGAGGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGT‑CAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTAA‑CCGGTGGCGGTTTCaa                                                                                                        <  1:1708350/238‑1 (MQ=18)

                                                                                                                                             cccggggcactCTGGTTTTCTCTCTGCTCCACGGTCAGGTCGCAGGGGCACGGCCTGCGGTGGCTTTCCTCA‑CTCACCCGGT‑CGGCCGGGATTGCTCTTT‑CCCGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCCCACTGCCCCCGT‑CTTCGCCGCCGGACTGCTGAGGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCC‑CCGGTGGCGGTTTCAACTCtgat                                                                                                 <  2:1994001/229‑1 (MQ=11)

                                                                                                                                                gggTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGT‑CCGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCA‑CCGGTGGCGGTTTCAACTCTGATGAg                                                                                              <  1:2051651/239‑1 (MQ=25)

                                                                                                                                                   tCCTTCGGGTTTGGCACCTGCGCCCGGTTCAGGCCGCACGTGGCCAGCCAGCGGCGGCTTTCCGCC‑CTCACCCGGT‑CAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACCCCACACTGCCACCGT‑CTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCA‑CCGGTGGCGGTTTCAACTCTGATGAGTg                                                                                            <  2:1368073/238‑1 (MQ=11)

                                                                                                                                                        gtGGTTTTGCTTCTGCTCGACGGCCAGGCCTCACGATCCCAGCCCGCGGCGGTTTTCCTCC‑CTCCCCCGGT‑CAGCCAGGATTTCACTTC‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCA‑CCGGTGGCGGTTTCAACTCTGATGAGTGG‑CGTGa                                                                                     <  2:1525438/238‑1 (MQ=11)

                                                                                                                                                        cTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCCCCCGGT‑CAGCCAGGATTTCCCTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCA‑CCGGTGGCGGTTTCAACTCTGATGAGTGG‑CGTGa                                                                                     <  1:859635/239‑1 (MQ=25)

                                                                                                                                                                  aTCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGT‑CAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCA‑CCGGTGGCGGTTTCAACTCTGATGAGTGG‑CGTGAAGCTCCCCgg                                                                           <  1:1616903/239‑1 (MQ=35)

                                                                                                                                                                  aTCTGCCCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCC‑CTCGCCCGGT‑CCGCCAGGATTTCCCTTT‑CCAGCCGGTCGAGCATCTGCGCCAGCACCGGCGGGAACACCACACTGCCACCGT‑CTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCA‑CCGGTGGCGGTTTCAACTCTGATGAGTGG‑CGTGAAGCTCCCCgg                                                                           <  1:969543/239‑1 (MQ=255)

                                                                                                                                                                          cccGGCCAGGCCGCAGCATGCCTGCCGGCGCCGGCTTGCCTCA‑CCGCCCCGGT‑CAGCCAGGCTTTCACTTT‑CCAGCCGGTCGGGCAGCTGCACCGGCACCGGCGGGAACACCACACTGCCACCGT‑CTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCA‑CCGGTGGCGGTTTCAACTCTGATGAGTGG‑CGTGAAGCTCCCCGGCTGATaaa                                                                   <  2:1632386/237‑1 (MQ=11)

                                                                                                                                                                                    ccGCACGATGCCAGCCGGCGGCGGCTTTCCTCA‑CTCCCCCGGT‑CAGCCGGGATTTCCCTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCA‑CCGGTGGCGGTTTCAACTCTGATGAGTGG‑CGTGAAGCTCCCCGGCTGATAAATCGTCTGGAt                                                         <  2:258178/239‑1 (MQ=16)

                                                                                                                                                                                    ccGCACGAGGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGT‑CAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCA‑CCGGTGGCGGTTTCAACTCTGATGAGTGG‑CGTGAAGCTCCCCGGCTGATAAATCGTCTg                                                            <  1:505017/236‑1 (MQ=18)

                                                                                                                                                                                      gCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGT‑CAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCCGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCA‑CCGGTGGCGGTTTCAACTCTGATGAGTGG‑CGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGc                                                       <  1:788141/239‑1 (MQ=18)

                                                                                                                                                                                         cGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCCCCCGGT‑CAGCCAGGATTTCCCTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCA‑CCGGTGGCGGTTTCAACTCTGATGAGTGG‑CGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCggg                                                    <  1:1258667/239‑1 (MQ=18)

                                                                                                                                                                                         cGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGT‑CAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCA‑CCGGTGGCGGTTTCAACTCTGATGAGTGG‑CGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCggg                                                    <  1:991497/239‑1 (MQ=33)

                                                                                                                                                                                         cGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGT‑CAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCA‑CCGGTGGCGGTTTCAACTCTGATGAGTGG‑CGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCggg                                                    <  1:220152/239‑1 (MQ=33)

                                                                                                                                                                                           aTGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGT‑CAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCGCCGT‑CTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCA‑CCGGTGGCGGTTTCAACTCTGATGAGTGG‑CGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGtt                                                  <  1:1143526/239‑1 (MQ=18)

                                                                                                                                                                                             gcgcgcgcgcgccgTCTTTCCTGC‑CACAGCCGGT‑CGGCCCGGATTCCCCATG‑CCAGCCGGTCGAGCATCCGCCCCAGCACCGGCCGGAACACCACGCTGCCGCCTT‑CTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGGGCTGCGTTTTCGCCAGCTCA‑CCGGTGGCGGTTTCAACTCTGATGAGTGG‑CGTGAAGCTCCCCGGCTGATAAAGCGTCTGGATGCGGGTTc                                                 <  2:1480626/230‑1 (MQ=11)

                                                                                                                                                                                              cctggcaggGGCGGCTTTCCTCC‑CACACCCGTC‑CAGCCAGGACTTCTCTTT‑CCATCCGGTCGGGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCTCCGT‑CTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCG‑CCGGTGGCGGTTTCAACTCTGATGAGTGG‑CGTGAAGCGCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTg                                               <  2:2038832/234‑1 (MQ=11)

                                                                                                                                                                                               cagccagcGGCGGCTTTCCTCA‑CTCACCCGGT‑CAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCA‑CCGGTGGCGGTTTCAACTCTGATGAGTGG‑CGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTg                                               <  1:418437/238‑1 (MQ=33)

                                                                                                                                                                                                   cggcgctggcTTTCCCCA‑CTC‑CCCGGTGCGGCCAGGATTTCACTTC‑CCAGCCGGTCGAGCATCTGCACCAGCCCCGGCGGGAACGCCACACTGCCACCGT‑CTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCA‑CCGGTGGCGGTTTCAACTCTGATGAGTGG‑CGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGtt                                          <  2:1429485/232‑1 (MQ=11)

                                                                                                                                                                                                   cagcGGCGGCTTTCCTCA‑CTCACCCGGT‑CAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCA‑CCGGTGGCGGTTTCAACTCTGATGAGTGG‑CGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTc                                             <  1:9411/236‑1 (MQ=33)

                                                                                                                                                                                                      cggcggCTTTCCTTA‑CTCACCCGGT‑CAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCA‑CCGGTGGCGGTTTCAACTCTGATGAGTGG‑CGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCt                                        <  1:1748743/238‑1 (MQ=18)

                                                                                                                                                                                                      cggcggCTTTCCTCA‑CTCACCCGGT‑CAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAAAACCCCCCCGCCCCCCG‑CTTTCCCCCCCGGCCGCCGGAACGCCTCCCCCAGGGCGCCGCGCTGCGGTTTCGCCAGCCCACCCGGGGCCGGTCTCACCCCCG‑TGAGTGGCCGTGAAGCTCCCCCGCCGGTAAATCGCCGGGGTGCGGGtgttgtgttcgt                                        >  2:392284/1‑228 (MQ=2)

                                                                                                                                                                                                         cggtTTCCCTCC‑CCCCCCCGGC‑CAGCCAGGATTTCCCTTT‑CCAGCCGGTCGGGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCA‑CCGGTGGCGGTTTCAACTCTGATGAGTGG‑CGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTa                                     <  2:609718/238‑1 (MQ=11)

                                                                                                                                                                                                         cggcTTTCCTCA‑CTCACCCGGT‑CAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCA‑CCGGTGGCGGTTTCAACTCTGATGAGTGG‑CGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTa                                     <  1:1866475/238‑1 (MQ=33)

                                                                                                                                                                                                          ggcTTTCCTCC‑CCCCCCCGGT‑CAGCCAGGATTTCCCTTT‑CCAGCCGGTCGGGCATCTGCACCAGCCCCGGCGGGAACACCACACTGCCACCGT‑CTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCA‑CCGGTGGCGGTTTCAACTCTGATGAGTGG‑CGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTAtt                                   <  2:1062731/239‑1 (MQ=255)

                                                                                                                                                                                                               tCCTCA‑CTCACCCGGT‑CAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCCCCGGCGGGAACACCACACTGCCACCGT‑CTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCA‑CCGGTGGCGGTTTCAACTCTGATGAGTGG‑CGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGt                              <  1:1226869/239‑1 (MQ=25)

                                                                                                                                                                                                                   ca‑ctcacCCGGT‑CAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCA‑CCGGTGGCGGTTTCAACTCTGATGAGTGG‑CGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGc                          <  1:7997/239‑1 (MQ=35)

                                                                                                                                                                                                                         ccccGGT‑C‑GCCAGGATTTGCCTGT‑CCAGCCGGTCGCGCATCTGCACCAGCACCGGCGGGAACACCACACCGCCACCGT‑CTTCCCCGCCGTACTGCTCAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCA‑CCGGTGGCGGTTTCAACTCTGATGAGTGG‑CGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCg                    <  2:1118838/238‑1 (MQ=11)

                                                                                                                                                                                                                          cccGGG‑CCGCCGGGCTTTCCCTTT‑CCAGCCGGTCGAGCCCCGGCCCCAGCACCGGCGGGAGCCCCACACTGCCACCGT‑CTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCA‑CCGGTGGCGGTTTCAACTCTGATGAGTGG‑CGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCg                    <  2:35803/239‑1 (MQ=11)

                                                                                                                                                                                                                          cccGCT‑CAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCCCACTGCCACCGT‑CTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCA‑CCGGTGGCGGTTTCAACTCTGATGAGTGG‑CGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCg                    <  2:888422/239‑1 (MQ=37)

                                                                                                                                                                                                                           ccGGT‑CAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCA‑CCGGTGGCGGTTTCAAct                                                                                                      <  2:468962/157‑1 (MQ=34)

                                                                                                                                                                                                                            cGGT‑CAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCA‑CCGGTGGCGGTTTCAACTCTGATGAGTGG‑CGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGcc                  <  1:258294/239‑1 (MQ=37)

                                                                                                                                                                                                                                 cAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCCCCAGCACCGGCGGGAACACCACACTGCCACCGT‑CTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGGGCGGGGTTTTCGCCAGCTCA‑CCGGGGGCGGTTTCAACTCTGATGAGTGG‑CGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGGCGCCGTc               <  2:307296/238‑1 (MQ=255)

                                                                                                                                                                                                                                    ccAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCA‑CCGGTGGCGGTTTCAACTCTGATGAGTGG‑CGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAg           <  1:1175254/239‑1 (MQ=37)

                                                                                                                                                                                                                                      ccgATTTTCCTTTCCCAGCGGGTCGGGCATCTGCCCCGGCACCGGCGGGACCCCCCCCCTGCCACCGT‑CTTCGCCGCCGGACTGCTGAGGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCA‑CCGGTGGCGGTTTCAACTCTGATGAGTGG‑CGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTccc             <  2:2100160/234‑1 (MQ=11)

                                                                                                                                                                                                                                         ttttCACTTC‑CCAGCCGGTCGGGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCA‑CCGGTGGCGGTTTCAACTCTGATGAGTGG‑CGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCg        <  1:1070616/236‑1 (MQ=37)

                                                                                                                                                                                                                                         aTTTCCCTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCA‑CCGGTGGCGGTTTCAACTCTGATGAGTGG‑CGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCg        <  1:458632/237‑1 (MQ=37)

                                                                                                                                                                                                                                          tttCACTTC‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCCCACTGCCACCGT‑CTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCA‑CCGGTGGCGGTTTCAACTCTGATGAGTGG‑CGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTAc     <  1:1116310/239‑1 (MQ=37)

                                                                                                                                                                                                                                           ttCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCA‑CCGGTGGCGGTTTCAACTCTGATGAGTGG‑CGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTAcc    <  1:1332595/239‑1 (MQ=37)

TGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCGTATTCTGCGCACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGT‑CAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCA‑CCGGTGGCGGTTTCAACTCTGATGAGTGG‑CGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAG  >  NC_007779/3874685‑3875158

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