BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_007779	3,874,971	T→C	6.7%	S1088G (AGT→GGT)	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,971	0	T	C	6.7%	605.4 / 2.6	263	S1088G (AGT→GGT)	Y75_RS19225	RHS element protein RhsA
Reads supporting (aligned to +/- strand): ref base T (15/218); new base C (21/2); total (37/226)
Fisher's exact test for biased strand distribution p-value = 3.94e-19
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 269 total aligned reads displayed.

ACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTT‑CACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTCCCGTTCA  >  NC_007779/3874734‑3875209

          tgtCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCGCCCGGTCAGCCAGGATTT‑CACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑Ct‑t                                                                                                                                                                                                                                     <  1:1571288/239‑1 (MQ=25)

          tgtCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTT‑CACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACCCTGccccactct                                                                                                                                                                                                                                        >  1:2033042/1‑231 (MQ=25)

          tgtCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTT‑CACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑Ct‑t                                                                                                                                                                                                                                     <  1:76226/239‑1 (MQ=35)

          tgtCCCTTCCTGGCTGATAGGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTT‑CACTTT‑CCAGCCGGTCGAGCATCGGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑Ct‑t                                                                                                                                                                                                                                     <  1:170513/239‑1 (MQ=25)

            tCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCCTTCCCCACCCACCCGGTCAGCCAGGATTT‑CACTTT‑CCAGGCGGTCGAGCATCCGCACCAGACCCGGCCGGAGCACCACACCGCCAcccaccgc                                                                                                                                                                                                                                      >  2:1594005/1‑232 (MQ=11)

              ccttccttGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTT‑CACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACAACCCACATCCACCgcccttcc                                                                                                                                                                                                                                    >  1:641763/1‑231 (MQ=255)

              ccttccttGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTT‑CACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAAAACCACCCTGCCCCCggcctccc                                                                                                                                                                                                                                    >  1:895782/1‑231 (MQ=255)

                 tccttGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTT‑CACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGGACACCACACTGCCCCCCT‑CT‑Tcccccc                                                                                                                                                                                                                               >  1:1376203/1‑236 (MQ=16)

                 tccttGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTT‑CACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCCGGAAAAACCCACCTCCAccccccttgcccc                                                                                                                                                                                                                                 >  1:5417/1‑227 (MQ=11)

                    ttGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTT‑CACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGa                                                                                                                                                                                                                           <  1:321245/239‑1 (MQ=33)

                     tGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTT‑CACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGGACACCCCACCGCCCCCCTCCT‑TCCCCCCCgg                                                                                                                                                                                                                            >  1:183982/1‑238 (MQ=11)

                          aTAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTT‑CACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGGACACCACACCGCCCCCCT‑CC‑TCCCCCCCGGActgctt                                                                                                                                                                                                                     >  1:714491/1‑238 (MQ=255)

                                 ccAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTT‑CACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAgcgc                                                                                                                                                                                                               <  1:1092681/238‑1 (MQ=31)

                                      ggCGGGCCGCGGTGGCTGCAGTGGTACGGGTGGATTTTTCCCGCCGGGGTGTCCACCGGGTCCCTCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTT‑CCCTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGc                                                                                                                                                                                                         <  2:148142/239‑1 (MQ=11)

                                      ggCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTGGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCCCCCGGTCAGCCAGGATTT‑CACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCg                                                                                                                                                                                                          <  1:1014886/238‑1 (MQ=14)

                                        cgggCCGGGTTGGTCGCGGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTCCCCCGGGTCCATCTGGTTTTGCATCTGCCCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTT‑CCCTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGcc                                                                                                                                                                                                        <  2:647164/236‑1 (MQ=255)

                                           gCCGGGAGGGTCGCAGGGTTTCAGGTGGATTTTTCGCGCCGGCGTGTCCCCCGGGTCCATCTGGTTTTGCATCGGCTCCACGGTCGGGCCGCGCGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTT‑CACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAgg                                                                                                                                                                                                     <  2:348730/238‑1 (MQ=11)

                                                gATGGTCGCAGTGGTGCAGGTGGATTTTTCCGGCGCGCGGGTCCACCGGGTCCATTTGGTTTTGCTTCTGCTCCACGGCCAGGCCGCCCGATGCCAGCCAGCGGCGGCTTTCCTCACTCCCCCGGTCAGCCAGGATTT‑CACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTgcg                                                                                                                                                                                                <  2:364773/238‑1 (MQ=11)

                                                       gCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTT‑CACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCg                                                                                                                                                                                        <  1:660881/239‑1 (MQ=31)

                                                       gCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTT‑CACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACCGCCCCCGG‑CT‑TCGCCCCCGGACCGCTGCAGCGCAACCCCCAGGCTGCGGcgccgcc                                                                                                                                                                                        >  1:113972/1‑238 (MQ=11)

                                                            ggTACAGGTGGTTTTTTCCCGCCGGGGGGTACCCCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCGGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTT‑CACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGtt                                                                                                                                                                                      <  2:182248/236‑1 (MQ=9)

                                                              tACAGGTGGATTTTTCGCGCCGGCGTGTGCACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTT‑CACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGc                                                                                                                                                                                 <  1:1594005/239‑1 (MQ=14)

                                                              tACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTT‑CACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTCCCACCGT‑CC‑TCGCCCCCGGGCTCCTGAAGCGCCTCCGCCAGGCCCCGGCGCCGCGTTTTcc                                                                                                                                                                                  >  1:1641542/1‑237 (MQ=11)

                                                               aCAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTT‑CACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGc                                                                                                                                                                                 <  2:1102589/238‑1 (MQ=31)

                                                                             cgcgCCGGGGTGTACACCGGTGCCATCTGGTTTTGCATCTGCCCCGCGGTCAGGCCGCCCGCTGCCGGCCAGCGGCGGTTTTCCGCACTCCCCGGGTCAGCCAGGATTT‑CACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACGCCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCGGCGTTTTCGCCAGCTCACCGGTGGc                                                                                                                                                                  <  2:421652/239‑1 (MQ=11)

                                                                             cgcgCCGGCGTTCACCCCGGGTCCCTCTGTTTTTCCATCTGCTCCACGGTCAGGCCGCACGGTGCCAGCCAGCGGCGGCTTTCCTCCCTCACCCGGTCAGCCAGGATTT‑CACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGc                                                                                                                                                                  <  2:113972/239‑1 (MQ=11)

                                                                             cgcgCCGGCGTGTCCACCGGGTCCATCGGGTTTTGCCTCTCCTCCGCGGTCAGGCCGCCGTATGCCAGCCAGCGGCGGCTTTCCTCACCCACCCTGTCAGCCAGGATTT‑CCCTTC‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGc                                                                                                                                                                  <  2:1275174/239‑1 (MQ=11)

                                                                             cgcgCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCCCGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTT‑CACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGc                                                                                                                                                                  <  1:1249862/239‑1 (MQ=14)

                                                                                gcCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTT‑CACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGt                                                                                                                                                               <  1:194110/239‑1 (MQ=32)

                                                                                  cgtggtgtACCCCGGGTCCATCTGGTTTTGCATCTGCTCCGCGGTCAGGCCGCACGATGCCAGCCGGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTT‑CGCTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGttt                                                                                                                                                             <  2:705580/235‑1 (MQ=12)

                                                                                  cgggatggacacGGGGCCCATCTGGTGTGCCATCTCCCCGACGGTAAGGCCCGACGATGCCAGCCAGCGGCGGCTTTCCTCCCTCTCCCGGTCAGCCAGAATTT‑CACTTC‑CCAGCCGGTCGAGCATCTGCACCAGCCCCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGCCTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGttt                                                                                                                                                             <  2:717586/231‑1 (MQ=11)

                                                                                  cGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCCCCCGGTCAGCCAGGATTT‑CCCTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGttt                                                                                                                                                             <  1:2209727/239‑1 (MQ=17)

                                                                                  cGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTT‑CACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGttt                                                                                                                                                             <  1:2063375/239‑1 (MQ=32)

                                                                                      gtgtACCCCGGGCCCATCTGTTTTTGCGTTTGTCCGCCGGGCAGGCCGCCCGTTGCCAGCCAGCGGCGGTTTTCCTCACTCACCCGGTCAGCCGGGATTT‑CGCGTT‑CCAGCCGGTCGTGCATCTGCACCAGCACCGGCGGGAACACCGCACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCaa                                                                                                                                                          <  2:1325542/238‑1 (MQ=11)

                                                                                      gtgtACCCCGGGCCCATCTGTTTTTGCATCTGCTCCACGGTCAGGCCGCACGAGGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTT‑CCCTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCaa                                                                                                                                                          <  2:842588/238‑1 (MQ=16)

                                                                                      gtgtACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCCCGAGGCCAGCCAGCGGCGGCTTTCCTCCCTCACCCGGTCAGCCAGGATTT‑CCCTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCaa                                                                                                                                                          <  2:1380740/238‑1 (MQ=16)

                                                                                            cccggggcactCTGGTTTTCTCTCTGCTCCACGGTCAGGTCGCAGGGGCACGGCCTGCGGTGGCTTTCCTCACTCACCCGGTCGGCCGGGATTG‑CTCTTT‑CCCGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCCCACTGCCCCCGT‑CT‑TCGCCGCCGGACTGCTGAGGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCCCCGGTGGCGGTTTCAACTCtgat                                                                                                                                                   <  2:1994001/229‑1 (MQ=11)

                                                                                               gggTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCCGCCAGGATTT‑CACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAg                                                                                                                                                <  1:2051651/239‑1 (MQ=25)

                                                                                                ggTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTT‑CACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGG‑CT‑TCGCCCCCCGACTGCTGAAGCGCATCCGCCAGGCCGCGGCGCGGCGGTTTCGCCCGCTCACCCGTGGCGGGTTCAACTCTGTTGAgg                                                                                                                                               >  1:1343362/1‑238 (MQ=11)

                                                                                                  tCCATCTGGTTTTGCATCTGCTCCCCGGTCAGGCCGCACGATTCCAGCCAGCGGCGGCTTTCCTCACTCCCCCGGTCAGCCAGGATTT‑CACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCTGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTg                                                                                                                                              <  1:1718919/238‑1 (MQ=255)

                                                                                                       gtGGTTTTGCTTCTGCTCGACGGCCAGGCCTCACGATCCCAGCCCGCGGCGGTTTTCCTCCCTCCCCCGGTCAGCCAGGATTT‑CACTTC‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGa                                                                                                                                        <  2:1525438/238‑1 (MQ=11)

                                                                                                              tGCATCTGCCCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTT‑CACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTccc                                                                                                                                 <  1:524741/239‑1 (MQ=25)

                                                                                                                    tGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCCCCCGGTCAGCCAGGATTT‑CACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCctgctg                                                                                                                           <  1:586192/239‑1 (MQ=25)

                                                                                                                         cccGGCCAGGCCGCAGCATGCCTGCCGGCGCCGGCTTGCCTCACCGCCCCGGTCAGCCAGGCTTT‑CACTTT‑CCAGCCGGTCGGGCAGCTGCACCGGCACCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATaaa                                                                                                                      <  2:1632386/237‑1 (MQ=11)

                                                                                                                               ctggCCTCACGATGCCAGCCAGCGGCGGCTTTCCCCACTCACCCGGTCAGCCAGGATTT‑CACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACGCCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTc                                                                                                                 <  1:1193351/236‑1 (MQ=18)

                                                                                                                                   ccGCCCGGTGCGCGCCGGCGCCGGCTTTCCTCCCTCACCCGGTCAGCCAGGATTG‑CACTTT‑CCGGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGAt                                                                                                            <  2:1476582/239‑1 (MQ=11)

                                                                                                                                   ccGCACGAGGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTT‑CACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTg                                                                                                               <  1:505017/236‑1 (MQ=18)

                                                                                                                                     ggccgagcccggccggcGGCGCCGTTCCTCCCCCCCCCGGTCAGCCCGGGTTT‑CACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGc                                                                                                          <  2:2066437/228‑1 (MQ=11)

                                                                                                                                     gCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTT‑CACTTT‑CCAGCCGGTCGAGCATCCGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGc                                                                                                          <  1:788141/239‑1 (MQ=18)

                                                                                                                                        ggtacccggcctgcGGCGGGTTACCTCCCCACCCCGGTCAGCCAGTATTT‑CCCTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCggg                                                                                                       <  2:1963785/231‑1 (MQ=11)

                                                                                                                                        cGATGCCAGCCAGCGGCGGCTTTCCTCACTCCCCCGGTCAGCCAGGATTT‑CCCTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCggg                                                                                                       <  1:1258667/239‑1 (MQ=18)

                                                                                                                                        cGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTT‑CACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCggg                                                                                                       <  1:991497/239‑1 (MQ=33)

                                                                                                                                        cGATGCCAGCCAGCGGCGGCTTTCCCCACTCCCCCGGTCAGCCAGGATTT‑CACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCggg                                                                                                       <  1:2199824/239‑1 (MQ=18)

                                                                                                                                           tGCCAGCCAGCGGCGGCGTTCCCCACTCACCCGGTCAGCCAGGATTT‑CACTTT‑CCAGCCGGTCTAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGTCTGATAAATCGTCTGGATGCGGGTTc                                                                                                    <  2:1764283/239‑1 (MQ=16)

                                                                                                                                              cagccagcGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTT‑CACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCCCACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTg                                                                                                  <  1:978389/238‑1 (MQ=18)

                                                                                                                                              cagccagcGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTT‑CACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTTCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTg                                                                                                  <  1:1114452/238‑1 (MQ=18)

                                                                                                                                              cagccagcGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTT‑CACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTg                                                                                                  <  1:418437/238‑1 (MQ=33)

                                                                                                                                              cagccagcGGCAGCTTTCCTCACCCACCCGGTCAGCCAGGATTT‑CACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTg                                                                                                  <  1:1500384/238‑1 (MQ=18)

                                                                                                                                                  cggccgcggcTTTCCCCCCCCACCCGGTCAGCCAGGATTT‑CCCTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACCCCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTc                                                                                                <  2:77026/231‑1 (MQ=255)

                                                                                                                                                  cagcGGCGGCTTTCCTCCCTCACCCGGTCAGCCAGGATTT‑CCCTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGtt                                                                                             <  1:1175083/239‑1 (MQ=18)

                                                                                                                                                  cagcGGCGGCTTTCCTCACTCCCCCGGTCAGCCAGGATTT‑CACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTc                                                                                                <  1:2156532/236‑1 (MQ=18)

                                                                                                                                                  cagcGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTT‑CACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACTCTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGtt                                                                                             <  1:1550397/239‑1 (MQ=18)

                                                                                                                                                  cagcGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTT‑CACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTc                                                                                                <  1:9411/236‑1 (MQ=33)

                                                                                                                                                  cagcGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTT‑CACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGtt                                                                                             <  1:1476431/239‑1 (MQ=33)

                                                                                                                                                     cggcggCTTTCCTTACTCACCCGGTCAGCCAGGATTT‑CACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCt                                                                                           <  1:1748743/238‑1 (MQ=18)

                                                                                                                                                     cggcggCTTTCCTCACTCACCCGGTCAGCCAGGATTT‑CACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGttt                                                                                                                                                             >  1:103578/1‑172 (MQ=33)

                                                                                                                                                      ggcggcTTTCCTCACTCACCCGGTCAGCCAGGATTT‑CACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTc                                                                                            <  1:958221/236‑1 (MQ=33)

                                                                                                                                                      ggcggcTTTCCTCACTCACCCGGTCAGCCAGGATTT‑CACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCt                                                                                           <  1:1905956/237‑1 (MQ=33)

                                                                                                                                                        cggcTTTCCTCACTCACCCGGTCAGCCAGGATTT‑CACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTa                                                                                        <  1:1866475/238‑1 (MQ=33)

                                                                                                                                                         ggcTTTCCTCCCCCCCCCGGTCAGCCAGGATTT‑CCCTTT‑CCAGCCGGTCGGGCATCTGCACCAGCCCCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTAtt                                                                                      <  2:1062731/239‑1 (MQ=255)

                                                                                                                                                         ggcTTTCCTCACTCACCCGGTCAGCCAGGATTT‑CACTTT‑CCAGCCGGTCGGGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTAtt                                                                                      <  1:2020777/239‑1 (MQ=21)

                                                                                                                                                              tCCTCACTCACCCGGTCAGCCAGGATTT‑CACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCCCCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGt                                                                                 <  1:1226869/239‑1 (MQ=25)

                                                                                                                                                                  cactcacCCGGTCAGCCAGGATTT‑CACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGc                                                                             <  1:7997/239‑1 (MQ=35)

                                                                                                                                                                       acCCGGTCAGCCAGGATTT‑CACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTc                                                                        <  1:2089223/239‑1 (MQ=37)

                                                                                                                                                                        cccgggcgcccAGGTTTC‑CACTTT‑CCAGCCGGCCGGGCATCTGCACCAGCACCGGCGGGAACCCCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCg                                                                       <  2:1357387/230‑1 (MQ=11)

                                                                                                                                                                        cccGGTCAGCCAGGATTT‑CACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCg                                                                       <  1:493356/239‑1 (MQ=37)

                                                                                                                                                                        cccGGTCAGCCAGGATTT‑CACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCg                                                                       <  1:235806/239‑1 (MQ=37)

                                                                                                                                                                         ccGGTCAGCCAGGATTT‑CACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAAct                                                                                                                                                        <  2:468962/157‑1 (MQ=34)

                                                                                                                                                                            gTCAGCCAGGATTT‑CACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCCGCTCACCGGGGGCGGGTTCCACTCTGATGAGGGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGGTCTGTCGTTCTGTATTGTGGTCAGCCGGGCGCCg                                                                    >  1:270968/1‑238 (MQ=255)

                                                                                                                                                                              cAGCCCGGGTTTCCACTTT‑CCAGCCGGCCGAGCATCTGCCCCGGCCCCGGCGGGACCCCCCCCCGGCCACCGT‑CT‑TCGCCGCCGGACTGCGGAGGCGCATCCGCCAGGCTGCGGGGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTc                                                                  <  2:270968/239‑1 (MQ=11)

                                                                                                                                                                              cAGCCAGGATTT‑CCCTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACCCCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAGCTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTc                                                                  <  1:1821921/238‑1 (MQ=255)

                                                                                                                                                                              cAGCCAGGATTT‑CACTTT‑CCAGCCGGTCGAGCATCTGCCCCAGCACCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGGGCGGGGTTTTCGCCAGCTCACCGGGGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGGCGCCGTc                                                                  <  2:307296/238‑1 (MQ=255)

                                                                                                                                                                                 ccagctttt‑CACTTT‑CCGGCCGGTCGGGCACCTGCACCAGCACCGGCGGGACCACCCCACTGCCACCGT‑CT‑CCGCCGCCGGACGGCTGAGGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGAGGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAg                                                              <  2:411877/233‑1 (MQ=11)

                                                                                                                                                                                 ccAGGATTT‑CCCTTC‑CCCGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAg                                                              <  2:895782/239‑1 (MQ=255)

                                                                                                                                                                                   ccgATTT‑TCCTTTCCCAGCGGGTCGGGCATCTGCCCCGGCACCGGCGGGACCCCCCCCCTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAGGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTccc                                                                <  2:2100160/234‑1 (MQ=11)

                                                                                                                                                                                      tttt‑CACTTC‑CCAGCCGGTCGGGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCg                                                           <  1:1070616/236‑1 (MQ=37)

                                                                                                                                                                                       ttt‑CCCTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTAc                                                        <  1:1684194/239‑1 (MQ=37)

                                                                                                                                                                                           cccTTC‑CCGGCCGGCCGGCCATCGGCCCCAGCACCGGCGGGAACACCCCCCTGCCCCCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAg                                                     <  2:983711/237‑1 (MQ=11)

                                                                                                                                                                                           cACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCgg                                                                                                                              <  1:1315546/166‑1 (MQ=35)

                                                                                                                                                                                            ccttccccAGCCGGCCGAGCCTCTGCCCCGGCCGCGGCGGCACCGCCCCCGTGCC‑CCCT‑CTGTCGCCGCCGGATGGCTGAAGCTCGTCCGCCGGGCTGCGGCGCGGCGTTTTCGCCAGCTCACCGGTGGTGGTTTCAACTCTGATGAGTGGCGTGAAGCCCCCCGGCTGATAAGGCGGCGTGAGGCGGGTTTTATCGTTCTGTATTGTGGTCAGGCGCGCGCGGGCCCAGCCGTTCCAg                                                     <  2:1203616/233‑1 (MQ=255)

                                                                                                                                                                                              ttt‑CCAGCCGGTCGAGCATCTGCCCCAGCACCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTc                                                  <  1:1287482/239‑1 (MQ=37)

                                                                                                                                                                                                       cGGTCGAGCATCTGCACCAGCCCCGGCGGGAACCCCACCCTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTg                                             <  1:932133/236‑1 (MQ=255)

                                                                                                                                                                                                        ggTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGt                                         <  1:82158/239‑1 (MQ=37)

                                                                                                                                                                                                            gAGCATCTGCACCAGCCCCGGCGGGAACACCACACTGCCCCCGT‑CT‑TCGCCGCCGGCCTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTcc                                     <  1:727290/239‑1 (MQ=255)

                                                                                                                                                                                                                    gccccAGCCCCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGc                             <  1:2110923/236‑1 (MQ=37)

                                                                                                                                                                                                                      aCCAGCACCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCg                            <  1:1423813/238‑1 (MQ=37)

                                                                                                                                                                                                                              cggcggGAACACCACTCTGCCCCCGT‑CG‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCGGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGc                   <  1:673634/239‑1 (MQ=255)

                                                                                                                                                                                                                                  gggAACACCCCACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGt               <  1:463664/239‑1 (MQ=37)

                                                                                                                                                                                                                                        cccaccgtGCCCCCGT‑AT‑TCCCCCCCGGACGGCTAATCCACGTCCGCCGGCCGCCGGCGCTGCGTTTTCGCCGGCTCACCGGTGCCGTTTTCAACTCGGATGAGTGCCGGGAACCTCCCCGGTCGATAAAGCGGCGGGATCCGGGTTCTGCCGTTCTGTATTGTGGTCAGCCGTTCGCCGCCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTc         <  2:5417/232‑1 (MQ=2)

                                                                                                                                                                                                                                               tGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGAGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTGTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGCGACAGCGACATCCAGCCCGTCAGGTCCCGTTc   >  1:282375/1‑238 (MQ=255)

ACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTT‑CACTTT‑CCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGT‑CT‑TCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCCGGTCGCCGTCCCAGCCGTACCAGGTCACTTGCGGTTTCCGTGACAGCGACATCCAGCCCGTCAGGTCCCGTTCA  >  NC_007779/3874734‑3875209

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