BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AE015451	1,499,508	+C	intergenic (+85/‑79)	trpS → / → zapE	Tryptophan‑‑tRNA ligase/putative enzyme with nucleoside triphosphate hydrolase domain

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	1,499,506	1	.	C	99.4%	316.7 / ‑0.9	317	intergenic (+83/‑81)	trpS/zapE	Tryptophan‑‑tRNA ligase/putative enzyme with nucleoside triphosphate hydrolase domain
Reads supporting (aligned to +/- strand): ref base . (1/1); new base C (128/186); total (129/188)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 336 total aligned reads displayed.

TTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGTT‑T‑CGGGCAC‑G‑C‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGT  >  AE015451/1499383‑1499625

  cAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCAc‑g‑c‑c                                                                                                                       <  2:2057644/125‑1 (MQ=25)

    gagCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCAc‑g‑c‑cc                                                                                                                      >  1:2058227/1‑123 (MQ=21)

      gCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCAc‑g‑c‑c                                                                                                                       <  2:2057649/121‑1 (MQ=21)

       ccTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGACCGCGAAACAGGCGCCGCGGTGCATGTGACCGGCTTTGCCGGTGTTCG‑CGGGCAc‑g‑c                                                                                                                         <  1:2057650/119‑1 (MQ=21)

             cGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCAC‑GCC‑CGCTCCccca                                                                                                              >  2:2057825/1‑124 (MQ=18)

             cGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCAC‑GCC‑CGCTCCcac                                                                                                               >  2:2057914/1‑123 (MQ=18)

              gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CAGGCAC‑GCC‑CGCTCCcac                                                                                                               >  1:2058402/1‑122 (MQ=18)

               aaGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTT‑CGCGGGCAC‑GCC‑CGCTCCCACAGAc                                                                                                           <  2:2057655/125‑1 (MQ=18)

                aGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGATGCATGGCACCGGCTTTGCCGGTGTT‑CGCGGGCAC‑GCC‑CGCTCCCACAg                                                                                                             <  1:2057656/122‑1 (MQ=18)

                 gCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCAC‑GCC‑CGCTCCcac                                                                                                               >  2:2058146/1‑119 (MQ=17)

                 gCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCAC‑GCC‑CGCTCCcac                                                                                                               >  1:2058248/1‑119 (MQ=17)

                  ccTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTT‑CGCGGGCAC‑GCC‑CGCTCCCACAGACTGa                                                                                                        <  2:2057660/125‑1 (MQ=18)

                      ggCACCACAGCAGGACTGAGGCTGCCAAGCGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTT‑CGCGGGCAC‑GCC‑CGCTCCCACAGACTGAt                                                                                                       <  1:2057667/122‑1 (MQ=17)

                      ggCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGGGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTT‑CGCGGGCAC‑GCC‑CGCTccc                                                                                                                 <  2:2057665/112‑1 (MQ=17)

                          ccacAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTT‑CGCGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCAtt                                                                                                <  1:2057673/125‑1 (MQ=18)

                           cacaGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTa                                                                                               >  2:2058055/1‑125 (MQ=18)

                           cacaGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTT‑CGCGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCa                                                                                                  <  1:2057678/122‑1 (MQ=18)

                           cacaGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTT‑CGCGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCAt                                                                                                 <  2:2057676/123‑1 (MQ=18)

                           cacaGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGCGTTCG‑CGGGCAC‑GCC‑CGCTCCCACAGACTg                                                                                                         >  1:2058129/1‑115 (MQ=17)

                             cagcagGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCa                                                                                                  >  2:2058195/1‑120 (MQ=18)

                               gcagGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGGACGG‑C‑CGCTCACACAGACTGATTg                                                                                                     >  1:2058297/1‑115 (MQ=17)

                                 agGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCAC‑GCC‑CGCTCCCACAGACTGATTg                                                                                                     >  2:2058255/1‑113 (MQ=17)

                                  ggACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTT‑CGCGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGc                                                                                        <  1:2057682/125‑1 (MQ=18)

                                   taCTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                                       >  1:2057888/2‑125 (MQ=18)

                                   gACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCAC‑GCC‑CGCTCCCACAGACTGa                                                                                                        >  2:2058159/1‑108 (MQ=14)

                                   gACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                                       >  1:2058024/1‑125 (MQ=18)

                                   gACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTT‑CGCGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                                       <  2:2057684/125‑1 (MQ=18)

                                    aCTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCt                                                                                      >  1:2058059/1‑125 (MQ=18)

                                       gAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATc                                                                                   <  2:2057693/125‑1 (MQ=18)

                                       gAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTT‑CGCGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTAt                                                                                    <  1:2057689/124‑1 (MQ=18)

                                       gAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTT‑CGCGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTAt                                                                                    <  1:2057694/124‑1 (MQ=18)

                                        aGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCg                                                                                  >  1:2058003/1‑125 (MQ=18)

                                         ggCTGCCAAACGGTGGGGGCGGGCGTGCCCGCGAGACAGGCGCCGCCGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGCCAC‑GCC‑CGCTCTCACAGACTGATTGCCATTAAGCGGGCCTATCGc                                                                                 <  1:2057697/125‑1 (MQ=18)

                                         ggCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGc                                                                                 <  2:2057699/125‑1 (MQ=18)

                                          gCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGc                                                                                 >  2:2057858/1‑124 (MQ=18)

                                          gCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTT‑CGCGGGCAC‑GCC‑CGCTCCCACAGACTGAt                                                                                                       <  1:2057704/102‑1 (MQ=14)

                                            tGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCt                                                                                >  2:2058023/1‑123 (MQ=18)

                                              ccAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCt                                                                                >  2:2058016/1‑121 (MQ=17)

                                              ccAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTT‑CGCGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAg                                                                            <  2:2057709/125‑1 (MQ=18)

                                               cAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTAc                                                                              >  2:2057996/1‑122 (MQ=18)

                                               cAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTT‑CGCGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGt                                                                           <  1:2057710/125‑1 (MQ=18)

                                                aaaCGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGt                                                                           >  1:2058541/1‑124 (MQ=18)

                                                 aaCGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGa                                                                         <  1:2057719/125‑1 (MQ=18)

                                                 aaCGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTT‑CGCGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGa                                                                         <  2:2057715/125‑1 (MQ=18)

                                                 aaCGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTT‑CGCGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGa                                                                         <  2:2057718/125‑1 (MQ=18)

                                                    ggTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACgg                                                                      <  2:2057727/125‑1 (MQ=18)

                                                    ggTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACgg                                                                      <  1:2057725/125‑1 (MQ=18)

                                                      tGGGACCGGTCGTGCCCGCGAAACAGGCGCCACGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCAC‑GCC‑CGATCCCGCAGACTGATTGCCATTAAGCg                                                                                           >  2:2058436/1‑102 (MQ=14)

                                                       gggAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTT‑CGCGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGGCTATCGCTACAGTGACGGcc                                                                    <  1:2057731/124‑1 (MQ=17)

                                                       gggAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTT‑CGCGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGccc                                                                   <  2:2057736/125‑1 (MQ=18)

                                                       gggAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTT‑CGCGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGcc                                                                    <  1:2057737/124‑1 (MQ=18)

                                                         gAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTT‑CGCGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCgg                                                                 <  2:2057740/125‑2 (MQ=18)

                                                           gCGGGCGTGCCCGCGAAACAGGTGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTT‑CGCGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGtt                                                                <  1:2057744/124‑1 (MQ=18)

                                                            cGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTT‑CGCGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGtttt                                                              <  1:2057748/125‑1 (MQ=21)

                                                             gggCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTT‑CGCGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGttttt                                                             <  2:2057751/125‑1 (MQ=21)

                                                              ggCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGtttttt                                                            <  2:2057754/125‑1 (MQ=21)

                                                              ggCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTT‑CGCGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGtttttt                                                            >  2:2058135/1‑125 (MQ=21)

                                                                  tGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTT‑CGCGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgttg                                                        <  2:2057760/125‑1 (MQ=21)

                                                                  tGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTT‑CGCGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgtt                                                         >  2:2058042/1‑124 (MQ=21)

                                                                   gCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgttg                                                        <  2:2057765/124‑1 (MQ=21)

                                                                   gCCCGCGAAACAGGCGCCGCCGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGc                                                       <  1:2057763/125‑1 (MQ=21)

                                                                     ccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTT‑CGCGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGcc                                                      <  2:2057770/124‑1 (MQ=21)

                                                                        cgAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGc                                                  >  1:2058076/1‑125 (MQ=25)

                                                                        cgAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGc                                                  <  2:2057774/125‑1 (MQ=25)

                                                                        cgAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTT‑CGCGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCg                                                   <  1:2057775/124‑1 (MQ=25)

                                                                          aaaCAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTCTTTTTTTGCCTCGCTa                                                >  1:2058141/1‑125 (MQ=25)

                                                                           aaCAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTaa                                               >  1:2058107/1‑125 (MQ=25)

                                                                           aaCAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTaa                                               >  1:2058454/1‑125 (MQ=25)

                                                                           aaCAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTaa                                               <  1:2057780/125‑1 (MQ=25)

                                                                            aCAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAAc                                              >  2:2058119/1‑125 (MQ=25)

                                                                             cAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTT‑CGCGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTaa                                               <  2:2057782/123‑1 (MQ=25)

                                                                               ggCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGaa                                           >  1:2058165/1‑125 (MQ=25)

                                                                                  gccgcGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                        <  2:2057787/125‑1 (MQ=25)

                                                                                  gccgcGGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                        >  1:2058322/1‑125 (MQ=25)

                                                                                  gccgcGGTGCATGGCACCGGCTTTGCCGGTGTT‑CGCGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                        <  1:2057788/125‑1 (MQ=25)

                                                                                     gcgGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGa                                     >  2:2058208/1‑125 (MQ=25)

                                                                                      cgGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGa                                     <  2:2057793/124‑1 (MQ=25)

                                                                                      cgGTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGAc                                    >  2:2058099/1‑125 (MQ=25)

                                                                                        gTGCATGGCACCGGCTTTGCCGGTGTTCG‑CGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTc                                  <  2:2057796/125‑1 (MQ=37)

                                                                                          gCATGGCACCGGCTTTGCCGGTGTT‑CGCGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTccc                                <  1:2057798/125‑1 (MQ=37)

                                                                                             tggcaccGGCTTTGCCGGTGTTCG‑CGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa                             <  2:2057804/119‑1 (MQ=37)

                                                                                             tggcaccGGCTTTGCCGGTGTTCG‑CGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa                             >  1:2057938/7‑125 (MQ=37)

                                                                                              ggcaccGGCTTTGCCGGTGTT‑CGCGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAg                            <  2:2057810/120‑1 (MQ=37)

                                                                                                  ccGGCTTTGCCGGTGTTCG‑CGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACg                        <  1:2057805/124‑1 (MQ=38)

                                                                                                   cGGCTTTGCCGGTGTTCG‑CGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGc                       <  2:2057815/125‑1 (MQ=38)

                                                                                                    ggCTTTGCCGGTGTT‑CGCGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGaa                          <  1:2057817/121‑1 (MQ=38)

                                                                                                     gCTTTGCCGGTGTTCG‑CGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACg                        >  2:2058177/1‑122 (MQ=38)

                                                                                                     gCTTTGCCGGTGTTCG‑CGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCt                      >  2:2058282/1‑124 (MQ=38)

                                                                                                      cTTTGCCGGTGTT‑CGCGGCCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGACCCGTTTTTTGTTGCCTCGCTAACGAATCACGACTCCCTTGGAACGCTAt                    <  1:2057821/125‑1 (MQ=37)

                                                                                                       tttGCCGGTGTTCG‑CGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATc                   >  2:2058331/1‑125 (MQ=38)

                                                                                                       tttGCCGGTGTTCG‑CGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATc                   >  1:2058365/1‑125 (MQ=38)

                                                                                                         tGCCGGTGTTCG‑CGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCaa                 >  2:2058084/1‑125 (MQ=39)

                                                                                                         tGCCGGTGTTCG‑CGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGATACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCaa                 >  2:2058237/1‑125 (MQ=38)

                                                                                                                   tcgcGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTg        >  1:2058093/4‑125 (MQ=255)

                                                                                                                    cgcGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATc          <  2:2057835/120‑1 (MQ=255)

                                                                                                                    cgcGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGa       <  2:2057840/123‑1 (MQ=255)

                                                                                                                      cGGGCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaaa     <  1:2057834/125‑1 (MQ=255)

                                                                                                                       gggCAC‑G‑CTCGCTCCCACAGACGGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaa      <  1:2057845/123‑1 (MQ=255)

                                                                                                                        ggCAC‑GCC‑CGCTCGCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACg   >  1:2058263/1‑125 (MQ=255)

                                                                                                                         gCAC‑GCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGt  <  2:2057849/125‑1 (MQ=255)

TTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGTT‑T‑CGGGCAC‑G‑C‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGT  >  AE015451/1499383‑1499625

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 15 ≤` `ATCG/ATCG` `< 27 ≤` `ATCG/ATCG` `< 34 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`