BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_014034	1,396,946	+A	4.9%	coding (828/1578 nt)	RCAP_rcc01307 →	hypothetical protein

Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_014034	1,396,946	1	.→A	4.9%	63.5 ‑2.9 ‑0.2	368	coding (828/1578 nt)	RCAP_rcc01307	hypothetical protein
Reads supporting (aligned to +/- strand): new base (5/13): ref base (180/170): total (185/183)
Fisher's exact test strand distribution p-value = 5.62e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.22e-01

Only 100 of 385 total aligned reads displayed.


CCGCGGGGGCGTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCG.C.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCG  >  NC_014034/1396846‑1397046


ccGCGGGGGCGTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.Gc                                                                                                              >  2:1720113/1‑101 (MQ=255)


   cgGGGGCGTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGc                                                                                                          >  2:1837274/1‑101 (MQ=255)


    gggggCGTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCa                                                                                                         <  2:1961925/101‑1 (MQ=255)


       ggCGTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCg                                                                                                     >  2:3784419/1‑101 (MQ=255)


         cGTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGc                                                                                                   <  2:1631884/101‑1 (MQ=255)


         cGTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGc                                                                                                   >  1:14634/1‑101 (MQ=255)


          gTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCt                                                                                                  <  1:1880600/101‑1 (MQ=255)


           tGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCt.t                                                                                                >  2:3586799/1‑101 (MQ=255)


            ggAGCGAAACCGGTCCAGAGATCGACCGCGCGGACCTCCTCAGGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.Tc                                                                                               <  2:1749438/101‑1 (MQ=255)


                 gAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCG.C.Ct                                                                                        <  1:3012471/101‑1 (MQ=255)


                  aaaCCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCG.c.ccc                                                                                       >  1:1561340/1‑99 (MQ=255)


                  aaaCCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCG.C.CTg                                                                                       >  2:318769/1‑101 (MQ=255)


                   aaCCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCG.C.CTG..c                                                                                    >  1:521643/1‑101 (MQ=255)


                   aaCCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCG.C.CTG..c                                                                                    >  1:2806305/1‑101 (MQ=255)


                   aaCCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCG.C.CTG..c                                                                                    >  2:285277/1‑101 (MQ=255)


                    aCCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCG.C.CTG..Cg                                                                                   <  1:1370727/101‑1 (MQ=255)


                    aCCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCG.C.CTG..Cg                                                                                   <  2:3943762/101‑1 (MQ=255)


                        gTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGCGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCG.C.CTG..CGCg.gg                                                                              >  1:3764269/1‑100 (MQ=255)


                         tCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCG.C.CTG..CGTG.GTc                                                                             <  2:953218/101‑1 (MQ=255)


                            agagATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCG.C.CTG..CGTG.GTCGCt                                                                          >  2:588571/1‑101 (MQ=255)


                              agaTAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCG.C.CTG..CGTG.GTCGCTcc                                                                        >  1:1769928/1‑101 (MQ=255)


                              agaTAGACCGCGCGGACCTCCTCAGACGTCGCGCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCG.C.CTG..CGTG.GTCGCTcc                                                                        <  1:2754403/101‑1 (MQ=255)


                                    aCCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCCAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCG.C.CTG..CGTG.GTGGCTGCCCCcgc                                                                  >  1:247800/1‑101 (MQ=255)


                                        cgcgGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCG.C.CTG..CGTG.GTCGCTCCCCCCGCgggg                                                              <  1:3551321/101‑1 (MQ=255)


                                         gcgGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCG.C.CTG..CGTG.GTCGCTCCCCCCGCggggg                                                             >  1:1047280/1‑101 (MQ=255)


                                           ggACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCG.C.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCg                                                           <  2:54571/101‑1 (MQ=255)


                                             aCCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCG.C.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTg                                                         >  2:74058/1‑101 (MQ=255)


                                              ctcgcggttgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T.TT.C.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTg                                                         <  1:1816521/91‑1 (MQ=0)


                                               attggccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATG.CGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGa                                                       >  1:2879910/11‑101 (MQ=0)


                                                gtcaacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTG.T.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGa                                                       >  2:2125338/10‑101 (MQ=9)


                                                gtcaacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTG.T.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGa                                                       >  1:3687396/10‑101 (MQ=9)


                                                gtcaacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTG.T.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGa                                                       >  1:100646/10‑101 (MQ=9)


                                                 tggccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATG.CGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATc                                                     >  1:2324378/9‑101 (MQ=2)


                                                 tcaacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTG.T.CTT..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGAt                                                      <  1:3971736/93‑1 (MQ=2)


                                                 tcaacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTG.T.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGAt                                                      <  1:3612837/93‑1 (MQ=9)


                                                 tcaacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTG.T.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGAt                                                      <  1:982790/93‑1 (MQ=9)


                                                 cctcAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCG.C.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATc                                                     <  2:2408563/101‑1 (MQ=255)


                                                 acccgcccGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AGCG.C.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTgg                                                        <  1:457664/94‑1 (MQ=0)


                                                  ctcAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCG.C.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCg                                                    <  2:1561235/101‑1 (MQ=255)


                                                  cccgcccGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AGCG.C.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGa                                                       <  1:2967672/95‑1 (MQ=2)


                                                   tcAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCG.C.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                   >  1:938767/1‑101 (MQ=255)


                                                   gccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATG.CGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                   <  1:3118427/95‑1 (MQ=2)


                                                   aacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTG.T.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCATGGATCg                                                    <  1:3871980/95‑1 (MQ=9)


                                                    cAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCG.C.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                  >  2:2737221/1‑101 (MQ=255)


                                                    acactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTG.T.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                   >  1:331521/6‑101 (MQ=9)


                                                     cccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATG.CGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                 >  2:3522513/5‑101 (MQ=2)


                                                     aGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCG.C.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                 <  1:318804/101‑1 (MQ=255)


                                                      gACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCG.C.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                >  1:590814/1‑101 (MQ=255)


                                                      ccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATG.CGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                <  2:4063411/98‑1 (MQ=2)


                                                       cggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATGAG.CTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAct                                               <  1:835095/99‑2 (MQ=2)


                                                       cggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATG.CGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAct                                               >  1:1689437/3‑100 (MQ=2)


                                                       cggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATG.CGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAct                                               <  1:4081188/99‑2 (MQ=2)


                                                        tgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTG.T.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAct                                               >  2:2499445/2‑100 (MQ=9)


                                                        tgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTG.T.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAct                                               >  1:4316482/2‑100 (MQ=9)


                                                        cGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCG.C.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAt                                              <  2:3509873/101‑1 (MQ=255)


                                                         gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.ACCT.GCGCG.GC.CGCA.TCGGCT.TCCG.C.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATc                                             <  1:3795696/101‑1 (MQ=255)


                                                          tCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTG.T.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcg                                             <  1:3408607/101‑4 (MQ=9)


                                                           cGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AGCG.C.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACac                                              <  2:2341368/101‑2 (MQ=2)


                                                            gCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T.TT.C.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcccgcc                                           <  1:2622314/101‑6 (MQ=2)


                                                            gCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TTCG.C.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGGGGGTCGAAACCTCAGc                                          >  1:2409259/1‑101 (MQ=37)


                                                            gCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATG.CGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActacacc                                          <  1:1925187/101‑7 (MQ=2)


                                                            gCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTG.T.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggc                                           <  1:215334/101‑6 (MQ=9)


                                                            ccTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCG.C.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGc                                          <  2:2041455/100‑1 (MQ=255)


                                                             cTCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AGCG.C.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCa                                            >  1:1117362/1‑101 (MQ=2)


                                                             cTCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AGCG.C.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTCGATCGAAACACCa                                            >  2:3707640/1‑101 (MQ=2)


                                                             cTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATG.CGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActacaccg                                         >  1:19824/1‑94 (MQ=2)


                                                             cTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTG.T.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcc                                          >  2:3795821/1‑95 (MQ=9)


                                                              tCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTG.T.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTCGATCGAAActcggcct                                         >  2:671292/1‑94 (MQ=9)


                                                               ccccccGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T.TT.C.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGGGGATCGAAAcccgccacg                                        >  2:4148970/1‑93 (MQ=0)


                                                               ccccccGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCG.C.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCtcg                                       >  2:1450816/1‑101 (MQ=255)


                                                                cccccGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AGCG.C.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCaatg                                         <  1:2297508/101‑4 (MQ=2)


                                                                cccccGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AGCG.C.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCaatg                                         <  2:1328175/101‑4 (MQ=2)


                                                                cccccGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T.TT.C.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcccgccccgt                                       >  2:1128457/1‑92 (MQ=2)


                                                                 ccccGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCG.C.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCtcgtc                                     >  1:4250044/1‑101 (MQ=255)


                                                                 ccccGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTG.T.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcctcga                                      >  2:3952928/1‑91 (MQ=2)


                                                                 ccccGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTG.T.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcctcga                                      >  1:3585290/1‑91 (MQ=2)


                                                                  cccGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AGCG.C.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCaatgcc                                       <  1:1610809/101‑6 (MQ=0)


                                                                   ccGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AGCG.C.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCaatgccg                                      <  2:1079132/101‑7 (MQ=0)


                                                                      tgggggCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCG.C.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTgg                                >  1:684610/2‑101 (MQ=255)


                                                                      cgGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.TGCA.TCGGCT.TCCC.C.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGGGGATCGAAACATCAGCTCGTCGTTgg                                >  1:3950927/1‑101 (MQ=255)


                                                                       gggggCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCG.C.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCCCGTCGTTGGt                               >  1:4154273/1‑101 (MQ=255)


                                                                       gggggCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCG.C.CTG..CGTG.GTCGCTCCCCCCCCGGGGGCGTGGATCGAAACATCAGCTCGTCTTTGGt                               >  2:4393506/1‑101 (MQ=255)


                                                                        ggggCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCG.C.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTc                              <  2:2995285/101‑1 (MQ=255)


                                                                         gggGGTGGATGGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCG.C.CTG..CGTG.GTCGCGCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCa                             <  1:2083712/101‑1 (MQ=255)


                                                                         gggCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCG.C.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCa                             <  2:852037/101‑1 (MQ=255)


                                                                         gggCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCG.C.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCa                             >  2:3873928/1‑101 (MQ=255)


                                                                            cGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCG.C.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGcc                          >  2:3425690/1‑101 (MQ=255)


                                                                             gTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCG.C.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCg                         <  2:3333458/101‑1 (MQ=255)


                                                                              tGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCG.C.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCgg                        >  1:1326538/1‑101 (MQ=255)


                                                                                   cGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCG.C.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGAt                   >  2:451913/1‑101 (MQ=255)


                                                                                      aaCA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCG.C.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTc                <  2:2391997/101‑1 (MQ=255)


                                                                                      aaCA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCG.C.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTc                >  2:21219/1‑101 (MQ=255)


                                                                                             cT.GCGCG.GC.CGCA.TCGGCT.TCCG.C.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGc          >  1:1664569/1‑101 (MQ=255)


                                                                                             cT.GCGCG.GC.CGCA.TCGGCT.TCCG.C.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGc          <  1:4439897/101‑1 (MQ=255)


                                                                                              t.GCGCG.GC.CGCA.TCGGCT.TCCG.C.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCt         >  2:2231599/1‑101 (MQ=255)


                                                                                              t.GCGCG.GC.CGCA.TCGGCT.TCCG.C.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCt         <  1:2021877/101‑1 (MQ=255)


                                                                                                  gcg.GC.CGCA.TCGGCT.TCCG.C.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTccc      <  1:1186559/101‑1 (MQ=255)


                                                                                                   cg.GC.CGCA.TCGGCT.TCCG.C.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTcccc     <  2:4271283/101‑1 (MQ=255)


                                                                                                      gC.CGCA.TCGGCT.TCCG.C.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTcccccc   <  2:2645353/101‑1 (MQ=255)


                                                                                                       c.cGCA.TCGGCT.TCCG.C.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTccccccc  >  2:1408120/1‑100 (MQ=255)


CCGCGGGGGCGTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCG.C.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCG  >  NC_014034/1396846‑1397046

Base quality scores:


ATCG

< 3 ≤


ATCG

< 18 ≤


ATCG

< 26 ≤


ATCG

< 33 ≤


ATCG

< 39 ≤


ATCG