BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_014034	448,780	A→C	2.3%	intergenic (+73/+147)	RCAP_rcc00419 → / ← pykA1	NADH:flavin oxidoreductase/NADH oxidase/pyruvate kinase

Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_014034	448,780	0	A→C	2.3%	3.0 ‑0.4 ‑0.2	181	intergenic (+73/+147)	RCAP_rcc00419/pykA1	NADH:flavin oxidoreductase/NADH oxidase/pyruvate kinase
Reads supporting (aligned to +/- strand): new base (4/2): ref base (131/44): total (135/46)
Fisher's exact test strand distribution p-value = 6.45e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.87e-01

Only 100 of 247 total aligned reads displayed.


TGCGCTGCGGCTGATGAAGGACGTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCCAATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAAGGCACTCG.GCCT.GCGCCCGCCCCGCCATCGGCGGGC.GCAGTCG.C.G.CCGCCGGGGCGGAGCGCAGGCCTCTGTTGTGCTTGGGGCAGGCGGTCGTGCCGCAGCGC  >  NC_014034/448680‑448880


tGCGCTGCGGCTGATGAATGACGTTTGAGGCTGACTTTGCCTCGGCCATTGCCGAGGCCAATCGCCACAGAGGTCCTGTGACCGCCGACGGCCCAAACGCa                                                                                                            >  1:6615004/1‑101 (MQ=38)


