BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_014034	1,867,833	+C	89.2%	intergenic (‑235/+80)	RCAP_rcc01729 ← / ← RCAP_rcc01730	pyridine nucleotide‑disulfide oxidoreductase/U32 family peptidase

Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_014034	1,867,833	1	.→C	89.2%	27.2 ‑4.4 ‑0.6	93	intergenic (‑235/+80)	RCAP_rcc01729/RCAP_rcc01730	pyridine nucleotide‑disulfide oxidoreductase/U32 family peptidase
Reads supporting (aligned to +/- strand): new base (34/49): ref base (0/10): total (34/59)
Fisher's exact test strand distribution p-value = 1.20e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.54e-01

Only 100 of 161 total aligned reads displayed.


GAAAAGCACTGCTTTTCCCCGCCCGTGTCGCGACAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAGCGCCCAGCCCCGGCGGCGCGAAA.G.C.G.CCGCCATGGGCGGGGCGGGCG.CTGGCC.G.A.GGGCCTTTGGGGCCCCCGGCGGTCAGGGGGATCTGTTGCGCATGGCCGCGTCAGCGCCGCGGCCA  >  NC_014034/1867733‑1867927


gAAAAGCACTGCTTTTCCCCGCCCGTGTCGCGACAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAAAGCCCAGCCCCCGCGGCGCGGAA.a.c.g                                                                                                     >  1:3874400/1‑101 (MQ=38)


  aaaGCACTGCTTTTCCCCGCCCGTGTCGCGACAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAGCGCCCAGCCCCGGCGGCGCGAAA.G.c.g.cc                                                                                                  <  2:867126/101‑1 (MQ=255)


  aaaGCACTGCTTTTCCCCGCCCGTGTCGCGACAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAGCGCCCAGCCCCGGCGGCGCGAAA.G.c.g.cc                                                                                                  <  2:1152811/101‑1 (MQ=255)


   aaGCACTGCTTTTCCCCGCCCGTGTCGCGACAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAACGCACCGCCCCCGCGGCGCGAAA.a.c.g.cgc                                                                                                 >  2:2956909/1‑99 (MQ=37)


      cACTGCTTTTCCCCGCCCGTGTCGCGACAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAACGCCCAGCCCCCGCGGCGCGGAA.A.c.g.cccgcc                                                                                              >  1:459391/1‑97 (MQ=37)


         tGCTTTCCCCCGCCCGTGTCGCGACAGCGTTGCAAAGCCCGCACCCCCAAGCGCGCGAGGGGCCAGCGCCCACCCCCGGCGGCGCGAAA.G.C.GCCCGCCATg                                                                                            <  2:3416217/101‑1 (MQ=25)


          gCTTTCCCCCGCCCGTGTCGCGACAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAGCGCCCAGCCCCGGCGGCGCGAAA.G.C.GCCCGCCATgg                                                                                           <  2:712709/101‑1 (MQ=255)


           cTTTTCCCCGCCGGTGTCGCGACAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAACGCCCCGCCCCCGCGGCGCGGAA.G.c.g.cccgccccggg                                                                                         >  1:2423619/1‑92 (MQ=255)


            ttttCCCCGCCCGTGTCGCGACAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAGCGCCCAGCCCCGGCGGCGCGAAA.G.C.GCCCGCCATgggc                                                                                         <  2:1872053/101‑1 (MQ=255)


               tCCCCGCCCGTGTCGCGACAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAGCGCCCAGCCCCGGCGGCGCGAAA.G.C.GCCCGCCATgggcggg                                                                                      <  2:4382347/101‑1 (MQ=255)


                ccccGCCGGTGTCGCGACAGCGTTGCAAAGACCGCATCCCCAAGCCCGCGAGAGGCAAGCGCCCCGCCCCGGCGGCGCGAAA.G.C.GCCCGCCATgggcgggg                                                                                     <  2:2896259/101‑1 (MQ=18)


                ccccGCCCGTGTCGCGACAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAACGCCCAGCCCCGGCGGCGCGAAA.G.C.GCCCGCCCTgggccggg                                                                                     >  2:1122110/1‑101 (MQ=25)


                 cccgcccgTGTCGCGACAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAACGCCCAGCCCCGGCGGCGCGAAA.G.C.GCCCGCCCTggggggggc                                                                                    >  1:399204/1‑101 (MQ=21)


                 cccgcccgTGTCGCGACAGCGTTGCAAAGACCGCACCCCCAAGCGCGCGAGAGGCCAGCGCCCAGCCCCGGCGGCGCGAAA.G.C.GCCCGCCATgggcggggc                                                                                    <  1:4054189/101‑1 (MQ=255)


                   cgcccgTGTCGCGACAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAACGCCCAGCCCCGGGGGCGCGCAA.G.C.GACCGCCCTGggcggggcgg                                                                                  >  2:2693265/1‑101 (MQ=25)


                    gcccgTGTCGCGACAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAGCGCCCAGCCCCCGCGGGGCGAAA.G.A.GCCCGCCACGGggggggcggg                                                                                 >  2:1124554/1‑101 (MQ=25)


                    gcccgTGTCGCGACAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAACGCCCCGCCCCGGCGGCGCGGAA.G.C.GCCCGCCCCGGgcggggcggg                                                                                 >  2:2920905/1‑101 (MQ=18)


                    ccccgTGTCGCGACAGCGTTGAAAAGCCCCCACCCCCAGGCGCGCGAGAGGCCAGCGCCCAGCCCCGGCGGCGCGAAA.G.C.GCCCGCCATGGgcggggcggg                                                                                 <  2:1577350/100‑1 (MQ=18)


                     cccgTGTCGCGACAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAGCGCCCAGCCCCGGCGGCGCGAAA.G.C.GCCCGCCATGGGCGgggcgggc                                                                                <  2:3978820/101‑1 (MQ=255)


                      cgttgtCGCGACAGCGTGCCAAAGCCCGCACCCCCAAGCGCGCGAGAGGCCAGCGCCCAGCCCCGGGGGCGCGAAA.G.C.GCCCGCCATGGGCGGgtcgggcg                                                                               <  2:130439/98‑1 (MQ=17)


                      ccgTGTCGCGACAGCGTTGCAAAGACCGCATCCCCAATCGCGGGAGAGGCCAGCGCCCAGCCCCGGCGGCGCGAAA.G.C.GCCCGCCATGGGCGGggcgggcg                                                                               <  2:119901/101‑1 (MQ=18)


                      ccgTGTCGCGACAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAGCCCCCAGCCCCGGCGGCGCGAAA.G.C.GCCCGCCATGGGCGGggcgggcg                                                                               <  2:2112957/101‑1 (MQ=255)


                      ccgTGTCGCGACAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAACGCCCCGCCCCCGCGGCGCGAAA.G.C.GCCCGCCCTCGgcgggggggggg                                                                               >  2:4154080/1‑99 (MQ=18)


                       cgTGTCGCGACGGCGTTGCAAAGCCCGCATCCCCAAGCGCGCGAGAGGCCAGCGCCCAGCCCCGGCGGCGCGAAA.G.C.CCCCGCCATGGGCGGGGCGGgcg.c                                                                             <  1:4309788/101‑1 (MQ=18)


                       cgTGTCGCGACAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAGCGCCCACCCCCCGCCGCGCGAAA.G.C.GCCCCCCCAGGGGGGGGGGGgcg.c                                                                             >  1:1514467/1‑101 (MQ=25)


                       cgTGTCGCGACAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAACGCCCAGCCCCGGCGGCGCGCAA.G.C.GCCCGCCCCGGGGGGGGGGGggg.c                                                                             >  1:2062653/1‑101 (MQ=25)


                       cgTGTCGCGACAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAACGCCCACCCCCGGCGGCGCGAAA.A.C.GCCCGCCCTGGGCGGGGGGGggg.c                                                                             >  1:66619/1‑101 (MQ=21)


                        gtgtCGCGACAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAACGCCCCGGCCCGGCGGCGCGGAA.G.CGG.CCGCCCTGGgcggggcggggggg                                                                             >  2:4199314/1‑97 (MQ=18)


                        gtgtCGCGACAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAACGCCCAGCCCCGGGGGCGCGAAA.A.C.GCCCGCCACGGGCGGGGCGGggg.cc                                                                            >  1:2831902/1‑100 (MQ=18)


                         tgtCGCGACAGCGTTGCAAAGCCCGCATCCCCAAGCGCGCGAGAGGCCAGCCCCCAGCCCCGGCGGCGCGAAA.G.C.GCCCGCCATGGGCGGGGCGGGCG.CTg                                                                           <  2:3384249/101‑1 (MQ=18)


                         tgccgcgACAGCGTTGCAAAGCCCCCATCCCCAAGCGCGCGAGAGGCCAGCGCCCAGCCCCGGCGGCGCGAAA.G.C.GCCCGCCATGGGCGGGGCGGGCG.CTg                                                                           <  2:2543825/98‑1 (MQ=18)


                           tCGCGACAGCGTTGCAAAGCCCGCACCCCCAAGCGCGCGAGAGGCCAGCGCCCAGCCCCGGCGGCGCGAAA.G.C.GCCCGCCATGGGCGGGGCGGGCG.CTGGc                                                                         <  2:3819974/101‑1 (MQ=18)


                            cgcgACAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAGCCCCCAGCCCCGGCGGCGCGAAA.G.C.GCCCGCCATGGGCGGGGCGGGCG.CTGGcc                                                                        <  1:3555853/101‑1 (MQ=255)


                            cgcgACAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAGAGCCCAGCCCCGGCGGCGCGAAA.A.C.GCCCGCCCTGGGCGGGGCGGGCG.CCgggc                                                                        >  2:3294023/1‑99 (MQ=18)


                              cgACAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAGCGCCCAGCCCCGGCGTCGCGAAA.G.C.GCCCGCCATGGGCGGGGCGGGCG.CCGGCC.G.a                                                                    >  1:1349332/1‑101 (MQ=18)


                                aCAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAACGCCCAGCCCCGGCGGCGCGAAA.G.C.GCCCGCCATGGGCGGGGCGGGGG.CCGGGC.G.G.gg                                                                 >  1:949847/1‑101 (MQ=18)


                                 cAGGGTTGCAAAGACCGTATCCCCAACCGCGCGAGAGGCCACCGCCCAGCCCCGGCGGCGCGAAA.G.C.GCCCGCCATGGGCGGGGCGGGCG.CTGGCC.G.A.ggg                                                                <  1:3694959/101‑1 (MQ=17)


                                 cAGCGTTGCAAAGACCGCACCCCCAAGCGCGCGAGAGGCCAGCGCCCAGCCCCGGCGGCGCGAAA.G.C.GCCCGCCATGGGCGGGGCGGGCG.CTGGCC.G.A.ggg                                                                <  2:321140/101‑1 (MQ=255)


                                    cGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAGCGCCCAGCCCCGGCGGCGCGAAA.G.C.GCCCCCCCTGGGCGGGGCGGGGG.CCGGGC.G.G.ggggcc                                                             >  1:270516/1‑98 (MQ=17)


                                     gtggCAAAGCCCGCATCCCCAAGCGGGGGGAAGGCAAGCGCCCAGCCCCGGCGGCCCGAAA.G.C.GCCCGCCATGGGCGGGGCGGGCG.CTGGCC.G.A.GGGCCtt                                                            <  1:4189876/98‑1 (MQ=14)


                                     gTTGCACAGCCCGCATCCCCAAGCGCGCGAGAGCCCAGCGCCCAGCCCCGGCGGCGCGAAA.G.C.GCCCGCCATGGGCGGGGCGGGCG.CTGGCC.G.A.GGGCCtt                                                            <  2:535185/101‑1 (MQ=17)


                                     gTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAGCGCCCAGCCCCGGCGGCGCGAAA.G.C.GCCCGCCCTGGGCGGGGCGGGCG.CCGGCC.G.A.GGGccgt                                                            >  1:2678557/1‑99 (MQ=18)


                                       tGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAACGCCCAGCCCCCGCGGCGCGAAA.A.C.GCCCGCCCATGGCGGGGGGGGGG.CCTGCC.G.ACGGGCCttt                                                           >  2:2803609/1‑101 (MQ=12)


                                        gCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAGCGCCCAGCCCCGGCGGCGCGAAA.G.G.GCCCCCCATGGGCGGGGCGGGGG.CGCGCC.G.C.GGGcccctgg                                                         >  2:4089849/1‑96 (MQ=17)


                                        gCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAACGCCCAGCCCCCGCGGCGCGAAA.G.C.GCCCCCCCTGGGCGGGGGGGGGG.CGGGGC.GCG.GGGCCTTTg                                                          >  1:52126/1‑101 (MQ=12)


                                         cAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAGCGCCCAGCCCCCGCGGCGCGAAA.C.A.GCCCGCCATGGGCGGGGGGGGGG.CTGGCC.G.A.GGGCCCCTggg                                                        >  2:307907/1‑101 (MQ=17)


                                            aGACCGCATCCCCAAGCGCGCGAGAGGCCAACGCCCAGCCCCGGCGGCGCGAAA.G.CGG.CCGCCCTGGGCGGGGCGGGGGCCGGGCG.G.G.GGGCCTTTGGGGc                                                      >  2:1568199/1‑101 (MQ=17)


                                             gcccGCACCCCCAAGCGCGCGAGAGGCCAGCGCCCAGCCCGGGCGGCGCGAAA.GCC.C.CCGCCATGGGCGGGGCGGGCG.CTGGCC.G.A.GGGCCTTTGGGGccc                                                    <  1:3360931/99‑1 (MQ=17)


                                             gACCGCATCCCCAAGCGCGCGAGAGGCCAACGCCCAGCCCCCGCGGCGCGAAA.G.C.GCCCGCCATGGGCGGGGCGGGCGCCGGGCC.C.G.GGGCCTTTGGGGcc                                                     >  2:3668871/1‑101 (MQ=17)


                                              cccGCATCCCCAAGCGCGCGAGAGGCCGGCGCCCAGCCCCGGCGGCGCGAAA.G.C.GCCCGCCATGGGCGGGGCGGGCG.CTGGCC.G.A.GGGCCTTTGGGGcccc                                                   <  2:399165/100‑1 (MQ=255)


                                               cctcacccccATGCGCGCGAGAGCCCCGCCCCCACCCCGGGCGTCGCGAAA.GCC.C.CCGCCATGGGCGGGGCGGGCG.CTGGCC.G.A.GGGCCTTTGGGGccccc                                                  <  2:3127838/95‑1 (MQ=14)


                                               ccccacccccAGGGGCGGGGGAGGCCGGCCCCCAGCCCGGGCGGCGCGAAACC.C.C.CCGCCATGGGCGGGGCGGGCG.CTGGCC.G.A.GGGCCTTTGGGGccccc                                                  <  1:906599/95‑1 (MQ=9)


                                                 ccacccccAAGCGCGCGAGAGGCCAGCGCCCAGCCCCGGCGGCGCGAAA.G.C.GCCCGCCATGGGCGGGGCGGGCG.CTGGCC.G.A.GGGCCTTTGGGGCCCCcgg                                                <  2:1172855/97‑1 (MQ=255)


                                                  cATCCCCAAGCGCGCGAGAGGCCAGCGCCCAGCCCCGGGGGCGCGAAA.G.C.GCCCGCCATGGGCGGGGCGGGCG.CTGGCC.G.A.GGGCCTTTGGGGCCCCcggc                                               <  1:562002/101‑1 (MQ=255)


                                                   cccccggaGCGGGCGAGAGCCCGGCCCCCCGCCCCGGCGGCGCGAAA.GCC.C.CCGCCATGGGCGGGGCGGGCG.CTGGCC.G.A.GGGCCTTTGGGGCCCCcggcg                                              <  1:111301/94‑1 (MQ=11)


                                                     ccccAAGCGCGCGAGAGGCCGGCGCCCAGCGCGGGCGGCTTGAAA.G.C.GCCCGCCATGGGCGGGGCGGGCG.CTGGCC.G.A.GGGCCTTTGGGGCCCCCGGCGGt                                            <  2:2868735/101‑1 (MQ=17)


                                                     ccccAAGCGCGCGAGAGGCCGGCCCCCAGCCCCGGCGGCGCGAAA.GCC.C.CCGCCATGGGCGGGGCGGGCG.CTGGCC.G.A.GGGCCTTTGGGGCCCCCGGCGGt                                            <  1:1394382/101‑1 (MQ=17)


                                                     ccccAAGCGCGCGAGAGGCCAGCGCCCAGCCCCGGCGGCGCGAAA.G.C.GCCCGCCATGGGCGGGGCGGGCG.CCGGGC.G.G.GGGCCCTTGGGGCCCCcggcggg                                            >  1:299475/1‑100 (MQ=17)


                                                       ccAGGCGCGCGAGGGCCCGGCCCCCGGCCCCGGCGGCGAAAAA.GCC.C.CCGCCATGGGCGGGGCGGGCG.CTGGCC.G.A.GGGCCTTTGGGGCCCCCGGCGGTCa                                          <  1:1627128/101‑1 (MQ=14)


                                                       ccAAGCGCGGGAGGGGCCAGCGCCCAGCCCCGGCGGCGCGAAA.G.C.GCCCGCCATGGGCGGGGCGGGCG.CTGGCC.G.A.GGGCCTTTGGGGCCCCCGGCGGTCa                                          <  1:1306825/101‑1 (MQ=14)


                                                       ccAAGCGCGCGGGAGGCCTGCGTCCCGCCCCGGCGGCGCGAAT.G.C.GCCCGCCATGGGCGGGGCGGGCG.CCGGCC.G.G.GGGGCTCCGGGTccctggggggccc                                          >  1:2342450/1‑91 (MQ=17)


                                                       ccAAGCGCGCGGGAAGCGAGCCCCCAGCCCCGGAGGCGCGAAA.GCC.C.CCGCCATGGGCGGGGCGGGCG.CTGGCC.G.A.GGGCCTTTGGGGCCCCCGGCGGTCa                                          <  1:4421062/101‑1 (MQ=17)


                                                       ccAAGCGCGCGAGGGGCCGGCGCCCAGCCCCGGCGGCGCGAAA.G.C.GCCCGCCATGGGCGGGGCGGGCG.CTGGCC.G.A.GGGCCTTTGGGGCCCCCGGCGGTCa                                          <  1:2222683/101‑1 (MQ=17)


                                                        cacgcgcgcGAGAGGCCAGCCCCCAGCCCCGGCGGCGCGAAA.G.C.GCCCGCCATGGGCGGGGCGGGCG.CTGGCC.G.A.GGGCCTTTGGGGCCCCCGGCGGTCAg                                         <  2:66606/98‑1 (MQ=17)


                                                        cAAGCGCGCGAGAGGCCAGCGCCCAGCCCCGGCGGCGCGAAA.G.C.GCCCGCCATGGGCGGGGCGGGCG.CCGGCC.G.A.GGGCCTTTGGGGCCCCCGGCGGTCAg                                         >  1:2996507/1‑101 (MQ=255)


                                                         aaGCGCGCGAGAGGCCAGCGCCCAGCCCCGGCGGCGCGAAA.G.C.GCCCGCCATGGGCGGGGCGGGGG.CCGGCC.G.G.GGGCCCTTGGGGGCCCccgcgggccgg                                        >  1:3240162/1‑96 (MQ=14)


                                                         aaGCGCGCGAGAGGCCAGCGCCCAGCCCCGGCGGCGCGAAA.G.C.GCCCGCCATGGGCGGGGCGGGCG.CCGGCC.G.C.GGGGCCTTTGGGCCCCcgggggggcgg                                        >  2:1406344/1‑96 (MQ=14)


                                                         aaGCGCGCGAGAGGCCAGCGCCCAGCCCCGGCGGCGCGAAA.G.C.GCCCGCCATGGGCGGGGCGGGCG.CCGGCC.G.C.GGGGCCTTGGGGCCCCCGGGGGGCAgg                                        >  2:3490078/1‑101 (MQ=17)


                                                             gcgcgAGAGGCCGGCGCCCAGCCCCGGCGGCGCGAAA.GCC.C.CCGCCATGGGCGGGGCGGGCG.CTGGCC.G.A.GGGCCTTTGGGGCCCCCGGCGGTCAGGGGGa                                    <  2:2699499/101‑1 (MQ=17)


                                                             gcgcgAGAGGCCAGCGCCCAGCCCCGGCGGCGCGAAA.G.C.GCCCGCCCTGGGCGGGGCGGGGG.CTGGCC.G.A.GGGCCCTTGGGGGCCCCCGCGGGCCgggggg                                    >  1:4213999/1‑100 (MQ=14)


                                                             cgcggacagGCCAGCACCCGCCCCGGGCGTCGCGAAA.G.C.GCCCGCCATGGGCGGGGCGGGCG.CTGGCC.G.A.GGGCCTTTGGGGCCCCCGGCGGTCAGGGGGa                                    <  2:1349235/94‑1 (MQ=14)


                                                              ggcaagagGCGAGCGCCCGGCCGGGGCGGCGAGAAG.G.C.GCCCGCCATGGGCGGGGCGGGCG.CTGGCC.G.A.GGGCCTTTGGGGCCCCCGGCGGTCAGGGGGAt                                   <  1:2368623/97‑1 (MQ=14)


                                                              cgcgAGAGGCCAGCGCCCAGCCCCGGCGGCGCGAAA.G.C.GCCCGCCATGGGCGGGGCGGGCG.CTGGGC.G.A.GGGCCTTTGGGGCCCCCGGCGGTCAggggggt                                   >  2:4170409/1‑99 (MQ=255)


                                                               gggaggGGCCAGCGCCCAGCCCCGGCGGCGCGAAA.G.C.GCCCGCCATGGGCGGGGCGGGCG.CTGGCC.G.A.GGGCCTTTGGGGCCCCCGGCGGTCAGGGGGATc                                  <  1:1156617/99‑1 (MQ=255)


                                                               gcgAGAGGCCAGCGCCCAGCCCCGGCGGCGCGAAA.G.C.GCCCGCCATGGGCGGGGCGGGCG.CTGGGC.G.G.GGGCCTTTGGGGCCCCCGGCGGGCCGGGGGGTc                                  >  1:1883413/1‑101 (MQ=17)


                                                                cgAGAGGCCAGCGCCCAGCCCCGGCGGCGCGAAA.G.C.GCCCGCCATGGGCGGGGCGGGCG.CTGGCC.G.A.GGGGCTTTGGGGCCCCCCGCGGGCAGGGGGGTCt                                 >  2:198578/1‑101 (MQ=17)


                                                                 gagaGGCCAGCGCCCAGCCCCGGCGGCGCGAAA.G.C.GCCCGCCATGGGGGGGGCGGGCG.CTGGGC.G.A.GGGCCCTTGGGGGCCCCCGCGGGCCAGGGGGTCTg                                >  1:189188/1‑101 (MQ=14)


                                                                 gagaGGCCAGCGCCCAGCCCCGGCGGCGCGAAA.G.C.GCCCGCCATGGGCGGGGCGGGCG.CTGGCC.G.A.GGGCCTTTGGGGGCCCCCGCGGGCAgggggggctg                                >  2:3127535/1‑96 (MQ=14)


                                                                  agagGCCAGCGCCCAGCCCCGGCGGCGCGAAA.G.C.GCCCGCCCTGGGCGGGGCGGGGG.GCGGCCGG.G.GGGCCTTTGGGG.CCCCGGCGGGCCGGGGGGTCTGt                               >  1:21763/1‑101 (MQ=2)


                                                                   gagGCCAGCGCCCAGCCCCGGCGGCGCGAAA.G.C.GCCCGCCATGGGCGGGGCGGGCG.CTGGCC.G.A.GGGCCCTTGGGGCCCCCGGCGGTCAGGGGGATCTGtt                              >  2:3568901/1‑101 (MQ=255)


                                                                   gagGCCAACGCCCAGCCCCGGCGGCGCGAAA.G.C.GCCCGCCCTGGGCCGGGCGGGAG.CTGGCC.C.A.CGGCCGTTGCGGCCCCCGGCGGGCAGGGGGAACTGtt                              >  2:3478536/1‑101 (MQ=14)


                                                                    agGCCAGCGCCCAGCCCCGGCGGCGCGAAA.G.C.GCCCGCCATGGGCGGGGCGGGCG.CTGGCC.G.A.GGGCCTTTGGGGCCCCCGGCGGTCCGGGGGGTCTGTTg                             >  2:3016215/1‑101 (MQ=17)


                                                                        cAGCGCCCAGCCCCGGCGGCGCGAAA.G.C.GCCCGCCATGGGCGGGGCGGGCG.CTGGCC.G.A.GGGCCTTTGGGGCCCCCGGCGGTCAGGGGGATCTGTTGCGCa                         >  1:4192310/1‑101 (MQ=255)


                                                                        cAGCCCCCACCCCCGGCGGCGCGAAG.GCC.C.CCGCCATGGGCGGGGCGGGCG.CTGGCC.G.A.GGGCCTTTGGGGCCCCCGGCGGTCAGGGGGATCTGTTGCGCa                         <  1:4154341/101‑1 (MQ=17)


                                                                         ggcgcCCGGCCCGGGGGGCGAGAAG.G.C.GCCCGCCATGGGCGGGGCGGGCG.CTGGCC.G.A.GGGCCTTTGGGGCCCCCGGCGGTCAGGGGGATCTGTTACGCAt                        <  2:1646133/100‑1 (MQ=14)


                                                                         ggcccccCGCCCCGGCGGCGGGAAA.GCC.C.CCGCCATGGGCGGGGCGGGCG.CTGGCC.G.A.GGGCCTTTGGGGCCCCCGGCGGTCAGGGGGATCTGTTGCGCAt                        <  1:1706524/97‑1 (MQ=17)


                                                                         aGCGCCCAGCCCCGGCGGCGCGAAA.G.C.GCCCGCCATGGGCGGGGCGGGCG.CTGGCC.G.A.GGGCCTTTGGGGCCCCCGGCGGTCAGGGGGGTCTGTTGGGCAt                        >  1:1901358/1‑101 (MQ=17)


                                                                          gcgcCCAGCCCCGGCGGCGCGAAA.G.C.GCCCGCCATGGGCGGGGCGGGCG.CTGGCC.G.A.GGGCCTTTGGGGCCCCCGGCGGTCAGGGGGATCTGTTGCGCATg                       <  1:1414710/101‑1 (MQ=255)


                                                                            gcCCAGCCCCGGCGGCGCGAAA.G.C.GCCCGCCATGGGCGGGGCGGGCG.CTGGCC.G.A.GGGCCTTTGGGGCCCCCGGCGGTCAGGGGGATCTGTTGCGCATGGc                     >  2:2311863/1‑101 (MQ=255)


                                                                             cccccccccGGCGGCGCAAAA.G.C.GCCCGCCATGGGCGGGGCGGGCG.CTGGCC.G.A.GGGCCTTTGGGGCCCCCGGCGGTCAGGGGGATCTGTTGCGCATGGcc                    <  1:3301704/96‑1 (MQ=18)


                                                                               cAGCCCCGGCGGCGCGAAA.G.C.GCCCGCCATGGGCGGGGCGGGCG.CTGGCC.G.A.GGGCCATTGGGGCCCCCCGCGGGCAGGGGGATCTGGTGCGCAGAGCcgc                  >  2:1505817/1‑101 (MQ=17)


                                                                                 cccccGGCGGCGCGAAA.G.C.GCCCGCCATGGGCGGGGCGGGCG.CTGGCC.G.A.GGGCCTTTGGGGCCCCCGGCGGTCAGGGGGATCTGTTGCGCATGGCCGCGt                <  1:1937212/100‑1 (MQ=255)


                                                                                  cccgGGGGGCGCGAAG.C.C.GCCCGCCATGGGCGGGGCGGGCG.CTGGCC.G.A.GGGCCTTTGGGGCCCCCGGCGGTCAGGGGGATCTGTTGCGCATGGCCGCGTc               <  1:3412860/101‑1 (MQ=18)


                                                                                   ccggggggcgCGAAA.GCC.C.CCGCCATGGGCGGGGCGGGCG.CTGGCC.G.A.GGGCCTTTGGGGCCCCCGGCGGTCAGGGGCATCTGTTGCGCATGGCCGCGTCa              <  2:1442844/95‑1 (MQ=18)


                                                                                   cccGGCGGCGCCAAA.GCC.C.CCGCCATGGGCGGGGCGGGCG.CTGGCC.G.A.GGGCCTTTGGGGCCCCCGGCGGTCAGGGGGATCTGTTGCGCATGGCCGCGTCa              <  1:2010630/101‑1 (MQ=21)


                                                                                     cggcggCGCGAAA.G.C.GCCCGCCATGGGCGGGGCGGGCG.CTGGCC.G.A.GGGCCTTTGGGGCCCCCGGCGGTCAGGGGGATCTGTTGCGCATGGCCGCGTCAgc            <  2:3981423/101‑1 (MQ=255)


                                                                                     cggcggCGCGAAA.G.C.GCCCGCCATGGGCGGGGCGGGCG.CTGGCC.G.A.GGGCCTTTGGGGCCCCCGGCGGGCAGGGGGGTCTGTTGCGCATGGCCGCGGCAgc            >  1:3066013/1‑101 (MQ=18)


                                                                                      ggcggcGCGAAA.G.C.GCCCGCCATGGGCGGGGCGGGCG.CTGGCC.G.A.GGGCCTTTGGGGCCCCCGGCGGTCAGGGGGATCTGTTGCGCATGGCCGCGTCAgcg           >  2:438547/1‑101 (MQ=255)


                                                                                      ggcggcGCGAAA.G.C.GCCCGCCATGGGCGGGGCGGGCG.CTGGCC.G.A.GGGCCTTTGGGGCCCCCGGCGGTCAGGGGGATCTGTTGCGCATGGCCGCGTCAgcg           <  1:4163684/101‑1 (MQ=255)


                                                                                               aaa.G.C.GCCCGCCATGGGCGGGGCGGGCG.CTGGCC.G.A.GGGCCTTTGGGGCCCCCGGCGGTCAGGGGGATCTGTTGCGCATGGCCGCGTCAGCGCCGCGGCCa  >  2:2407113/1‑101 (MQ=255)


GAAAAGCACTGCTTTTCCCCGCCCGTGTCGCGACAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAGCGCCCAGCCCCGGCGGCGCGAAA.G.C.G.CCGCCATGGGCGGGGCGGGCG.CTGGCC.G.A.GGGCCTTTGGGGCCCCCGGCGGTCAGGGGGATCTGTTGCGCATGGCCGCGTCAGCGCCGCGGCCA  >  NC_014034/1867733‑1867927

Base quality scores:


ATCG

< 3 ≤


ATCG

< 12 ≤


ATCG

< 20 ≤


ATCG

< 31 ≤


ATCG

< 40 ≤


ATCG