BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_007779	3,953,876	A→G	7.3%	S879G (AGC→GGC)	bcsC →	cellulose biosynthesis protein BcsC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007779	3,953,876	0	A	G	7.3%	513.3 / 5.2	200	S879G (AGC→GGC)	bcsC	cellulose biosynthesis protein BcsC
Reads supporting (aligned to +/- strand): ref base A (78/102); new base G (14/4); total (94/106)
Fisher's exact test for biased strand distribution p-value = 6.13e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 200 total aligned reads displayed.

GGTTACTCCGATCTGAAAGCGCACACTACCATGTTGCAGGTGGATGCGCCGTATTCTGACGGGCGGATGTTCTTTCGCAGTGATTTCGTCAATATGAACGTCGGCAGTTTCTCCACTAATGCCGATGGCAAATGGGATGACAACTGGGGCACCTGTACATTACAGGACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGATATCGGTACC‑ACGCCGATGGGCTTCAACGTGGT‑GGATGTGGTCGGCGGCATCAGTTACAGCG‑ATGAT‑ATCGGGCCGCTGGGTTACACCGTTAACGCCCACCGTCGGCCCATCTCCAG‑TTCTTTGCTGGCCTTTGGTGGGC‑AAAAAGAC‑TCCCCGAGCAATACCGGGAAAAAATGGGGTGGCGTACGTGCCGACGGTGTGGGGCTAAGTCTGAGCTACGATAAAGGT  >  NC_007779/3953639‑3954109

