BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_014034	1,395,666	T→C	5.1%	intergenic (+478/‑453)	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,666	0	T→C	5.1%	19.4 ‑2.0 ‑0.2	238	intergenic (+478/‑453)	RCAP_rcc01306/RCAP_rcc01307	Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand): new base (4/8): ref base (134/91): total (138/100)
Fisher's exact test strand distribution p-value = 1.30e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.49e-01

Only 100 of 260 total aligned reads displayed.


CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTC  >  NC_014034/1395566‑1395766


cgtcgCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGTTCGAAACAt                                                                                                          >  1:4076265/1‑101 (MQ=255)


      tCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCAt                                                                                                    >  2:2706700/1‑101 (MQ=255)


      tCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCAt                                                                                                    >  1:2334190/1‑101 (MQ=255)


       ccccccGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATg                                                                                                   <  1:4189995/101‑1 (MQ=255)


       ccccccGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATg                                                                                                   <  2:4209063/101‑1 (MQ=255)


        cccccGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGt                                                                                                  <  1:1019182/101‑1 (MQ=255)


          cccGCGGGGGCGTGGATCGAAACCCCCTTGCGGCTTGTCACTTCCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAg                                                                                                <  2:593523/101‑1 (MQ=255)


          cccGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAg                                                                                                <  1:3178283/101‑1 (MQ=255)


           ccGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGt                                                                                               >  1:3488789/1‑101 (MQ=255)


           ccGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGt                                                                                               <  2:3507559/101‑1 (MQ=255)


           ccGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGAGTCGCTCCCCCCGCGGGGGGGGGGGTCGAAACACCTGCAGTTAGt                                                                                               >  2:3689276/1‑101 (MQ=38)


               gggggCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTcggc                                                                                           <  1:2171506/101‑1 (MQ=255)


               gggggCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTcggc                                                                                           <  2:3168764/101‑1 (MQ=255)


                ggggCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTcggcg                                                                                          <  2:3389475/101‑1 (MQ=255)


                 gggCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTcggcg..g                                                                                       >  1:2272101/1‑101 (MQ=255)


                    cGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCa                                                                                   <  2:2249878/101‑1 (MQ=255)


                      tGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGt                                                                                 <  1:796365/101‑1 (MQ=255)


                      tGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGt                                                                                 >  1:2102832/1‑101 (MQ=255)


                       ggATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAGCATCTGCATGTAGTCCGCG..G.GCAGtt                                                                                >  1:1446933/1‑101 (MQ=255)


                       ggATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGtt                                                                                <  1:1249426/101‑1 (MQ=255)


                        gATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTc                                                                               <  2:3716027/101‑1 (MQ=255)


                        gATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTc                                                                               >  2:2078314/1‑101 (MQ=255)


                        gATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTc                                                                               <  1:3398163/101‑1 (MQ=255)


                         aTCGAAACCACCTTGCGGATTGTCACTTTCGGGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCg                                                                              >  1:4130328/1‑101 (MQ=255)


                         aTCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCg                                                                              >  1:2719858/1‑101 (MQ=255)


                         aTCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCg                                                                              <  1:1759300/101‑1 (MQ=255)


                         aTCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCg                                                                              >  1:1378243/1‑101 (MQ=255)


                         aTCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCg                                                                              >  1:4392892/1‑101 (MQ=255)


                         aTCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCg                                                                              >  2:4313303/1‑101 (MQ=255)


                         aTCGAAACCACCTTGCGGATTGTCACGTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCg                                                                              <  2:2991352/101‑1 (MQ=255)


                          tCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGAGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGCAGACGGCG..G.GCAGTTCGa                                                                             >  2:1339793/1‑101 (MQ=38)


                            gAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.Gc                                                                          >  2:2749207/1‑101 (MQ=255)


                            gAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.Gc                                                                          <  1:2759516/101‑1 (MQ=255)


                            gAAACCACCTTGCGGATTGTCACGTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.Gc                                                                          <  2:655753/101‑1 (MQ=255)


                             aaaCCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCACTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.Gcc                                                                         >  2:943537/1‑101 (MQ=255)


                                ccaccTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCa                                                                      <  2:4232643/101‑1 (MQ=255)


                                 caccTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAg                                                                     >  2:1056063/1‑101 (MQ=255)


                                    cTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCg                                                                  >  2:2829198/1‑101 (MQ=255)


                                    cTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCg                                                                  <  2:3550826/101‑1 (MQ=255)


                                      tGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCt                                                                <  1:4413393/101‑1 (MQ=255)


                                        cGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTcc                                                              >  2:1391325/1‑101 (MQ=255)


                                          gATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTcccc                                                            >  2:4133052/1‑101 (MQ=255)


                                          gATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTcccc                                                            >  1:3339632/1‑101 (MQ=255)


                                          gATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTcccc                                                            >  2:465538/1‑101 (MQ=255)


                                           gttGTCACTTTCGCGGGGCGAAGCGGCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTccccc                                                           <  2:874945/100‑1 (MQ=255)


                                           aTTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTccccc                                                           <  1:1519990/101‑1 (MQ=255)


                                             tGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTAGCTCCCCCcg                                                         >  2:3978852/1‑101 (MQ=255)


                                              gTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCcgc                                                        >  1:1189058/1‑101 (MQ=255)


                                                cACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCgg                                                      >  1:525253/1‑101 (MQ=255)


                                                 ccTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCggg                                                     <  1:1624763/100‑1 (MQ=255)


                                                 aCTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCggg                                                     >  1:4258449/1‑101 (MQ=255)


                                                   tttCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCggggg                                                   >  1:3239532/1‑101 (MQ=255)


                                                    ttCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGc                                                  >  2:3062841/1‑101 (MQ=255)


                                                    ttCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGc                                                  <  1:2187180/101‑1 (MQ=255)


                                                      cgcgGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGt                                                >  1:2213866/1‑101 (MQ=255)


                                                      cgcgGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGt                                                >  1:2566944/1‑101 (MQ=255)


                                                       gcgGGGCGAAGCGTCCCCCCCCCCGCGGGGGCGTGGATCCAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTg                                               <  2:2365053/101‑1 (MQ=255)


                                                        cgGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTgg                                              >  1:1708800/1‑101 (MQ=255)


                                                             cGAAGCGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                         <  1:374188/101‑1 (MQ=255)


                                                              gAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                        >  1:3120485/1‑101 (MQ=255)


                                                               ctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCG..GCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                        <  1:4626361/99‑1 (MQ=2)


                                                               aaGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                       >  1:3595838/1‑101 (MQ=255)


                                                                tgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCG..GCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                       >  2:2043302/2‑101 (MQ=2)


                                                                aGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGagatcgga                                      >  2:4429769/1‑94 (MQ=255)


                                                                aGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                      <  1:280870/101‑1 (MQ=255)


                                                                aGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                      >  2:3471186/1‑101 (MQ=255)


                                                                 gCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCGGCAGTGTCA.GGTTTTCCA.GCGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                      >  2:171963/1‑101 (MQ=2)


                                                                 gCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCG..GCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                      <  2:1202685/101‑1 (MQ=2)


                                                                 ccggccccccccccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACt                                     <  2:1983125/93‑1 (MQ=255)


                                                                  cGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTg                                    >  2:489493/1‑101 (MQ=255)


                                                                   gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTgg                                   >  2:1738643/1‑101 (MQ=255)


                                                                       ccccccccGGGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATg                               <  2:2462964/99‑1 (MQ=255)


                                                                        tCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATgg                              >  2:1477728/1‑101 (MQ=255)


                                                                         ccccccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGa                             >  2:2570285/1‑101 (MQ=255)


                                                                            cccGGGGGGGCGGGGATCGAAACATCTGCGGGTAGTCGGCG..G.GCAGTGGGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCt                          <  2:1140907/101‑1 (MQ=37)


                                                                             ccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTg                         >  1:2046872/1‑101 (MQ=255)


                                                                               gcgGGGGCGTGGATCGAAACATGTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGgc                       >  2:4599516/1‑101 (MQ=255)


                                                                               gcgGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGgc                       <  1:3068577/101‑1 (MQ=255)


                                                                               gcgGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGGGGATCGAAACTGGTAGGGATCTGgc                       >  2:374939/1‑101 (MQ=255)


                                                                                ggggggCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGgcg                      <  2:221146/100‑1 (MQ=255)


                                                                                 gggggCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGGGGATCGAAACTGGTATGGATCTGgcgc                     >  1:908596/1‑101 (MQ=255)


                                                                                  ggggCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCa                    >  2:2883898/1‑101 (MQ=255)


                                                                                    ggCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCaaa                  >  2:2667554/1‑101 (MQ=255)


                                                                                     gCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAAc                 <  1:3781842/101‑1 (MQ=255)


                                                                                      cGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCTCAAACg                >  1:761529/1‑101 (MQ=255)


                                                                                        tGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGa              >  2:3619179/1‑101 (MQ=255)


                                                                                        tGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGa              <  1:353206/101‑1 (MQ=255)


                                                                                         ggATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGaa             >  1:4435724/1‑101 (MQ=255)


                                                                                          gATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGaaa            >  1:196784/1‑101 (MQ=255)


                                                                                          gATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGaaa            <  2:389278/101‑1 (MQ=255)


                                                                                          gATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGaaa            >  2:3274515/1‑101 (MQ=255)


                                                                                           aTCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAg           <  1:3653671/101‑1 (MQ=255)


                                                                                            tCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGt          >  1:259671/1‑101 (MQ=255)


                                                                                              gAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCt        <  1:1056151/101‑1 (MQ=255)


                                                                                                  cATCTGCATGTGGGCGGGG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgtcg    <  1:397770/101‑1 (MQ=255)


                                                                                                  cATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTcagatc    >  1:2410716/1‑96 (MQ=255)


                                                                                                  cATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgtcg    >  1:1927152/1‑101 (MQ=255)


                                                                                                   ctCTGGATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgtcgt   <  2:761454/100‑1 (MQ=255)


                                                                                                    tCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgtcgtc  <  2:4617803/101‑1 (MQ=255)


                                                                                                    tCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgtcgtc  >  2:1550351/1‑101 (MQ=255)


CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTC  >  NC_014034/1395566‑1395766

Base quality scores:


ATCG

< 3 ≤


ATCG

< 18 ≤


ATCG

< 26 ≤


ATCG

< 33 ≤


ATCG

< 40 ≤


ATCG