BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_007779	3,874,897	A→C	9.3%	S1112R (AGT→AGG)	Y75_RS19225 ←	RHS element protein RhsA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007779	3,874,897	0	A	C	9.3%	639.2 / 6.0	319	S1112R (AGT→AGG)	Y75_RS19225	RHS element protein RhsA
Reads supporting (aligned to +/- strand): ref base A (74/198); new base C (18/28); total (92/227)
Fisher's exact test for biased strand distribution p-value = 1.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 328 total aligned reads displayed.

CAGGCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCGTATTCTGCGCACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCC  >  NC_007779/3874660‑3875135

         aaGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCGTATTCTGCGCACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGt                                                                                                                                                                                                                                      <  2:2212061/239‑1 (MQ=37)

               ctgcAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCGTATTCTGCGCACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACCGGGGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCGGGTTTTGCATCTGCTCCCCGGTCAGGCCGCACGATTCCAGCCCCAGGCGGCGTCCCTCC‑CTCACCCCGCCACCCa                                                                                                                                                                                                                                >  2:1645754/1‑239 (MQ=11)

                 gcAGCTGATGCGTGTGCTCGTCATTCAGCAGGGTGCCCCATTCATCGTATTCTGCGCACCATTCTGTTGTCCTTTCCTTGCTTATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGTCAGCCAgg                                                                                                                                                                                                                              <  2:1077096/239‑1 (MQ=255)

                 gcAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCGTATTCTGCGCACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGTCAGCCAgg                                                                                                                                                                                                                              <  1:2226431/239‑1 (MQ=37)

                      tGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCGTATTCTGCGCACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGTCAGCCAGGAttt                                                                                                                                                                                                                          <  1:728808/238‑1 (MQ=37)

                            gggTTCTCTTCATTCAGCAGGTTGCCCCATTCATCGTATTCTGCGCACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCCCCCGGCCAGCCAGGATTTCACttt                                                                                                                                                                                                                    >  1:1830170/1‑238 (MQ=37)

                               ttCTCTTCATTCAGCAGGTTGCCCCATTCATCGTATTCTGCGCACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTAC‑CTCACCCCGTCAGCCAGGCGTTCACTTTCCa                                                                                                                                                                                                                 >  1:644152/1‑238 (MQ=255)

                                 ctattccttcaGCATGTTGCCCCATTCATCGTATTCTGCGCGCCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCCCCCGGTCAGCCAGGATTTCACTTTCCAGcc                                                                                                                                                                                                              <  2:1812234/236‑1 (MQ=255)

                                  tctTCATTCAGCAGGTTGCCCCATTCATCGTATTCTGCGCACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACCGGTGGATTTTTCGCGCCGGCGTTTACACCAGGTCAATCTGGTTTTGCCTCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCC‑CCCACCCGGTCAGCCCGGGCTTCACTTTTCAgtcc                                                                                                                                                                                                             >  2:1476431/1‑236 (MQ=11)

                                    ttcattcaGCAGGTTGCCCCATTCATCGTATTCTGCGCACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCCGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGTCAGCCAGGCTTTCACTTTCCAGCCgg                                                                                                                                                                                                            >  1:11034/1‑238 (MQ=37)

                                    ttcattcaGCAGGTTGCCCCATTCATCGTATTCTGCGCACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCCCCCGGTCAGCCAGGGTTTCCCTTTCCAGCCgg                                                                                                                                                                                                            >  1:1062731/1‑238 (MQ=255)

                                        ttcaGCAGGTTGCCCCATTCATCGTATTCTGCGCACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCC‑CTCCCCCGGTCCGCCAGGCTTTCACTTTCCAGCCGGTCgcg                                                                                                                                                                                                       >  1:1325542/1‑237 (MQ=255)

                                         tcaGCAGGTTGCCCCATTCATCGTATTCTGCGCACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGACGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGTCCGCCAGGATTTCACTTTCCAGCCGGTCGAg                                                                                                                                                                                                       >  1:1764283/1‑238 (MQ=37)

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

                                            gcagGTTGCCCCATTCATCGTATTCTGCGCACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGCGCa                                                                                                                                                                                                     >  1:1242624/1‑237 (MQ=37)

                                            gcagGTTGCCCCATTCATCGTATTCTGCGCACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCCTTCCTCA‑CTCCACCCGTCAGCCAGGGTTTCACTTTCCACCCGGTCCGGCAt                                                                                                                                                                                                    >  2:194110/1‑238 (MQ=11)

                                             cagGTGGCCCCATTCATCGTTGTGTGCACACCAGTCGTGTGTCCTTGCCTGCCGGATAAGGGCCAGCGGCAGGTCGAGAGGGTCGCAGTGGGACAGGTGGATTTTTCGGGCCGTCGTGTACACCGGGTCCATCTGGTTGTGCATCTTCTCCGCGGTCAGGCGTCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGTTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCAta                                                                                                                                                                                                   <  2:2069811/238‑2 (MQ=11)

                                                gTTGCCCCATTCATCGTATTCTGCGCACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCCCCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGc                                                                                                                                                                                                >  1:541972/1‑238 (MQ=37)

                                                 ttGCCCCATTCATCGTATTCTGCGCACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGc                                                                                                                                                                                                <  1:1821335/237‑1 (MQ=37)

                                                     cccATTCATCGTATTCTGCGCACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAg                                                                                                                                                                                           <  1:2173077/238‑1 (MQ=37)

                                                      ccATGCATCGTCTTCTGCGCCCCATTCTGTTGTCCCTGCCTTGCTGATAAGGGCCGGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGCTTTTTCGCGCCGGCGTGTACCCCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCCCCCGGTCAGCCAGGAATTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCa                                                                                                                                                                                         <  2:541972/239‑1 (MQ=11)

                                                           cATCGTATTCTGCGCACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACcggc                                                                                                                                                                                    <  1:1375626/239‑1 (MQ=37)

                                                             tctatttCTGCGCCCCATTTTGTTGTCCTTTCCTTGCGGATGAGGGCCTGCGGCAGGCCGCAATGGTGGTAGTGGTTCAGGTGGATTTTTCGCGCCGGCGTTTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCGGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGTCAGCCAGGATTTCCCTTTCCAGCCGGTCGAGCATCTGCACCAGCACcggcgg                                                                                                                                                                                  <  2:1717084/234‑1 (MQ=11)

                                                             tCGTTTTGTGCGCACCATTCTGTTGTCCCTTCCTTGTGGATAAGGGCCCGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGGGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACcggcgg                                                                                                                                                                                  <  2:1836354/239‑1 (MQ=255)

                                                             tCGTATTCTGCCCCCCATTCTTTTGTCCCTTCCGTGCTGATAAGGGCCAGCGGCGGGCCGCGATGGTCGCAGTGTTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACcggcgg                                                                                                                                                                                  <  2:679289/239‑1 (MQ=255)

                                                             tCGTATTCTGCCCCCCATTCCGTTGTCCCTTCCTTTCTGATAAGGGCCAGCTGCAGGCCGCGATGGTCGCATTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGGTTTGCATCTGCTCCACGGTCAGGCCGCACTATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCTCCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACcggcgg                                                                                                                                                                                  <  1:1984494/239‑1 (MQ=11)

                                                                 tttgtGCGCACCATTCTGTTGGCCCTTCCTGGTTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTCCAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGaa                                                                                                                                                                               <  2:1775017/234‑1 (MQ=255)

                                                                 aTTCTGCGCCCCTTGCTGTTGTCCCTCCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGAATTTTCCGCGCCGGCGTGTACCCCGGGTCCAGCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCCCGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAac                                                                                                                                                                              <  2:1590935/239‑1 (MQ=11)

                                                                 aTTCTGCGCACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGaa                                                                                                                                                                               <  1:1563501/238‑1 (MQ=37)

                                                                 aTTCTGCGCACCATTCTGGTGTCCCTTCCTTGGTGATAGGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAac                                                                                                                                                                              <  2:1544635/239‑1 (MQ=255)

                                                                   tCTGCGCGCCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAacac                                                                                                                                                                            <  1:451428/239‑1 (MQ=37)

                                                                   tCTGCGCACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCctca‑ct                                                                                                                                                                                                                                              <  2:2180663/173‑1 (MQ=39)

                                                                   tCTGCGCACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGGaca                                                                                                                                                                             >  1:386381/1‑238 (MQ=37)

                                                                   tCTGCGCACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGGaca                                                                                                                                                                             >  1:909088/1‑238 (MQ=37)

                                                                   tCTGCGCACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAacac                                                                                                                                                                            <  1:489902/239‑1 (MQ=37)

                                                                       cgcACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGGAaaccac                                                                                                                                                                         >  1:411877/1‑238 (MQ=37)

                                                                       cgcACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAAaacccc                                                                                                                                                                         >  1:1476582/1‑236 (MQ=37)

                                                                       cgcACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCCTTCCTCA‑CTCACCCGGGCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGGAcacccc                                                                                                                                                                         >  1:296840/1‑236 (MQ=255)

                                                                          aCCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCCCACGATTCCAGCCAGCGGCGGCTTTCCCCC‑CTCACCCGGTCAGCCAGGCATTCACCTTCCCGCCGGCCGAGCATCCCCACCCACACCCGCgggtcacacacccc                                                                                                                                                                      >  2:1132632/1‑227 (MQ=11)

                                                                            cATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAAcaccaccccgc                                                                                                                                                                    >  1:1357387/1‑233 (MQ=37)

                                                                             aTTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCa                                                                                                                                                                  <  1:1670742/239‑1 (MQ=37)

                                                                               tCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCC‑CTCACCCGGTCAGCCAGGATTTCACTTTTCAACCGGTCGCGCATCTGCCCCAGCACCCGCGGGAGAACCacaccacccc                                                                                                                                                                 >  2:1028292/1‑232 (MQ=11)

                                                                               tCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGGCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCCCGATGCCAGCCAGCGGCGGCTTTCCCCC‑CCCACCCGGTCAGCCCGGATTTCACCTTCCAGCCCGTCGCGCATCCGCACCCGCACCCGCCGGCACACCccaccacccc                                                                                                                                                                 >  2:1423813/1‑232 (MQ=11)

                                                                                    tgtCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCtt                                                                                                                                                           <  1:76226/239‑1 (MQ=35)

                                                                                    tgtCCCTTCCTGGCTGATAGGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCGGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCtt                                                                                                                                                           <  1:170513/239‑1 (MQ=25)

                                                                                      tCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGGACACCcccccgcccccgccctc                                                                                                                                                          >  1:27411/1‑233 (MQ=255)

                                                                                              ttGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGa                                                                                                                                                 <  1:1631215/239‑1 (MQ=33)

                                                                                                   ggtcaGGGCCAGCGGCAGGCCGAGATGGTGGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGActgct                                                                                                                                            <  2:1320671/235‑1 (MQ=14)

                                                                                                    aTAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGGACACCACACCGCCCCCCTCCTCCCCCCCGGActgctt                                                                                                                                           >  1:714491/1‑238 (MQ=255)

                                                                                                      aaGGCCCAGCGGAAGGCCGCGATGGTCGCGTTGGGACAGGTGGATTTGTCGCGCCGTCGTGTACGCCGGGTCCATCTGGTTTTGCATCTGCTCCCCGGCCAGGCCTCCGGATGCCCACCAGCGGCGGCTTTCCTCA‑CTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCATCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGaa                                                                                                                                         <  2:1981538/239‑1 (MQ=11)

                                                                                                                ggCGGGCCGCGGTGGCTGCAGTGGTACGGGTGGATTTTTCCCGCCGGGGTGTCCACCGGGTCCCTCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGTCAGCCAGGATTTCCCTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGc                                                                                                                               <  2:148142/239‑1 (MQ=11)

                                                                                                                ggCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGc                                                                                                                               <  1:1285102/239‑1 (MQ=31)

                                                                                                                ggCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTGGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCCCCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCg                                                                                                                                <  1:1014886/238‑1 (MQ=14)

                                                                                                                   aGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGCAACCCACACCGCCCCCGGCCTCGCCGCCGGCCTGCTGAAgcgccaccaccac                                                                                                                            >  1:1844953/1‑230 (MQ=11)

                                                                                                                   aGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCa                                                                                                                             <  1:1872903/238‑1 (MQ=31)

                                                                                                                          gATGGTCGCAGTGGTGCAGGTGGATTTTTCCGGCGCGCGGGTCCACCGGGTCCATTTGGTTTTGCTTCTGCTCCACGGCCAGGCCGCCCGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCCCCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTgcg                                                                                                                      <  2:364773/238‑1 (MQ=11)

                                                                                                                               tCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCt                                                                                                                 <  1:1091970/238‑1 (MQ=31)

                                                                                                                                 gCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCg                                                                                                              <  1:660881/239‑1 (MQ=31)

                                                                                                                                                 aTTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGGAAACCACACCGCCCCCCTCTTCCCCCCCCGGCCTCTGGAGCGCCCCCGCCAGGCTCCGGCCCTGCGCTTTccccccgccac                                                                                               >  1:1733571/1‑231 (MQ=11)

                                                                                                                                                     tccgggcCGGAGTGTCCACCGGTTCCATCCGGTTTGGCATCTGCTCCACGGTCGGGCCGCAGGAGGCCACCCGGCGGCGGCTTTCCTCC‑CTCCCCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTg                                                                                          <  2:1179364/234‑1 (MQ=11)

                                                                                                                                                      tCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTgg                                                                                         <  1:569034/239‑1 (MQ=31)

                                                                                                                                                       cgcgCCGGGGTGTACACCGGTGCCATCTGGTTTTGCATCTGCCCCGCGGTCAGGCCGCCCGCTGCCGGCCAGCGGCGGTTTTCCGCA‑CTCCCCGGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACGCCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCGGCGTTTTCGCCAGCTCACCGGTGGc                                                                                        <  2:421652/239‑1 (MQ=11)

                                                                                                                                                       cgcgCCGGCGTGTCCACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGc                                                                                        <  1:617895/239‑1 (MQ=14)

                                                                                                                                                          gcCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGt                                                                                     <  1:194110/239‑1 (MQ=32)

                                                                                                                                                            cgggatggacacGGGGCCCATCTGGTGTGCCATCTCCCCGACGGTAAGGCCCGACGATGCCAGCCAGCGGCGGCTTTCCTCC‑CTCTCCCGGTCAGCCAGAATTTCACTTCCCAGCCGGTCGAGCATCTGCACCAGCCCCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGCCTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGttt                                                                                   <  2:717586/231‑1 (MQ=11)

                                                                                                                                                            cGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGttt                                                                                   <  1:2041135/239‑1 (MQ=32)

                                                                                                                                                                gtgtACCCCGGGCCCATCTGTTTTTGCATCTGCTCCACGGTCAGGCCGCACGAGGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGTCAGCCAGGATTTCCCTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCaa                                                                                <  2:842588/238‑1 (MQ=16)

                                                                                                                                                                gtgtACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCCCGAGGCCAGCCAGCGGCGGCTTTCCTCC‑CTCACCCGGTCAGCCAGGATTTCCCTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCaa                                                                                <  2:1380740/238‑1 (MQ=16)

                                                                                                                                                                        cGGGTCCATCTGGTTTGGCATCGGCTCCACGGTCAGCCCCCGCGATGCGCGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGTCAGCCAGGATTTCCCTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACGCCACACTGCCACCGTCTTCGCCGCCGGGCTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCtgatga                                                                       <  2:396620/239‑1 (MQ=11)

                                                                                                                                                                          ggTCCATCTGGTTTTGCTTCTGCCCCACGTTCAGGCCGCACGCTGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAg                                                                      <  2:1608182/238‑1 (MQ=255)

                                                                                                                                                                            tCCATCTGGTTTTGCATCTGCTCCCCGGTCAGGCCGCACGATTCCAGCCAGCGGCGGCTTTCCTCA‑CTCCCCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCTGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTg                                                                    <  1:1718919/238‑1 (MQ=255)

                                                                                                                                                                                        tGCATCTGCCCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTccc                                                       <  1:524741/239‑1 (MQ=25)

                                                                                                                                                                                         gCATCTTCTCCACGGTCCGGGCGCGCGTTGCCAGCCAGCGGCGGCTTTCCTCA‑CTCGCCCGGTCAGCCAGGATTTCCCTCGCCAGCCGGTCGAGCATCCGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTcccc                                                      <  2:1354446/239‑1 (MQ=11)

                                                                                                                                                                                                   cACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATaaa                                            <  1:862953/239‑1 (MQ=35)

                                                                                                                                                                                                        tCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCGCCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCa                                                                                 <  1:812127/197‑1 (MQ=18)

                                                                                                                                                                                                         ctggCCTCACGATGCCAGCCAGCGGCGGCTTTCCCCA‑CTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACGCCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTc                                       <  1:1193351/236‑1 (MQ=18)

                                                                                                                                                                                                             ccGCACGATGCCAGCCGGCGGCGGCTTTCCTCA‑CTCCCCCGGTCAGCCGGGATTTCCCTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGAt                                  <  2:258178/239‑1 (MQ=16)

                                                                                                                                                                                                               gCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCCGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGc                                <  1:788141/239‑1 (MQ=18)

                                                                                                                                                                                                                  cGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCCCCCGGTCAGCCAGGATTTCCCTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCggg                             <  1:1258667/239‑1 (MQ=18)

                                                                                                                                                                                                                  cGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCggg                             <  1:220152/239‑1 (MQ=33)

                                                                                                                                                                                                                       cctggcaggGGCGGCTTTCCTCC‑CACACCCGTCCAGCCAGGACTTCTCTTTCCATCCGGTCGGGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCTCCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCGCCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCGCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTg                        <  2:2038832/234‑1 (MQ=11)

                                                                                                                                                                                                                        cagccagcGGCGGCTTTCCTCC‑CTCCCCCGGGCAGCCCGGATTTCACTTTCCAGCCCGTCCAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTg                        <  1:2058134/238‑1 (MQ=255)

                                                                                                                                                                                                                        cagccagcGGCGGCTTTCCTCA‑CTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCCCACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTg                        <  1:978389/238‑1 (MQ=18)

                                                                                                                                                                                                                        cagccagcGGCAGCTTTCCTCA‑CCCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTg                        <  1:1500384/238‑1 (MQ=18)

                                                                                                                                                                                                                            cggcggcgGCTTTCCTCA‑CCCCCCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGtt                   <  2:2098146/237‑1 (MQ=18)

                                                                                                                                                                                                                            cagcGGCGGCTTTCCTCC‑CTCACCCGGTCAGCCAGGATTTCCCTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGtt                   <  1:1175083/239‑1 (MQ=18)

                                                                                                                                                                                                                            cagcGGCGGCTTTCCTCC‑CTCACCCGGTCAGCCAGGATTTCACTTCCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGtt                   <  1:1592661/239‑1 (MQ=18)

                                                                                                                                                                                                                            cagcGGCGGCTTTCCTCC‑CCCACCCGGCCAGCCAGGATTTCCCTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGtt                   <  2:757414/239‑1 (MQ=16)

                                                                                                                                                                                                                            cagcGGCGGCTTTCCTCA‑CTCCCCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTc                      <  1:2156532/236‑1 (MQ=18)

                                                                                                                                                                                                                            cagcGGCGGCTTTCCTCA‑CTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTc                      <  1:1552071/236‑1 (MQ=33)

                                                                                                                                                                                                                            cagcGGCGGCTTTCCTCA‑CTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTc                      <  1:1212007/236‑1 (MQ=33)

                                                                                                                                                                                                                            cagcGGCGGCTTTCCTCA‑CTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTc                      <  1:9411/236‑1 (MQ=33)

                                                                                                                                                                                                                            cagcGGCGGCTTTCCTCA‑CTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGtt                   <  1:1456358/239‑1 (MQ=33)

                                                                                                                                                                                                                            cagcGGCGGCTTTCCTCA‑CTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGt                    <  1:2094536/238‑1 (MQ=33)

                                                                                                                                                                                                                                ggcggcTTTCCTCA‑CTCCCCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCt                 <  1:1414285/237‑1 (MQ=18)

                                                                                                                                                                                                                                ggcggcTTTCCTCA‑CTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTc                  <  1:958221/236‑1 (MQ=33)

                                                                                                                                                                                                                                  cggtTTCCCTCC‑CCCCCCCGGCCAGCCAGGATTTCCCTTTCCAGCCGGTCGGGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTa              <  2:609718/238‑1 (MQ=11)

                                                                                                                                                                                                                                   ggcTTTCCTCA‑CTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGGGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTAtt            <  1:2020777/239‑1 (MQ=21)

                                                                                                                                                                                                                                   ggcTTTCCTCA‑CTCACCCGGTCAGCCAGGATTTCACTCTCCAGCCGGTCGAGCCTCTGCACCAGCACCGGCGGGAACACCCCACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTAt             <  1:1141383/238‑1 (MQ=18)

                                                                                                                                                                                                                                             a‑ctcacCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGcc  <  1:673120/239‑1 (MQ=37)

CAGGCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCGTATTCTGCGCACCATTCTGTTGTCCCTTCCTTGCTGATAAGGGCCAGCGGCAGGCCGCGATGGTCGCAGTGGTACAGGTGGATTTTTCGCGCCGGCGTGTACACCGGGTCCATCTGGTTTTGCATCTGCTCCACGGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCA‑CTCACCCGGTCAGCCAGGATTTCACTTTCCAGCCGGTCGAGCATCTGCACCAGCACCGGCGGGAACACCACACTGCCACCGTCTTCGCCGCCGGACTGCTGAAGCGCATCCGCCAGGCTGCGGCGCTGCGTTTTCGCCAGCTCACCGGTGGCGGTTTCAACTCTGATGAGTGGCGTGAAGCTCCCCGGCTGATAAATCGTCTGGATGCGGGTTCTGTCGTTCTGTATTGTGGTCAGCC  >  NC_007779/3874660‑3875135

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