BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,880,847	A→G	9.6%	L125P (CTC→CCC)	recF ←	gap repair protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,880,847	0	A	G	9.6%	423.9 / 2.4	156	L125P (CTC→CCC)	recF	gap repair protein
Reads supporting (aligned to +/- strand): ref base A (61/73); new base G (17/2); total (80/76)
Fisher's exact test for biased strand distribution p-value = 3.47e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 164 total aligned reads displayed.

CATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCAC‑GGGCG‑TAGCTGTTCGTAACG‑TGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGGTGCGCCAGTTCCGCGACCTTATGCCCGTCTGTACCGTCGATGCGGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAG  >  NC_000913/3880612‑3881085

cATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCAC‑GGGCG‑TAGCTGTTCGTAACG‑TGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGag                                                                                                                                                                                                                                                <  2:989540/237‑1 (MQ=255)

  tatCGGCCGCGATACCGTCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCAC‑GGCCG‑TAGCTGTTCGCACCG‑TGTCACCTGGCGCAGCGCCGCATTGCCCTGCTTGAGCAACCGCTTGAGATTGCTCCAGGCGGCGACAAATCCGGGTGCCTTGAGAAAGGATCCCCAGACGAGGAGGGGCCTTCTCTATTTGGGGCCGCCGCTGcataaa                                                                                                                                                                                                                                            >  1:4368/1‑239 (MQ=255)

        caggtagccgggCGCTAATCCCCGCGCCCCGGGTGTGAATTTGGCCCGCCGGCGGGATCAGCTCTTTATCCCAC‑GGGCG‑TAGCTGGTGGTACCG‑TGGCACTTGGCGCGGCGCCGCATTGCGTTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAgtaaagtaaa                                                                                                                                                                                                                                       <  1:48712/228‑1 (MQ=255)

         cgcgATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCAC‑GGGCG‑TAGCTGTTCGTAACG‑TGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTTTGAAAGCATCCCCAGGCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAAc                                                                                                                                                                                                                                      >  2:1322408/1‑238 (MQ=255)

         cgcgATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCAC‑GGGCG‑TAGCTGTTCGTAACG‑TGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAAGGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGGAAAc                                                                                                                                                                                                                                      >  2:1268428/1‑238 (MQ=255)

         cgcgATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCAC‑GGGCG‑TAGCTGTTCGTAACG‑TGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCACTCCCAGTCGAGGAATGCTCTTCTGGATTTGTGGCCGCCGTTGGGGAAAGCAAAc                                                                                                                                                                                                                                      >  2:668989/1‑238 (MQ=255)

            gATACCGGGTCTATCCTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCAC‑GGGCG‑TAGCTGTTCGTAACG‑TGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTCCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCtt                                                                                                                                                                                                                                  <  2:699277/239‑1 (MQ=255)

                     gctttccTCCGCGCGCCAGGTGCTTATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCAC‑GGGCG‑TAGCTGTTCGTAACG‑TGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGtt                                                                                                                                                                                                                          <  2:498953/232‑1 (MQ=255)

                         taCTCCGCGCGCCAGGTGCTGATTTGCTCCGCCATCGGGATCAGCTCTTTATCCCAC‑GGGCG‑TAGCTGTTCGTACCG‑TGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGCAAAACCCGGGTTCCGTATGACAGACTCCCCACTCCAGGAAGCGTCTTCTTGTTTTGGGCTCGCCGTTGGGGAAAGGAAAccccttccgccctacta                                                                                                                                                                                                                      >  1:315592/1‑224 (MQ=255)

                              cgcggcaacggtggtGCTTTTCCCCCCCAGCGGTATCAGCTCTGTACCCAAG‑GGGCG‑GAGCGTTACGTAACT‑GGTCCCCTGGCTCAGCGCGGCATTGCGTTGCTTGAGCATCCGCTTGAGAGTGCTCCAGGCGGTGAAAAATCCGTGTTCGTTGTGAAAGCATCCCCAGGCGAGGAATGCTCTTCTGGATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCa                                                                                                                                                                                                               <  1:208071/230‑1 (MQ=255)

                              cgcgcgCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCGC‑GGGCG‑TAGCTGTTCGTAACG‑TGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCa                                                                                                                                                                                                               <  2:1033647/239‑1 (MQ=255)

                              cgcgcgCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCAC‑GGGCG‑TAGCTGTTCGTAACG‑TGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGCAAGCATCCCCAATCGGGGAATGGCCTTCCGGTTTTGGGGCCGCCGCTGGGTAAATGAAACCCTCCTGGGGTT‑ATCaaagaca                                                                                                                                                                                                               >  1:53984/1‑234 (MQ=255)

                              cgcgcgCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCAC‑GGGCG‑TAGCTGTTCGTAACG‑TGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGGGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGGGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTgaa                                                                                                                                                                                                               >  2:951529/1‑237 (MQ=255)

                                cgcgcCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCAC‑GGGCG‑TAGCTGTTCGTAACG‑TGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGGAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGCCGAGGAAGGCTCTTCTGTCTTTTGGACCGCCGTTTGGTAAAGTCAACCCTTcctcggttttaaacagaatt                                                                                                                                                                                                              >  1:785706/1‑220 (MQ=255)

                                cgcgcCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCAC‑GGGCG‑TAGCTGTTCGTAACG‑TGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTTTGAAAGCATCCCCAGTCGAGGGATGCTCTTCTGTATTTGGGGCCGCCGTTGGGTAAAGTAAACCCTTCTTGCCGT‑TTCCACTGCAtt                                                                                                                                                                                                             >  2:286613/1‑239 (MQ=255)

                                cgcgcCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCAC‑GGGCG‑TAGCTGTTCGTAACG‑TGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGGCGAGGAATGCTCTTCTGTTTTTGGGGCCGCCGTTTGTTAAAGTAAACCCTTCCGGCGTT‑ATCaaactcatt                                                                                                                                                                                                             >  2:1050081/1‑232 (MQ=255)

                                   gcCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCAC‑GGGCG‑TAGCTGTTCGTAACG‑TGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGCGGAATGCCCTTCTTTATTTGGGGCCGCCGTTGGGGAAAGAAAAACCTTCTGGGGTT‑ATAAAACGCATggg                                                                                                                                                                                                           >  1:518052/1‑238 (MQ=255)

                                     cAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCAC‑GGGCG‑TAGCTGTTCGTAACG‑TGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGCGCAAGCCTCCCCAGGCGGGGAATGCTCTTCTTCATTTTGGGCAGCCGTTGAGTACACGAAACCCCTCTGGCGTT‑ATGCCCTAGATTGGCAt                                                                                                                                                                                                        >  1:709529/1‑239 (MQ=255)

                                                      ccccAGGGGGATCAGTTTTTTTTCCCAA‑GGGCG‑GTGCTTTTCGACCCGTTGTC‑CCTGGCGCCGGCCCGCTTCTCGCTGCTTGGCCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTGCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGGTGCGCCAGTTCcg                                                                                                                                                                                       <  1:412943/237‑1 (MQ=255)

                                                      cGCCAGCGGGATCAGCTTTTTATCCCAC‑GGGCG‑TAGCTGTTCGTAACG‑TGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGGTGCGCCAGTTCcg                                                                                                                                                                                       <  2:724748/239‑1 (MQ=255)

                                                       gCCCGCGGGTTCAGCTCTTT‑TCCCCC‑GGGGGTTGGTTGTTTGTAACT‑TGTCCCCGTGCCAAGCGCCGTATTGCGTGTCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGGTGCGCCAGTTCcgc                                                                                                                                                                                      <  1:1127816/239‑1 (MQ=255)

                                                       gCCAGCGGGATCAGCTCTTTATCCCAC‑GGGCG‑TAGCTGTTCGTAACG‑TGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGGTGCGCCAGTTCcgc                                                                                                                                                                                      <  2:390820/239‑1 (MQ=255)

                                                       gCCAGCGGGATCAGCTCTTTATCCCAC‑GGGCG‑TAGCTGTTCGTAACG‑TGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGGTGCGCCAGTTCcgc                                                                                                                                                                                      <  2:681359/239‑1 (MQ=255)

                                                         cAGCGGGAACAGCTCTTTATCCCAC‑GGGCG‑TAGCTGTTCGTAACG‑TGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGGTGCGCCAGTTCCGCGa                                                                                                                                                                                    <  2:1060260/239‑1 (MQ=255)

                                                            gggggtccgtttttttTCCCCC‑GGGGG‑TGTGTGTGCGTAACG‑TGTCCCCTGGGGCGGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCCCCAGGGGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGGTGCGCCAGTTCCGCGAcc                                                                                                                                                                                  <  1:605643/226‑1 (MQ=255)

                                                                aTCAGCTCTTTACCCCAC‑GGGGG‑TAGCTGTTCGGAACG‑TGTCCCCTGGCGCAGCGCCGCATTGCGCTGCTTGAGC‑ATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGTTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGGTGCGCCAGTTCCGCGACCTTATg                                                                                                                                                                             <  2:1235855/238‑1 (MQ=255)

                                                                 tcgccgttttATCCCAA‑GGGCG‑TGTTTTTTCGTAACT‑TGTCACCTGGCGCAGGGCAGCATTGCCCTGCTGGGGCAATCGCTTGAAATTGCCCCAGGCGGTTAAAAATCCGGGTTCGTGGTGAAAGCATCCCCAGTTGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGGTGCACCAGTTCCGCGACCTTATg                                                                                                                                                                             <  1:1173047/231‑1 (MQ=255)

                                                                 tCAGGTCTTTACCCCAC‑GGGCG‑TAGTGGTTCGTAACG‑GTTCCCCTGGCGCAGCCCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGGTGCGCCAGTTCCGCGACCTTATg                                                                                                                                                                             <  2:1148304/238‑1 (MQ=255)

                                                                 tCAGCTTTTTATCCCTG‑GGGTG‑GAGCTTGTCGGAACG‑TGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGGTGCGCCAGTTCCGCGACCTTATg                                                                                                                                                                             <  2:125760/238‑1 (MQ=255)

                                                                 tCAGCCTTTTATCCCAC‑GGGCG‑TAGCTGTTCGTAACG‑TGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGGTGCGCCAGTTCCGCGACCTTATg                                                                                                                                                                             <  2:664765/238‑1 (MQ=255)

                                                                 ctgcctttttATTCCCC‑GGGCG‑TGGCTGTTCGTAAGG‑TGTACCCTGGGGCTGGGCTGCATTGTGGTGGTTGAGCAATCGCTTGGGATTGCTCCAGGGGGTTAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGCAGTTGGGGCCGCCGTTGAGTAAAGTAACCCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGGTGCGCCAGTTCCGCGACCTTATg                                                                                                                                                                             <  1:1247645/231‑1 (MQ=255)

                                                                  ttcgctttttTCCCACGGGGCT‑TGGGTTTTCGATACG‑TGGCCCTTGGCGCCGGGCCGTCTTGGGCTGCTTGGACAATTGCTTGAGCTTGCTCCAGGCGGTGAAAAATCCGGGTTTGTTTTGAAAGCCTCCCCAGTCGAGGAATGCTCTTCTGTATTGGGGGCCGCCGTTGAGTAAAGTAAGCCCTTCTGGCGTT‑ATCAACTGCATTGGCTTCACGTGCGCCAgccgcttggacttatg                                                                                                                                                                             <  1:772843/232‑16 (MQ=255)

                                                                     tctttttcccccC‑GGGCT‑TAGCTTTCTGTATCG‑TGTCACGTGGCCATGCGCCGAATTGAGCGGTGTGAGCGATGGCTTGAGATTGCTCCGGGCGGTGAGAAACCGGGGTTCGTGGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCGGTATTTGGGGCGGCCGTTGAGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGGTGCGCCAGTTCCGCGACCTTATGCCCGt                                                                                                                                                                        <  1:169003/231‑1 (MQ=255)

                                                                        tttATCCCAC‑GGGCG‑TAGCTGTTCGTAACG‑TGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGCTTGCTCCAGGCGGCGCAAACTCCGGGTTAGCTGTGAAAACATCCCCCGTCGAGGAGAGCTCTTCTGTCTTTGGGTCCGCCGTTGGGTAAAGTAAACCCTTCCGGCGTTAATAAAGTCCA‑TAGCATCTGTTGCGCCAGTCCCCGGACCtttaccgcc                                                                                                                                                                         >  1:850952/1‑228 (MQ=255)

                                                                          tATCCCAC‑GGGCG‑TAGCTGTTCGTAACG‑TGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGGGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGGTGCGCCAGTTCCGCGACCTTATGCCCGTCTGt                                                                                                                                                                    >  2:288580/1‑238 (MQ=255)

                                                                             cccAG‑GGGCG‑TAGGTGTTCGTAACG‑TGTCCCCTGGCGCAGCGCCGCATGGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGGTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGGTGCGCCAGTTCCGCGACCTTATGCCCGTCTGTACCg                                                                                                                                                                <  2:284938/239‑1 (MQ=255)

                                                                                aC‑GGGCG‑TAGCTGTTCGTAACG‑TGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGGTGGGCCAGTTCCGCGACCTTATGCCCGCCTGTACCGTCg                                                                                                                                                             >  2:318828/1‑239 (MQ=255)

                                                                                   gggCG‑TAGCTGTTCGTAACG‑TGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGGTGCGCCAGTTCCGCGGCCTTATGCCCGTCTGTACCGTCCAt                                                                                                                                                           >  2:1133182/1‑239 (MQ=255)

                                                                                   gggCG‑TAGCTGTTCGTAACG‑TGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGGGAAACATCCGGGTTCGTTTTGTAAGCTTCCCCAGCCGAGCGATGCTCTTCAGAATTTGGGGCACCCGTTGTGTAGAGAAAACCCTTCTGGCGTT‑ATAAACTGCATTGGCACCAGGGGCGCCAATTCCGCGACCTTCTGCCCCTCCGTAcccactag                                                                                                                                                           >  1:1283503/1‑233 (MQ=255)

                                                                                             tGTTCGGAACG‑TGTCACCTGGCGCGGGGCCGCATGGCGCGGTTTGAGCAATCGGTTGAGATGGCTCCGGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATTCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGGTGCGCCAGTTCCGCGACCTTATGCCCGTCTGTACCGTCGATGCGGACtt                                                                                                                                                   <  1:288580/238‑1 (MQ=255)

                                                                                                tCGTAACG‑TGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGGGTAAAGTAAACCCTTCTGGCGTT‑TTCAACTGCATTTGCATCCGGTGGGCCCGTTCCGCGACCTTATTCCCTTCTTTACCCTCGATGCGGACttt                                                                                                                                                  >  2:1196325/1‑236 (MQ=255)

                                                                                                  gTAACG‑TGTCCCCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCCCCAGGCGGTGAAAAATCCGGGTTCGGTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGGTGCGCCAGTTCCGCGACCTTATGCCCGTCTGTACCGTCGATGCGGACTTTGCTGt                                                                                                                                             <  2:597911/239‑1 (MQ=255)

                                                                                                  gTAACG‑TGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGGTGCGCCAGTTCCGCGACCTTATGCCCGTCTGTACCGTCGATGCGGACTTTGCTGt                                                                                                                                             <  2:1103661/239‑1 (MQ=255)

                                                                                                          gtCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGGTGCGCCAGTTCCGCGACCTTATGCCCGTCTGTACCGTCGATGCGGACTTTGCTGTCGccccg                                                                                                                                      >  2:1112317/1‑237 (MQ=255)

                                                                                                          gtCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGGATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGGTGCGCCAGTTCcg                                                                                                                                                                                       >  1:914996/1‑190 (MQ=255)

                                                                                                           tccccTGGGGCAGGGCCGCATTGCGTTGTTTTAGCAATCGCTTGGGATTCCTCCAGGCGGGGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGGTGCGCCAGTTCCGCGACCTTATGCCCGTCTGTACCGTCGATGCGGACTTTGCTGTCGCCCTGt                                                                                                                                     <  1:486550/236‑1 (MQ=255)

                                                                                                           tCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGGTGCGCCAGTTCCGCGACCTTATGCCCGTCTGTACCGTCGATGCGGACTTTGCTGTCGCCCTGt                                                                                                                                     <  2:106733/239‑1 (MQ=255)

                                                                                                           gccccttgcggcgggttgtcttGCGCTGGTTGGGCCATCGGTTTAGTTTGCATCGGGCGGTGAAAAATCCGGGTTCGTTTTGAAAGAATCCACAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAACCCCTTCTGGCGTT‑ATCACCTGCATTGGCATCAGGTGCGCCAGTTCCGCGACCTTATGCCCGTCTGTACCGTCGATGCGGACTTTGCTGTCGCCCTGt                                                                                                                                     <  1:1268428/219‑1 (MQ=255)

                                                                                                            cgccTGGCGCAGCGCCGCATGGCGCTGCTTGAGCATTCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGGCGTTGAGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGGTGCGCCAGTTCCGCGACCTTATGCCCGTCTGTACCGTCGATGCGGTCTTTGCTGTCGCCCTGt                                                                                                                                     <  2:327358/236‑1 (MQ=255)

                                                                                                            cccctggtggcgggccgcATTGCGCTTTTGGACCAAGCGTATGAGATTGCTCCGGGCGGTGAAAAGTCCGGGTTCGTTGCGAAAGCATCCCCAGGCGAGGAATGCTCTTCTTTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTGCTGGCGTT‑ATCAACTGCATTGGCATCAGGTGCGCCAGTTCCGCGACCTTATGCCCGTCTGTACCGTCGATGCGGACTTTGCTGTCGCCCTGt                                                                                                                                     <  1:519977/225‑1 (MQ=255)

                                                                                                                 ggCCCAGCGCCGCTTTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGGTGCGCCAGTTCCGCGACCTTATGCCCGTCTGTACCGTCGATGCGGACTTTGCTGTCGCCCTGTTTGTCt                                                                                                                               <  2:1283503/239‑1 (MQ=255)

                                                                                                                  gcgcAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGACAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGTTGCGCCCGTTCCGCGACCTTTTGCCCGTCTGTACCGTCGATGCGGACTTTGCTGTCGCCCCGTTTGTCtt                                                                                                                              >  2:75380/1‑239 (MQ=255)

                                                                                                                       gcgcCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGGTGCGCCAGTTCCGCGACCTTATGCCCGTCTGTACCGTCGCTGCGGACTTTGCTGGCGCCCTGTTTGTCTTTGGtt                                                                                                                         >  2:1182136/1‑239 (MQ=255)

                                                                                                                        cgccgcATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTTAAAAATCCGGGTTCGTTGTGGAAGCATGCCCAGTAGAGGAATGCTCTTTTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTG‑CTCAACTGCATGGTCATCAGGTGCGCCAGTTCACCGACCTTATGCCCGTCTTATACGTCGCTGCGGCCTTTGCTGTCCCCCTGTTTTCCTTTGGGTa                                                                                                                        >  1:1131734/1‑239 (MQ=255)

                                                                                                                        cgccgcATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGGTCGTTGTGAAAGCATCCCCAGTCGCGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCAATGGCCTCAGGTGCGCCAGTTCCGCGACCTTATGCCCGTCTTTCCCTTCGATGCGGACTTTCCTGTCGCCCTGTTTGCCTTTGGTTa                                                                                                                        >  2:134771/1‑239 (MQ=255)

                                                                                                                            ggaTTGCGCTGGTTGAACAATCGCTTGAAATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAACGATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGGTGCGCCAGTTCCGCGACCTTATGCCCGTCTGTACCGTCGATGCGGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGcc                                                                                                                    <  2:349973/237‑1 (MQ=255)

                                                                                                                                  cgcTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGGATGCGCTTCTGTATTTTGGGCCGCCGGTGAGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCATTGGCATCCGGTGGGCCAGTTCCACGACCCTAAGCCCGTCTTGACCCGCGATGCCGGCCTTGCTGTCGCCCCTTTTTTCTTTTGTTAAGACAAACGc                                                                                                              >  1:1055879/1‑239 (MQ=255)

                                                                                                                                       ctTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGGAGCGCCAGTTCCGCGACCTTATGCCCGTCTGTACCGTCGATCCGGACTTTGCTTTCGCCCTGTTTGTCTTTGGTTAAGCCAATcccggtt                                                                                                          >  2:93110/1‑232 (MQ=255)

                                                                                                                                         tGAGAAATCGCTTGAGTTTGGTCCAGGGGGTGAAAAATCCGGGTTGGTGGGGAAAGCATCCCCAGGCGAGGAATGCTTTTCTGAATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCGGGCGTT‑ATCAACTGCATTGGCATCAGGTGCGCCAGGTCCGCGACCTTATGCCCGTCTGTACCGTCGATGCGGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTcg                                                                                                       <  1:75380/239‑1 (MQ=255)

                                                                                                                                            gCAATCGCTTGAGCTTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGGTGCGCCAGTTCCGCGACCTTATGCCCGTCTGTACCGTCGATGCGGACTTTGCTGTCGCCCTGTTTGTCTTTTGTTAAGCCAATCGCCGGCTcccg                                                                                                     >  2:1044141/1‑238 (MQ=255)

                                                                                                                                             cAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGGTGCGCCAGTTCCGCGACCTTATg                                                                                                                                                                             <  2:211088/165‑1 (MQ=255)

                                                                                                                                                tgccTTGGGATTCCTCGAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCACCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGGTGCGCCAGTTCCGCGACCTTATGCCCGTCTGTACCGTCGATGCGGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCtt                                                                                                <  1:1017009/236‑1 (MQ=255)

                                                                                                                                                tCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGGTGCGCCAGTTCCGCGACCTTATGCCCGTCTGTACCGTCGATGCGGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCtt                                                                                                <  2:518052/239‑1 (MQ=255)

                                                                                                                                                tCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCCTCTGGCGTT‑ATCAACTGCATTGGCATCAGGTGGACCCGTTCCCCGACCTTATGCCCGTCAGTACCGTCGATGCGGACCTTGCTGTGGCCCCGTTTTTCTTTGGTTAGACCaaatggtgttcgcccccctc                                                                                                >  1:1025505/1‑221 (MQ=255)

                                                                                                                                                tCGCTTGAGATTGCTCCAGGCGGTGAAAAAGCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGGTGCGCCAGTTCCGCGACCTTATGCCCGTCTGTACCGTCGATGCGGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCtt                                                                                                <  2:3246/239‑1 (MQ=255)

                                                                                                                                                 cGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGCGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTTGCGTT‑ATCAACTGCATTGGCATCAGGTGCGCCAGTTCCGCGACCTTATGCCCGTCTGTACCGTCGATGCGGACTTTGCTGTCGCCCTGTTTTTCTTTGGTTAAGCCAATCGCTGTTTCCCGCCCtt                                                                                                >  2:11689/1‑238 (MQ=255)

                                                                                                                                                 cGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGGTGCGCCAGTTCCGCGACCTTATGCCCGTCTGTACCGTCGATGCGGACTTTGCTGTCGCCCTGttt                                                                                                                                   >  2:928043/1‑203 (MQ=255)

                                                                                                                                                   cTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCATT‑ATCAACTGCATTGGCATCAGGTGCGCCAGTTCCGCGACCTTATGCCCGTCTGTACCGTCGATGCGGACTTTGCTGTCGCCCTGt                                                                                                                                     <  2:553649/199‑1 (MQ=255)

                                                                                                                                                             cTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGGTGCGCCAGTTCCGCGACCTTATGCCCGTCTGTACCGTCGATGCGGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCCCTTCCCCCCGTAATc                                                                                    >  2:435492/1‑238 (MQ=255)

                                                                                                                                                              tgtaGGGGGGGAAAAACCCGGTTTCTTGGGGAATGCATCCCCAGTCAAGGAATGCTTTTCTGTA‑TTGGCGCCCGCGTGGATTAAAGTAAACCCTTCTGGCTTT‑ATCAACTGCATTGGCATCAGGTGCGCCAGTTCCGCGACCTTATGCCCGTCTGTACCGTCGATGCGGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCGCGCGCTCTTCGCCCTGTAATCGc                                                                                  <  1:182718/235‑1 (MQ=255)

                                                                                                                                                                  ggCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGGTGCGCCAGTTCCGCGACCTTATGCCCGTCTGTACCGTCGATGCGGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTTTCTCGCCCCCTTCGCCCTGAAATTGccccg                                                                              >  2:99622/1‑237 (MQ=255)

                                                                                                                                                                  ggCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGGTGCGCCAGTTCCGCGACCTTATGCCCGTCTGTACCGTCGATGCGGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCgg                                                                              >  2:1149848/1‑238 (MQ=255)

                                                                                                                                                                    ggtggAAAAATCTGGGTCTGGTGGGAAACCGTCCCCAGTCAAGGAATGCTTTTCTGTATTGGGGGCCGCCGTTGGGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGGTGCGCCAGTTCCGCGACCTTATGCCCGTCTGTACCGTCGATGCGGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCg                                                                               <  1:266439/232‑1 (MQ=255)

                                                                                                                                                                     ggggAAAATTCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGGTGCGCCAGTTCCGCGACCTTATGCCCGTCTGTACCGTCGATGCGGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTgg                                                                            <  2:839479/235‑1 (MQ=255)

                                                                                                                                                                             atttggtttCGTTGTGAAAGCATCCCCGGTCGGGGAATGCTTTTTTGTATTGGGGGCCGGCGTTTAGTAAAGTAAACCCCTCTGGCGTT‑ATCAACTGCATGGGCATAAGGTGCGCCAGTTCCGCGACCTTATGCCCGTCTGTACCGTCGATGCGGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAAc                                                                   <  1:1044141/232‑1 (MQ=255)

                                                                                                                                                                                 gggTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGGTGCGCCAGTTCCGCGACCTTATGCCCGTCTGTACCGTCGATGCGGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGCCTCGCGCTCTTCGCCCGGTAATCGCCCGTGGGGAAAAAACGcc                                                                >  1:155903/1‑238 (MQ=255)

                                                                                                                                                                                   tttCGTTGGGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGGTGCGCCAGTTCCGCGACCTTATGCCCGTCTGTACCGTCGATGCGGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGcctc                                                              <  2:630147/237‑1 (MQ=255)

                                                                                                                                                                                     tCGTTGTGAAAACATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGGTGCGCCAGTTCCGCGACCTTATGCCCGTCTGTACCGTCGCTGCGGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGGGGAGAACAAACGcctcct                                                            >  2:709739/1‑238 (MQ=255)

                                                                                                                                                                                                   tCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGGAGCGCAAGTTCCGCGACCTTATGCCCGTCTGTACCGTCGATGCGGACTTTGCTGTCGCCCTGTTTGTCTTTGTTTAAGCCAATCGCTCTCTCGCGCTCTTCGCCCTGTACTCGCCCGTGGAGAACAAACCACTCCTGCTCATGGCGAATg                                              >  2:141363/1‑238 (MQ=255)

                                                                                                                                                                                                        atttggggATTGTTCTTTTGTTGTTGGGGCCGCCGTTGCGTAAAGTACACCCTTCTTGCGTT‑AGCACCTGCATGGGCATCAGGTGCGCCAGTTCCGCGACCTTATGGCCGTCTGTACCGTCGATGCGGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGa                                        <  1:1182136/233‑1 (MQ=255)

                                                                                                                                                                                                           cGAGGAATGCTCTTTTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGGTGCGCCAGTTCCGCGACCTTATGCCCGTCTGTACCGTCGATGCGGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGa                                        <  2:443503/236‑1 (MQ=255)

                                                                                                                                                                                                           cGAGGAATGCTCTTCTGTATGTGGGGCCGCCGTTGAGTAACGTCATCCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGGTGCGCCAGTTCCGCGACCTTATGCCCGTCTGTACCGTCGATGCGGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGAcc                                      <  2:785706/238‑1 (MQ=255)

                                                                                                                                                                                                             aGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGGTGCGCCAGTTCCCCGACCTTATGCCCGTCTGTACCGTCGATGCGGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTTTCCCGCTCTTCCCCCTGTCATCGCCCGTGGTGAAAAAACGCCCCCCGCTCTTGCCGaaacacccacccaa                                    >  1:432360/1‑226 (MQ=255)

                                                                                                                                                                                                                  ttcttttttttttttGGGCCCGCCGTTGAGTAAAGTAACCCCTTCTGCGGTT‑ATCAACTGCATTGGCACCAGGGGCGCCAGTTCCGCGCCCTTATGCCCGGCTGTCCCGGCGATGCGGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCa                              <  1:1322408/227‑1 (MQ=255)

                                                                                                                                                                                                                    ctctTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGGTGCGCCAGTTCCGCGACCTTATGCCCGTCTGTACCGTCGATGCGGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCaa                             >  2:877703/1‑238 (MQ=255)

                                                                                                                                                                                                                    ctctTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGGTGCGCCAGTTCCGCGACCTTATGCCCGTCTGTACCGTCGATGCGGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGCATCGCGCTTGTCGCCCTGTAATTGCCCTGGGAGAACAAACGCCTCCTGCTCTTGGCGatagcaccgaccacactccaac                            >  1:38494/1‑218 (MQ=255)

                                                                                                                                                                                                                          tGTTTTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTT‑AACAATTGCATTGGCATCAGGGGGGCCAGTTCCGCGACCTTATGCCCGTCTGTACCGTCGATGCGGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGa                      <  1:1149848/239‑1 (MQ=255)

                                                                                                                                                                                                                          tGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGGTGCGCCAGTTCCGCGACCTTATGCCCGTCTGTACCGTCGATGCGGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGa                      >  2:583236/1‑239 (MQ=255)

                                                                                                                                                                                                                            tATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGGTGCGCCAGTTCCGCGACCTTATGCCCGTCTGTACCGTCGATGCGGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCCCTCCGCCCCGTAATCGCCCGTGGAGTACAAACGCCTCCTGCCCAGGGCGAATGAAGCACCCAAACTGaaaaatccgca                     >  1:1079473/1‑231 (MQ=255)

                                                                                                                                                                                                                               tttgggCCCTCGTTGGTTAAAATAAACCCCTCTGGCGTT‑ATCAACTGCCTTGGCATGAGTTGCGCCAGTTCCGCGGCCTTATGCCCGTCTGTACCGGCGATGCGGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCGCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGc                 <  1:668989/236‑1 (MQ=255)

                                                                                                                                                                                                                                 ggggCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGGTGCGCCAGTTCCGCGACCTTATGCCCGTCTGTACCGTCGATGCGGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCTTGGACAACAACCGCCTCCTGCTCCTGGCGAATAGCGGGCCCAACATGCACACTTCGaaaaacc                >  1:604520/1‑234 (MQ=255)

                                                                                                                                                                                                                                 ggggCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGGTGCGCCAGTTCCGCGACCTTATGCCCGTCTGTACCGTCGATGCGGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGGGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGccc               >  2:737303/1‑239 (MQ=255)

                                                                                                                                                                                                                                 ggggCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGGTGCGCCAGTTCCGCGACCTTATGCCCGTCTGTACCGTCGATGCGGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAAACTGCAAACTGCGAAACAcc                >  2:970437/1‑238 (MQ=255)

                                                                                                                                                                                                                                   ggCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGGTGCGCCAGTTCCGCGACCTTATGCCCGTCTGTACCGTCGATGCGGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCg                                                                                              <  1:970045/158‑1 (MQ=255)

                                                                                                                                                                                                                                      cgccgTTGAGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGGTGCGCCAGTTCcgc                                                                                                                                                                                      >  2:1111241/1‑67 (MQ=255)

                                                                                                                                                                                                                                      cgccgTTGAGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGGTGCGCCAGTTCCGCGACCTTATGCCCGTCTGTACCGTCGATGCGGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCa          >  2:1288160/1‑239 (MQ=255)

                                                                                                                                                                                                                                          gTTTAGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGGTGCGCCAGTTCCGCGACCTTATGCCCGTCTGTACCGTCGATGCGGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGc      <  2:458993/239‑1 (MQ=255)

                                                                                                                                                                                                                                           tttagtaaagtaaCCCCTTCGGGCGTT‑ATCAACTGCATTGGCATAAGGTGCGCCAGTTCCGCGCCCTGATGCCCGTCTGTACCGTCGATGCGGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGcc     <  1:70807/236‑1 (MQ=255)

                                                                                                                                                                                                                                           ttGAGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGGTGCGCCAGTTCCGCGACCTTATGCCCGTCTGTACCGTCGATGCGGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGcc     <  2:47611/239‑1 (MQ=255)

                                                                                                                                                                                                                                              agtaaagtaaACCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGGTGCGCCAGTTCCGCGACCTTATGCCCGTCTGTACCGTCGATGCGGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAg  <  2:155903/239‑1 (MQ=255)

CATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCAC‑GGGCG‑TAGCTGTTCGTAACG‑TGTCACCTGGCGCAGCGCCGCATTGCGCTGCTTGAGCAATCGCTTGAGATTGCTCCAGGCGGTGAAAAATCCGGGTTCGTTGTGAAAGCATCCCCAGTCGAGGAATGCTCTTCTGTATTTGGGGCCGCCGTTGAGTAAAGTAAACCCTTCTGGCGTT‑ATCAACTGCATTGGCATCAGGTGCGCCAGTTCCGCGACCTTATGCCCGTCTGTACCGTCGATGCGGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAG  >  NC_000913/3880612‑3881085

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