BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_014034	1,867,831	+C	6.2%	intergenic (‑233/+82)	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,831	1	.→C	6.2%	5.9 ‑2.1 ‑0.0	64	intergenic (‑233/+82)	RCAP_rcc01729/RCAP_rcc01730	pyridine nucleotide‑disulfide oxidoreductase/U32 family peptidase
Reads supporting (aligned to +/- strand): new base (0/4): ref base (28/32): total (28/36)
Fisher's exact test strand distribution p-value = 1.25e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.03e-01

Only 100 of 163 total aligned reads displayed.


AGGAAAAGCACTGCTTTTCCCCGCCCGTGTCGCGACAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAGCGCCCAGCCCCGGCGGCGC.GAAA.G.CG.CCG.CCATGGGCGGGGCGGGCG.CT.GGCCGA.GGGCCTTTGGGGCCCCCGGCGGTCAGGGGGATCTGTTGCGCATGGCCGCGTCAGCGCCGCGGCCAA  >  NC_014034/1867731‑1867928


 ggAAAAGCACTGCTTTTCCCCGCCCGTGTCGCGACAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAGCGCCCAGCCCCGGCGGCGC.GAAA.g.c                                                                                                       <  2:1315535/101‑1 (MQ=255)


 ggAAAAGCACTGCTTTTCCCCGCCCGTGTCGCGACAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAGCGCCCAGCCCCGGCGGCGC.GAAA.g.c                                                                                                       <  1:1242046/101‑1 (MQ=255)


 ggAAAAGCACTGCTTTTCCCCGCCCGTGTCGCGACAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAACGCCCCGCCCCCGCGGCGC.GGAA.g.c                                                                                                       >  1:2281421/1‑101 (MQ=37)


  gAAAAGCACTGCTTTTCCCCGCCCGTGTCGCGACAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAGCGCCCCGCCCCCGCGGGGC.GAAA.g.cg                                                                                                      >  1:4374093/1‑101 (MQ=255)


  gAAAAGCACTGCTTTTCCCCGCCCGTGTCGCGACAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAGCGCCCACCCCCGGCGGCGC.CAAA.a.cg                                                                                                      >  2:1058702/1‑101 (MQ=255)


   aaaaGCACTGCTTTTCCCCGCCCGTGTCGCGACAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAACGCCCCGCCCCCGCGGCGC.GAAA.g.cg.c                                                                                                    >  1:4546400/1‑101 (MQ=255)


    aaaGCACTGCTTTTCCCCGCCCGTGTCGCGACAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAGCGCCCAGCCCCGGCGGCGC.GAAA.G.cg.cc                                                                                                   <  1:754343/101‑1 (MQ=255)


     aaGCCCTGCTTTTCCCCGCCCGTGTCGCGACAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAGCGCCCAGCCCCGGCGGCGC.GAAA.G.cg.ccc                                                                                                  <  2:3056886/101‑2 (MQ=255)


     aaGCACTGCTTTTCCCCGCCCGTGTCGCGACAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAGAGCCCCGCCCCCGCGGCGC.GAAA.G.cg.ccc                                                                                                  >  1:4543477/1‑100 (MQ=255)


     aaGCACTGCTTTTCCCCGCCCGTGTCGCGACAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAACGCCCCGCCCCCGCGGCGC.GCAA.G.ag.ccc                                                                                                  >  2:3686823/1‑100 (MQ=37)


      aGCACTGCTTTTCCCCGCCCGTGTCGCGACAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAACGCCCCGCCCCCGCGGCGC.GAAA.A.cg.cccc                                                                                                 >  1:3230933/1‑99 (MQ=37)


         aCTGCTTTTCCCCGCCCGTGTCGCGACAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAACGCCCAGCCCCCCCGGCGC.CAAA.G.Cg.ccc.cccc                                                                                             >  2:2059460/1‑99 (MQ=37)


          cgccttccccccGCCGGTCTCGCGACGGTTTTGCAAACACCGCATCCCCAAGCGGGAGAGGGGCCAGCGCCCACCCCCGGCGGCGC.GAAA.G.CGCCCG.CCAt                                                                                             <  2:4691494/93‑1 (MQ=12)


          cTGCTTCCCCCCCCCGGTGTCGCGAGAGCGTTGCAAAGCCCGCCCCCCCAGGCGCGCGAGGGGCCGGCCCCCCCCCCCGGCGGCGC.GAAA.GCCC.CCG.CCAt                                                                                             <  2:4270436/101‑1 (MQ=12)


           ttcTTTCCCCCGCCCGTGTCGCGACAGCGTTGCAAAGCCCGCACCCCCAAGCGCGCGAGGGCCCAGCGCCCAGCCCGGGCGGCGC.AAAA.G.CGCCCG.CCATg                                                                                            <  2:1192849/99‑1 (MQ=25)


             cTTTTCCCCGCCCGTGTCGCGACAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAGCGCCCAGCCCCGGCGGCGC.GAAA.A.CGCCCG.CCATggg                                                                                          >  1:2934060/1‑101 (MQ=255)


               tttCCCCGCCCGTGTCGCGACAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAACGCCCCGCCCCCGCGGCGC.GAAA.A.AGCCCG.CCCTgggcg                                                                                        >  2:142137/1‑101 (MQ=21)


               tttCCCCGCCCGTGTCGCGACAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAACGCCCCCCCCCCGCGGCGC.GCAA.A.Cg.ccc.ccccccggcg                                                                                       >  1:1080668/1‑93 (MQ=25)


                  ccccGCCCGTGTCGCGACAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAGCGCCCAGCCCCGGCGGCGC.GAAA.G.CGCCCG.CCATgggcgggg                                                                                     <  2:4366644/101‑1 (MQ=255)


                  ccccGCCCGTGTCGCGACAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAGCGCCCAGCCCCGGCGGCGC.GAAA.G.CGCCCG.CCATgggcgggg                                                                                     <  1:141836/101‑1 (MQ=255)


                  ccccGCCCGTGTCGCGACAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAACGGCCCGCCCCCGCGGCGCGGAAA.C.GG.CCC.CCCCgggcgggg                                                                                     >  1:2525728/1‑101 (MQ=21)


                     cgcccgTGTCGCGACAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAGCGCCCAGCCCCGGCGGCGC.GAAA.G.CGCCCG.CCATGggcggggcgg                                                                                  <  2:4899570/101‑1 (MQ=255)


                      gcccgTGTCGCGACAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAGCGCCCAGCCCCGGCGGCGC.GGAA.G.CGCCCG.CCCAGGgcggggcggg                                                                                 >  2:5051536/1‑101 (MQ=25)


                      gcccgTGTCGCGACAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAAAGCCCAGCCCCCGCGGCGC.GAAA.G.CGCCCG.CCCAGGgcgggggggg                                                                                 >  1:2781908/1‑101 (MQ=21)


                          gtgtCGCGACAGCGTTGAAAAGACCGCATCCCCAGGCGCGCGAGAGGCCAGCGCCCAGCCCCGGCGGCGC.GAAA.G.CGCCCG.CCATGGGCGGGGCGGGCG.Ct                                                                            <  2:671474/101‑1 (MQ=18)


                           tgtCGCGACAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAGCGCCCAGCCCCGGCGGCGC.GAAA.G.CGCCCG.CCATGGGCGGGGCGGGCG.CT.g                                                                          <  2:246599/101‑1 (MQ=255)


                           tgtCGCGACAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAAGGCCCCGCCCCCGCGGCGC.GAAA.A.CGCCCG.CCACGGGCGGGGCCGggg.ccg                                                                           >  2:4605384/1‑99 (MQ=18)


                           tgtCGCGACAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAACGCCCAGCCCCGGCGGGGCGGAAA.C.CG.CCG.CCCAGGGCCGGGCGGggg.ccc                                                                           >  2:2422043/1‑99 (MQ=25)


                           ggcggccagggCGTTGCAAACGCCCCATCCCCAAGCGCGCGAGAGGCCAGCGCCCAGCCCCGGCGCCGC.CAAA.GCCC.CCG.CCATGGGCGGGGCGGGCG.CT.g                                                                          <  2:4971786/91‑1 (MQ=14)


                            gtCGCGACAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAACGGCCACGCCCGGCGGGGC.GAAA.A.CGCCCG.CCGGGGGAGGGGGGGGGG.GT.gg                                                                         >  1:2299784/1‑101 (MQ=255)


                             tCGCGACAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAGCGCCCAGCCCCGGCGGCGC.GAAA.G.CGCCCG.CCATGGGCGGGGCGGGCG.CT.GGc                                                                        <  2:3009565/101‑1 (MQ=255)


                              cgcgACAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAGCGCCCAGCCCCGGCGGCGC.GAAA.G.CGCCCG.CCCTGGGCGGGGCGGGGG.CC.gggc                                                                       >  2:3596226/1‑99 (MQ=18)


                              cgcgACAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAAAGCCCAGCCCCCGCGGCCC.GAAA.G.CGCCCG.CCCTGGGCGGGGGGGGGG.CG.GGcc                                                                       >  2:3272478/1‑101 (MQ=18)


                               gcgACAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAAAGCCCAGCCCCCGCGGCGC.GAAA.G.CGCCCG.CCCTGGGCGGGGGGGGGG.CC.GGCCg                                                                      >  1:3891053/1‑101 (MQ=18)


                               gcgACAGCGTTGCAAAGACCCCATCCCAAAGCGCGCGAGAGGCCAGCGCCCAGCCCCGGCGGCGC.GAAA.G.CGCCCG.CCATGGGCGGGGCGGGCG.CT.GGCCg                                                                      <  2:3131633/101‑1 (MQ=18)


                                cgACAGCTTTCCAAAGCCCCCACCCCCAAGCGCGAGAGAGCCCAGCCCCCAGCCCCGGGGGCGC.GAAA.GCCC.CCG.CCATGGGCGGGGCGGGCG.CT.GGCCGa                                                                     <  1:3166436/101‑1 (MQ=14)


                                cgACAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCGAACGCCCAGCCCCGGCGGCGC.GAAA.G.CGCCCG.CCCTGGGCGGGGCGGGGG.CG.GGCCgg                                                                     >  2:2580107/1‑100 (MQ=18)


                                cgACAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAGCGCCCAGCCCCGGCGGCGC.GAAA.A.CGCCCG.CCCTGGGCGGGGCGGGGG.CC.GGGCgg                                                                     >  1:2049256/1‑100 (MQ=18)


                                 gACAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAGCGCCCAGCCCCGGCGGCGC.GAAA.G.CGCCCG.CCATGGGCGGGGCGGGCG.CT.GGCCGA.g                                                                   <  2:4546080/101‑1 (MQ=255)


                                 gACAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAACGCCCAGCCCCGGCGGCGC.GAAA.G.CGCCCG.CCACGGGCGGGGCGGGGG.CC.GGCCggg                                                                    >  2:178042/1‑99 (MQ=18)


                                 gACAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAACGCCCAGCCCCCGCGGCGC.GAAA.A.CGCCCG.CCCTCGGCGGGGCGGGGG.CC.gggccag                                                                    >  2:1182092/1‑96 (MQ=18)


                                  aCAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAGGGCCCAGCCCCGGCGGCGC.GAAA.G.CGCCCG.CCATGGGCGGGGCGGGCG.CG.GGGCGC.gg                                                                  >  1:4752450/1‑101 (MQ=18)


                                  aCAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAGCCCCCAGCCCCGGCGGCGC.GAAA.G.CGCCCG.CCATGGGCGGGGCGGGCG.CT.GGCCGA.gg                                                                  <  2:451707/101‑1 (MQ=255)


                                  aCAGCGTTGAAAAGACCGCACCCCCAAGCGCGCGAGAGGCCAGCGCCCAGCCCCGGCGGCGC.GAAA.G.CGCCCG.CCATGGGCGGGGCGGGCG.CT.GGCCGA.gg                                                                  <  1:3571973/101‑1 (MQ=17)


                                    agggTTGCAAACCCCCCCCCCCCAGGCGCGCGAGAGGCCGGCCCCCAGCCCCGGCGGCGC.AAAA.GCCC.CCG.CCATGGGCGGGGCGGGCG.CT.GGCCGA.GGGc                                                                <  2:942503/98‑1 (MQ=9)


                                     tctgcggcaggACCGCATCCCCAAGCACGGCAGAGGCCAGCGTCCCGCCCCGGCGGCGC.GAAA.G.CGCCCG.CCATGGGCGGGGCGGGCGCCG.GGCCGG.GGGc                                                                >  1:4734197/11‑101 (MQ=21)


                                     gCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAGCGCCCAGCCCCGGCGGCGC.GAAA.G.CGCCCG.CCACGGGCGGGGCGGGGG.CC.GGCCCA.ggggc                                                               >  1:5022/1‑99 (MQ=18)


                                     gCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAACGCCCCGCCCCCGCGGCGC.GAAA.G.CGCCCG.CCCTGGGCGGGGGGGGCG.CC.GGccccggggc                                                                >  1:2021896/1‑94 (MQ=17)


                                      cGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAACGCCCAGCCCCGGCGGCGC.GAAA.A.CG.CCGCCCATGGGCGGGGGGGGGG.CC.GGCCGG.ggggcc                                                              >  1:4301455/1‑98 (MQ=18)


                                        ttGCAAAGACCGCATCCCAAAGCGCGCGAGAGGCCAGCGCCCAGCCCCGGCGGCGC.GAAA.G.CGCCCG.CCATGGGCGGGGCGGGCG.CT.GGCCGA.GGGCCttt                                                            <  2:851386/101‑1 (MQ=255)


                                        gtGCAAGGCCCCCCTCCCCAGGCGGGCGAGGGGCGGGCGCCCGGCCCGGGGGGCCC.GAAA.GCCC.CCG.GCATGGGCGGGGCGGGCG.CT.GGCCGA.GGGCCttt                                                            <  1:1871745/100‑1 (MQ=2)


                                          gCAAAGCCCGCACCCCAAAGCGCGCGAGAGGCCAGCGCCCAGCCCCGGCGGCGC.GAAA.G.CGCCCG.CCATGGGCGGGGCGGGCG.CT.GGCCGA.GGGCCTTTgg                                                          <  2:3913892/101‑1 (MQ=17)


                                          gCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAGCCCCCAGCCCCCGCGGCGC.GAAA.A.CGCCCG.CCATGGGCGGGGGGGGCG.GC.GGCCCACGGGCCttgg                                                           >  2:3416817/1‑99 (MQ=17)


                                          gCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAACGCCCAGCCCCGGCGGCGC.GAAA.A.CGCCCC.CCCTGGGCGGGGCGGGGG.CC.GGCCGC.GGGGCCTTgg                                                          >  1:4021497/1‑101 (MQ=17)


                                          gCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAACGCCCAGCCCCCGCGGCGC.GAAA.G.CGCCCG.CCATGGGCGGGGCGGGCG.CT.GGCCGA.GGGCCctttg                                                          >  2:1668401/1‑99 (MQ=17)


                                          gCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAAAGGCCACCCCCGGCGGGGC.GAAA.A.CGCCCG.CCATGGGCGGGGCGGGGG.CG.GGCCGC.GGGCCttgtg                                                          >  2:2515222/1‑98 (MQ=14)


                                               gcccGCACCCCCAAGCGCGCGAGAGGCCAGCGCCCAGCCCCGGCGGCGC.GAAA.G.CGCCCG.CCATGGGCGGGGCGGGCG.CT.GGCCGA.GGGCCTTTGGGGccc                                                     <  1:1240356/99‑1 (MQ=255)


                                                cccccccccccAGGCCCGGGGGGGGCCGGCCCCCCCCCCCGGCGGCGC.GAAG.G.CGCCCG.CCATGGGCGGGGCGGGCG.CT.GGCCGA.GGGCCTTTGGGGcccc                                                    <  2:1785655/94‑1 (MQ=0)


                                                aCCGCATCCCCAAGCGCGCGAGAGGCCAACGCCCCGCCCCCGCGGCGC.GAAA.G.CGCCCG.CCCTGGGCGGGGCGGGGC.CC.GGCCGC.GGGGCCTTGGGGgccc                                                    >  1:914988/1‑101 (MQ=9)


                                                   ccgtCCCCAAGCGCGCGGGGGGCGGGCGGCGAGCCCCGGCGGCGC.GAAG.G.CGCCCG.CCATGGGCGGGGCGGGCG.CT.GGCCGA.GGGCCTTTGGGGCCCCcgg                                                 <  2:3401403/98‑1 (MQ=14)


                                                   ccacccccAGGCGCGCGAGAGGCCAGCGCCCAGCCCCGGCGGCGC.GAAA.G.CGCCCG.CCATGGGCGGGGCGGGCG.CT.GGCCGA.GGGCCTTTGGGGCCCCcgg                                                 <  1:87440/97‑1 (MQ=17)


                                                    cATCCCCAAGCGCGCGAGAGGCCAGCGCCCAGCCCCGGCGGCGC.GAAA.G.CGCCCG.CCATGGGGGGGGCGGGCG.CTCGGCGCG.GGGCCTTGGGGGccccccc                                                 >  2:2180418/1‑99 (MQ=17)


                                                    cATCCCCAAGCGCGCGAGAGGCCAGCGCCCAGCCCCGGCGGCGC.GAAA.G.CGCCCG.CCATGGGCGGGGCGGGGG.CT.GGCCGA.GGGCCCTTGGGGGCCCccgc                                                >  1:2840026/1‑101 (MQ=17)


                                                     aTCCCCAAGCGCGCGAGAGGCCAGGCCCCCGCCCCGGCGGCGC.GGAA.G.CGCCCG.CCCTGGGCGGGGCGGGCG.CT.GGGCGAGGGGCCTTTGGGGcccctgcg                                                >  2:420537/1‑97 (MQ=14)


                                                        cccAGGCGCGGGGGGGCCCCGCCCCCCCCCCGGGCGGCGC.AAAC.C.CCCCCG.CCATGGGCGGGGCGGGCG.CT.GGCCGA.GGGCCTTTGGGGCCCCCGGCGGTc                                            <  1:4240491/101‑1 (MQ=2)


                                                          cacgcgcgcGAGAGCCCAGCGCCCAGCCCCGGGGGCGC.GAAG.G.CGCCCG.CCATGGGCGGGGCGGGCG.CT.GGCCGA.GGGCCTTTGGGGCCCCCGGCGGTCAg                                          <  1:3637388/98‑1 (MQ=17)


                                                           aaGCGCGCGAGAGGCCAGCGCCCAGCCCCGGCGGCGC.GAAA.G.CGCCCG.CCATGGGCGGGGCGGGCG.CT.GGCCGG.GGGCCCTTGGGGccccccgggggccgg                                         >  1:600910/1‑91 (MQ=17)


                                                            aGCGCGCGAGAGGCCAGCGCCCAGCCCCGGCGGCGC.GAAA.G.CGCCCG.CCCTGGGCGGGGCGGGCG.CC.CGCCGG.GGGGCTTTGGGGGCCCCCGCGGTcgcgg                                        >  1:4657737/1‑97 (MQ=14)


                                                              tccgatctagGCCAGCGCCCAGCCCCGGCGGCGC.GAAA.G.CGCCCG.CCATGGGCGGGGCGGGCG.CT.GGCCGA.GGGCCTTTGGGGCCCCCGGCGGTCAggggg                                      <  2:4377322/93‑1 (MQ=17)


                                                                cgcgAGAGGCCAGCGCCCAGCCCCGGCGGCGC.GAAA.G.CGCCCG.CCATGGGCGGGGCGGGCG.CT.GGCCGA.GGGCCTTTGGGGCCCCCGGCGGTCAGGGGGAt                                    <  2:3051123/101‑1 (MQ=255)


                                                                cgcgAGAGGCCAGCCCCCCGCCCCGGCGGCGC.GAAA.GCCC.CCG.CCATGGGCGGGGCGGGCG.CT.GGCCGA.GGGCCTTTGGGGCCCCCGGCGGTCAGGGGGAt                                    <  1:70717/101‑1 (MQ=14)


                                                                cccgAGGGCCCGGCCCCCGGCCCCGGCGGCGG.AAAA.G.CGCCCG.CCATGGGCGGGGCGGGCG.CT.GGCCGA.GGGCCTTTGGGGCCCCCGGCGGTCAGGGGGAt                                    <  1:3932548/99‑1 (MQ=14)


                                                                  cgAGAGGCCAGCGCCCAGCCCCGGCGGCGC.GAAA.G.CGCCCG.CCATGGGCGGGGCGGGCG.CT.GGCCGA.GGGCCCTTGGGGCCCCCGGCGGGCAGGGGGGTCt                                  >  2:1372610/1‑101 (MQ=14)


                                                                  cgAGAGGCCAGCGCCCAGCCCCGGCGGCGC.GAAA.G.CGCCCG.CCATGGGCGGGGCGGGCG.CC.GGCCGA.GGGCCTTTGGGGCCCCCCGCGGGCCgggggggct                                  >  2:2685671/1‑97 (MQ=14)


                                                                   gagaGGCCAGCGCCCAGCCCCGGCGGCGC.GAAA.G.CGCCCG.CCCTGGGCGGGGCGGGCG.CT.GGCCGC.GGGCCTTTGGGGCCCCCGGCGGGCAGGGGGGTCTg                                 >  1:4200571/1‑101 (MQ=17)


                                                                   gagaGCCCGGCGCCCAGCCCCGGCGGCGC.GAAA.GCCC.CCG.CCATGGGCGGGGCGGGCG.CT.GGCCGA.GGGCCTTTGGGGCCCCCGGCGGTCAGGGGGATCTg                                 <  2:610406/101‑1 (MQ=17)


                                                                     ggggCCGGCGCCCGGCCCGGGCGGCGC.GAAA.G.CGCCCG.CCATGGGCGGGGCGGGCG.CT.GGCCGA.GGGCCTTTGGGGCCCCCGGCGGTCAGGGGGATCTGtt                               <  2:1876301/99‑1 (MQ=17)


                                                                     gagGCCAGCGCCCAGCCCCGGCGGCGC.GAAA.G.CGCCCG.CCATGGGCGGGGCGGGCG.CT.GGCCGA.GGGGCTTTGGGGCCCCCGGCGGGCAGGGGGGGCTGtt                               >  1:1015198/1‑101 (MQ=17)


                                                                     gagGCCAGCGCCCAGCCCCGGCGGCGC.GAAA.G.CGCCCG.CCATGGGCGGGGCGGGCG.CT.GGCCGA.GGGCCTTTGGGGCCCCCGGCGGTCAGGGGGATCTGtt                               <  1:2515236/101‑1 (MQ=255)


                                                                      agGCCAGCGCCCAGCCCCGGCGGCGC.GAAA.G.CGCCCG.CCATGGGCGGGGCGGGCG.CT.GGCCGA.GGGCCTTTGGGGCCCCCGGCGGTCAGGGGGGTCTGTTg                              >  2:550014/1‑101 (MQ=255)


                                                                       ggCCAGCCCCCGGCCCCGGCGGCGC.GAAG.G.CGCCCG.CCATGGGCGGGGCGGGCG.CT.GGCCGA.GGGCCTTTGGGGCCCCCGGCGGTCAGGGGGATCTGTTgc                             <  2:2781706/101‑1 (MQ=17)


                                                                        gCCAGCGCCCAGCCCCGGCGGCGC.GAAA.G.CGCCCG.CCATGGGCGGGGCGGGCG.CT.GGCCGA.GGGCCTTTGGGGCCCCCGGCGGTCAGGGGGATCTGTTgcg                            >  2:798255/1‑101 (MQ=255)


                                                                         ccAGCCCCCCGCCCGGGCGGCGC.GAAA.G.CGCCCG.CCATGGGCGGGGCGGGCG.CT.GGCCGA.GGGCCTTTGGGGCCCCCGGCGGTCAGGGGGATCTGTTgcgc                           <  1:4977916/101‑1 (MQ=17)


                                                                         ccAGCCCCCAGCCCCGGCGGCGC.GAAA.GCCC.CCG.CCATGGGCGGGGCGGGCG.CT.GGCCGA.GGGCCTTTGGGGCCCCCGGCGGTCAGGGGGATCTGTTgcgc                           <  2:1641295/101‑1 (MQ=18)


                                                                          cggcccccccccccGGCGGCGC.GAAA.G.CGCCCG.CCATGGGCGGGGCGGGCG.CT.GGCCGA.GGGCCTTTGGGGCCCCCGGCGGTCAGGGGGATCTGTTGCGCa                          <  2:1495926/91‑1 (MQ=17)


                                                                               cccAGCCCCGGCGGCGC.GAAA.G.CGCCCG.CCATGGGCGGGGCGGGCG.CT.GGCCGA.GGGCCTTTGGGGCCCCCGGCGGTCAGGGGGGTCTGTTGCGCATGGcc                     >  2:4373538/1‑101 (MQ=255)


                                                                                   gCCCGGGGGGCGC.GAAC.G.CGCCCG.CCATGGGCGGGGCGGGCG.CT.GGCCGA.GGGCCTTTGGGGCCCCCGGCGGTCAGGGGGATCTGTTGCGCATGGCCGCGt                 <  1:35887/101‑1 (MQ=18)


                                                                                    ccccGGCGGCGC.GAAG.GCCC.CCG.CCATGGGCGGGGCGGGCG.CT.GGCCGA.GGGCCTTTGGGGCCCCCGGCGGTCAGGGGGATCTGTTGCGCATGGCCGCGTc                <  2:4096455/101‑1 (MQ=18)


                                                                                     cgggggggcgC.GAAG.G.CGCCCG.CCATGGGCGGGGCGGGCG.CT.GGCCGA.GGGCCTTTGGGGCCCCCGGCGGTCAGGGGGATCTGTTGCGCATGGCCGCGTCa               <  1:3408281/95‑1 (MQ=18)


                                                                                      ccGGCGCCGC.GAAA.G.CGCCCG.CCATGGGCGGGGCGGGCG.CT.GGCCGA.GGGCCTTTGGGGCCCCCGGCGGTCAGGGGGATCTGTTGCGCATGGCCGCGTCAg              <  1:413411/101‑1 (MQ=255)


                                                                                       ggggggcgC.GAAG.G.CGCCCG.CCATGGGCGGGGCGGGCG.CT.GGCCGA.GGGCCTTTGGGGCCCCCGGCGGTCAGGGGGATCTGTTGCGCATGGCCGCGTCAgc             <  1:2111469/97‑1 (MQ=18)


                                                                                        ggcggcGC.GAAA.G.CGCCCG.CCATGGGCGGGGCGGGCG.CT.GGCCGA.GGGCCTTTGGGGCCCCCGGCGGTCAGGGGGATCTGTTGCGCATGGCCGCGTCAgcg            >  2:2862945/1‑101 (MQ=255)


                                                                                          gggcgC.GAAG.GCCC.CCG.CCATGGGCGGGGCGGGCG.CT.GGCCGA.GGGCCTTTGGGGCCCCCGGCGGTCAGGGGGATCTGTTGCGCATGGCCGCGTCAGcgcc          <  1:1628950/100‑1 (MQ=21)


                                                                                          cggcgC.GAAACC.CC.CCG.CCATGGGCGGGGCGGGCG.CT.GGCCGA.GGGCCTTTGGGGCCCCCGGCGGTCAGGGGGATCTGTTGCGCATGGCCGCGTCAGcgcc          <  2:446374/101‑1 (MQ=21)


                                                                                                 gAAA.G.CGCCCG.CCATGGGCGGGGCGGGCG.CT.GGCCGA.GGGCCTTTGGGGCCCCCGGCGGTCAGGGGGATCTGTTGCGCATGGCCGCGTCAGCGCCGCgggc    >  2:280947/1‑99 (MQ=255)


                                                                                                 gAAA.G.CGCCCG.CCATGGGCGGGGCGGGCG.CT.GGCCGA.GGGCCTTTGGGGCCCCCGGCGGTCAGGGGGATCTGTTGCGCATGGCCGCGTCAGCGCCGCGGcc    >  2:3323375/1‑101 (MQ=255)


                                                                                                 gAAA.G.CGCCCG.CCATGGGCGGGGCGGGCG.CT.GGCCGA.GGGCCTTTGGGGCCCCCGGCGGTCAGGGGGATCTGTTGCGCATGGCCGCGTCAGCGCCGCGGcc    <  1:1668530/101‑1 (MQ=255)


                                                                                                 gAAA.G.CGCCCG.CCATGGGCGGGGCGGGCG.CT.GGCCGA.GGGCCTTTGGGGCCCCCGGCGGTCAGGGGGATCTGTTGCGCATGGCCGCGTCAGCGCCGCGGcc    >  1:630929/1‑101 (MQ=255)


                                                                                                  aaa.G.CGCCCG.CCATGGGCGGGGCGGGCG.CT.GGCCGA.GGGCCTTTGGGGCCCCCGGCGGTCAGGGGGGTCTGTTGCGCATGGCCGCGTCAGCGCCGCGGCCa   >  1:4638509/1‑101 (MQ=255)


                                                                                                  aaa.G.CGCCCG.CCATGGGCGGGGCGGGCG.CT.GGCCCA.GGGCCTTTGGGGCCCCCGGCGGGCAGGGGGATCTGTTGGGCATGGCCGCGTCCGCGCCGCGGCCa   >  1:888866/1‑101 (MQ=21)


AGGAAAAGCACTGCTTTTCCCCGCCCGTGTCGCGACAGCGTTGCAAAGACCGCATCCCCAAGCGCGCGAGAGGCCAGCGCCCAGCCCCGGCGGCGC.GAAA.G.CG.CCG.CCATGGGCGGGGCGGGCG.CT.GGCCGA.GGGCCTTTGGGGCCCCCGGCGGTCAGGGGGATCTGTTGCGCATGGCCGCGTCAGCGCCGCGGCCAA  >  NC_014034/1867731‑1867928

Base quality scores:


ATCG

< 3 ≤


ATCG

< 15 ≤


ATCG

< 22 ≤


ATCG

< 33 ≤


ATCG

< 41 ≤


ATCG