BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_007779	3,875,144	T→C	8.9%	D1030G (GAC→GGC)	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,875,144	0	T	C	8.9%	574.0 / 5.3	243	D1030G (GAC→GGC)	Y75_RS19225	RHS element protein RhsA
Reads supporting (aligned to +/- strand): ref base T (97/117); new base C (2/25); total (100/143)
Fisher's exact test for biased strand distribution p-value = 1.02e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 244 total aligned reads displayed.

TCAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTCCCGTTCACGCCGCCACACCCGTTTTGCCACCCTGCGGCCCAGCGGGTCGTAAAGATAGCGACTTTCGACCAGCGGCTCTTCATATTGTGTCCGCGTGTAGTGCACCAGCCGGTGCTGACTGTCGTAATGGTACCGGTGAGTCCGCTCATCATCCGTGCGG  >  NC_007779/3874907‑3875362

tCAGCCAGGATTTCCCTTT‑CCAGCCGGTCGAGCATCTGCACCAGCTCCGGCGGGAACACCCCCCTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTc                                                                                                                                                                                                                           <  1:366734/239‑1 (MQ=255)

tCAGCCAGGATTTCCCTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGt                                                                                                                                                                                                                            <  1:1244158/238‑1 (MQ=37)

 cAGCCAGGATTTCACTTC‑CCCGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCCCACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTcc                                                                                                                                                                                                                          <  1:2020994/239‑1 (MQ=255)

    ccgggttctcccTTC‑CCGGCCGGTCGGGCATCTGCACCAGCACCGGCGGGAACCCCCCACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAg                                                                                                                                                                                                                       <  2:710959/228‑1 (MQ=11)

    ccAGGATTTCCCTTC‑CCCGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAg                                                                                                                                                                                                                       <  2:895782/239‑1 (MQ=255)

     cggtgtttCCCTTC‑CCCCCCGGTCGCGCATCTGCACCAGCACCGGCGGCAACCCCGCACGGCCACCTTCTTC‑GCCGCCGGACGGCTGAAGCGCAGCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGc                                                                                                                                                                                                                      <  2:27411/234‑1 (MQ=11)

     cctgtgttCCCTTTCCCAGCCGGTCGGGCATCCGCCCCAGCCCCGGCGGGAACCCCCCACTGCCACCGTCTCC‑GCCGCCGGACTGCTGAGGCGCATCCGCCCGGCTGGGGCGCTGCGTTTGCGCCAGCTAACCGGTGGCGGTTTCAACTCTGATGGGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAg                                                                                                                                                                                                                       <  2:2033042/233‑1 (MQ=11)

     ccggATTTCACTGT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGc                                                                                                                                                                                                                      <  1:1063459/237‑1 (MQ=37)

     cAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAg                                                                                                                                                                                                                       <  1:217034/238‑1 (MQ=37)

           ttCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTAcc                                                                                                                                                                                                                <  1:1332595/239‑1 (MQ=37)

             cACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAg                                                                                                                                                                                                              <  1:1959824/239‑1 (MQ=37)

                ttt‑CCAGCCGGTCGAGCATCTGCCCCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTc                                                                                                                                                                                                           <  1:1287482/239‑1 (MQ=37)

                    cccgCCGGTCGCGCATCTGCCCCGGCACCGGCGGAAACCCCACACTGCCACCGTCGTC‑GCCGCCGGACTGCTGAGGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGGGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACt                                                                                                                                                                                                        <  2:1080612/236‑1 (MQ=11)

                    ccAGCCGGCCCAGCTTCTGCACCACGACCGGGAGTTTCACCCCGTGGCCCCCGGTTTC‑GCCGCCTGCCTCGCGAAGCGTGTCCCCCCGGTTGTGGCATTGCGTTGTCGCCGGCTTCCCGGGGGCGGTTTCAACTCTGATGGGTGGCGTGAAGCTCCCCGTCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTTTGTATTGTGTTCAGCCGGTCCCCGTCCCAGCCGTTCCAGGTCACt                                                                                                                                                                                                        <  2:1988577/239‑1 (MQ=2)

                      acccGGTCGAGCATCTGCACCAGCACCGGCTGGAACGCCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGGGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCACCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTg                                                                                                                                                                                                      <  1:771494/237‑1 (MQ=255)

                                     tGCCCCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAg                                                                                                                                                                                       <  1:995756/239‑1 (MQ=37)

                                      gccccAGCCCCGGCGGAAACCCCACACTGCCCCCTTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGGGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGc                                                                                                                                                                                      <  2:913820/236‑1 (MQ=255)

                                      gCACCAGCACCGGCGGGAACCCCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGc                                                                                                                                                                                      <  1:1511195/239‑1 (MQ=37)

                                      cccccggccccGGCGGGAACCCCCCCCAGCCACCG‑CTTCTGCCGCCGGACTGCTGAGGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCCGCTCACCGGTGGCGGTTTCAACTCTGATGAGGGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGCCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGc                                                                                                                                                                                      <  2:1186400/230‑1 (MQ=11)

                                                cggcggGAACACCTCCCTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGc                                                                                                                                                                            <  1:206643/239‑1 (MQ=37)

                                                cggcggGAACACCACTCTGCCCCCGTCGTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCGGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGc                                                                                                                                                                            <  1:673634/239‑1 (MQ=255)

                                                      gAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTc                                                                                                                                                                       <  1:519258/238‑1 (MQ=37)

                                                         cccccctccgCCCCCGTCTTG‑GCCGCCGGCCGGCGGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGGGGGGCGGGAGGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTCCCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGt                                                                                                                                                                   <  2:390195/230‑1 (MQ=11)

                                                          accacACTGCCACCGTCTTC‑GCCGCCGGACTGCTTAAGCGCATCCGCCAGGCTGCGGCGCGGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTc                                                                                                                                                                  <  1:1300748/239‑1 (MQ=37)

                                                                cTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCCGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGCCATCCCGCCCGGCAGGTCCCGtt                                                                                                                                                             >  2:1454763/1‑238 (MQ=255)

                                                                 tGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGCCCCAGCCGTACCAGGTCACt                                                                                                                                                                                                        <  2:292353/194‑1 (MQ=38)

                                                                 tGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGAGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTGTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGCGACAGCGACATCCAGCCCGTCAGGTCCCGTTc                                                                                                                                                            >  1:282375/1‑238 (MQ=255)

                                                                      ccGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCCGCCCGTCCGGTCCCGTTTAcgcc                                                                                                                                                       >  2:1288058/1‑238 (MQ=255)

                                                                      ccGTCTCC‑GCCCCCGGCCTGTTGACGCGCTCCCGCAAGGCTGCGGGGCTGCTTTTTCGCCGGCTCACGGGTGGCGGTTTCACCCCTGATGAGTGGCGTGAGGCCCCCGGGCTGATAATTCGTCTGGTTGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCGGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTCCCGTTCAcgccg                                                                                                                                                      <  2:355979/239‑1 (MQ=11)

                                                                          cTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTCCCGTTCACGCCGCca                                                                                                                                                   <  1:1209877/238‑1 (MQ=38)

                                                                           ttC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCCGGTCCCGGTCACGCCGCcac                                                                                                                                                  >  2:677816/1‑238 (MQ=255)

                                                                             c‑gccgcCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCCACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCCGGATGCGGGTTCTGGCGTTCTGGATTGTGGGCAGCCGGGCGCCGTCCCGGCCGTACCAGGTCCCTTGCGGTTTCCGTGGCCGGCGCATCCCGCCCCCCAGGTCCCGGTCCCCCCCCCccac                                                                                                                                                >  2:150275/1‑238 (MQ=11)

                                                                                   ccgGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTTACAGCGACCTCCAGCCCCCCAGGTCCCGTTCACTCCGCCCCACCCGttt                                                                                                                                          >  2:1301959/1‑239 (MQ=255)

                                                                                    cgGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTCCCGTTCACGCCGCCACACCCGtttt                                                                                                                                         >  1:2029481/1‑239 (MQ=38)

                                                                                           ctgAAGCGCATCCGCCAGGCTGCGGCGCTGCGGTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGCGTGGCGTTCAGCTCCCCGGCTGATAAATCGTCTGGACGCCGGTTCTGTCGTTCTGGATTCGGAGCTGTCGGCGGCCGTGCCACCCGTACCAGGTCAGTGGCGGTTTCCGTGACAGCCTCCTCCAGCGCGTCAAGGCCAGTTCACGCCGCCGCGCCCGGTTTcccccc                                                                                                                                   >  2:1873145/1‑238 (MQ=2)

                                                                                                gcgcATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCg                                                                                                                                                                                                    <  2:938492/168‑1 (MQ=39)

                                                                                                  gcAGCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTCCCGTTCACGCCGCCACACCCGTTTTGCCACCCTGCGGcc                                                                                                                           <  1:118082/239‑1 (MQ=25)

                                                                                                     tCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCg                                                                                                                                                                                                    >  1:1872822/1‑163 (MQ=35)

                                                                                                     tCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTc                                                                                                                                                                       >  1:662366/1‑192 (MQ=34)

                                                                                                     tCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTCCCGTTCACGCCGCCACACCCGTTTTGCCACCCTGCGGCCCa                                                                                                                         >  1:843884/1‑238 (MQ=35)

                                                                                                     tCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTCCCGTTCACGCCGCCACACCCGTTTTGCCACCCTGCGGCCCa                                                                                                                         >  1:186584/1‑238 (MQ=35)

                                                                                                     tCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTCCCGTTCACGCCGCCACACCCGTTTTGCCACCCTGCGGCCCa                                                                                                                         >  2:2153069/1‑238 (MQ=35)

                                                                                                      ccgccgGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTCCCGTTCACGCCGCCACACCCGTTTTGCCACCCTGCGGCCCAGc                                                                                                                       <  1:1544609/239‑1 (MQ=25)

                                                                                                      ccGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCCCCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTCCCGTTCACGCCGCCACACCCGTTTTGCCACCCTGCGGCCCAGc                                                                                                                       <  1:1174699/239‑1 (MQ=25)

                                                                                                      ccGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTCCCGTTCACGCCGCCACACCCGTTTTGCCACCCTGCGGCCCAGc                                                                                                                       <  1:1551885/239‑1 (MQ=35)

                                                                                                                 cggcgCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTCCCGTTCACCCCGCCACCCCCGTTTTGCCACCCCGCGGCCCAGCGGGTCGTAAAg                                                                                                            >  2:162311/1‑239 (MQ=25)

                                                                                                                      cTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCACCTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAACCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGGGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTCCCGTTCACGCCGCCACACCCGTTTTGCCACCCTGCGGCCCAGCGGGTCGTAAAGATAg                                                                                                        <  1:560396/238‑1 (MQ=25)

                                                                                                                      cTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTCCCGTTCACGCCGCCACACCCGTTTTGCCACCCTGCGGCCCAGCGGGTCGTAAAGATAg                                                                                                        >  1:1478639/1‑238 (MQ=35)

                                                                                                                                      tCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTCCCGTTCACGCCGCCACACCCGTTTTGCCACCCTGCGGCCCAGCGGGTCGTAAAGATAGCGACTTTCGACCAg                                                                                          <  1:930212/236‑1 (MQ=35)

                                                                                                                                      tCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTCCCGTTCACGCCGCCACACCCGTTTTGCCACCCTGCGGCCCAGCGGGTCGTAAAGATAGCGACTTTCGACCAGCg                                                                                        <  1:1454763/238‑1 (MQ=35)

                                                                                                                                       cgccGGGGGGGGTTTCAACTCGGTGGAGGGGCGTGAGGCGCCCCGGCTGATAAATCGTCTGGGTGCGGGTTCGGTCGTTCTGTTTTGTGGTCAGCCGGTCGCCGCCCCAGCCGTACCAGGTCACTTGCGGTTTCCGGGACAGCGACATCCAGCCCGTCAGGTCCCGTTCACGCCGCCACACCCGTTTTGCCACCCTGCGGCCCAGCGGGTCGTAAAGATAGCGACTTTCGACCAGCGGc                                                                                      <  2:852195/237‑1 (MQ=11)

                                                                                                                                       ccccGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTCCCGTTCACGCCGCCACACCCGTTTTGCCACCCTGCGGCCCAGCGGGTCGTAAAGATAGCGACTTTCGACCAGCGGc                                                                                      <  1:1758540/237‑1 (MQ=25)

                                                                                                                                       cACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTTAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTCCCGTTCACGCCGCCACACCCGTTTTGCCACCCTGCGGCCCAGCGGGTCGTAAAGATAGCGACTTTCGACCAGCGGc                                                                                      <  1:946773/239‑1 (MQ=25)

                                                                                                                                       cACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGTTGATAAATCGTTTGGATGCGGATTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTCCCGTTCACGCCGCCACACCCGTTTTGCCACCCTGCGGCCCAGCGGGTCGTAAAGATAGCGACTTTCGACCAGCGGc                                                                                      <  1:127103/239‑1 (MQ=25)

                                                                                                                                       cACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTCCCGTTCACGCCGCCACACCCGTTTTGCCACCCTGCGGCCCAGCGGGTCGTAAAGATAGCGACTTTCGACCAGCGGc                                                                                      <  1:449241/239‑1 (MQ=35)

                                                                                                                                       cACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTCCCGTTCACGCCGCCACACCCGTTTTGCCACCCTGCGGCCCAGCGGGTCGTAAAGATAGCGACTTTCGACCAGCGGc                                                                                      <  1:257135/239‑1 (MQ=35)

                                                                                                                                       cACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATGGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTCCCGTTCACGCCGCCACACCCGTTTTGCCACCCTGCGGCCCAGCGGGTCGTAAAGATAGCGACTTTCGACCAGCGGc                                                                                      <  1:974004/239‑1 (MQ=18)

                                                                                                                                       cACCGGTGGCGGTTTCAACTCGGATGGGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGCCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTCCCGTTCACGCCGCCACACCCGTTTTGCCACCCTGCGGCCCAGCGGGTCGTAAAGATAGCGACTTTCGACCAGCGGc                                                                                      <  1:2117180/239‑1 (MQ=25)

                                                                                                                                       cACCGGTGGCGGTTTCAACTCGGATGATTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGTTTCTTTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTCCCGTTCACGCCGCCACACCCGTTTTGCCACCCTGCGGCCCAGCGGGTCGTAAAGATAGCGACTTTCGACCAGCGGc                                                                                      <  1:1288167/239‑1 (MQ=255)

                                                                                                                                            gTGGCGGTTCCAACTCGGATGAGGGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCGGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTCCCGTTCACGCCGCCACACCCGTTTTGCCACCCTGCGGCCCAGCGGGTCGTAAAGATAGCGACTTTCGACCAGCGGCTCtt                                                                                  <  1:1980916/238‑1 (MQ=255)

                                                                                                                                                gggTTTCAACTCTGTTGAGTGGCTTGAATCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTTTTGTGGTCAGCCGGCCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTCCCGTTCACGCCGCCACACCCGTTTTGCCACCCTGCGGCCCAGCGGGTCGTAAAGATAGCGACTTTCGACCAGCGGCTCTTCata                                                                              <  1:490589/237‑1 (MQ=255)

                                                                                                                                                ctttgtCAAATCTGGTGCGTCGCGTGAAGCCCCCCGGCTGATAAATCTTCGGGTCGCGGGTTCTGCTTTTCTGTATTGGGGTGAGCCGGCCGCCTCCCCAGCCGTCCCAGGTCACTGGCGGTTTCCGTGACTGCGACATCCAGCCCGTCAGGTCCCTTTCACGCCGCCACACCCGTTTTGCCACCCTGCGGCCCAGCGGGTCGTAAAGATAGCGACTTTCGACCAGCGGCTCTTCata                                                                              <  2:1581789/233‑1 (MQ=11)

                                                                                                                                                cGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTCCCGTTCACGCCGCCACACCCGTTTTGCCACCCTGCGGCCCAGCGGGTCGTAAAGATAGCGACTTTCGACCAGCGGCTCTTCata                                                                              <  1:173078/238‑1 (MQ=37)

                                                                                                                                                    ttCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCCGGTCACTTGCGGGTTCCGTGACAGCGACATCCACCCCCTCAAGTCCCCTTCCCGCCGCCACCCCCCCTTTGCCACCCTGCGGCCCAGCGGGTCGTAAAGATAGCGACATTTGACCAGCGGCCCTTCAAATtgtt                                                                         >  2:766798/1‑238 (MQ=11)

                                                                                                                                                      cacctctGAGGGGGGGTGTGAAGCCCCCGTGTTGATAATTCTTCTGGATGGGGTTTCTGTCTTTCTGTTTGGTGCTAAGCCGGCGCCCGCCCCAGCCGTACCAGGTCCCTTGCGGTTGCCGTGACAGCGACATCCAGCCCGTCAGGTCCCGTTCACGCCGCCACACCCGTTTTGCCACCCTGCGGCCCAGCGGGTCGTAAAGATAGCGACTTTCGACCAGCGGCTCTTCATATTGTGTc                                                                       <  2:645379/236‑1 (MQ=11)

                                                                                                                                                          tctGATGTGTGGCGTGAAGCTCCCCGGCTGGTAAATCGTCTGGATGCGGGTTCTGTCGTTCGGTATTGGGGTCAGCCGGTCGCCTTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTCCCGTTCACGCCGCCACACCCTTTTTGCCACCCTGCGGCCCAGCGGGTCGTAAAGATAGCGACTTTCGACCAGCGGCTCTTCATATTGTGTCcgcg                                                                   <  1:710763/239‑1 (MQ=255)

                                                                                                                                                             gctgtggggCTTGAAGCTCCCGGGTTGATAACCCGCCTGGAGCCGGGTTCTGTCGTTGTGTATTGGGGTCGGCCGGTCGCCCTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTCCCGTTCACGCCGCCACACCCGTTTTGCCACCCTGCGGCCCAGCGGGTCGTAAAGATAGCGACTTTCGACCAGCGGCTCTTCATATTGTGTCCGCgtg                                                                 <  2:1129499/231‑1 (MQ=11)

                                                                                                                                                             gatgaGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTCCCGTTCACGCCGCCACACCCGTTTTGCTACCCTGCGGCCCAGCGGGTCGTAAAGATAGCGACTTTCGACCAGCGGCTCTTCATATTGTGTCCGCgtgt                                                                >  1:597845/1‑239 (MQ=18)

                                                                                                                                                              atgaGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTCCCGTTCACGCCGCCACACCCGTTTTGCCACCCTGCGGCCCAGCGGGTCGTAAAGATAGCg                                                                                                      >  1:1564551/1‑200 (MQ=34)

                                                                                                                                                              atgaGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTCCCGTTCACGCCGCCACACCCGTTTTGCCACCCTGCGGCCCAGCGGGTCGTAAAGATAGCGACTTTCGACCAGCg                                                                                        <  1:818664/214‑1 (MQ=34)

                                                                                                                                                                   tGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTCCCGTTCACGCCGCCACACCCGTTTTGCCACCCTGCGGCCCAGCGGGTCGTAAAGATAGCGACTTTCGACCAGCGGCTCTTCATATTgtgt                                                                        >  1:1858787/1‑225 (MQ=34)

                                                                                                                                                                   tGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTCCCGTTCACGCCGCCACACCCGTTTTGCCACCCTGCGGCCCAGCGGGTCGTAAAGATAGCGACTTTCGACCAGCGGCTCTTCATATTGTGTCCGCGTGTAGTGc                                                           >  1:1930110/1‑238 (MQ=33)

                                                                                                                                                                      cGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTCCCGTTCACGCCGCCACACCCGTTTTGCCACCCTGCGGCCCAGCGGGTCGTAAAGATAGCGACTTTCGACCAGCGGCTCTTCATATTGTGTCCGCGTGTAGTGCACCa                                                       >  1:1631443/1‑239 (MQ=33)

                                                                                                                                                                      cGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTCCCGTTCACGCCGCCACACCCGTTTTGCCACCCTGCGGCCCAGCGGGTCGTAAAGATAGCGACTTTCGACCAGCGGCTCTTCATATTGTGTCCGCGTGTAGTGCACCa                                                       >  1:1597965/1‑239 (MQ=33)

                                                                                                                                                                        tGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTCCCGTTCACGCCGCCACACCCGTTTTGCCACCCTGCGGCCCAGCGGGTCGTAAAGATAGCGCCTTTCGACCAGCGGCTCCTCACATTTGGTCCGCGTGTAGTGTACCAg                                                      >  2:1390109/1‑238 (MQ=255)

                                                                                                                                                                        tGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTCCCGTTCACGCCGCCACACCCGTTTTGCCACCCTGCGGCCCAGCGGGTCGTAAAGATAGCGACTTTCGACCAGCgg                                                                                       >  2:580004/1‑205 (MQ=34)

                                                                                                                                                                        tGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTCCCGTTCACGCCGCCACACCCGTTTTGCCACCCTGCGGCCCAGCGGGTCGTAAAGATAGCGACTTTCGACCAGCgg                                                                                       <  1:580004/205‑1 (MQ=34)

                                                                                                                                                                         gAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATGGTGGTCAGCCGGTCGCCGCCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTCCCGTTCACGCCGCCACACCCGTTTTGCCACCCTGCGGCCCAGCGGGTCGTAAAGATAGCGACTTTCGACCAGCGGCTCTTCATATTGTGTCCGCGTGTAGTGCACCAGc                                                     <  1:1102373/238‑1 (MQ=17)

                                                                                                                                                                               ccccGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTTTTTTGTGGTCAGCCGGTCGCCGCCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTCCCGTTCACGCCGCCACACCCGTTTTGCCACCCTGCGGCCCAGCGGGTCGTAAAGATAGCGACTTTCGACCAGCGGCTCTTCATATTGTGTCCGCGTGTAGTGCACCAGCCGGTGCt                                              <  1:1130952/239‑1 (MQ=18)

                                                                                                                                                                               ccccGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTCCCGTTCACGCCGCCACACCCGTTTTGCCACCCTGCGGCCCAGCGGGTCGTAAAGATAGCGACTTTCGACCAGCGGCTCTTCCTATTGTGTCCGCGTGTATTGCACCCGCCGGt                                                 >  2:2106033/1‑236 (MQ=18)

                                                                                                                                                                               ccccGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTCCCGTTCACGCCGCCACACCCGTTTTGCCACCCTGCGGCCCAGCGGGTCGTAAAGATAGCGACTTTCGACCAGCGGCTCTTCATCTTGTGTCCGCGGGTAGTGCACCAGCCggc                                                 >  2:260185/1‑235 (MQ=18)

                                                                                                                                                                               ccccGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCACCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTCCCGTTCACGCCGCCACACCCGTTTTGCCACCCTGCGGCCCAGCGGGTCGTAAAGATAGCGACTTTCGACCAGCGGCTCTTCATATTGTGTCCGCGTGTAGTGCACCAGCCGGTGCt                                              <  1:1373002/239‑1 (MQ=18)

                                                                                                                                                                                cccGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTCCCGGTTCACTTGCGGTTTCCGTGACAGCGCCCTCCAGCCCGTCAGGTCCCGTTCACGCCGCCACCCCCGTTTTGCCACCCTGCGGCCCAGCGGGTCGCAACGTGCGCGACTTTCGACCAGCGGCTCTTCATATGTTGACCCCGTGTAGTGCACCAGCCGGTGCTg                                             >  2:476287/1‑239 (MQ=11)

                                                                                                                                                                                cccGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTCCCGTTCACGCCGCCACACCCGTTTTGCCACCCTGCGGCCCAGCGGGTCGTAAAGATAGCGCCTTTCGACCAGCGGCTCTTCATATTGTGTCCGCGTGTAGTGCACCAGCCGGTGCTg                                             >  2:1837771/1‑239 (MQ=18)

                                                                                                                                                                                 ccGGCTGATAAATCGTCTGGATGCGGGTGCTGTCGTTCTGTATTGTGTTCAGCGGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTCCCGTTCACGCCGCCACACCCGTTTTGCCACCCTGCGGCCCAGCGGGTCGTAAAGATAGCGACTTTCGACCAGCGGCTCTTCATATTGTGTCCGCGTGTAGTGCACCAGCCGGTGCTGa                                            <  1:340249/239‑1 (MQ=17)

                                                                                                                                                                                 ccGGCTGATAAATCGTCGGGTTGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTCCCGTTCACGCCGCCACACCCGTTTTGCCACCCTGCGGCCCAGCGGGTCGTAAAGATAGCGACTTTCGACCAGCGGCTCTTCATATTGTGTCCGCGTGTAGTGCACCAGCCGGTGCTGa                                            <  1:422/239‑1 (MQ=18)

                                                                                                                                                                                  cGGCTGATAAAGCGTCTGGGTGCGGGTTCTGTCGTTCTTTTTGGTGGTCAGCCGGTCGCCGCCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTCCCGTTCACGCCGCCACACCCGTTTTGCCACCCTGCGGCCCAGCGGGTCGTAAAGATAGCGACTTTCGACCAGCGGCTCTTCATATTGTGTCCGCGTGTAGTGCACCAGCCGGTGCTGAc                                           <  1:1014289/239‑1 (MQ=255)

                                                                                                                                                                                    gCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTCCCGTTCACGCCGCCACACCCGTTTTGCCACCCTGCGGCCCAGCGGGTCGTAAAGATAGCGACTTTCGACCAGCGGCTCTTCATATTGTGTCCGCGTGTAGTGCACCAGCCGGTGCTg                                             >  1:713421/1‑235 (MQ=33)

                                                                                                                                                                                      tGAAAAATCGTCTGGATGCGGTTTCTGTCGTTCTGTTTTGTGGTCAGCCGGTCGCCGCCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTCCCGTTCACGCCGCCACACCCGTTTTGCCACCCTGCGGCCCAGCGGGTCGTAAAGATAGCGACTTTCGACCAGCGGCTCTTCATATTGTGTCCGCGTGTAGTGCACCAGCCGGTGCTGACTGt                                        <  1:1511652/238‑1 (MQ=12)

                                                                                                                                                                                          aaaTCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTCCCGTTCACGCCGCCACACCCGTTTTGCCACCCTGCGGCCCAGCGGGTCGTAAAGATAGCGACTTTCGACCAGCGGCTCTTCATATTGTGTCCGCGTGTAGTGCACCAGCCGGTGCTGACTGTCGTaa                                   <  1:1662066/239‑1 (MQ=33)

                                                                                                                                                                                              cGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTCCCGTTCACGCCGCCACACCCGTTTTGCCACCCTGCGGCCCAGCGGGTCGTAAAGATAGCGACTTTCGACCACCGGCTCTTCATAtt                                                                            >  1:18812/1‑194 (MQ=21)

                                                                                                                                                                                                  tGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTCCCGTTCACGCCGCCACACCCGTTTTGCCACCCTGCGGCCCAGCGGGTCGTAAAGATAGCGACTTTCGACCAGCGGCTCTTCATATTGTGTCCGCGTGTAGTGCACCAGCCGGTGCTGACTGTCGTAATGGTa                              >  1:1304305/1‑236 (MQ=33)

                                                                                                                                                                                                  tGGATGCGCGTCCTGTCATTTTGTATTGTGGTCAGCCGGTCGCCGCCCCCGCCCTACCAGGTCACTTGCGGTTTCCGTGACCGCGACATCTACCCCGCCAGCTCCCGTTCACGCCGCCACACCCGTTTTGCCACCCTGCGCCCCACCGGGTCGTAACGATAGCGACCTTCGTCCAGCGGCCCTTCATCTTGTTCCCGTGCGTAGTTCACCAGCCCGTCCTGCCTGCCGTAATGGGACCg                           >  2:1354025/1‑239 (MQ=11)

                                                                                                                                                                                                          gtttCTGTCGTTCGGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTCCCAGGTCACTTGCGGTTTCCGTGCCAGCGACATCCAGCCCGTCAGGTCCCGTTCACGCCGCCACACCCGTTTTGCCACCCTGCGGCCCAGCGGGTCGTAAAGATAGCGACTTTCGACCAGCGGCTCTTCATATTGTGTCCGCGTGTAGTGCACCAGCCGGTGCTGACTGTCGTAATGGTACCGGTGAGTc                    <  1:738253/236‑1 (MQ=17)

                                                                                                                                                                                                          gggtCTGTCATTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGCACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGACCCGTTCACGCCGCCACACCCGTTTTGCCACCCTGCGGCCCAGCGGGTCGAAAAGATAGCGACTTTCCACCAGCGGCCCTTCCTATTGTGTCCGCGAGTAGTGCACCAGCCGGTGCTGACTGTCGTAACAGTACCGGGGAGTca                   >  2:456396/4‑238 (MQ=11)

                                                                                                                                                                                                          ggTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTCCCGTTCACGCCGCCACACCCGTTTTGCCACCCTGCGGCCCAGCGGGTCGTAAAGATAGCGACTTTCGACCAGc                                                                                         >  2:1547388/1‑169 (MQ=34)

                                                                                                                                                                                                                      cTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTCCCGTTCACGCCGCCACACCCGTTTTGCCACCCTg                                                                                                                                >  2:1043865/1‑118 (MQ=32)

                                                                                                                                                                                                                      cTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTCCCGTTCACGCCGCCACACCCGTTTTGCCACCCTg                                                                                                                                <  1:1043865/118‑1 (MQ=32)

                                                                                                                                                                                                                         tATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTCCCGTTCACGCCGCCACACCCGTTTTGCCACCCTGCGGCCCAGCGGGTCGTAAAGATAGCGACTTTCGACCAGCGGCTCTTCATATTGTGTCCGCGTGTAGTGCACCAGCCGGTGCTGACTGTCGTAATGGTACCGGTGAGTCCGCTCATCATcc        >  1:1580115/1‑235 (MQ=31)

                                                                                                                                                                                                                           ttGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTCCCGTTCACGCCGCCACACCCGTTTTGCCACCCTGCGGCCCAGCGGGTCGTAAAGATAGCGACTTTCGACCAGCGGCTCTTCATATTGTGTCCGCGTGTAGTGCACCAGCCGGTGCTGACTGTCGTAATGGTACCGGTGAGTCCGCTCATCATCCGTGCgg  >  1:1463385/1‑239 (MQ=11)

TCAGCCAGGATTTCACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTCCCGTTCACGCCGCCACACCCGTTTTGCCACCCTGCGGCCCAGCGGGTCGTAAAGATAGCGACTTTCGACCAGCGGCTCTTCATATTGTGTCCGCGTGTAGTGCACCAGCCGGTGCTGACTGTCGTAATGGTACCGGTGAGTCCGCTCATCATCCGTGCGG  >  NC_007779/3874907‑3875362

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: `‑`