BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_007779	3,874,967	A→C	9.4%	V1089G (GTG→GGG)	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,967	0	A	C	9.4%	508.1 / 7.7	259	V1089G (GTG→GGG)	Y75_RS19225	RHS element protein RhsA
Reads supporting (aligned to +/- strand): ref base A (27/196); new base C (9/23); total (36/223)
Fisher's exact test for biased strand distribution p-value = 2.61e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 275 total aligned reads displayed.

TGCGCACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAG‑CCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGT‑CGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTC  >  NC_007779/3874729‑3875202

ttcgcgCCATTCTTTTGTCCCTTCTTGGCGGATAAGGGCCAGCGGCAGGCCGCGAGGGTGGCAGTGGTACAGGTGGATTTTTCGCGCAGGCGTGTACCCCGGGTCCATCTGGTTTTGCATCTGCTCAACGGTCAGGCCGCCCGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAG‑CCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAAcacca                                                                                                                                                                                                                                               <  2:1419057/237‑1 (MQ=11)

       cATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAG‑CCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAAcaccaccccgc                                                                                                                                                                                                                                         >  1:1357387/1‑233 (MQ=37)

         ttCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAG‑CCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGGAaaccaccaccccc                                                                                                                                                                                                                                       >  1:1203616/1‑231 (MQ=255)

          tCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCCCTCACCCGGTCAG‑CCAGGATTTCACTTTTCAACCGGTCGCGCATCTGCCCCAGCACCCGCGGGAGAACCacaccacccc                                                                                                                                                                                                                                      >  2:1028292/1‑232 (MQ=11)

          tCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGCTACAGGTGGATTTTTCGCGCCGGCGTGTACAACGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTTCTCACCCACCCGGGCAG‑CCCGGGATTCCATTTCCCGCCGGTCGCGCACTTGCACACGCACCGGGGGGGACACCacacccgcac                                                                                                                                                                                                                                      >  2:366734/1‑232 (MQ=11)

               tgtCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCGCCCGGTCAG‑CCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCtt                                                                                                                                                                                                                                <  1:1571288/239‑1 (MQ=25)

               tgtCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAG‑CCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCtt                                                                                                                                                                                                                                <  1:76226/239‑1 (MQ=35)

                 tCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCCCCCCGTCAG‑CCCGGATTTCACTTTCCACCCGGTCGAGCATCTGCACCAGCACCGGCGGGCACACCCCACCGcccaccccctc                                                                                                                                                                                                                               >  2:991497/1‑229 (MQ=11)

                 tCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCCTTCCCCACCCACCCGGTCAG‑CCAGGATTTCACTTTCCAGGCGGTCGAGCATCCGCACCAGACCCGGCCGGAGCACCACACCGCCAcccaccgc                                                                                                                                                                                                                               >  2:1594005/1‑232 (MQ=11)

                      tccttGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAG‑CCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGGACACCACACTGCCCCCCTCTTc‑ccccc                                                                                                                                                                                                                         >  1:1376203/1‑236 (MQ=16)

                      tccttGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAG‑CCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCCGGAAAAACCCACCTCCAccccccttgcccc                                                                                                                                                                                                                          >  1:5417/1‑227 (MQ=11)

                       ccttGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTGCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAG‑CCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑Gccgccg                                                                                                                                                                                                                       <  1:2150997/239‑1 (MQ=18)

                         ttGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAG‑CCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGa                                                                                                                                                                                                                     <  1:321245/239‑1 (MQ=33)

                              ggtcaGGGCCAGCGGCAGGCCGAGATGGTGGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAG‑CCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGActgct                                                                                                                                                                                                                <  2:1320671/235‑1 (MQ=14)

                              gATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAG‑CCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGActgct                                                                                                                                                                                                                <  1:932450/239‑1 (MQ=32)

                                      cgaGCGGCAGCCAGTGCGGGTCGCAGTGGTCCAGGTGGATTTTGCGCGCCGGCGTTTACCCCGGGTCCATCTGGTTTTGCATCTGCTCCCCGGTCAGGCCGCCCGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGTTCAG‑CCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACCCCACACTGCCACCGTCTTC‑GCCGCCGGACTGCAGAAgcgc                                                                                                                                                                                                         <  2:807933/236‑1 (MQ=11)

                                      ccAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAG‑CCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAgcgc                                                                                                                                                                                                         <  1:1823006/238‑1 (MQ=31)

                                      ccAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAG‑CCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAgcgc                                                                                                                                                                                                         <  1:1092681/238‑1 (MQ=31)

                                           ggCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTGGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCCCCCGGTCAG‑CCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCg                                                                                                                                                                                                    <  1:1014886/238‑1 (MQ=14)

                                            gcgggCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCGGGCGGTTACCCCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAG‑CCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGc                                                                                                                                                                                                   <  2:1376203/235‑1 (MQ=9)

                                             cgggCCGGGTTGGTCGCGGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTCCCCCGGGTCCATCTGGTTTTGCATCTGCCCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAG‑CCAGGATTTCCCTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGcc                                                                                                                                                                                                  <  2:647164/236‑1 (MQ=255)

                                                gCCGGGAGGGTCGCAGGGTTTCAGGTGGATTTTTCGCGCCGGCGTGTCCCCCGGGTCCATCTGGTTTTGCATCGGCTCCACGGTCGGGCCGCGCGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAG‑CCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAgg                                                                                                                                                                                               <  2:348730/238‑1 (MQ=11)

                                                 ccGCGATGGTCGCAGTGGTACAGGTGGA‑TTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAG‑CCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCGCCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAgg                                                                                                                                                                                               <  1:2219384/236‑1 (MQ=14)

                                                            gCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAG‑CCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCg                                                                                                                                                                                  <  1:660881/239‑1 (MQ=31)

                                                            gCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAG‑CCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACCGCCCCCGGCTTC‑GCCCCCGGACCGCTGCAGCGCAACCCCCAGGCTGCGGcgccgcc                                                                                                                                                                                  >  1:113972/1‑238 (MQ=11)

                                                                 ggTACAGGTGGTTTTTTCCCGCCGGGGGGTACCCCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCGGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAG‑CCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGtt                                                                                                                                                                                <  2:182248/236‑1 (MQ=9)

                                                                   tACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAG‑CCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTCCCACCGTCCTC‑GCCCCCGGGCTCCTGAAGCGCCTCCGCCAGGCCCCGGCGCCGCGTTTTcc                                                                                                                                                                            >  1:1641542/1‑237 (MQ=11)

                                                                   tACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCGCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAG‑CCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCg                                                                                                                                                                            <  1:1729849/238‑1 (MQ=14)

                                                                    aCAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAG‑CCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGc                                                                                                                                                                           <  1:93108/238‑1 (MQ=31)

                                                                    aCAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAG‑CCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGc                                                                                                                                                                           <  1:184443/238‑1 (MQ=31)

                                                                        gTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAG‑CCAGGATTTCACTTTCCAGCCGGTCGCGCATCTGCACCAGCACCGGCGGGAACACCCCACTGCCCCCGCCCTC‑CCCGCCGGACTGCTGAAGCGCAACCCCCACGCCGCGGCGCCGCGTTTTCCCCCGc                                                                                                                                                                       >  1:396620/1‑238 (MQ=11)

                                                                        gTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAG‑CCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGGACACCACACCGCCCCCGGCCTC‑TCCCCCGGCGCGCCGGAGCGCCGCCCCCACGCGGCGGCGCCGCGGTTTCCCCAGc                                                                                                                                                                       >  2:418437/1‑238 (MQ=11)

                                                                                 tCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAG‑CCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTgg                                                                                                                                                             <  1:569034/239‑1 (MQ=31)

                                                                                  cgcgCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAG‑CCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGc                                                                                                                                                            <  1:79631/239‑1 (MQ=31)

                                                                                  cacgcCGGGGTGTACACCGGGTCCGTCGGGTTTTGCCTCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAG‑CCAGGATTTCCCTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGc                                                                                                                                                            <  2:1228014/237‑1 (MQ=9)

                                                                                       cGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAG‑CCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGttt                                                                                                                                                       <  1:2041135/239‑1 (MQ=32)

                                                                                       cGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAG‑CCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGttt                                                                                                                                                       <  1:2063375/239‑1 (MQ=32)

                                                                                       cGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAG‑CCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGttt                                                                                                                                                       <  1:237254/239‑1 (MQ=32)

                                                                                           gtgtACCCCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCCCTCCCCCGGTCAG‑CCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCaa                                                                                                                                                    <  1:507940/238‑1 (MQ=18)

                                                                                           gtgtACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGAGGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAG‑CCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTAACCGGTGGCGGTTTCaa                                                                                                                                                    <  1:1708350/238‑1 (MQ=18)

                                                                                            tgtACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAG‑CCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCaa                                                                                                                                                    <  1:49580/237‑1 (MQ=33)

                                                                                                 cccggggcactCTGGTTTTCTCTCTGCTCCACGGTCAGGTCGCAGGGGCACGGCCTGCGGTGGCTTTCCTCACTCACCCGGTCGG‑CCGGGATTGCTCTTTCCCGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCCCACTGCCCCCGTCTTC‑GCCGCCGGACTGCTGAGGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCCCCGGTGGCGGTTTCAACTCtgat                                                                                                                                             <  2:1994001/229‑1 (MQ=11)

                                                                                                   cGGGTCCATCTGGTTTGGCATCGGCTCCACGGTCAGCCCCCGCGATGCGCGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAG‑CCAGGATTTCCCTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACGCCACACTGCCACCGTCTTC‑GCCGCCGGGCTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCtgatga                                                                                                                                           <  2:396620/239‑1 (MQ=11)

                                                                                                    gggTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCCG‑CCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAg                                                                                                                                          <  1:2051651/239‑1 (MQ=25)

                                                                                                            cTGGTTTACCCTCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCCGCGGCTTTCCGCACTCCCCCGGTCAG‑CCGGGATTTCACTTGCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGa                                                                                                                                  <  2:745744/239‑1 (MQ=11)

                                                                                                                    gCATCTTCTCCACGGTCCGGGCGCGCGTTGCCAGCCAGCGGCGGCTTTCCTCACTCGCCCGGTCAG‑CCAGGATTTCCCTCGCCAGCCGGTCGAGCATCCGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTcccc                                                                                                                          <  2:1354446/239‑1 (MQ=11)

                                                                                                                      aTCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAG‑CCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCgg                                                                                                                        <  1:1616903/239‑1 (MQ=35)

                                                                                                                      aTCTGCCCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCCCTCGCCCGGTCCG‑CCAGGATTTCCCTTTCCAGCCGGTCGAGCATCTGCGCCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCgg                                                                                                                        <  1:969543/239‑1 (MQ=255)

                                                                                                                              cccGGCCAGGCCGCAGCATGCCTGCCGGCGCCGGCTTGCCTCACCGCCCCGGTCAG‑CCAGGCTTTCACTTTCCAGCCGGTCGGGCAGCTGCACCGGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATaaa                                                                                                                <  2:1632386/237‑1 (MQ=11)

                                                                                                                              cACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAG‑CCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATaaa                                                                                                                <  1:862953/239‑1 (MQ=35)

                                                                                                                                        ccGCCCGGTGCGCGCCGGCGCCGGCTTTCCTCCCTCACCCGGTCAG‑CCAGGATTGCACTTTCCGGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGAt                                                                                                      <  2:1476582/239‑1 (MQ=11)

                                                                                                                                        ccGCACGATGCCAGCCGGCGGCGGCTTTCCTCACTCCCCCGGTCAG‑CCGGGATTTCCCTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGAt                                                                                                      <  2:258178/239‑1 (MQ=16)

                                                                                                                                         cGCACGATGCCAGCCAGCGGCGGCTTTCCTCACCCACCCGGTCAG‑CCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATg                                                                                                     <  1:457124/239‑1 (MQ=18)

                                                                                                                                          ggccgagcccggccggcGGCGCCGTTCCTCCCCCCCCCGGTCAG‑CCCGGGTTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGc                                                                                                    <  2:2066437/228‑1 (MQ=11)

                                                                                                                                          gCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAG‑CCAGGATTTCACTTTCCAGCCGGTCGAGCATCCGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGc                                                                                                    <  1:788141/239‑1 (MQ=18)

                                                                                                                                             ggtacccggcctgcGGCGGGTTACCTCCCCACCCCGGTCAG‑CCAGTATTTCCCTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCggg                                                                                                 <  2:1963785/231‑1 (MQ=11)

                                                                                                                                             cGATGCCAGCCAGCGGCGGCTTTCCTCACTCCCCCGGTCAG‑CCAGGATTTCCCTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCggg                                                                                                 <  1:1258667/239‑1 (MQ=18)

                                                                                                                                             cGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAG‑CCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCggg                                                                                                 <  1:991497/239‑1 (MQ=33)

                                                                                                                                             cGATGCCAGCCAGCGGCGGCTTTCCCCACTCCCCCGGTCAG‑CCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCggg                                                                                                 <  1:2199824/239‑1 (MQ=18)

                                                                                                                                                tGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAG‑CCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGtt                                                                                               <  1:633745/238‑1 (MQ=33)

                                                                                                                                                tGCCAGCCAGCGGCGGCGTTCCCCACTCACCCGGTCAG‑CCAGGATTTCACTTTCCAGCCGGTCTAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGTCTGATAAATCGTCTGGATGCGGGTTc                                                                                              <  2:1764283/239‑1 (MQ=16)

                                                                                                                                                   cagccagcGGCAGCTTTCCTCACCCACCCGGTCAG‑CCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTg                                                                                            <  1:1500384/238‑1 (MQ=18)

                                                                                                                                                       cagcGGCGGCTTTCCTCCCTCCCCCGGTCAG‑CCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGtt                                                                                       <  2:1576199/239‑1 (MQ=18)

                                                                                                                                                       cagcGGCGGCTTTCCTCCCTCACCCGGTCAG‑CCAGGATTTCCCTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGtt                                                                                       <  1:1175083/239‑1 (MQ=18)

                                                                                                                                                       cagcGGCGGCTTTCCTCACTCCCCCGGTCAG‑CCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTc                                                                                          <  1:2156532/236‑1 (MQ=18)

                                                                                                                                                       cagcGGCGGCTTTCCTCACTCACCCGGTCAG‑CCAGGATTTCCCTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGtt                                                                                       <  1:957634/239‑1 (MQ=18)

                                                                                                                                                       cagcGGCGGCTTTCCTCACTCACCCGGTCAG‑CCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACTCTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGtt                                                                                       <  1:1550397/239‑1 (MQ=18)

                                                                                                                                                       cagcGGCGGCTTTCCTCACTCACCCGGTCAG‑CCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGtt                                                                                       <  1:1476431/239‑1 (MQ=33)

                                                                                                                                                          cggcggCTTTCCTTACTCACCCGGTCAG‑CCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCt                                                                                     <  1:1748743/238‑1 (MQ=18)

                                                                                                                                                          cggcggCTTTCCTCACTCACCCGGTCAG‑CCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGttt                                                                                                                                                       >  1:103578/1‑172 (MQ=33)

                                                                                                                                                           ggcggcTTTCCTCACTCACCCGGTCAG‑CCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCt                                                                                     <  1:1905956/237‑1 (MQ=33)

                                                                                                                                                                   tCCTCACTCACCCGGTCAG‑CCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCCCCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGt                                                                           <  1:1226869/239‑1 (MQ=25)

                                                                                                                                                                         cctcccccGTCAGCCCAGGATTTCCCTTCCCCGCCGGACGGGCATCCGCACCAGCACCGGCGGGAACACCACCCTGCCGCCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGcc                                                                      <  2:1495475/235‑1 (MQ=11)

                                                                                                                                                                          tcacCCGGTCAG‑CCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCCTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGGGGCGTGAAGCTCCCCCGCTGATAAAACTTCTGGATGCGGGTTCTGTCGTTCTGTTTTGTGGGCAGCCgg                                                                    >  1:710959/1‑239 (MQ=255)

                                                                                                                                                                            ccccGGTC‑G‑CCAGGATTTGCCTGTCCAGCCGGTCGCGCATCTGCACCAGCACCGGCGGGAACACCACACCGCCACCGTCTTC‑CCCGCCGTACTGCTCAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGT‑Cg                                                                <  2:1118838/238‑1 (MQ=11)

                                                                                                                                                                            acCCGGTCAG‑CCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGT‑c                                                                 <  1:2089223/239‑1 (MQ=37)

                                                                                                                                                                              ccGGTCAG‑CCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAAct                                                                                                                                                  <  2:468962/157‑1 (MQ=34)

                                                                                                                                                                                 gTCAG‑CCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCCGCTCACCGGGGGCGGGTTCCACTCTGATGAGGGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGGTCTGTCGTTCTGTATTGTGGTCAGCCGGG‑CGCCg                                                             >  1:270968/1‑238 (MQ=255)

                                                                                                                                                                                  tCAT‑CCTGGATTGCACTTTCCCGCCGCTCCAGGACCCGCACCAGCCCCGGCGGGCTCCCCCCACGGCCCCCGTCTCC‑TCCTCCGGACTGCGGAAGCGCATCCGCCGGGCTGCGGCGCTGCGTTTTCGCCAGCTCCCCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCGGGATGCGGGTACTGTCTTTCTGTATTGTGGTCAGCCGGT‑CGCCGTc                                                           <  2:544750/239‑1 (MQ=11)

                                                                                                                                                                                  tCAG‑CCAGGATTTCCCTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGT‑CGCCGt                                                            <  1:1244158/238‑1 (MQ=37)

                                                                                                                                                                                   cAG‑CCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCCCCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGGGCGGGGTTTTCGCCAGCTCACCGGGGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGG‑CGCCGTc                                                           <  2:307296/238‑1 (MQ=255)

                                                                                                                                                                                   cAG‑CCAGGATTTCACTTCCCCGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCCCACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGT‑CGCCGTcc                                                          <  1:2020994/239‑1 (MQ=255)

                                                                                                                                                                                       ccgggttctcccTTCCCGGCCGGTCGGGCATCTGCACCAGCACCGGCGGGAACCCCCCACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGT‑CGCCGTCCCAg                                                       <  2:710959/228‑1 (MQ=11)

                                                                                                                                                                                       ccAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACCGCCACCGTCTTC‑GCCGCCCGACCGCTGAAGCGCCTCCGCCCGGCTGCGGCGCTGCGTTTTCGCCAGCTCCCCCGCGGCGGGTTCAACTCCGGTGAGTGGCGTGAAGGCCCCCCGCCGGTAAATCGTCTGGATGCGGGTTTTGTCGTTCTGTCTTTTTGGCTGCCGGTCCGCCCTccc                                                         >  2:869236/1‑238 (MQ=11)

                                                                                                                                                                                            ttttCACTTCCCAGCCGGTCGGGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGT‑CGCCGTCCCAGCCg                                                    <  1:1070616/236‑1 (MQ=37)

                                                                                                                                                                                             tttCCCTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGT‑CGCCGTCCCAGCCGTAc                                                 <  1:1684194/239‑1 (MQ=37)

                                                                                                                                                                                             tttCACTTCCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCCCACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGT‑CGCCGTCCCAGCCGTAc                                                 <  1:1116310/239‑1 (MQ=37)

                                                                                                                                                                                              ttCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGT‑CGCCGTCCCAGCCGTAcc                                                <  1:1332595/239‑1 (MQ=37)

                                                                                                                                                                                                ctcccacggcccGGTCGAGCTTGCGCGCCACCCCCGGGGTGACAGCCACGCTGCCACCGTGTTC‑GCCGGCGCCCTCCGTATGCGCATCGGCCGGGCTGCGGCCGTCCTTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATTAGTGGCGTGAAGTTCCCCGGCGGATAAGTCCTCTGGATGCGGTTCCTGTCGTTCTGTATTGTGGTCCGCCGGT‑CGCCGTCCCAGCCGTACCAg                                              <  2:2110142/229‑1 (MQ=2)

                                                                                                                                                                                                cACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGT‑CGCCGTCCCAGCCGTACCAg                                              <  1:1959824/239‑1 (MQ=37)

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

                                                                                                                                                                                                            ggTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGT‑CGCCGTCCCAGCCGTACCAGGTCACTTGCGGt                                  <  1:82158/239‑1 (MQ=37)

                                                                                                                                                                                                                gAGCATCTGCACCAGCCCCGGCGGGAACACCACACTGCCCCCGTCTTC‑GCCGCCGGCCTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGT‑CGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTcc                              <  1:727290/239‑1 (MQ=255)

                                                                                                                                                                                                                 ggCAGCTGCTCCAGCCCCGGCGGGAACACCCCCCTGCCACCGTCTTC‑GCCGCCTGACTTCTGAAGCGCAGCCGCCAGGCTGGGGGGCGGCGTTTTCGCCAGCTTAGGGGGGGTGGGTTAAGTATTGATGAGTGGCGTGAAGCTCCCAGGCGGGTAACTCGTCTGGATGCGGGTCCTTTCGTTCGGTATTGTGGTCAGCCGGT‑CGCCGCCCCAGCCGTACCAGGTCACTTGCGGTTTcc                              <  2:909088/237‑1 (MQ=11)

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

                                                                                                                                                                                                                          aCCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGT‑CGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCg                     <  1:1423813/238‑1 (MQ=37)

                                                                                                                                                                                                                                  cggcggGACCACCACACGGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCCCCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGGCAGCCGGT‑CGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGc            <  1:1871480/239‑1 (MQ=255)

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

                                                                                                                                                                                                                                      gggAACACCCCACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGT‑CGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGt        <  1:463664/239‑1 (MQ=37)

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

TGCGCACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAG‑CCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTC‑GCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGT‑CGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTC  >  NC_007779/3874729‑3875202

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