BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_014034	1,965,279	+C	100%	intergenic (+47/+45)	bdhA → / ← RCAP_rcc01818	3‑hydroxybutyrate dehydrogenase/hypothetical protein

Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_014034	1,965,279	1	.→C	100.0%	22.4	108	intergenic (+47/+45)	bdhA/RCAP_rcc01818	3‑hydroxybutyrate dehydrogenase/hypothetical protein
Reads supporting (aligned to +/- strand): new base (40/58): ref base (9/1): total (49/59)
Fisher's exact test strand distribution p-value = 4.97e-03
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.60e-07

Only 100 of 154 total aligned reads displayed.


CGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTGCGCCTC.GGGGCAAGCCC.GAGG.CG.ACCCG.CC..GCGGG.CT.C.CG.CCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGTCGGC  >  NC_014034/1965180‑1965374


cgcAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTC.GGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.CG.C.CG.c                                                                                                >  1:7583846/1‑101 (MQ=255)


  cAGGTCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTC.GGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.CG.C.CC.ccc                                                                                              >  1:7229250/1‑101 (MQ=38)


   aGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTC.GGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.GG.C.CC.cccc                                                                                             >  1:4735328/1‑101 (MQ=38)


   aGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTC.GGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.CT.C.CG.cccc                                                                                             >  2:3442388/1‑101 (MQ=255)


       cACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTC.GGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.GG.C.CC.cccccccg                                                                                         >  2:410114/1‑97 (MQ=37)


        aCCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTC.GGGGCAAGCCC.GAGG.CG.ACCCG.CC..CGGGG.GG.C.CG.cccccccgc                                                                                        >  2:2143392/1‑96 (MQ=37)


            gCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTC.GGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.GG.C.CC.ccccccggccccc                                                                                    >  2:6821550/1‑92 (MQ=25)


               cgacgaTTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTC.GGGGCAAGCCC.GAGG.CG.ACCCG.CC..CGGGG.CT.CGGG.CCCCGGTCCCCCGAg                                                                                  >  2:501321/1‑101 (MQ=37)


                gacgaTTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTC.GGGGCAAGCCC.GAGG.CG.ACCCG.CC..CGGGG.GT.C.CGCCCCCGGACCCCCGAgg                                                                                 >  2:5963352/1‑101 (MQ=37)


                   gaTTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTTGCCTC.GGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.GT.C.CGCCCCCGGGGCCCCGAGGGAt                                                                              >  1:2989998/1‑101 (MQ=25)


                    aTTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTC.GGGGCAAGCCC.GAGG.CG.ACCCG.CC..CGGGG.CT.C.CGCCCCCCGACCCCCGAGGGAtt                                                                             >  1:2802488/1‑101 (MQ=25)


                     tttCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTC.GGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.GTCC.CC.CCCCCGGCCCCCCAGGGAttt                                                                            >  1:5936307/1‑101 (MQ=21)


                      ttCTGCGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTC.GGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.CTCC.CC.CCCCGGACCCCCGAgggagtgc                                                                           >  2:8970865/1‑95 (MQ=21)


                      tggtggggggggcggcGGGCGGGCGTTGGGAGAGATTCCGCCGG.GGGGCACCCCC.GGGG.CC.CCCCG.CC..GGGGG.CT.C.CGCCCCCGGACCCCCGAGGTATTTc                                                                           <  2:7221878/91‑1 (MQ=9)


                       tCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTC.GGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.GT.C.CGCCCCCCGGCCCCCCAGGGATTTCa                                                                          >  2:3288021/1‑101 (MQ=21)


                        gtggggggggcggcGGCGCGGCTCGGTGAGAGTTCGCGCCGG.GGGCCAGCCCC.GGGC.CC.CCCCG.CC..GGGGG.CT.C.CGCCCCCGGACCCCCGAGGTATTTCAg                                                                         <  2:5936072/93‑1 (MQ=2)


                        cTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTC.GGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.CTCC.CC.CCCCGGACCCCCGAGGGATTTccg                                                                         >  2:4280292/1‑99 (MQ=25)


                         tgtgGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTC.GGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.GT.T.CGCCCCCCGGCCCCCCCGGTATTTTAgg                                                                        >  2:6616805/1‑101 (MQ=18)


                          gtgGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTC.GGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.GTCC.CC.CCCCCGGCCCCCCAGGGGTTTCAGGc                                                                       >  1:5802370/1‑101 (MQ=21)


                          gtgGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTC.GGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.CTCC.CC.CCCCCGACCCCCGCGGGATTTCCGGc                                                                       >  1:2249429/1‑101 (MQ=21)


                           tgGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTC.GGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.GT.C.CGCCCCCGGACCCCCCGGGGATTTCAGGCa                                                                      >  2:5636575/1‑101 (MQ=21)


                           tgGATGGCGGCTGGACGGCGCTGTGAGAGATTTCGCCTC.GGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.CC.C.CGCCCCCCGGCCCCCCGGGGATTTCCGGCa                                                                      >  1:7108094/1‑101 (MQ=18)


                             gggggcggcGGGCGGGCGTGGTGAGGGTTTCCGCCGG.GGGGCACCCCC.GGGC.CC.CCCCG.CC..GGGGG.CT.C.TGCCCCCGGACCCCCGAGGTATTTCAGGCAAg                                                                    <  1:6107488/98‑1 (MQ=2)


                               gggcggcGGGCGGGGGTTGTGAGAGTTTCCCCCGG.GGGGCACGCCG.GGGG.CC.CCCCG.CC..GGGGG.CT.C.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGAt                                                                  <  2:5539946/100‑1 (MQ=9)


                                ggcggcTGGACGGCGCTGTGAGAGATTTCGCCTC.GGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.CT.C.CGCCCCCGGACCCCCGAGGGATTTCAGGGCAGATg                                                                 >  2:6675437/1‑101 (MQ=25)


                                ggcggcGGGGGGGCGGGGTGAGAGTTCCCGCCGG.GGGCCAGCCCC.GAGG.CG.CCCCG.CC..GGGGG.CT.C.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATg                                                                 <  2:7690983/101‑1 (MQ=2)


                                   ggcTGGACGGCGCTGTGAGAGATTTCGCCTC.GGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.CT.C.CGCCCCCGGACCCCCGAGGGATTTCAGGCAAGATGaaa                                                              >  1:1299692/1‑101 (MQ=21)


                                    gggggcgggCGGTGGGAGAGTTTCCGCCGG.GGGCCAACCCG.GGGC.CC.CCCCG.CC..GGGGG.CT.C.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAag                                                             <  2:4718566/94‑1 (MQ=2)


                                    gcgggccGGCGCTGTGAGAGATCCCGCCTC.GGGGCACGCCC.GGGC.CC.CCCCG.CC..GGGTG.CT.C.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAag                                                             <  2:3196919/95‑1 (MQ=16)


                                    gcTGGACGGCGCTGTGAGAGATTTCGCCTC.GGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGGCCC.C.CC.CCCCCGGCCCCCCCGGGATTTTCGGCAAGATGaaaa                                                             >  1:621864/1‑100 (MQ=12)


                                     cTGGACGGCGCTGTGAGAGATTTCGCCTC.GGGGCAAGCCC.GAGG.CG.ACCCG.CCGGGGGGG.CC.C.CC.CCCCGGCCCCCCCGGGGTTTTCAGCGAAGAGGaaaa                                                             >  1:2513049/1‑100 (MQ=12)


                                      tGGACGGCGCTGTGAGAGATTTCGCCTC.GGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.CT.C.CGCCCCCGGACCCCCGAGGGATTTCAGGCAAGATGAAagag                                                           >  2:5093515/1‑101 (MQ=21)


                                        gggggCGGTGTGAGAG.TTCCCCCGCGGGGGCACCCCC.GAGC.CC.CCCCG.CC..GGGGG.CT.C.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAggg                                                         <  2:3614019/98‑1 (MQ=9)


                                            gatcTGTGAGAGATTTCGCCTC.GGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.GTCC.CC.CCCCCGACCCCCGAGGGATTTCAGGCCAGATTGAAGAGGGgggc                                                     >  1:8543887/4‑101 (MQ=16)


                                              gcTGTGAGAGATTTCGCCTC.GGGGCAAGCCC.GAGG.CG.ACCCG.CCGGGGGGG.GC.C.CC.CCCCGCGCCCCCGCGGGGTTTCAGCCAGGAGGAAAAAGGGgggcc                                                    >  2:2069025/1‑100 (MQ=9)


                                              gcTGTGAGAGATTTCGCCTC.GGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.GTCC.CC.CCCCCGGCCCCCGCGGGATTTCAGGCAAGAAGGAAGAGGGgcgggc                                                   >  2:3458240/1‑98 (MQ=16)


                                                tgtgAGAGATTTCGCCTC.GGGGCAAGCCC.GAGG.CG.ACCCGCCC..GGGGG.CC.C.CC.CCCCCGGCCCCCCCGGGGTTTTAAGGGAGAAGGAAGAGGGGGGCTCAg                                                 >  1:5986250/1‑101 (MQ=9)


                                                 gtgAGAGATTTCGCCTC.GGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.GT.C.CGCCCCCGGGCCCCCCAGGGATTTTAGGCAAGATGAAAGAggggggggcagg                                                >  2:7787975/1‑93 (MQ=18)


                                                  tgAGAGATTTCGCCTC.GGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.CT.C.CGCCCCCGGACCCCCGCGGGATTTCAGGCAAGATGAAAGAGGGGCGCTCAggg                                               >  2:8295888/1‑100 (MQ=21)


                                                  tgAGAGATTTCGCCTC.GGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.CT.C.CGCCCCCGGACCCCCGAGGGATTTCAGGCAAGATGAAAGAGGGGCGCTCCggg                                               >  2:5277745/1‑100 (MQ=21)


                                                    agagaTTTCGCCTC.GGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.GG.C.CGCCCCCCGCCCCCCGAGGAATTTCAGGCAAGATGAAAGAGGGGCGGGCAggggg                                             >  2:115514/1‑101 (MQ=18)


                                                    agagaTTTCGCCTC.GGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.CT.C.CGCCCCCGGACCCCCGAGGGATTTCAGGCAAGATGAAAGAGGGGCGCTCAggcgg                                             >  2:1545964/1‑101 (MQ=21)


                                                    agagaTTTCGCCTC.GGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.CT.C.CGCCCCCGGACCCCCCAGGGATTTCAGGCAAGATGAAAGAGGGGCGCTCCggcgg                                             >  2:8234648/1‑101 (MQ=21)


                                                      agaTTTCGCCTC.GGGGCAAGCCC.GAGG.CG.ACCCG.CCG.GGGGG.CC.C.CC.CCCCCGGCCCCCCCGGGATTTTAGGCAAAATGAAAGAGGGGGGCTCAGGCGGcc                                           >  1:963026/1‑101 (MQ=16)


                                                        aTTTCGCCTC.GGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.CT.C.CGCCCCCGGACCCCCCAGGGATTTCAGGCAAGATGAAAGAGGGGCGCCCAGGCGGCctt                                         >  1:5798850/1‑101 (MQ=21)


                                                        aTTTCGCCTC.GGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.CT.C.CGCCCCCCGGCCCCCCAGGGCTTTCGGGCAAGATGAAAAAGGGGGGCTCAGGCGGcccc                                         >  1:8385315/1‑99 (MQ=16)


                                                         tctccgcgc.GGGGCAACCCCGGCGG.CC.CCCCG.CC..GGGGG.CT.C.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCctt                                         <  2:8669000/93‑1 (MQ=18)


                                                           tcgcgcggggggCCGGCCC.GGGG.CGCCCCCG.CC..GGGGG.CT.C.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCcttc                                        <  2:5476407/93‑1 (MQ=18)


                                                            cgcGTC.GGGGCAAGCGC.GAGG.CG.ACCCG.CC..GGGGG.CT.C.CGCCCCCGGACCCCCGAGGGATTTCAGGCAAGATGAAAGAGGAGCGCTCAGGCGGCCttcttc                                     >  1:8825620/1‑101 (MQ=21)


                                                             gcCTC.GGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.CT.C.CGCCCCCGGACCCCCGACGGATTTCCGGCAAGAAGAAAGAGGGGCGCTCAGGCGGGCTTCTTcg                                    >  1:8972582/1‑101 (MQ=21)


                                                              gccc.GGGGCACCCCCGGACG.CC.ACCCG.CC..GGGGG.CT.C.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTcg                                    <  1:5419641/98‑1 (MQ=18)


                                                              ccgc.GGGGCCACCCCGGAGC.CG.CCCCG.CC..GGGGG.CT.C.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTcg                                    <  2:8094216/98‑1 (MQ=18)


                                                                ccgggggCCAGCCC.GGGGCCC.CCCCG.CC..GGGGG.CT.C.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTcgc                                   <  2:5360321/98‑1 (MQ=21)


                                                                 c.GGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.GT.C.CGCCCCCCGACCCCCCAGGGATTTCCGGCAAGATGGAAGAGGGGGGCTCAGGCGGCCTTTCTcgccgc                                >  2:8142677/1‑101 (MQ=21)


                                                                 c.GGGGCAAGCCC.GAGG.CG.ACCCG.CC..GGGGG.CT.C.CGCCCCCCGGCCCCCCCGGGATTTCAGGCCAGATGAAAGAGGGGCGCTCAGGCGGCCCTCTTcgccgc                                >  1:1493935/1‑101 (MQ=21)


                                                                        aaGCCC.GAGG.CG.ACCCG.CC..GGGGG.CTCC.CC.CCCCCGGCCCCCGAGGGATTTCAGGCAAGATGAAAGAGGGGCGCTCCGGCGGCCTTCTTCGCCGCCCTGcc                          >  2:8364642/1‑101 (MQ=21)


                                                                        aaGCCC.GAGG.CG.ACCCG.CC..GGGGG.CT.C.CGCCCCCCGAGCCCCCCGGTATTTCAGGCAAGATGAAAGAGGGGCGGGCGGGCGGGCTTCCTCcccccccgtcc                          >  2:2834626/1‑96 (MQ=9)


                                                                           ccc.GAGG.CG.ACCCG.CC..GGGGG.CTCC.CC.CCCCCGACCCCCCAGGGATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCCTGccccg                       >  2:7044733/1‑99 (MQ=21)


                                                                             c.GAGG.CG.ACCCG.CC..GGGGG.GTCC.CC.CCCCCGACCCCCCAGGGAATTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCCTGCCCgggg                     >  2:5625296/1‑99 (MQ=21)


                                                                             c.GAGC.CG.CCCCG.CC..GGGGG.CT.C.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGgtg                     <  2:3202607/101‑1 (MQ=25)


                                                                               gAGG.CG.ACCCG.CC..GGGGG.CT.C.CGCCCCCGGACCCCCGAGGGATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGgtgt                    >  2:7238006/1‑101 (MQ=25)


                                                                               cggggcccccccG.CC..GGGGG.CT.C.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGgt                      <  2:1627820/92‑1 (MQ=25)


                                                                               ccggccccccccG.CC..GGGGG.CT.C.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGgt                      <  1:1943294/92‑1 (MQ=25)


                                                                               ccgcggccccccG.CC..GGGGG.CT.C.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGgt                      <  1:3687327/92‑1 (MQ=25)


                                                                                aGG.CG.ACCCG.CC..GGGGG.GT.C.CGCCCCCCGACCCCCGAGGGATTTCCGGCAAGATGAAAGAGGGGCGCCCAGGCGGCCCTCTTCGCCGCCCTGCCCCgtggc                   >  2:2894662/1‑99 (MQ=18)


                                                                                  gggccccccG.CC..GGGGG.CT.C.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGtc                   <  2:4123805/95‑1 (MQ=25)


                                                                                   gccgccccG.CC..GGGGG.CT.C.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGtct                  <  2:8972161/96‑1 (MQ=25)


                                                                                   gcccccccG.CC..GGGGG.CT.C.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGtct                  <  2:8754701/96‑1 (MQ=25)


                                                                                   ccccccccG.CC..GGGGG.CT.C.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGtct                  <  2:1995825/96‑1 (MQ=25)


                                                                                    cccccccG.CC..GGGGG.CT.C.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGtctc                 <  1:7277124/97‑1 (MQ=25)


                                                                                    cccccccG.CC..GGGGG.CT.C.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGtctc                 <  1:5406111/97‑1 (MQ=25)


                                                                                    cccccccG.CC..GGGGG.CT.C.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGtctc                 <  1:5277929/97‑1 (MQ=25)


                                                                                     ccccccG.CC..GGGGG.CT.C.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCg                <  1:8056243/98‑1 (MQ=37)


                                                                                     ccccccG.CC..GGGGG.CT.C.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCg                <  2:7107790/98‑1 (MQ=37)


                                                                                     ccccccG.CC..GGGGG.CT.C.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCg                <  2:8221958/98‑1 (MQ=37)


                                                                                        cccG.CC..GGGGG.CT.C.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGc             <  1:1662958/101‑1 (MQ=255)


                                                                                         ccg.cc..gGGGG.CT.C.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGct            <  2:7401024/101‑1 (MQ=255)


                                                                                          cg.cc..ggGGG.CT.C.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGctt           <  2:3738524/101‑1 (MQ=255)


                                                                                          cg.cc..ggGGG.CT.C.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGctt           <  1:8752861/101‑1 (MQ=255)


                                                                                           g.cc..ggGGG.CT.C.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGcttc          <  2:3831917/101‑1 (MQ=255)


                                                                                           g.cc..ggGGG.CC.C.CGCCCCCGGACCCCCGAGGGATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCCTCTTCGCCGCCTTGCCCGGGGTCTCGGGgttc          >  2:6491933/1‑101 (MQ=25)


                                                                                                cggggg.CT.C.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGcttctt        <  1:2894760/98‑1 (MQ=255)


                                                                                                cggggg.CT.C.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGcttctt        <  2:5798615/98‑1 (MQ=255)


                                                                                                cggggg.CT.C.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGcttctt        <  1:1546011/98‑1 (MQ=255)


                                                                                                cggggg.CT.C.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGcttctt        <  2:7356430/98‑1 (MQ=255)


                                                                                                cggggg.CT.C.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGGCTCGGGcttctt        <  1:8096862/98‑1 (MQ=37)


                                                                                                 ggggg.CT.C.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTg       >  1:7268734/3‑101 (MQ=255)


                                                                                                 ggggg.CT.C.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTg       <  1:1135582/99‑1 (MQ=255)


                                                                                                 ggggg.CT.C.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCCCCTTGCCCGGGGCCTCGGGCTTCTTg       >  1:5169890/3‑101 (MQ=37)


                                                                                                  gggg.CT.C.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGt      >  2:1165664/2‑101 (MQ=255)


                                                                                                  gggg.CT.C.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGt      >  1:4636400/2‑101 (MQ=255)


                                                                                                  gggg.CT.C.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGt      >  1:3924748/2‑101 (MQ=255)


                                                                                                   ggg.CT.C.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGTc     <  2:4593225/101‑1 (MQ=255)


                                                                                                   ggg.CT.C.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGTc     <  2:7706394/101‑1 (MQ=255)


                                                                                                    gg.CT.C.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGTCg    >  1:6171899/1‑101 (MQ=255)


                                                                                                     g.CT.C.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGTCgg   >  1:937916/1‑101 (MQ=255)


                                                                                                     g.CT.C.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGTCgg   >  1:2512145/1‑101 (MQ=255)


                                                                                                     g.CT.C.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGTCgg   <  2:904870/101‑1 (MQ=255)


                                                                                                       cT.C.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGTCgtc  >  1:7758973/1‑99 (MQ=255)


                                                                                                       cT.C.CGCCCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGTCGGc  >  1:2111976/1‑101 (MQ=255)


CGCAGATCACCGGCACGACGATTTCTGTGGATGGCGGCTGGACGGCGCTGTGAGAGATTGCGCCTC.GGGGCAAGCCC.GAGG.CG.ACCCG.CC..GCGGG.CT.C.CG.CCCCGGACCCCCGAGGTATTTCAGGCAAGATGAAAGAGGGGCGCTCAGGCGGCCTTCTTCGCCGCCTTGCCCGGTGTCTCGGGCTTCTTGTCGGC  >  NC_014034/1965180‑1965374

Base quality scores:


ATCG

< 3 ≤


ATCG

< 15 ≤


ATCG

< 27 ≤


ATCG

< 34 ≤


ATCG

< 41 ≤


ATCG