BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_014034	1,396,956	+C	11.9%	coding (838/1578 nt)	RCAP_rcc01307 →	hypothetical protein

Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_014034	1,396,956	1	.→C	11.9%	152.1 ‑0.8 ‑0.4	387	coding (838/1578 nt)	RCAP_rcc01307	hypothetical protein
Reads supporting (aligned to +/- strand): new base (26/20): ref base (171/170): total (197/190)
Fisher's exact test strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 4.13e-01

Only 100 of 413 total aligned reads displayed.


GTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCGGGGGCGTG  >  NC_014034/1396856‑1397056


gTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGCCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCt                                                                                                           >  1:1701842/1‑101 (MQ=255)


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


   gATCGAAAGCGGTCCAGAGATAGACCGCGCGGAGCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.Tcc                                                                                                       <  1:4098130/101‑1 (MQ=255)


   gATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.Tcc                                                                                                       <  1:2345289/101‑1 (MQ=255)


    aTCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCg                                                                                                      <  1:3467721/101‑1 (MQ=255)


       gAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.Ct                                                                                                  >  2:3071471/1‑101 (MQ=255)


       gAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCCCG.GC.CGCA.TCGGCT.TCCGC.Ct                                                                                                  >  2:1891541/1‑101 (MQ=255)


        aaaCCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTg                                                                                                 <  2:3983125/101‑1 (MQ=255)


        aaaCCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTg                                                                                                 >  1:1997387/1‑101 (MQ=255)


          aCCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..Cg                                                                                             >  2:3236514/1‑101 (MQ=255)


           ccGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGt                                                                                            <  2:699647/101‑1 (MQ=255)


            cGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTg                                                                                           >  1:2696866/1‑101 (MQ=255)


            cGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTg                                                                                           >  2:2828122/1‑101 (MQ=255)


              gTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCC.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.Gt                                                                                        <  1:2685880/101‑1 (MQ=255)


                  agagATGGGCCGCGCGGCCCTCCTCAGACGTCTGGCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCt                                                                                    <  2:1701737/101‑1 (MQ=38)


                  agagATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCt                                                                                    <  1:857447/101‑1 (MQ=255)


                  agagATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCt                                                                                    >  1:2086580/1‑101 (MQ=255)


                   gacaTAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTc                                                                                   >  2:3798167/4‑101 (MQ=255)


                     gaTAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTccc                                                                                 >  1:1625200/1‑101 (MQ=255)


                       tAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GGCGCTccccc                                                                               >  1:1381005/1‑101 (MQ=255)


                        aGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTcccccc                                                                              >  1:4188940/1‑101 (MQ=255)


                           ccGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCcccg                                                                           >  1:4044900/1‑101 (MQ=255)


                               gcgGACCCCCTCAGACGGCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCggggg                                                                       <  1:1858512/101‑1 (MQ=255)


                                    cctcctCAGACGTCGCTCCCCCCGCGGGGTCGTTGATCGATACA.GCCT.GCGCT.GC.CGCA.TCGGCT.TCCGC.CTG..CGTT.GTCGCCCCCCCCGCGGGGGCGTgg                                                                  <  1:4712738/101‑1 (MQ=38)


                                    cctcctCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCgcat                                                                  <  2:4388010/101‑4 (MQ=255)


                                     tcgcggttgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTgg                                                                  <  2:1250707/92‑1 (MQ=0)


                                     tcgcggttgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTgg                                                                  <  1:2243023/92‑1 (MQ=0)


                                     attggccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGa                                                                 >  2:1594426/11‑101 (MQ=0)


                                      tcctcAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTgggt                                                                >  2:4038780/1‑99 (MQ=255)


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


                                       cctcAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATc                                                               >  1:1556395/1‑101 (MQ=255)


                                       acccgcccGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTgg                                                                  >  1:3949070/8‑101 (MQ=0)


                                        ctcaagcggcGCGCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCg                                                              <  2:266986/92‑1 (MQ=37)


                                        caacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATc                                                               <  2:1669496/94‑1 (MQ=9)


                                         gccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                             >  2:2799351/7‑101 (MQ=2)


                                         aacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGGGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCACGGGCGCGCGGacca                                                              >  2:2367170/7‑98 (MQ=2)


                                         aacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCg                                                              >  2:1380891/7‑101 (MQ=9)


                                          cgcccGTCCCTCCCCCCCCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATc                                                               <  2:2263893/97‑1 (MQ=2)


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


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


                                          cAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                            <  2:102252/101‑1 (MQ=255)


                                          cAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                            <  2:1997270/101‑1 (MQ=255)


                                           ttgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                            <  1:389107/98‑1 (MQ=2)


                                           cactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                            >  1:1234552/5‑101 (MQ=9)


                                            tgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                           <  2:148623/99‑1 (MQ=2)


                                            gACGTGGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                          <  1:2857765/101‑1 (MQ=255)


                                            gACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGGTCGAAAc                                                          >  1:4618125/1‑101 (MQ=255)


                                            ccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                          <  1:4547005/98‑1 (MQ=2)


                                            ccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                          <  1:332461/98‑1 (MQ=2)


                                             gcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                          >  2:83022/2‑101 (MQ=2)


                                             ctgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                          >  2:3294759/3‑101 (MQ=9)


                                             cggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAct                                                         <  1:783361/99‑2 (MQ=2)


                                               gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                          >  1:2782668/1‑101 (MQ=2)


                                               gTCGCGCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActac                                                       <  1:444856/101‑4 (MQ=2)


                                               gTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAccc                                                        <  2:4067860/101‑3 (MQ=2)


                                                tCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACa                                                         <  1:3492632/101‑1 (MQ=2)


                                                tCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGGGGATCGAAACATCa                                                      >  1:614202/1‑101 (MQ=255)


                                                tCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcg                                                       >  2:2521528/1‑98 (MQ=9)


                                                 cGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACac                                                        >  1:439379/1‑100 (MQ=2)


                                                 cGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActacac                                                     <  2:2162605/101‑6 (MQ=2)


                                                  gcgccccccGCGGGGGCGTGGATCGAAAC..GCTG.GCCCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggc                                                     <  1:1324092/98‑6 (MQ=9)


                                                  gCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActacacc                                                    >  2:540992/1‑95 (MQ=2)


                                                  gCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggc                                                     <  1:3065801/101‑6 (MQ=9)


                                                  gCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggc                                                     >  2:3613733/1‑96 (MQ=9)


                                                  gCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggc                                                     >  1:3969387/1‑96 (MQ=9)


                                                   ccccccccGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcc                                                    <  2:3862502/99‑7 (MQ=9)


                                                   cTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcc                                                    <  2:3733892/101‑7 (MQ=9)


                                                    tCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCaa                                                     >  2:1767536/1‑100 (MQ=2)


                                                    tCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCaa                                                     >  1:79411/1‑100 (MQ=2)


                                                    tCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCtc                                                  >  1:363386/1‑101 (MQ=255)


                                                    tCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcct                                                   >  2:1726268/1‑94 (MQ=9)


                                                    cccccccGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcct                                                   <  1:1000110/100‑8 (MQ=9)


                                                     ccccccGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCaat                                                    <  1:2406552/101‑3 (MQ=2)


                                                     ccccccGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCtcg                                                 <  2:2047798/101‑1 (MQ=255)


                                                     ccccccGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcctc                                                  <  1:2159914/101‑9 (MQ=9)


                                                      cccccGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcctcg                                                 <  1:2505007/101‑10 (MQ=9)


                                                        cccGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCaatgcc                                                 >  1:1463097/1‑96 (MQ=0)


                                                        cccGCGGGGGCGTGGATCGAAACA.GCCG.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCtcgtcg                                              <  2:3455038/101‑1 (MQ=255)


                                                          cgcgGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCaatgccga                                               <  1:1083774/101‑8 (MQ=0)


                                                           gcgGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGGCGTTg                                           >  2:4195614/1‑101 (MQ=255)


                                                               gggCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCa                                       <  1:1143073/101‑1 (MQ=255)


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


                                                               gggCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGGCa                                       >  2:779461/1‑101 (MQ=255)


                                                                  cGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCATCGTTGGTCAGcc                                    >  1:3184106/1‑101 (MQ=255)


                                                                    tGGAGCGAAACA.GCCT.GCGCG.GC.CGCA.GCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCgg                                  <  1:2582197/101‑1 (MQ=255)


                                                                     ggATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGa                                 <  2:873478/101‑1 (MQ=255)


                                                                      gATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTTGTCCGCCGGAc                                >  1:712573/1‑101 (MQ=255)


                                                                      gATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGAc                                >  1:3277656/1‑101 (MQ=255)


                                                                           aaaCA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGt                           <  2:2717255/101‑1 (MQ=255)


                                                                                ggCCC.GCGCG.GC.AGCA.CCGTCT.CCCGC.TTG..CGTG.GTCACTCCCCCCGCTGTGTCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGt                       <  1:898098/100‑1 (MQ=255)


                                                                                      gcgcG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTc                  <  2:4337311/101‑1 (MQ=255)


                                                                                       cgcg.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTcc                 <  2:642784/101‑1 (MQ=255)


                                                                                       cgcg.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTcc                 <  1:3798402/101‑1 (MQ=255)


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


                                                                                            gC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTcccccc             <  2:1469522/101‑1 (MQ=255)


                                                                                            gC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTcccccc             <  1:3019088/101‑1 (MQ=255)


                                                                                                gCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCcgcg          >  2:484125/1‑101 (MQ=255)


                                                                                                 cA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCgg         <  2:2433963/101‑1 (MQ=255)


                                                                                                        cT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCGGGGGCGt   >  2:3385689/1‑101 (MQ=255)


                                                                                                        cT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCGGGGGCGt   >  1:2621771/1‑101 (MQ=255)


GTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCGGGGGCGTG  >  NC_014034/1396856‑1397056

Base quality scores:


ATCG

< 3 ≤


ATCG

< 15 ≤


ATCG

< 24 ≤


ATCG

< 33 ≤


ATCG

< 39 ≤


ATCG