ggtTACTCCGATCTGAAAGCGCACACTACCATGTTGCAGGTGGATGCGCCGTATTCTGACGGGCGGATGTTCTTTCGCAGTGATTTCGTCAATATGAACGTCGGCAGTTTCTCCACTAATGCCGATGGCAAATGGGATGACAACTGGGGCACCTGTACATTACAGGACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAg                                                                                                                                                                                                                                                 <  1:25704/239‑1 (MQ=255)

 gtTACTCCGATCTGAAAGCGCACACTACCATGTTGCAGGTGGATGCGCCGTATTCTGACGGGCGGATGTTCTTTCGCAGTGATTTCGTCAATATGAACGTCGGCAGTTTCTCCACTAATGCCGATGGCAAATGGGATGACAACTGGGGCACCTGTACATTACAGGACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAg                                                                                                                                                                                                                                                 <  2:1765952/238‑1 (MQ=255)

  ttgcTCCGATCTGAAAGCGCACACTACCATGTTGCAGGTGGATGCGCCGTATTCTGACGGGTGGATGTTCTTGCGCAGTGATTTCGTCAATATGAACGTCTGCAGTTTCTCCACTAATGCCGATGGCAAATGGGATGACAACTGGGGCACCTGTACATTACAGGACTGTGGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCt                                                                                                                                                                                                                                               <  2:199885/236‑1 (MQ=255)

     cTCCGAGCTGAAACCGCACGCTACAAAGTTACAGGTGGAGGCGCCGTATGCTGACGGAGGGGTGTTCGTTCGCAGTGATTTCGTTAATATGAACGTCGGCAGTTACTCCACTAATGAGGATGTCAAATGGGATGACAGCTGGGGCTCCTGTACATTACAGGATTGTAGCGGCAACCGCAGCCAGTCGTATTGCGGTGCAAGCGTTGCGGTCGGCTGGCGAAATGACGTCTGGAGCTggg                                                                                                                                                                                                                                            <  2:1308259/239‑1 (MQ=255)

         gATCTGAAAGCGCACACTACCATGTTGCAGGTGGATGCGCCGTATTCTGACGGGCGGATGTTCTTTCGCAGTGATTTCGTCAATATGAACGTCGGCAGTTTCTCCACTAATGCCGATGGCAAATGGGATGACAACTGGGGCACCTGTACATTACAGGACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGatat                                                                                                                                                                                                                                        <  1:1266647/239‑1 (MQ=255)

         gATCTGAAAGCGCACACTACCATGTTGCAGGTGGATGCGCCGTATTCTGACGGGCGGATGTTCTTTCGCAGTGATTTCGTCAATATGAACGTCGGCAGTTTCTCCACTAATGCCGATGGCAAATGGGATGACAACTGGGGCACCTGTACATTACAGGACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGata                                                                                                                                                                                                                                         >  1:2128307/1‑238 (MQ=255)

                    gcACACTACCATGTTGCAGGTGGATGCGCCGTATTCTGACGGGCGGATGTTCTTTCGCAGTGATTTCGTCAATATGAACGTCGGCAGTTTCTCCACTAATGCCGATGGCAAATGGGATGACAACTGGGGCACCTGTACATTACAGGACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGATATCGGTACC‑ACg                                                                                                                                                                                                                             >  1:1050813/1‑238 (MQ=255)

                     cacaCTACCATGTTGCAGGTGGATGCGCCGTATTCTGACGGGCGGATGTTCTTTCGCAGTGATTTCGTCAATATGAACGTCGGCAGTTTCTCCACTAATGCCGATGGCAAATGGGATGACAACTGGGGCACCTGTACATTACAGGACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGATATCGGTACC‑AAGcc                                                                                                                                                                                                                           >  1:1941560/1‑239 (MQ=255)

                                  tGCAGGTGGATGCGCCGTATTCTGACGGGCGGATGATCTTTCGCAGTGATTTCGTCAATATGAACGTCGGCAGTTTCTCCACTAATGCCGATGGCAACTGGGATGACAACTGGGGCACCTGTACATTACAGGACTGTAGCGGCAACCGCAGCCAGTCGGCTTCCGGTGCCAGCGTGGCGGACGGCCGGCGAAATGACGTGTGGAGCGGGGGTCTCGGTAAC‑ACGCCGCTGGGCCTCCAc                                                                                                                                                                                                              >  1:1474732/1‑239 (MQ=255)

                                   gCAGGTGGGTGCGCCGTATGATGACGGGCGGATGTTCTTTCGCAGTGATTTCGTCAATATGAACGTCGGCAGTTTCTCCACTAATGCCGATGGCAAATGGGATGACAACTGGGGCACCTGTACATTACAGGACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGATATCGGTACC‑ACGCCGATGGGCTTCAACg                                                                                                                                                                                                             <  2:1941560/239‑1 (MQ=255)

                                    cAGGTGGATGCGCCGTATTCTGACGGGCGGATGTTCTTTCGCAGTGATTTCGTCAATATGAACGTCGGCAGTTTCTCCACTAATGCCGATGGCAAATGGGATGACAACTGGGGCACCTGTACATTACAGGACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGatat                                                                                                                                                                                                                                        <  2:615412/212‑1 (MQ=255)

                                    cAGGTGGATGCGCCGTATTCTGACGGGCGGATGTTCTTTCGCAGTGATTTCGTCAATATGAACGTCGGCAGTTTCTCCACTAATGCCGATGGCAAATGGGATGACAACTGGGGCACCTGTACATTACAGGACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGATATCGGTACC‑ACGCCGATGGGCTTCAACgt                                                                                                                                                                                                            >  1:2078701/1‑239 (MQ=255)

                                    cAGGTGGATGCGCCGTATTCTGACGGGCGGATGTTCTTTCGCAGTGATTTCGTCAATATGAACGTCGGCAGTTTCTCCACTAATGCCGATGGCAAATGGGATGACAACTGGGGCACCTGTACATTACAGGACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGATATCGGTACC‑ACGCCGATGGGCTTCAACgt                                                                                                                                                                                                            >  1:1547298/1‑239 (MQ=255)

                                    cAGGTGGATGCGCCGTATTCTGACGGGCGGATGTTCTTTCGCAGTGATTTCGTCAATATGAACGTCGGCAGTTTCTCCACTAATGCCGATGGCAAATGGGATGACAACTGGGGCACCTGTACATTACAGGACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGATATCGGTACC‑ACGCCGATGGGCTTCAACgt                                                                                                                                                                                                            >  1:1352148/1‑239 (MQ=255)

                                      ggTGGATGCGCCGTATTCTGACGGGCGGATGTTCTTTCGCAGTGATTTCGTCAATATGAACGTCGGCAGTTTCTCCACTAATGCCGATGGCAAATGGGATGACAACTGGGGCACCTGTACATTACAGGACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAAAGAGCACCGGGGCTGGGGTAGCGGAACC‑ACGACGATAGGGTTCAAccggg                                                                                                                                                                                                          >  2:2131879/1‑235 (MQ=255)

                                        tGGATGCGCCGTATTCTGACGGGCGGATGTTCTTTCGCAGTGATTTCGTCAATATGAACGTCGGCAGTTTCTCCACTAATGCCGATGGCAAATGGGATGACAACTGGGGCACCTGTACATTACAGGACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGATATCGGTACC‑ACGCCGATGGGCTTCAACgtggt                                                                                                                                                                                                         <  1:594518/238‑1 (MQ=255)

                                           aTGCGCCGTATTCTGACGGGCGGATGTTCTTTCGCAGTGATTTCGTCAATATGAACGTCGGCAGTTTCTCCACTAATGCCGATGGCAAATGGGATGACAACTGGGGCACCTGTACATTACAGGACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGGCTGGAGCTGGGGCATCGGGAGC‑ACGCCGATGGGCGTCACAGgggt‑gggg                                                                                                                                                                                                    >  2:744981/1‑237 (MQ=255)

                                              cgcCGTATTCTGACGGGCGGATGTTCTTTCGCAGTGATTTCGTCAATATGAACGTCGGCAGTTTCTCCACTAATGCCGATGGCAAATGGGATGACAACTGGGGCACCTGTACATTACAGGACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCCGGCGGAAAGAACGCCGGGGCTGGGGTACCGGTTCC‑ACGCCGGCGGGCGTCCACACGGG‑GGGtgtg                                                                                                                                                                                                 >  2:1539603/1‑239 (MQ=255)

                                                 cGTATTCTGACGGGCGGATGTTCTTTCGCAGTGATTTCGTCAATATGAACGTCGGCAGTTTCTCCACTAATGCCGATGGCAAATGGGATGACAACTGGGGCACCTGTACATTACAGGACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGGCTGGAGCTGGGATATCGGTACC‑ACGCCGATGGGCTTCAACGTGGT‑GGATGTGGTc                                                                                                                                                                                              >  1:2058852/1‑239 (MQ=255)

                                                    aTTCTGACGGGCGGATGTTCTTTCGCAGTGATTTCGTCAATATGAACGTCGGCAGTTTCTCCACTAATGCCAATGGCAAATGGGATGACAACTGGGGCACCTGTACATTACAGGACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGATATCGGTACC‑ACGCCGATGGGCTTCAACGTGGT‑GGATGTGGTcgg                                                                                                                                                                                            <  2:383501/238‑1 (MQ=255)

                                                     ttCTGACGGGCGGATGTTCTTTCGCAGTGATTTCGTCAATATGAACGTCGGCAGTTTCTCCACTAATGCCGATGGCAAATGGGATGACAACTGGGGCACCTGTACATTACAGGACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCGGGCGGAATGACGGTTGGGGGCGGGGTATCGGGACC‑CCGCCCACGGGCTTCAACGTGGG‑GGACGTGGTcgggg                                                                                                                                                                                          >  2:2126465/1‑237 (MQ=255)

                                                       cTGACGGGCGGATGTTCTTTCGCAGTGATTTCGTCAATATGAACGTCGGCAGTTTCTCCACTAATGCCGATGGCAAATGGGATGACAACTGGGGCACCTGTACATTACAGGACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGATATCGGTACC‑ACGCCGATGGGCTTCAACGTGGT‑GGATGTGGTCggcggc                                                                                                                                                                                        <  2:374732/239‑1 (MQ=255)

                                                        tgccGGGCGTGTGTTCTTTGGCGGTGATGTCGTCAATATGAACGACGGAGGTTTCTCCACTAAGGCCGAGGGCAAATGGGCTGACAACGGGAGCCCCTGTTCATTACAGGACTGTAGCGGCAACCGCAGCCAGCCGGATTCCGGTGCCAGCGTGGGGGTGGGGTGGGGAAATGACGTCTGGAGCTGGGATATGGGTACC‑ACGCCGATGGGCTTCAACGTGGG‑GGATGTGGTCGGCGGCa                                                                                                                                                                                       <  2:2066831/236‑1 (MQ=255)

                                                         gACGGGCGGATGTTCTTTCGCAGTGATTTCGTCAATATGAACGTCGGCAGTTTCTCCACTAATGCCGATGGCAAATGGGATGACAACTGGGGCACCTGTACATTACAGGACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGATATCGGTACC‑ACGCCGATGGGCTTCAACGTGGT‑GGATGTGGTCGGCGGCAt                                                                                                                                                                                      <  2:2085810/239‑1 (MQ=255)

                                                             ggCGGATGTTCTTTCGCAGTGATTTCGTCAATATGAACGTCGGCAGTTTCTCCACTAATGCCGATGGCAAATGGGATGACAACTGGGGCACCTGTACATTACAGGACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGATATCGGTACC‑CCGCCGATGGGCTTCAACGTGGT‑GGATGGGGGCGGGGGCATCAGt                                                                                                                                                                                  >  1:1025663/1‑239 (MQ=255)

                                                             ggCGGATGTTCTTTCGCAGTGATTTCGTCAATATGAACGTCGGCAGTTTCTCCACTAATGCCGATGGCAAATGGGATGACAACTGGGGCACCTGTACATTACAGGACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGATATCGGTACC‑ACGCCGATGGGCTTCAACGTGGG‑GGATGTGGTCGGCGGCATCAGt                                                                                                                                                                                  >  1:1997271/1‑239 (MQ=255)

                                                                     ttCTTTCGCAGTGATTTCGTCAATATGAACGTCGGCAGTTTCTCCACTAATGCCGATGGCAAATGGGATGACAACTGGGGCACCTGTACATTACAGGACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGATATCGGTACC‑ACGCCGATGGGCTTCAACGTGGT‑GGATGTGGTCGGCGGCATCAGTTACAGCg‑a                                                                                                                                                                         <  2:1050813/239‑1 (MQ=255)

                                                                          tCGCAGTGATTTCGTCAATATGAACGTCGGCAGTTTCTCCACTAATGCCGATGGCAAATGGGATGACAACTGGGGCACCTGTACATTACAGGACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGATATCGGTACC‑ACGCCGATGGGCTTCAACGTGGT‑GGATGTGGTCGGCGGCATCAGTTACAGCG‑ATGat‑a                                                                                                                                                                   <  1:83376/239‑1 (MQ=255)

                                                                              aGTGATTTCGTCAATATGAACGTCGGCAGTTTCTCCACTAATGCCGATGGCAAATGGGATGACAACTGGGGCACCTGTACATTACAGGACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGATATCGGTACC‑ACGCCGATGGGCTTCAACGTGGT‑GGATGTGGTCGGCGGCATCAGTTACAGCG‑ATGAT‑ATCg                                                                                                                                                                <  1:768694/238‑1 (MQ=255)

                                                                                  aTTTCGTCAATATGAACGTCGGCAGTTTCTCCACTAATGCCGATGGCAAATGGGATGACAACTGGGGCACCTGTACATTACAGGACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGCAATGACGTCTGGAGGTGGGATATCGGGACC‑ACGCCGCGGGGCTTCAACGGGGT‑GGGGGGGGTCGGCGGGATCAGTTACAGCG‑CTGAT‑ATACGGcc                                                                                                                                                            >  2:1646519/1‑238 (MQ=255)

                                                                                  aTTTCGTCAATATGAACGTCGGCAGTTTCTCCACTAATGCCGATGGCAAATGGGATGACAACTGGGGCACCTGTACATTACAGGACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGATATCGGTACC‑ACGCCGATGGGCTTCaa                                                                                                                                                                                                               <  2:3408/190‑1 (MQ=255)

                                                                                  aTTTCGTCAATATGAACGTCGGCAGTTTCTCCACTAATGCCGATGGCAAATGGGATGACAACTGGGGCACCTGTACATTACAGGACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGATATCGGTACC‑ACGCCGATGGGCTTCaa                                                                                                                                                                                                               >  1:3408/1‑190 (MQ=255)

                                                                                  aTTTCGTCAATATGAACGTCGGCAGTTTCTCCACTAATGCCGATGGCAAATGGGATGACAACTGGGGCACCTGTACATTACAGGACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGATATCGGTACC‑ACGCCGATGGGCTTCAACgtg                                                                                                                                                                                                           <  1:870858/194‑1 (MQ=255)

                                                                                   tttCGTCAATATGAACGTCGGCAGTTTCTCCACTAATGCCGATGGCAAATGGGATGACAACTGGGGCACCTGTACATTACAGGACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGATATCGGTACC‑ACGCCGATGGGCTTCAACGTGGT‑GGATGTGGTCGGCGGCATCAGTTACAGCG‑ATGAT‑ATCGGGCCg                                                                                                                                                           <  2:219451/238‑1 (MQ=255)

                                                                                   tttCGTAAATATGACCGTCGGCAGTTTCTCCCCTAATGCCGATGGCAAATGGGATGACAACTGGGGCACCTGTACATTACAGGACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGATATCGGTACC‑ACGCCGATGGGCTTCAACGTGGT‑GGATGTGATCGGCGGCATCAGTTACAGCG‑ATGAT‑ATCGGGCCg                                                                                                                                                           <  2:407286/238‑1 (MQ=255)

                                                                                       gTCAATATGAACGTCGGCAGTTTCTCCACTAATGCCGATGGCAAATGGGATGACAACTGGGGCACCTGTACATTACAGGACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGATATCGGTACC‑ACGCCGATg                                                                                                                                                                                                                       >  1:373152/1‑177 (MQ=255)

                                                                                       gTCAATATGAACGTCGGCAGTTTCTCCACTAATGCCGATGGCAAATGGGATGACAACTGGGGCACCTGTACATTACAGGACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGATATCGGTACC‑ACGCCGATGGGCTTCAACGTGGT‑GGATGTGGTCggcggc                                                                                                                                                                                        <  2:1427070/207‑1 (MQ=255)

                                                                                       gTCAATATGAACGTCGGCAGTTTCTCCACTAATGCCGATGGCAAATGGGATGACAACTGGGGCACCTGTACATTACAGGACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGATATCGGTACC‑ACGCCGATGGGCTTCAACGTGGT‑GGATGTGGTCggcggc                                                                                                                                                                                        >  1:1427070/1‑207 (MQ=255)

                                                                                             tgggaCGTGGGCAGTTTCTCCACTAATGCCGATGGCAAATGGGATGACAACTGGGGCACCTGTACATTACAGGACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGGGCTGGGATATCGGTACC‑ACGCCGATGGGCTTCAACGTGGT‑GGATGTGGTCGGCGGCATCAGTTACAGCG‑ATGAT‑ATCGGGCCGCTGGGTTacac                                                                                                                                                <  2:1692571/235‑1 (MQ=255)

                                                                                             atGAACGTCGGCAGTTTCTCCACTAATGCCGATGGCAAATGGGATGACAACTGGGGCACCTGTACATTACAGGACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGCCGTCTGGAGCTGGGATATCGGTACC‑ACGCCGATGGGCTTCAACGTGGT‑GGATGTGGTCGGCGGCATCAGTTACAGCG‑ATGAT‑ATCGGGCCGCTGGGTTacac                                                                                                                                                <  2:881816/239‑1 (MQ=255)

                                                                                                 aCGTCGGCAGTTTCTCCACTAATGCCGATGGCAAATGGGATGACAACTGGGGCACCTGTACATTACAGGACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCCGGCGAAATGACGTCCGGAGCTGGGGGTTCGGGACC‑ACGCCGGCGGGCTTTCACAGGGG‑GGGTGCGGGCGGCGGCCTCAATTACCGGG‑ATGGT‑TTCGGGCGCCCGGGGTACACCGtt                                                                                                                                            >  2:1288336/1‑239 (MQ=255)

                                                                                                      ggCAGTTTCTCCACTAATGCCGATGGCAAATGGGATGACAACTGGGGCACCTGTACATTACAGGACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGCTATCGGTGCC‑ACGGCGATGGGCTTCAACGTGGT‑GGCAGTGGGCTGGGGCAGCCGTTCGCGCG‑GTGAT‑ATCGGGGTGCGGGGGCACACCGTTAACg                                                                                                                                        >  2:1235974/1‑238 (MQ=255)

                                                                                                      ggCAGTTTCTCCACTAATGCCGATGGCAAATGGGATGACAACTGGGGCACCTGTACATTACAGGACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGATATCGGTACC‑ACGCCGATGGGCTTCAACGTGGT‑GGATGTGGTCGGCGGCATCAGTTACAGCG‑ATGAT‑ATCGGGCCGCTGGGTTACACCGTTAACg                                                                                                                                        >  1:503949/1‑238 (MQ=255)

                                                                                                      ggCAGTTTCTCCACTAATGCCGATGGCAAATGGGATGACAACTGGGGCACCTGTACATTACAGGACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGATATCGGTACC‑ACGCCGATGGGCTTCAACGTGGT‑GGATGTGGTCGGCGGCATCAGTTACAGCG‑ATGAT‑ATCGGGCCGCTGGGTTACACCGTTAACg                                                                                                                                        >  1:1336085/1‑238 (MQ=255)

                                                                                                                 cACTAATGCCGATGGCAAATGGGATGACAACTGGGCCGGCTGTACATTACAGGACTGTAGCGGCAACCGCGGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGATATCGGTACC‑ACGCCGATGGGCTTCAACGTGGT‑GGATGTGGTCGGCGGCATCAGTTACAGCG‑ATGAT‑ATCGGGCCGCTGGGTTACACCGTTAACGCCCACCGTCGGc                                                                                                                            <  1:608210/239‑1 (MQ=255)

                                                                                                                   cTAATGCCGATGGCAAATGGGATGACAACTGGGGCACCTGTACATTACAGGACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGATATCGGTACC‑ACGCCGATGGGCTTCAACGTGGT‑GGATGTGGTCGGCGGCATCAGTTACAGCG‑ATGAT‑ATCGGGCCGCTGGGTTACACCGTTAACGCCCACCGTCGGccc                                                                                                                          <  1:210672/239‑1 (MQ=255)

                                                                                                                   cTAAGGCCGATGGTAAATGGGATGACAACTGGGGCACCTGGACATTGAAGGACTGTAGCGGCAGCCGGAGCCAGTCGGATTCGGGTGCCGGGGTGGCGGTCGGCTGGCGAAATGGCGTCTGGAGCTGGGATATCGGTACC‑ACGCCGATGGGCTTCAACGTGGT‑GGATGTGGTCGGCGGCATCAGTTACAGCG‑ATGAT‑ATCGGGCCGCTGGGTTACACCGTTAACGCCCACCGTCGGccc                                                                                                                          <  2:1000064/239‑1 (MQ=255)

                                                                                                                     aaTGCCGATGGCAAATGGGATGACAACTGGGGCACCTGTACATTACAGGACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGATATCGGTACC‑ACGCCGATGGGCTTCAACGTGGT‑GGATGTGGTCGGCGGCATCAGTTACAGCG‑ATGAT‑ATCGGGCCg                                                                                                                                                           >  1:966349/1‑204 (MQ=255)

                                                                                                                     aaTGCCGATGGAAAATGGGATGACAACTGGGGCACCTGTACATTACAGGACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGATATCGGTACC‑ACGCCGATGGGCTTCAACGTGGT‑GGATGTGGTCGGCGGCATCAGTTACAGCG‑ATGAT‑ATCGGGCCg                                                                                                                                                           <  2:966349/204‑1 (MQ=255)

                                                                                                                       tGCCGATGGCAAATGGGATGACAACTGGGGCACCTGTACATTACAGGACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGGAATGACGTCCGGAGCTGGGGTATCGGTACC‑ACGCCGGAGGGGCTGAAAGGGGT‑GGGTGTGGGCGGCGGGCTCCGTTTAAGCG‑ATGGT‑ACCGGGCCGCGGGGGTACACCGTTAAAGCCCACCGTCGGccactc                                                                                                                       >  2:1105221/1‑234 (MQ=255)

                                                                                                                           gATGGCAAATGGGATGACAACTGGGGCACCTGTACATTACAGGACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGATATCGGTACC‑ACGCCGATGGGCTTCAACGTGGT‑GGATGTGGTCGGCGGCATCAGTTACAGCG‑ATGAT‑ATCGGGCCGCTGGGTTACACCGTTAACGCCCACCGTCGGccc                                                                                                                          <  2:60178/231‑1 (MQ=255)

                                                                                                                           gATGGCAAATGGGATGACAACTGGGGCACCTGTACATTACAGGACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGATATCGGTACC‑ACGCCGATGGGATTCAACGTGGT‑GGATGTGGTCGGCGGCCTCAGATACAGCG‑ATGAG‑ATCGGGCCGCTGGGTTACACCGGTAACGCCCACCGTCGGCCCCTCTCCAg                                                                                                                  >  1:1000064/1‑239 (MQ=255)

                                                                                                                            aTGGCAAATGGGATGACAACTGGGGCACCTGTACATTACAGGACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGATATCGGTACC‑ACGCCGGTGGGCTTCCAAGCGGT‑GGGTGTGGTAGGCGGGATCGGTTACAGCG‑ATGGA‑ATCGGGCCGCTGGGTTACACCGTGAACGCGCACCGTCTGCCCATCCCCact                                                                                                                 >  2:608695/1‑237 (MQ=255)

                                                                                                                            aTGGCAAATGGGATGACAACTGGGGCACCTGTACATTACAGGACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGATATCGGTACC‑ACGCCGATGGGCTTCAACGTGGG‑GGGTGTGGTCGGCGGCATCAGTTACAGCG‑ATGAT‑ATCGGGCCGCGGGGTTACACCGTTAACGCCCACCGTCGGCCCATCTCCAG‑t                                                                                                                >  1:407163/1‑239 (MQ=255)

                                                                                                                             tGGCAAATGGGATGACAACTGGGGCACCTGTACATTACAGGACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGATATCGGTACC‑ACGCCGATGGGCTTCAACGTGGT‑GGATGTGGTCGGCGGCATCAGTTACAGCG‑ATGAT‑ATCGGGCCGCTGGGTTACACCGTTAACGCCCACCGTCGGCCCATCTCCAG‑tt                                                                                                               <  1:1832860/239‑1 (MQ=255)

                                                                                                                                       gATGACAACTGGGGCACCTGTACATTACAGGACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGATATCGGTACC‑ACGCCGATGGGCTTCaa                                                                                                                                                                                                               >  1:1458563/1‑137 (MQ=255)

                                                                                                                                        aTGACAACTGGGGCACCTGTACATTACAGGACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGATATCGGTACC‑ACGCCGATGGGCTTCAACGTGGT‑GGATGTGGTCGGCGGCATCAGTTACAGCG‑ATGAT‑ATCGGGCCGCTGGGTTACACCGTTAACGCCCACCGTCGGCCCAtctc                                                                                                                     <  1:1566400/223‑1 (MQ=255)

                                                                                                                                             aaCTGGGGCACCTGTACATTACAGGACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTACGGAGCTGGGCTATCGGTACC‑ACGGCGATGGGGTTCAACGGGGG‑GGGTGTGGGCGGTGGCATCAGTTACAGCG‑ATGAT‑ATCGGGGCGCCGGGGTACACCGTTAACGCCCACTGGCGGCCCATCTCCAG‑TTCGTTTGTGGCCTTtgg                                                                                               >  2:565946/1‑239 (MQ=255)

                                                                                                                                               cTGGGGCACCTGTACATTACAGGACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGATATCGGTACC‑ACGCCGATGGGCTTCAACGTGGT‑GGATGTGGTCGGCGGCATCAGTTACAGCG‑ATGAT‑ATCGGGCCGCTGGGTTACACCGTTAACGCCCACCGTCGGCCCATCTCCAG‑TTCTTTGCTGGCCTTtggtg                                                                                             <  1:1691168/239‑1 (MQ=255)

                                                                                                                                                   ggCACCTGTACATTACAGGACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCCGGCGGAATGACGGCTGGAGCCGGGGTATCGGGACC‑AAGCCCGTGGGGCTCAACGTGGG‑GGATGTGGGCGGCGGCCTCCGGTACCGGG‑GTGGT‑AACCGGCCGCCGGGTTACACCGCTTACCCCCACCCTCCGCCCCTCTCCAG‑TTCCTTGCTGGCCTTTTGggggtc                                                                                         >  2:2057162/1‑237 (MQ=255)

                                                                                                                                                          gTACATTACAGGACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGATATCGGTACC‑ACGCCGATGGGCTTCAACGTGGT‑GGATGTGGTCGGCGGCATCAGTTACAGCG‑ATGAT‑ATCGGGCCGCTGGGTTACACCGTTAACGCCCACCGTCGGCCCATCTCCAG‑TTCTTTGCCGGCCTTTGGTGGGC‑AAAAAGa                                                                                  >  1:1398900/1‑238 (MQ=255)

                                                                                                                                                           tACATTACAGGACTGTAGCGGCAACCGCGGCCAGTAGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGATATCGGTACC‑ACGCCGATGGGCTTCAACGTGGT‑GGATGTGGTCGGCGGCATCAGTTACAGCG‑ATGAT‑ATCGGGCCGCTGGGTTACACCGTTAACGCCCACCGTCGGCCCATCTCCAG‑TTCTTTGCt                                                                                                        <  1:1178775/216‑1 (MQ=255)

                                                                                                                                                           tACATTACAGGACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGATATCGGTACC‑ACGCCGATGGGCTTCAACGTGGT‑GGATGTGGTCGGCGGCATCAGTTACAGCG‑ATGAT‑ATCGGGCCGCTGGGTTACACCGTTAACGCCCACCGTCGGCCCATCTCCAG‑TTCTTTGCTGGCCTTtggt                                                                                              <  1:1489026/226‑1 (MQ=255)

                                                                                                                                                           tACATTACAGGACGGGAGGGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGGGTGGCGGTCGGCTGGGGAAATGACGTCTGGAGCTGGGATATCGGTACC‑ACGCCGATGGGTTTCAACGTGGT‑GGATGTGGTCGGCGGCATCAGTTACAGCg‑a                                                                                                                                                                         <  2:826614/153‑1 (MQ=255)

                                                                                                                                                             cATTACAGGAGTGTGGCGGCAACCGCAGCCGGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGATATCGGTACC‑ACGCCGATGGGCTTCAACGTGGT‑GGATGTGGTCGGCGGCATCAGTTACAGCG‑ATGAT‑ATCGGGCCGCTGGGTTACACCGTTAACGCCCACCGTCGGCCCATCTCCAG‑TTCTTTGCTGGCCTTTGGTGGGC‑AAAAAGAC‑Tcc                                                                             <  2:807500/239‑1 (MQ=255)

                                                                                                                                                                tACAGGACTGTAGCGGCAACCGCATCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCAGCCGGCGAAATGACGTCTGGGGCTGGGACATAGGTACC‑CCCCCCATGGGCTTCAACGTGGG‑GGCTTGGGGCGGCGGCGATAGGTACAGCG‑AAGATGATGGGCCCGC‑GGGTTACACCGGTTACGCCCACCTTCGCTCCAACTCCCA‑TTCTTTTCATGTCTTTAGAGGGG‑AACAAAAC‑TACCCg                                                                          >  2:1834769/1‑239 (MQ=255)

                                                                                                                                                                tACAGGACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGATATCGGTACC‑ACGCCGATGGGCTTCAACGTGGT‑GGATGTGGTCGGCGGCATCAGTTACAGCG‑ATGAT‑ATCGGGCCg                                                                                                                                                           <  1:811185/161‑1 (MQ=255)

                                                                                                                                                                tACAGGACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGATATCGGTACC‑ACGCCGATGGGCTTCAACGTGGT‑GGATGTGGTCGGCGGCATCAGTTACAGCG‑ATGAT‑ATCGGGCCGCTGGGTTacac                                                                                                                                                <  1:1247780/172‑1 (MQ=255)

                                                                                                                                                                  cAGGACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGTCAGCGTGGCGGTCGGCCGGCGAACAGGAGGCTGGGACTGGGATCTCGGGACCAACGCGGAGGGGC‑TCAACGCGGTGGGATCGGGAGGGCGCACACGGTTCAAGCG‑‑GGAT‑ATCGGGCCGATGGGGTTAACCGTTAACGCCCACCGCGGGCCCGAGCCCCG‑GTCTGGTCTGGCCTCTTGTGGGG‑AAAAAGAC‑CCCCCgc                                                                         >  2:1185224/1‑237 (MQ=255)

                                                                                                                                                                  cAGGACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGATATCGGTACC‑ACGCCGATGGGCTTCAACGTGGT‑GGATGTGGTCGGCGGCATCAGTTACAGCG‑ATGAT‑ATCg                                                                                                                                                                >  1:1377722/1‑154 (MQ=255)

                                                                                                                                                                  cAGGACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGCCTGGGGCTGGGATATCGGTACC‑ACGCCGATGGGCTTCCACGGGGG‑GGAAGTGGTCGGCGGCCTCAGTTACAGCG‑GTGAT‑AACGGGCCGCTGGGTTACACCGTTAACGCCCACCGTCGGCCCAACTCCAG‑GTCTTTGCTGGCCTTTGGGGGGC‑AAAAAGAA‑TCCCCgc                                                                         >  2:261976/1‑237 (MQ=255)

                                                                                                                                                                     gACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATTACGTCCGGAGCTGGGGTATCGGTACC‑ACGCCGATGGGCCTCCACGTGGG‑GGGTGTGGTCGGCGGCCTCAGTTACAGCG‑ACGAT‑AACGGGCCGCTGGGTTACACCGTTAACCCCCACCCGCGGCCCATCTCCCGTTTCTTGGCGGGCC‑TTGGGGGGG‑ACAAAGAC‑CTCCCGGGCaa                                                                     >  2:503413/1‑239 (MQ=255)

                                                                                                                                                                     gACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGATATCGGTACC‑ACGCCGATGGGCTTCAACGTGGT‑GGATGTGGTCGGCGGCATCAGTTACAGCG‑ATGAT‑ATCGGGCCGCTGGGTTACACCGTTAACGCCCACCGTCGGCCCATCTCCAG‑TTCTTTGCTGGCCTTTGGTGGGC‑AAAAAGAC‑TCCCCGAGCaa                                                                     >  1:1185539/1‑239 (MQ=255)

                                                                                                                                                                     gACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGATATCGGTACC‑ACGCCGATGGGCTTCAACGTGGT‑GGATGTGGTCGGCGGCATCAGGTACAGCG‑ATGAT‑ATCGGGCCGCTGGGTTACACCGTTAACGCCCACCGTCGGCCCATCTCCAG‑TTCTTTGCTGGCCTTTGGTGGGC‑AAAAAGAC‑TCCCCGAGcca                                                                     >  1:1224841/1‑237 (MQ=255)

                                                                                                                                                                     gACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGATATCGGTACC‑ACGCCGATGGGCTTCAACGTGGT‑GGATGTGGGCGGCGGCATCAGTTACAGCG‑ATGAT‑ATCggg                                                                                                                                                              >  1:826614/1‑153 (MQ=255)

                                                                                                                                                                     gACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCCGGGATATCGGTACC‑ACGCCGATGGGCTTCAACGTGGT‑GGATGTGGTCGGCGGCATCCGGTACAGCG‑ATGAT‑ATCGGGCCGCTGGGGTACACCGTTAACGCCCACCCGCGGCCCATCTCCAG‑TTCTTTTCTGGCCTTTGGGGGGC‑AAAAAGAC‑tgtcccttaaa                                                                     >  2:1489026/1‑229 (MQ=255)

                                                                                                                                                                      aCTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGATATCGGTACC‑ACGCCGATGGGCTTCAACGTGGT‑GGATGTGGTCGGCGGCATCAGTTACAGCG‑ATGAT‑ATCGGGCCGCTGGGTTACACCGTTAACGCCCACCGTCGGCCCATCTCCAG‑TTCTTTGCTGGCCTTTGGTGGGC‑AAAAAGAC‑TCCCCGAGCaa                                                                     <  1:2126465/238‑1 (MQ=255)

                                                                                                                                                                        tGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGATATCGGTACC‑ACGCCGATGGGCTTCAACGTGGT‑GGATGTGGTCGGCGGCATCAGTTACAGCG‑ATGAT‑ATCGGGCCGCTGGGTTACACCGTTAACGCCCACCGTCGGCCCATCTCCAG‑TTCTTTGCTGGCCTTTGGTGGGC‑AAAAAGAC‑TCCCCGCGCAATa                                                                   >  1:738916/1‑238 (MQ=255)

                                                                                                                                                                          tAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGGGCCGGGGTATCGGGACC‑AAGCCGATGGGCCTTAAAGTGGG‑GGATGGGGGCGGCCGCACCAGTTACCGCG‑GTGGT‑AACGGGCCGCTGGGGTGCCCCGTTAACGCCCACCGTCGGCCCCACCCCAG‑TTCTTTGCCGGCCCTTGGGTGGC‑CAAAAGAC‑CCCCCCAACCATAccc                                                                >  2:1742695/1‑238 (MQ=255)

                                                                                                                                                                          tAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGATATCGGTACC‑ACGCCGATGGGCTTCAACGTGGT‑GGATGTGGTCGGCGGCATCAGTTACAGCG‑ATGAT‑ATCGGGc                                                                                                                                                             <  2:916573/149‑1 (MQ=255)

                                                                                                                                                                          tAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGATATCGGTACC‑ACGCCGATGGGCTTCAACGTGGT‑GGATGTGGTCGGCGGCATCAGTTACAGCG‑ATGAT‑ATCGGGCCGCTGCGTTACACCGTTAACGCCCACCGTCGGCCCATCTCCAG‑TTCTTTGCTGGCCTTTGGTGGGC‑AAAAAGAG‑TCCCCGAGCAATACCg                                                                >  1:650868/1‑239 (MQ=255)

                                                                                                                                                                          tAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGATATCGGTACC‑ACGCCGATGGGCCTCAACGTGGG‑GGATGTGGGCGGCGGCATCAGTTACAGCG‑ATGAT‑ATCGGGCCGCTGGGTTACACCGTGAACGCCCACCGCCGGGCCCTTTCCAG‑TTTTTTGCTGGCCCTTGGTGGGC‑AAAAAGAC‑ACCCCGCGCAAtccag                                                                >  2:696162/1‑235 (MQ=255)

                                                                                                                                                                          tAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGATATCGGTACC‑ACCCCGATGGGCTTCAACGTGGG‑GGATGTGGTCGGCGGCATCAGTTACAGCG‑GTGAT‑ATCCGGCCGCCTGGTTACACCGATAACGTCCACCGGCGGCCCATCTCCAG‑TTCTTTTCTGGCCTTTGGGGGGC‑AAAAAGAC‑TCCCCGAGCAAAACCg                                                                >  2:196101/1‑239 (MQ=255)

                                                                                                                                                                          aaGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGGACTGTCCGCTGGGTCTGCGGAATCGGTACC‑ACGCCGAAGCGCTTCAACGCGGC‑GGAAGGGGACGGCGGACACCGCTACAGCA‑GTGAT‑GACAGGTCGCAGGGTTACACCGTTAACGCACACCGCCGTCCCAACCCCCG‑TTTTTTGCTGGCCTTTTGTGGGC‑CGAAAGGCTTCCCCAAGCacaaca                                                                 >  2:148759/2‑234 (MQ=255)

                                                                                                                                                                                  aCCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGGAATGACATCTTGGGCTGGGAAAACGGTACC‑ACGCCGTATGGCTTCAACGTGGG‑GGATGTGGGCGGCTGCATCTGTTACAGCGTAGGAT‑TCCGGGCCG‑TGGGGTACAACGTTAATGCCCACTCTCTGCCCCTCTCCAG‑ATATTTGATTGCCTTTGGGTGGC‑AAAAAAAC‑TCCCCGAGCAAAACCAGGtaaaaa                                                        >  2:1273870/1‑239 (MQ=255)

                                                                                                                                                                                    cGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTTGGCTGGCGAAAAGACGTCTGGAGCTGGGGCATCGGGGAC‑ACGCCGAGGGGCTTCAAAGACGT‑GGAGGGGGGCGGCGGCATCAGTTACAGCG‑AGGAA‑ATCGGGCCGCTGGGGTAAACCGTGAAAGCCCACCGTCGGCCCATCCCCAG‑GTCTTTGCTGGCCATTGTCGGGCAAAAAAGGC‑CCCCCCA‑CAATACCggggcaatatg                                                      >  2:449925/1‑231 (MQ=255)

                                                                                                                                                                                       aGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGATATCGGTACC‑ACGCCGATGGGCTTCAACGTGGT‑GGATGTGGTCGGCGGCATCAGTTACAGCG‑ATGAT‑ATCGGGCCGCTGGGTTACACCGTTAACGCCCACCGTCGGCCCATCTCCAG‑TTCTTTGCTGGCCTTTGGTGGGC‑AAAAAGAC‑TCCCCGAGCAATACCGGGAAAAAATggg                                                    <  1:696162/238‑1 (MQ=255)

                                                                                                                                                                                          cAGTCGGATTCCGGTGCCAGGGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGATATCGGTACC‑ACGCCGATGGGCTTCAACGTGGT‑GGATGTGGTCGGCGGCATCAGTTACAGCG‑ATGAT‑ATCGGGCCGCTGGGTTACACCGTTAACGCCCACCGTCGGCCCATCTCCAG‑TTCTTTGCTGGCCTTTGGTGGGC‑AAAAAGAC‑TCCCCGAGCAATACCGGGAAAAAATgggg                                                   <  1:2179956/236‑1 (MQ=255)

                                                                                                                                                                                             tCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGCCTGGAGCTGGGGTAACGGTACC‑CCGCCGATGGGGTTCAACGTGGG‑GGATGTGGGCGGCGGCATCAGTTACAGCG‑ATGAT‑ATCGGGCCGGTGGGTTACACCGGTAACGCCCACCGTCGGCCCATCTCCAG‑TTTTTTGCTGGCCTTTGGTGGGC‑AAAAAGAC‑TCCCCGTGCAATACCGGGGAAAAAAGGGGGGGGGt                                             >  2:925725/1‑239 (MQ=255)

                                                                                                                                                                                              cGGGTTCCGGTGCCAGGGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGATATCGGTACC‑ACGCCGATGGGCTTCAACGTGGT‑GGATGTGGTCGGCGGCATCAGTTACAGCG‑ATGAT‑ATCGGGCCGCTGGGTTACACCGTTAACGCCCACCGTCGGCCCATCTCCAG‑TTCTTTGCTGGCCTTTGGTGGGC‑GAAAAGAC‑TCCCCGAGCAATACCGGGAAAAAATGGGGTGGCg                                              <  1:924530/237‑1 (MQ=255)

                                                                                                                                                                                                   tCCGGGGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGATATCGGTACC‑ACGCCGATGGGCTTCAACGTGGT‑GGATGTGGTCGGCGGCATCAGTTACAGCG‑ATGAT‑ATCGGGCCGCTGGGTTACACCGTTAACGCCCACCGTCGGCCCATCTCCAG‑TTCTTTGCTGGCCTTTGGTGGGC‑AAAAAGAC‑TCCCCGAGCAATACCGGGAAAAAATGGGGTGGCGTACGTGc                                       <  1:1105221/239‑1 (MQ=255)

                                                                                                                                                                                                      ggTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGATATCGGTACC‑ACGCCGATGGGCTTCAACGTGGT‑GGATGTGGTCGGCGGCATCAGTTACAGCG‑ATGAT‑ATCGGGCCGCTGGGTTACACCGTTAACGCCCACCGTCGGCCCATCTCCAG‑TTCTTTGCTGGCCTTTGGTGGGC‑AAAAAGAC‑TCCCCGAGCAATACCGGGAAAAAATGGGGTGGCGTACGTGCCGa                                    <  1:1915486/239‑1 (MQ=255)

                                                                                                                                                                                                      ggTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGATATCGGTACC‑ACGCCGATGGGCTTCAACGTGGT‑GGATGTGGTCGGCGGCATCAGTTACAGCG‑ATGAT‑ATCGGGCCGCTGGGTTACACCGTTAACGCCCACCGTCGGCCCATCTCCAG‑TTCTTTGCTGGCCTTTGGTGGGC‑AAAAAGAC‑TCCCCGAGCAATACCGGGAAAAAATGGGGTGGCGTACGTGCCGa                                    <  1:284561/239‑1 (MQ=255)

                                                                                                                                                                                                        tGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGATATCGGTACC‑ACGCCGATGGGCTTCAACGTGGT‑GGATGTGGTCGGCGGCATCAGTTACAGCG‑ATGAT‑ATCGGGCCGCTGGGTTACACCGTTAACGCCCACCGTCGGCCTATCTCCAG‑TTCTTTGCTGGCCTTTGGTGGGC‑AAAAAGAC‑TCCCCGAGCAATACCGGGAAAAAATGGGGTGGCGTACGTGCCGAc                                   <  1:196101/238‑1 (MQ=255)

                                                                                                                                                                                                             gCGTGGCGGTCGGCTGGCGCAATGACGTCTGGAGCTGGGATATCGGTACC‑ACGCCGATGGGCTTCAACGTGGT‑GGATGTGGTCGGCGGCATCAGTTACAGCG‑ATGAT‑ATCGGGCCGCTGGGTTACACCGTTAACGCCCACCGTCGGCCCATCTCCAG‑TTCTTTGCTGGCCTTTGGTGGGC‑AAAAAGAC‑TCCCCGAGCAATACCGGGAAAAAATGGGGTGGCGTACGTGCCGACGgtgt                              <  2:1555235/238‑1 (MQ=255)

                                                                                                                                                                                                               gTGGCGGTCGGCTGGGGAAATGAGGTCTGGAGCTGGGATATCGGTCCC‑ACGCCGATGGGCTTCAACGTGGT‑GGATGTGGTCGGCGGCATCAGTTACAGCG‑ATGAT‑ATCGGGCCGCTGGGTTACACCGTTAACGCCCACCGTCGGCCCATCTCCAG‑TTCTTTGCTGGCCTTTGGTGGGC‑AAAAAGAC‑TCCCCGAGCAATACCGGGAAAAAATGGGGTGGCGTACGTGCCGACGgtgt                              <  1:884540/236‑1 (MQ=255)

                                                                                                                                                                                                                       cggcggGCGAAATGACGTCTGGAGGTGGGATATCGGTCCC‑ACGCCGATGGGCTTCAACGTGGT‑GGATGTGGTCGGCGGCATCAGTTACAGCG‑ATGAT‑ATCGGGCCGCTGGGTTACACCGTTAACGCCCACCGTCGGCCCATCTCCAG‑TTCTTTGCTGGCCTTTGGTGGGC‑AAAAAGAC‑TCCCCGAGCAATACCGGGAAAAAATGGGGTGGCGTACGTGCCGACGGTGTGGGGCTAAGt                    <  2:1352148/238‑1 (MQ=255)

                                                                                                                                                                                                                            ggCGAAATGACGTCTGGAGCTGGGATATCGGTACC‑ACGCCGATGGGCTTCAACGTGGT‑GGATGTGGTCGGCGGCATCAGTTACAGCG‑ATGAT‑ATCGGGCCGCTGGGTTACACCGTTAACGCCCACCGTCGGCCCATCTCCAG‑TTCTTTGCTGGCCTTTGGTGGGC‑AAAAAGAC‑TCCCCGAGCAATACCGGGAAAAAATGGGGTGGCGTACGTGCCGACGGTGTGGGGCTAAGTCTGAGc              <  1:565946/239‑1 (MQ=255)

                                                                                                                                                                                                                            ggCGAAAGGGGGTGTGGAGCTGGGGTATCGGGACC‑ACGCGGATGGGAGTCATCATGGT‑GGTTGTGGTCGGCGGCCTCACTTACAGCG‑ATGCT‑ATCGCGACGCTGGGTTTGGCCGTTCATGCCCAACGTCGGCCCAACTCCAG‑TTCTTTGCTGGCCGTTGGTGGGC‑AGCAAGAC‑TCCCCGAGCAATGCCGGGAAAAAATGGGGTGGCGTCCGTGCCGACGGTGTGGGGCTAAGTCTGAga              <  2:1474732/239‑2 (MQ=255)

                                                                                                                                                                                                                                   tgccGTCTGGGGCGTGGATATCGGTACC‑ACGCCGATGGGCTTCAACGTGGT‑GGATGTGGTCGGCGGCATCAGTTACAGCG‑ATGAT‑ATCGGGCCGCTGGGTTACACCGTTAACGCCCACCGTCGGCCCATCTCCAG‑TTCTTTGCTGGCCTTTGGTGGGC‑AAAAAGAC‑TCCCCGAGCAATACCGGGAAAAAATGGGGTGGCGTACGTGCCGACGGTGTGGGGCTAAGTCTGAGCTACGATa       <  2:322862/236‑1 (MQ=255)

GGTTACTCCGATCTGAAAGCGCACACTACCATGTTGCAGGTGGATGCGCCGTATTCTGACGGGCGGATGTTCTTTCGCAGTGATTTCGTCAATATGAACGTCGGCAGTTTCTCCACTAATGCCGATGGCAAATGGGATGACAACTGGGGCACCTGTACATTACAGGACTGTAGCGGCAACCGCAGCCAGTCGGATTCCGGTGCCAGCGTGGCGGTCGGCTGGCGAAATGACGTCTGGAGCTGGGATATCGGTACC‑ACGCCGATGGGCTTCAACGTGGT‑GGATGTGGTCGGCGGCATCAGTTACAGCG‑ATGAT‑ATCGGGCCGCTGGGTTACACCGTTAACGCCCACCGTCGGCCCATCTCCAG‑TTCTTTGCTGGCCTTTGGTGGGC‑AAAAAGAC‑TCCCCGAGCAATACCGGGAAAAAATGGGGTGGCGTACGTGCCGACGGTGTGGGGCTAAGTCTGAGCTACGATAAAGGT  >  NC_007779/3953639‑3954109

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