tGCGCTGCGGCTGATGAAGGACGTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCCAATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAAGGCa                                                                                                            <  1:5510622/101‑1 (MQ=255)


 gcgcTGCGGCTGATGAAGGACGTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCCAATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAAGGCAc                                                                                                           <  2:1008912/101‑1 (MQ=255)


 gcgcTGCGGCTGATGAAGGACGTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCCAATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAAGGCAc                                                                                                           <  1:215882/101‑1 (MQ=255)


   gcTGCGGCTGATGAAGGACGTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCCAATCGCAACAGCGGTCCTGTGACCGCCGAGGGCCCAAAGGCACTc                                                                                                         >  1:957557/1‑101 (MQ=255)


    cTGCGGCTGATGAAGGACGTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCCAATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAAGGCACTCg                                                                                                        >  2:826413/1‑101 (MQ=255)


    cTGCGGCTGATGAAGGACGTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCCAATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAAGGCACTCg                                                                                                        <  1:720078/101‑1 (MQ=255)


      gCGGCTGATGAAGGACGTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCCAATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCCAAGGCACTCG.Gc                                                                                                     >  1:3993250/1‑101 (MQ=255)


      gCGGCTGATGAAGGACGTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCCAATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAAGGCACTCG.Gc                                                                                                     >  2:1649514/1‑101 (MQ=255)


      gCGGCTGATGAAGGACGTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCCAATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAAGGCACTCG.Gc                                                                                                     <  2:5656201/101‑1 (MQ=255)


       cGGCTGATGAAGGACGTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCCAATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAAGGCACTCG.Gcc                                                                                                    <  1:478964/101‑1 (MQ=255)


        ggCTGATGAAGGACGTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCCAATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAAGGCACTCG.GCCt                                                                                                   <  2:5943504/101‑1 (MQ=255)


        ggCTGATGAAGGACGTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCCAATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAAGGCACTCG.GCCt                                                                                                   <  1:2280152/101‑1 (MQ=255)


          cTGATGAAGGACGTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCCAATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAAGGCACTCG.GCCtcc                                                                                                 >  1:2375042/1‑99 (MQ=255)


           tgatgaAGGACGTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCCAATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAAGGCACTCG.GCCT.gcg                                                                                               >  2:4229883/1‑101 (MQ=255)


           tgatgaAGGACGTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCCAATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAAGGCACTCG.GCCT.gcg                                                                                               >  2:4873250/1‑101 (MQ=255)


           tgatgaAGGACGTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCCAATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAAGGCACTCG.GCCT.gcg                                                                                               >  1:827229/1‑101 (MQ=255)


            gatgaAGGACGTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCCAATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAAGGCACTCG.GCCT.gcgc                                                                                              >  1:2527473/1‑101 (MQ=255)


             atgaAGGACGTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCCAATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAAGGCACTCG.GCCT.Gcgcc                                                                                             <  1:444460/101‑1 (MQ=255)


             atgaAGGACGTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCCAATCGCAACAGAGGCCCTGTGACCGCCGAGGGCCCAAAGGCACTCG.GCCT.gcggc                                                                                             >  2:483054/1‑99 (MQ=255)


              tgaAGGACGTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCCAATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAAGGCACTCG.GCCT.Gcgccc                                                                                            <  2:6004801/101‑1 (MQ=255)


              tgaAGGACGTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCCAATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAAGGCACTCG.GCCT.Gcgccc                                                                                            >  2:1515574/1‑101 (MQ=255)


               gaAGGACGTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCCAATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAAGGCACTCG.GCCT.Gcgcccg                                                                                           >  1:3021774/1‑101 (MQ=255)


               gaAGGACGTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCCAATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAAGGCACTCG.GCCT.Gcgcccg                                                                                           >  1:2072167/1‑101 (MQ=255)


                aaGGACGTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCCAATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAAGGCACTCG.GCCT.Gcgcccgc                                                                                          >  1:11493/1‑101 (MQ=255)


                 aGGCCGTTGGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCCAATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAAGGCACTCG.GCCT.Gcgcccgcc                                                                                         <  2:1744595/101‑1 (MQ=255)


                 aGGACGTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCCAATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAAGGCACTCG.GCCT.Gcgcccgcc                                                                                         >  2:6092591/1‑101 (MQ=255)


                 aGGACGTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCCAATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAAGGCACTCG.GCCT.Gcgcccgcc                                                                                         >  1:3170127/1‑101 (MQ=255)


                 aGGACGTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCCAATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAAGGCACTCG.GCCT.Gcgcccgcc                                                                                         <  1:291983/101‑1 (MQ=255)


                   gACGTTTGAGGCTGACTTGGCCTAGGCCATTGCCGAGGCCAATCGCAAGAGAGGTCCTGTCACCGCCGAGGGCCCAAAGGCACTCG.GCCT.GCGcccgcccc                                                                                       >  1:5809923/1‑101 (MQ=255)


                      gTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCCAATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAAGGCACTCG.GCCT.GCGCccgccccgcc                                                                                    >  2:1840754/1‑101 (MQ=255)


                      gTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCCAATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAAGGCACTCG.GCCT.GCGCccgccccgcc                                                                                    >  2:5012338/1‑101 (MQ=255)


                        ttGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCCAATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAAGGCACTCG.GCCT.GCGCCCGCCCCGCCAt                                                                                  >  2:517504/1‑101 (MQ=255)


                         tGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCCAATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAAGGCACTCG.GCCT.GCGCCCGCCCCGCCATc                                                                                 >  2:7008806/1‑101 (MQ=255)


                          gAGGCTGACTTGGCCTCGGCCATTGCCGAGGCCAATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAAGGCACTCG.GCCT.GCGCCCCCCCCGCCCTcg                                                                                >  2:3299081/1‑101 (MQ=255)


                          gAGGCTGACTGGGCCTCGCCCATTGCCGAGGCCAATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAAGGCACTCG.GCCT.GCGCCCGCCCCGCCATcg                                                                                <  2:1260267/101‑1 (MQ=255)


                              cTGACTTGGCCTCGGCCATTGCCGAGGCCAATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAAGGCACTCG.GCCT.GCGCCCGCCCCGCCCTcggcgg                                                                            >  1:1387080/1‑101 (MQ=255)


                              cTGACTTGGCCTCGGCCATTGCCGAGGCCAATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAAGGCACTCC.GCCT.GCCCCCGCCCCGCCATcggcgg                                                                            >  2:1888571/1‑101 (MQ=255)


                               tgccTTGGCCCCGGCCATGGCCGGGGCCAAGCGCAACAGAGGTCCTGGGACCGCCGAGGGCCCAAAGGCACTCG.GCCT.GCGCCCGCCCCGCCATCggcggg                                                                           <  1:508275/98‑1 (MQ=18)


                               tGACTTGGCCTCGGCCATTGCCGAGGCCAATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAAGGCACTCG.GCCT.GCGCCCGCCCCGCCATCggcggg                                                                           >  1:5158487/1‑101 (MQ=255)


                                 aCTTGGCCTCGGCCATTGCCGAGGCCAATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAAGGCACTCG.GCCT.GCGCCCGCCCCGCCATCggcggggg                                                                         >  2:2898096/1‑99 (MQ=255)


                                 aCTTGGCCTCGGCCATTGCCGAGGCCAATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAAGGCACTCG.GCCT.GCGCCCGCCCCGCCATCggcggggg                                                                         >  2:3666947/1‑99 (MQ=255)


                                    tGGCCTCGGCCATTGCCGAGGCCAATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAAGGCACTCG.GCCT.GCGCCCGCCCCGCCATCggcgggggggg                                                                      >  1:3421339/1‑96 (MQ=255)


                                    tGGCCTCGGCCATTGCCGAGGCCAATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAAGGCACTCG.GCCT.GCGCCCGCCCCGCCATCggcgggccggg                                                                      >  2:3927301/1‑97 (MQ=255)


                                    tGGCCTCGGCCATTCCCGGGGCCAATCGCAACAGAGGTCCTGTGCCCGCCGAGGGCCCAAAGGCACTCG.GCCT.GCGCCCGCCCCGCCATCGGCGGGC.GCAg                                                                     <  1:5838501/101‑1 (MQ=255)


                                    gggCCTCGGCCATGGCCGAGGCCAAGCGCAACAGAGGTCCTGTGCCCGCCGAGGGCCCAAAGGCACTCG.GCCT.GCGCCCGCCCCGCCATCggcggggcggc                                                                      <  2:3577882/100‑6 (MQ=18)


                                     ggCCTCGGCCATTGCCGAGGCCAATCGCAACAGAGGTCCTGTGCCCGCCGAGGGCCCAAAGGCCCTCG.GCCT.GCGCCCGCCCCGCCATCGGCGGGC.GCAGt                                                                    <  2:4520577/101‑1 (MQ=255)


                                      gCCTGGGCCATTGCGGAGGCCAAGCGCAACAGAGGCCCTGTGCCCGCCGAGGGCCCAAAGGCCCTCG.GCCT.GCGCCCGCCCCGCCATCGGCGGGC.GCAGTc                                                                   <  1:6622549/101‑1 (MQ=17)


                                        cTCGGCCATTGCCGAGGCCAATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAAGGCACTCG.GCCT.GCGCCCGCCCCGCCATCggcgggggggggcgc                                                                  >  1:3209859/1‑92 (MQ=255)


                                         tCGGCCATTGCCGAGGCCAATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAAGGCACTCG.GCCT.GCGCCCGCCCCGCCATCGGCGGGG.GCGGGcg.c.g                                                              >  1:4164847/1‑101 (MQ=255)


                                           ggCCATTGCCGAGGCCAATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAAGGCACTCG.GCCT.GCGCCCGCCCCGCCATCGGCGGGG.GCCGGCG.g.g.cc                                                           >  1:679422/1‑101 (MQ=18)


                                            gccttgcccGGGGCAAATCGCAACGGAGGTCCTGTGACCGCCGAGGGCCCAAAGCCCTTCG.GCCT.GCGCCCGCCCCGCCATCGGCGGGC.GCAGTCG.c.g.ccg                                                          <  2:2105197/94‑1 (MQ=17)


                                            gCCATTGCCGGGCCCAATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAAGGCACTCG.GCCT.GCGCCCGCCCCGCCATCGGCGGGC.GCAGTCG.c.g.ccg                                                          <  2:6354399/101‑1 (MQ=255)


                                            gCCATGGCCGAGGCAAATCGCAACAGGGGTCCTGTGCCCGCCGAGGGCCCAAAGGCACTCG.GCCT.GCGCCCGCCCCGCCATCGGCGGGC.GCAGTCG.c.g.ccg                                                          <  1:598324/101‑1 (MQ=18)


                                            gCCATGGCCGAGCCCAATCGCAACAGAGGCCCGGTGCCCGCCGAGGGCCAAAAGGCACTCG.GCCT.GCGCCCGCCCCGCCATCGGCGGGC.GCAGTCG.c.g.ccg                                                          <  1:6423356/101‑1 (MQ=18)


                                              catgcccGGGGCCAATCGCAACGGAGGTCCTGTGCCCGCGGAGGGCCCAAAGGCACTCG.GCCT.GCGCCCGCCCCGCCATCGGCGGGC.GCAGTCG.C.g.ccgcc                                                        <  2:558110/96‑1 (MQ=18)


                                              cATTGCCGAGGCCAATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAAGGCACTCG.GCCT.GCGCCCGCCCCGCCATCGGCGGGG.GGGGGCG.C.g.ccccc                                                        >  1:1010171/1‑101 (MQ=21)


                                              cATGGCCGGGGCCAATCGCAACAGAGGCCCTGTGCCCGCCGAGGGCCAAAAGGCACGCG.GCCT.GCGCCCGCCCCGCCATCGGCGGGC.GCAGTCG.C.g.ccgcc                                                        <  1:1094135/101‑1 (MQ=18)


                                               aTTGCCGAGGCCAATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAAGGCACTCG.GCCT.GCGCCCGCCCCGCCATCGGCGGGG.GCGGGGG.C.g.ccgccc                                                       >  2:6342996/1‑100 (MQ=21)


                                                    cGAGGCCAATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAAGGCACTCG.GCCT.GCGCCCGCCCCGCCATCGGCGGGG.GGGGGCG.C.G.GCGCCcggggg                                                  >  1:1502883/1‑99 (MQ=21)


                                                    cGAGGCCAATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAAGGCACTCG.GCCT.GCGCCCGCCCCGCCATCGGCGGGG.GGGGGCG.C.G.CCCCCcggggg                                                  >  1:166362/1‑99 (MQ=21)


                                                    cGAGGCCAATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAAGGCACTCG.GCCT.GCGCCCGCCCCGCCATCGGCGGGG.GGCGTCG.G.G.GCCCCGCGGCg                                                  >  1:1702980/1‑101 (MQ=18)


                                                    cGAGGCCAATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAAGGCACTCG.GCCT.GCGCCCGCCCCGCCATCGGCGGGG.GCGGGGG.C.G.CCCCCcggggg                                                  >  2:6422755/1‑99 (MQ=21)


                                                    cGAGGCCAATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAAGGCACTCG.GCCT.GCGCCCGCCCCGCCATCGGCGGGG.GCGGGCG.C.G.CCGCCcggggg                                                  >  1:5090881/1‑99 (MQ=21)


                                                     gAGGCCAATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAAGGCACTCG.GCCT.GCGCCCGCCCCGCCATCGGCGGGG.GCGGGCG.G.G.GCGCCCGGGGgg                                                 >  2:3699642/1‑101 (MQ=18)


                                                     gAGGCCAATCGCAACAGAGGTCCTGTGACCCCCGAGGGCCAAAAGGCCCCCG.GCCT.GCGCCCGCCCCGCCATCGGCGGGC.GCAGTCG.C.G.CCGCCGGGGCgg                                                 <  1:6343320/101‑1 (MQ=18)


                                                        gccaacggaaaCAGGGGTCCGGTGCCCGCGGAGGGCCCAAGGGCCCTGG.GCCT.GCGCCCGCCCCGCCATCGGCGGGC.GCAGTCG.C.G.CCGCCGGGGCGGAgc                                              <  1:1440153/92‑1 (MQ=11)


                                                        gCCAATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAAGGCACTCG.GCCT.GCGCCCGCCCCGCCATCGGCGGGG.GGCGGCG.C.G.CCCCCCGGGCggggg                                              >  2:4776492/1‑98 (MQ=18)


                                                        ccaaatccaaaCAGGGGCCCTGTGCCCGCCGAGGCCCAAAAGGCCCTCG.GCCT.GCGCCCGCCCCGCCATCGGCGGGC.GCAGTCG.C.G.CCGCCGGGGCGGAgc                                              <  1:3337689/92‑1 (MQ=11)


                                                          cAATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAAGGCACTCG.GCCT.GCGCCCGCCCCGCCATCGGCGGGG.GGGGACC.C.G.CCCCCCGGGCGGGgcgg                                            >  2:2752372/1‑100 (MQ=18)


                                                           aaTCGCAACAGAGGTCCTGTGCCCGCCGAGGGCCCAAAGGCACTCG.GCCG.GCGCCCGCCCCGCCATCGGCGGGC.GCAGTCG.C.G.CCGCCGGGGCGGAGCGCa                                           <  1:2872563/101‑1 (MQ=255)


                                                           aaTCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAAGGCACTCG.GCCT.GCGCCCGCCCCGCCATCGGCGGGG.GGCGGCG.CGG.CCGCCGGGGGggggggg                                            >  2:3225343/1‑96 (MQ=18)


                                                           aaTCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAAGGCACTCG.GCCT.GCGCCCGCCCCGCCATCGGCGGGG.GCGGGCG.C.G.CCCCCCGGGCggggggcc                                           >  2:4296094/1‑95 (MQ=18)


                                                           aaTCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAAGGCACTCG.GCCT.GCGCCCGCCCCGCCATCGCCGGGG.GGGGGCG.C.G.CCGCCCGGGGggggggcg                                           >  1:4860038/1‑95 (MQ=18)


                                                              cgaccgggagGTCCTGGGCCCGCCGAGGGCCAAAGGGCCCGCG.GCGT.GCGCCCGCCCCGCCATCGGCGGGC.GCAGTCG.C.G.CCGCCGGGGCGGAGCGCAGGc                                        <  1:4195875/94‑1 (MQ=14)


                                                              cGCAACGGAGGTCCTGTGCCCGCGGGGGGCCCAAAGGCCCTCG.GCCT.GCGCCCGCCCCGCCATCGGCGGGC.GCAGTCG.C.G.CCGCCGGGGCGGAGCGCAGGc                                        <  2:351995/101‑1 (MQ=14)


                                                              cGCAACGGAGGCCCTGTGACCGCCGGGGGCCCAAAGGCCCTCG.GCCG.GCGCCCGCCCCGCCATCGGCGGGC.GCAGTCG.C.G.CCGCCGGGGCGGAGCGCAGGc                                        <  1:3820737/101‑1 (MQ=14)


                                                              cGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAAGGCACTCG.GCCT.GCGCCCGCCCCGCCATCGGCGGGG.GCGGCGG.G.G.CCGCCGGGGGgggggggcggc                                        >  1:6906907/1‑92 (MQ=17)


                                                                cAACAGAGGTCCTGTGACCGCCGAGGGCCCAAAGGCACTCG.GCCT.GCGCCCGCCCCGCCATCGGCGGGG.GGAGGCG.C.G.GCGCCGGGGCGGGGCGCGGGccc                                      >  1:5403039/1‑100 (MQ=17)


                                                                 aaCAGGGGCCCTGTGCCCGCGGAGGCCCAAAAGGCCCTGG.CCCT.GCGCCCGCCCCGCCATCGGCGGGC.GCAGTCG.C.G.CCGCCGGGGCGGAGCGCAGGCctc                                     <  2:1989114/101‑1 (MQ=11)


                                                                  acggagGTCCTGTGCCGGCGGGGGGCCCAAGGGCCCTCG.GCCG.GCGCCCGCCCCGCCATCGGCGGGC.GCAGTCG.C.G.CCGCCGGGGCGGAGCGCAGGCctct                                    <  2:3209743/98‑1 (MQ=11)


                                                                   cAGAGTCCCTGTGCCCGCGGGGGCCCCAAGGCCCCTCG.GCCT.GCGCCCGCCCCGCCATCGGCGGGC.GCAGTCG.C.G.CCGCCGGGGCGGAGCGCAGGCCTCtg                                   <  2:3170014/101‑1 (MQ=11)


                                                                      gggTCCTGTGCCCGCCGAGGGCCAAAAGGCACGCG.GCCT.GCGCCCGCCCCGCCATCGGCGGGC.GCAGTCG.C.G.CCGCCGGGGCGGAGCGCAGGCCTCtgttg                                <  1:5753269/100‑1 (MQ=255)


                                                                      agGTCCTGTGACCGCCGAGGGCCCAAAGGCACTCG.GCCT.GCGCCCGCCCCGCCATCGGCGGGG.GCGGGCG.C.G.CCGCCCGGGGGGGGGGCGGGcccccgcct                                >  1:4081457/1‑94 (MQ=14)


                                                                      aagTCCCGTGACCGCCGGGGGCCCAAAGGCCCTCG.GCCT.GCGCCCGCCCCGCCATCGGCGGGG.CCCGCCG.C.C.CCCCCCGGGCGGGGGGCCGGccgcccgcg                                >  1:1603405/3‑94 (MQ=9)


                                                                       ggTCCTGTGCCCGCCGGGGGCCAAAGGGCCCCCG.CCCT.GCCCCCGCCCCGCCATCGGCGGGC.GCAGTCG.C.G.CCGCCGGGGCGGAGCGCAGGCCTCtgttgt                               <  1:3299197/101‑1 (MQ=14)


                                                                       ggTCCGTTGCCCGCGGGGGGCCCAAAGGCACTCG.GCGT.GCGCCCGCCCCGCCATCGGCGGGC.GCAGTCG.C.G.CCGCCGGGGCAGAGCGCAGGCCTCtgttgt                               <  1:4282088/101‑1 (MQ=14)


                                                                         tCCGGTGCCCGCGGGGGCCCCAAGGGCCCTCG.CCCT.GCGCCCGCCCCGCCATCGGCGGGC.GCAGTCG.C.G.CCGCCGGGGCGGAGCGCAGGCCTCTGTTGTGc                             <  2:2527379/101‑1 (MQ=11)


                                                                         cccTGTGCCCGCCGAGGGCCCAAAGGCACTCG.CCCT.GCCCCCGCCCCGCCATCGGCGGGC.GCAGTCG.C.G.CCGCCGGGGCGGAGCGCAGGCCTCTGTTGTGc                             <  1:6080426/100‑1 (MQ=255)


                                                                          ccTGTGCCCGCCGAGGGCCAAAGGGCACTCG.GCCT.GCGCCCGCCCCGCCATCGGCGGGC.GCAGTCG.C.G.CCGCCGGGGCGGAGCGCAGGCCTCTGTTGTGCt                            <  2:2677271/101‑1 (MQ=255)


                                                                             gtgACCGCCGAGGGCCCAAAGGCACTCG.GCCT.GCGCCCGCCCCGCCATCGGCGGGCGGCCGCCC.C.C.CCCCCGGGG.GGGGCGCCCGCCTCCCGTGTGCTTgg                         >  2:3835230/1‑101 (MQ=2)


                                                                              ggcccgccgAGGCCCCAAGGGCCCTCG.GCCT.GCGCCCGCCCCGCCATCGGCGGGC.GCAGTCG.C.G.CCGCCGGGGCGGAGCGCAGGCCTCTGTTGTGCTTggg                        <  1:2625106/98‑1 (MQ=11)


                                                                                     cgAGGGCCCAAAGGCACTCG.GCCT.GCGCCCGCCCCGCCATCGGCGGGC.GGAGGCG.C.G.GCGCCCGGGGGGGGGGCAGGACGCCGGTGTGCttttgggggggg                 >  2:4576022/1‑91 (MQ=9)


                                                                                        gggCCCAAAGGCCCCCGCGCCT.GCCCCCGCCCCGCCATCGGCGGGC.GCAGTCG.C.G.CCGCCGGGGCGGAGCGCAGGCCTCTGTTGTGCTTGGGGCAGGCGGt               <  2:2563789/101‑1 (MQ=14)


                                                                                         ggCCCAAAGGCACTCG.GCCT.GCGCCCGCCCCGCCATCGGCGGGG.GCAGGCGCC.C.CCCCCGGGGGGGGGGGCGGGCCCCGGTT.TGCTTtggggggggggggg             >  1:5942324/1‑91 (MQ=2)


                                                                                           gccAAAGGCCCCCC.GCCCGGCGCCCGCCCCGCCATCGGCGGGC.GCAGTCG.C.G.CCGCCGGGGCGGAGCGCAGGCCTCTGTTGTGCTTGGGGCAGGCGGTCGt            <  2:6721342/100‑1 (MQ=14)


                                                                                            ccAAAGGCCCCCC.GGCTGGCGCCCGCCCCGCCATCGGCGGGC.GCAGTCG.C.G.CCGCCGGGGCGGAGCGCAGGCCTCTGTTGTGCTTGGGGCAGGCGGTCGTg           <  2:921180/101‑1 (MQ=14)


                                                                                              aaaGGCCCTCG.GCCT.GCGCCCGCCCCGCCATCGGCGGGG.GGAATCG.C.T.CCGCCGGGGCGGGGGGGAGGGCTCTGTCGCGCCTGGGGGAGACgggctgtgac        >  2:2035101/1‑93 (MQ=2)


                                                                                                aGGCACTCG.GCCT.GCGCCCGCCCCGCCATCGGCGGGG.GGCCGCG.C.G.CCCCCCGGGCGGGGGGGCGGGCCCTCTTGTGCTTTGGGGGGGCggggcgggcgcc      >  1:80036/1‑91 (MQ=2)


                                                                                                aGGCACTCG.GCCT.GCGCCCGCCCCGCCATCGGCGGGC.GCCGGCG.C.GCCCGCCGGGG.GGGGCGCGGGGCGCTGGTTTTCGTTGGGTAGGCGGGGGCGCCGcc      >  2:2216928/1‑100 (MQ=2)


TGCGCTGCGGCTGATGAAGGACGTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCCAATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAAGGCACTCG.GCCT.GCGCCCGCCCCGCCATCGGCGGGC.GCAGTCG.C.G.CCGCCGGGGCGGAGCGCAGGCCTCTGTTGTGCTTGGGGCAGGCGGTCGTGCCGCAGCGC  >  NC_014034/448680‑448880

Base quality scores:


ATCG

< 3 ≤


ATCG

< 18 ≤


ATCG

< 29 ≤


ATCG

< 34 ≤


ATCG

< 41 ≤


ATCG