BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_014034	1,395,865	G→C	7.8%	intergenic (+677/‑254)	RCAP_rcc01306 → / → RCAP_rcc01307	Cas2 family CRISPR‑associated protein/hypothetical protein

Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_014034	1,395,865	0	G→C	7.8%	72.3 ‑1.0 ‑1.1	449	intergenic (+677/‑254)	RCAP_rcc01306/RCAP_rcc01307	Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand): new base (14/21): ref base (202/212): total (216/233)
Fisher's exact test strand distribution p-value = 3.80e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.46e-02

Only 100 of 478 total aligned reads displayed.


CGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGC.CGCG.C.A.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGG  >  NC_014034/1395766‑1395965


cgCCCCCCCCGCGGGGGCGTGGATCGAAACCCGCCGCGTGAGTTGCCAGTGATAGACGGCGTGGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTgc                                                                                                            <  2:2286733/101‑1 (MQ=37)


  cTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTgcgc                                                                                                          >  2:1270132/1‑101 (MQ=255)


      ccccGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCACCCGCGGGGGGGTGGATCGAAACTGCgccgcc                                                                                                      >  1:5173274/1‑101 (MQ=255)


      ccccGCGGGGGCCTCGATCGAAACCTGCGGCGGGAGTCGCCAGTGATAGACGCCGTGGGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCgccgcc                                                                                                      <  2:4409119/101‑1 (MQ=37)


       cccGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCt                                                                                                     <  1:1451398/101‑1 (MQ=255)


        ccGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTa                                                                                                    <  1:590738/101‑1 (MQ=255)


         cgcgGGGGCGTGGAGCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTAt                                                                                                   <  2:4252732/101‑1 (MQ=255)


             ggggCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCa                                                                                               >  1:4007822/1‑101 (MQ=255)


              gggCGTTGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCTCTCCCCCCGCGGGGGCGTGGATCGAAGCTGCGCCGCCTATCGCAt                                                                                              <  1:6723444/101‑1 (MQ=255)


              gggCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCAt                                                                                              <  2:714684/101‑1 (MQ=255)


              gggCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGACGCCTATCGCAt                                                                                              <  1:4542491/101‑1 (MQ=255)


               ggCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATc                                                                                             <  2:534464/101‑1 (MQ=255)


                 cGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCg                                                                                           >  1:4431916/1‑101 (MQ=255)


                 cGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCg                                                                                           <  2:5069139/101‑1 (MQ=255)


                   tGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...t                                                                                      >  1:3878181/1‑101 (MQ=255)


                     gATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCg                                                                                    >  1:1954574/1‑101 (MQ=255)


                      aTCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGa                                                                                   <  1:4091486/101‑1 (MQ=255)


                      aTCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGa                                                                                   <  1:6006397/101‑1 (MQ=255)


                       tCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAg                                                                                  >  1:351732/1‑101 (MQ=255)


                         gAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAgcg                                                                                >  2:1450736/1‑101 (MQ=255)


                          aaaCCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAgcgc                                                                               >  2:6606208/1‑101 (MQ=255)


                          aaaCCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAgagc                                                                               >  2:939758/1‑101 (MQ=255)


                           aaCCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGcgc.c                                                                             >  2:3389142/1‑101 (MQ=255)


                             ccTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGcgc.cgc                                                                           >  2:1920034/1‑101 (MQ=255)


                             ccTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGcgc.cgc                                                                           <  2:5394707/101‑1 (MQ=255)


                               tGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGC.Cgcg.c                                                                        <  1:421279/101‑1 (MQ=255)


                                  ggcgTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGC.CGCG.C.A.Gt                                                                   >  2:5723030/1‑101 (MQ=255)


                                   gcgTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...CCGAGCGC.CGCG.C.A.Gtg                                                                  >  2:3603329/1‑100 (MQ=255)


                                           tGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGC.CGCG.C.A.GTCGCTccccc                                                          >  1:1153256/1‑101 (MQ=255)


                                            gCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGC.CGCG.C.A.GTCGCTcccccc                                                         <  2:2963274/101‑1 (MQ=255)


                                             ccAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGC.CGCG.C.A.GTCGCTCCCCCcg                                                        <  2:2349701/101‑1 (MQ=255)


                                             ccAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGC.CGCG.C.A.GTCGCTCCCCCcg                                                        <  2:557807/101‑1 (MQ=255)


                                              cAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGC.CGCG.C.A.GTCGCTCCCCCcgc                                                       >  2:3449845/1‑101 (MQ=255)


                                                    tAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGC.CGCG.C.A.GTCGCTCCCCCCGCTGGGGc                                                 <  1:5637625/101‑1 (MQ=255)


                                                    tAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGC.CGCG.C.A.GTCGCTCCCCCCGCGGGGGc                                                 >  2:2481263/1‑101 (MQ=255)


                                                      gACGCCGTGGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGC.CGCG.C.A.GTCGCTCCCCCCGCGGGGGCGt                                               >  2:1004433/1‑101 (MQ=255)


                                                       aCGCCGTGGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGC.CGCG.C.A.GTCGCTCCCCCCGCGGGGGCGTg                                              <  2:6837355/101‑1 (MQ=255)


                                                        gccaacctatgTCGCTCACCCCGCGGGGGCGTGGATCGAAACTGTGCCGGTTCCCGGAGACGC...TGGGGCCCTCGCG.C...GTCGCTCCCCCCGCGGGGGCGTgg                                             <  2:6316965/91‑1 (MQ=9)


                                                        cGCGGTGGACGTCCCTCCCCCCGTGGGGGCGGGGATCGAACCTGCGCCACCTATCGCATCCGG...TCGAGCGC.CGCG.C.A.GTCGCTCCCCCCGCGGGGGCGTgg                                             <  2:4070035/101‑1 (MQ=38)


                                                         gCCGTGGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGC.CGCG.C.A.GTCGCTCCCCCCGCGGGGGCGTGGa                                            >  1:2419733/1‑101 (MQ=255)


                                                         gCCGTGGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGC.CGCG.C.A.GTCGCTCCCCCCGCGGGGGCGTGGa                                            >  1:1726573/1‑101 (MQ=255)


                                                         gCCGTGGACGTCGCCCCCCCCGCGGGGGCGTGGATCGAGACTGCGCCGCCTATCGCATCCGG...TCGAGCGC.CGCG.C.A.GTCGCTCCCCCCGCGGGGGCGTGGa                                            <  1:5085521/101‑1 (MQ=255)


                                                         ccaacctatgTCGCTCACCCCGCGGGGGCGTGGATCGAAACTGTGCCGGTTCCCGGAGACGC...TGGGGC.C.CTCG.C.ACGTCGCTCCCCCCGCGGGGGCGTGGa                                            <  1:966489/92‑1 (MQ=9)


                                                           gcccggcGTCGCTCCCCCCGCGGGGTCGTGGATCGAAACCCCGCGGCGTCCCGAAGACGGTCTTCTGGTGC.CGCT.G.C.GTCGCTCCCCCCGCGGGGGCGTgg                                             <  2:1118135/95‑1 (MQ=2)


                                                           gcccggcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCCGCGGCGTCCCGAAGACGGTCTTCTGGTGC.CGCT.G.C.GTCGCTCCCCCCGCGGGGGCGTgg                                             <  2:4712033/95‑1 (MQ=2)


                                                           aacctatgTCGCTCACCCCGCGGGGGCGTGGATCGAAACTGTGCCGGTTCCCGGAGACGC...TGGGGC.C.CTCG.C.ACGTCGCTCCCCCCGCGGGGGCGTGGATc                                          >  2:1563305/8‑101 (MQ=9)


                                                            gTGGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGC.CGCG.C.A.GTCGCTCCCCCCGCGGGGGCGTGGATCg                                         >  2:4813314/1‑101 (MQ=255)


                                                            gTGGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGC.CCCG.C.A.GTCGCTCCCCCCGCGGGGGCGTGGATCg                                         >  1:3635253/1‑101 (MQ=255)


                                                             ccggcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCCGCGGCGTCCCGAAGACGGTCTTCTGGTGC.CGCT.G.C.GTCGCTCCCCCCGCGGGGGCGTGGAt                                           <  2:6999633/97‑1 (MQ=2)


                                                             ccggcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCCGCGGCGTCCCGAAGACGGTCTTCTGGTGC.CGCT.G.C.GTCGCTCCCCCCGCGGGGGCGTGGAt                                           >  2:5117400/5‑101 (MQ=2)


                                                              ggACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGC.CGCG.C.A.GTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                       <  2:3868151/101‑1 (MQ=255)


                                                              ctatgTCGCTCACCCCGCGGGGGCGTGGATCGAAACTGTGCCGGTTCCCGGAGACGC...TGGGGC.C.CTCG.C.ACGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                       >  2:1366227/5‑101 (MQ=9)


                                                              cggcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCCGCGGCGTCCCGAAGACGGTCTTCTGGTGC.CGCT.G.C.GTCGCTCCCCCCGCGGGGGCGTGGATc                                          <  1:4892822/98‑1 (MQ=2)


                                                               tatgTCGCTCACCCCGCGGGGGCGTGGATCGAAACTGTGCCGGTTCCCGGAGACGC...TGGGGC.C.CTCGCC.C.GTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                      >  1:1899076/4‑101 (MQ=2)


                                                               ggcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCCGCGGCGTCCCGAAGACGGTCTTCTGGTGC.CGCT.G.C.GTCGCTCCCCCCGCGGGGGCGTGGATCg                                         >  1:1100936/3‑101 (MQ=2)


                                                                 tgTCGCTCACCCCGCGGGGGCGTGGATCGAAACTGTGCCGGTTCCCGGAGACGC...TGGGGC.C.CTCG.C.ACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAca                                    >  2:5975960/2‑100 (MQ=9)


                                                                 tgTCGCTCACCCCGCGGGGGCGTGGATCGAAACTGTGCCGGTTCCCGGAGACGC...TGGGGC.C.CTCG.C.ACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAca                                    >  1:4688879/2‑100 (MQ=9)


                                                                 cGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGC.CGCG.C.A.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACg                                    <  1:1287205/101‑1 (MQ=255)


                                                                 cGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGC.CGCG.C.A.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACg                                    >  2:3250990/1‑101 (MQ=255)


                                                                  gTCGCTCCCCTCGCGGGGGCGTGGATCGAAACCCCGCGGCGTCCCGAAGACGGTCTTCTGGTGC.CGCT.G.C.GTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                      <  2:4146182/101‑1 (MQ=2)


                                                                  gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGC.CGCG.C.A.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGa                                   >  1:3103248/1‑101 (MQ=255)


                                                                   tCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGC.CGCG.C.A.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGAt                                  >  2:2632018/1‑101 (MQ=255)


                                                                   tCGCGCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGC.CGCG.C.A.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGAt                                  <  1:6178001/101‑1 (MQ=255)


                                                                    cGCTCACCCCGCGGGGGCGTGGATCGAAACTGTGCCGGTTCCCGGAGACGC...TGGGGC.C.CTCG.C.ACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATg                                 >  1:3744026/1‑101 (MQ=9)


                                                                    cGCCCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGC.CGCG.C.A.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATg                                 <  2:1659801/101‑1 (MQ=255)


                                                                     gCTCCCCCCGGGGGGGCGGGGATCGAAACGGCGCCGCCTATCGCATCCGG...TCGAGCGC.CGCG.C.A.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGc                                <  2:1792732/101‑1 (MQ=255)


                                                                      ctcgcgccGCGGGGGCGTGGATCGAAACTGTGCCGGTTCCCGGAGACGC...TGGGGC.C.CTCG.CTC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCg                               <  2:5855268/95‑1 (MQ=2)


                                                                      cTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGC.CGCG.C.A.GTCGCTCCCCCCCCGGGGGCGCGGATCGAAACGAGGCg                               >  1:1793627/1‑101 (MQ=255)


                                                                       tcaccccGCGGGGGCGTGGATCGAAACTGTGCCGGTTCCCGGAGACGC...TGGGGC.C.CTCG.C.ACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCgg                              <  2:1050617/98‑2 (MQ=9)


                                                                       tcaccccGCGGGGGCGTGGATCGAAACTGTGCCGGTTCCCGGAGACGC...TGGGGC.C.CTCG.C.ACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCgg                              <  1:322104/98‑2 (MQ=9)


                                                                       tCCCCCCGCGGGGGCGTGGATCGAAACTGTGCCGGTTCCCGGAGACGC...TGGGGC.C.CTCG.C.ACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCgg                              <  1:2533538/101‑2 (MQ=9)


                                                                       tCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGC.CGCG.C.A.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGt                              >  1:4816346/1‑101 (MQ=255)


                                                                         cccccGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGC.CGCG.C.A.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTc                            >  1:3796049/1‑101 (MQ=255)


                                                                         cccccGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGC.CGCG.C.A.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTc                            <  2:2067123/101‑1 (MQ=255)


                                                                         cccccGCGGGGGCGTGGATCGAAACTGCGCCCCCTATCGCATCCGG...TCGAGCGC.CGCG.C.A.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTc                            <  1:3309835/101‑1 (MQ=255)


                                                                          ccccGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGC.CGCG.C.A.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCa                           <  2:6411699/101‑1 (MQ=255)


                                                                            ccGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGC.CGCG.C.A.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATg                         >  1:1623534/1‑101 (MQ=255)


                                                                            ccGCGGGGGCGTGGATCGAAACCCCGCGGCGTCCCGAAGACGGTCTTCTGGTGC.CGCT.G.C.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcaccccgaac                            <  2:6021689/101‑10 (MQ=2)


                                                                             cgcgGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGC.CGCG.C.A.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATgg                        <  2:3635104/101‑1 (MQ=255)


                                                                              gcgGGGGCGTGGATCGAAACTGTGCCGGTTCCCGGAGACGC...TGGGGC.C.CTCG.C.ACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCggcgcgggc                       <  1:7043461/101‑9 (MQ=2)


                                                                               cgGGGGCGTGGATCGAAACTGTGCCGGTTCCCGGAGACGC...TGGGGC.C.CTCG.C.ACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCggcgcgggcc                      >  1:1394390/1‑92 (MQ=2)


                                                                                gggggCGTGGATCGAAACTGTGCCGGTTCCCGGAGACGC...TGGGGC.C.CTCG.C.ACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCggcgcgggcct                     <  2:6124289/101‑11 (MQ=2)


                                                                                gggggCGTGGATCGAAACTGTGCCGGTTCCCGGAGACGC...TGGGGC.C.CTCG.C.ACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAAAATGCggcgcgggcct                     >  2:4569629/1‑91 (MQ=2)


                                                                                gggggCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGC.CGCG.C.A.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCt                     >  1:5411603/1‑101 (MQ=255)


                                                                                 ggggCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGC.CGCG.C.A.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTa                    >  1:5566782/1‑101 (MQ=255)


                                                                                  gggCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGC.CGCG.C.A.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTAc                   <  1:2713429/101‑1 (MQ=255)


                                                                                  gggCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGC.CGCG.C.A.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTAc                   <  1:6521361/101‑1 (MQ=255)


                                                                                  gggCGCTGATCGAACCTGCCCCGCCTGCCGCATCCGG...TCGAGCGC.CGCG.C.A.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCttc                   <  2:6725519/101‑3 (MQ=25)


                                                                                   ggCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGC.CGCG.C.A.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACg                  <  1:344881/101‑1 (MQ=255)


                                                                                    gCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGC.CGCG.C.A.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGt                 >  1:2930011/1‑101 (MQ=255)


                                                                                    gCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGC.CGCG.C.A.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGt                 <  1:1002962/101‑1 (MQ=255)


                                                                                     cGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGC.CGCG.C.A.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGCTGCGTTCATGGCCTACGTc                >  1:6844405/1‑101 (MQ=255)


                                                                                      gTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGC.CGCG.C.A.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCg               >  2:2677663/1‑101 (MQ=255)


                                                                                       tGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGC.CGCG.C.A.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCgg              >  1:2071116/1‑101 (MQ=255)


                                                                                       tGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGC.CGCG.C.A.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCgg              >  1:5784231/1‑101 (MQ=255)


                                                                                        ggATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGC.CGCG.C.A.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCggg             >  1:6454286/1‑101 (MQ=255)


                                                                                              aaaCTGCTCCGCCTATCGCATCCGG...TCGAGCGC.CGCG.C.A.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTgg       <  2:5621019/101‑1 (MQ=255)


                                                                                                 cTGCGCCGCCTATCGCATCCGG...TCGAGCGC.CGCG.C.A.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGCCt    >  2:788534/1‑101 (MQ=255)


                                                                                                  tGCGCCGCCTATCGCATCCGG...TCGAGCGC.CGCG.C.A.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTg   >  2:5217067/1‑101 (MQ=255)


                                                                                                   gcgcCGCCTATCGCATCCGG...TCGAGCGC.CGCG.C.A.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTgg  >  2:4005149/1‑101 (MQ=255)


CGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGG...TCGAGCGC.CGCG.C.A.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGG  >  NC_014034/1395766‑1395965

Base quality scores:


ATCG

< 3 ≤


ATCG

< 20 ≤


ATCG

< 27 ≤


ATCG

< 33 ≤


ATCG

< 40 ≤


ATCG