BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_014034	586,917	+G	100%	intergenic (‑145/+37)	RCAP_rcc00548 ← / ← purL	BolA family protein/phosphoribosylformylglycinamidine synthase II

Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_014034	586,917	1	.→G	100.0%	42.4	154	intergenic (‑145/+37)	RCAP_rcc00548/purL	BolA family protein/phosphoribosylformylglycinamidine synthase II
Reads supporting (aligned to +/- strand): new base (77/62): ref base (1/14): total (78/76)
Fisher's exact test strand distribution p-value = 2.47e-04
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 4.76e-01

Only 100 of 196 total aligned reads displayed.


GGCGACAAGAGCCAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGA.GCT.C.GGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCTGAGCTCTGCCAGCG  >  NC_014034/586817‑587012


ggCGACAAGAGCCAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGA.GCT.c                                                                                                    >  1:1344839/1‑101 (MQ=255)


ggCGACAAGAGCCAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGA.GCT.c                                                                                                    <  1:4753134/101‑1 (MQ=255)


  cGACAAGAGCCAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGA.GCT.C.gg                                                                                                 >  1:3578756/1‑101 (MQ=255)


     cAAGAGCCAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGA.GCT.C.ggggg                                                                                              >  2:4481370/1‑101 (MQ=255)


     cAAGAGCCAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGA.GCT.C.ggggg                                                                                              <  2:156175/101‑1 (MQ=255)


     cAAGAGCCAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGA.GCT.C.ggggg                                                                                              >  2:4590200/1‑101 (MQ=255)


     cAAGAGCCAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTAGCCGA.GCT.C.ggggg                                                                                              >  2:441165/1‑101 (MQ=255)


      aaGAGCCAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGA.GCT.C.gggggg                                                                                             >  2:1071009/1‑101 (MQ=255)


       agagCCAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGA.GCT.C.gggggg.g                                                                                           >  2:159493/1‑101 (MQ=255)


       agagCCAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGA.GCT.C.gggggg.g                                                                                           >  2:123026/1‑101 (MQ=255)


       agagCCAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGA.GCT.C.gggggg.g                                                                                           <  1:726524/101‑1 (MQ=255)


          gCCAAGCACCGGTTGCGTGTGCCAAGGGACGCAGCGGGCAGGGGCCGCGGCGAGGAGGGGGGTCGGGGCGTGGGCAGGTGTTGACCA.GGC.C.gggggg.ggcc                                                                                        >  2:2128407/1‑98 (MQ=255)


             aaGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGA.GCT.CGGGGGGG.G.CGCtt                                                                                     >  1:1254510/1‑101 (MQ=255)


             aaGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGA.GCT.C.gggggg.gggggtt                                                                                     >  1:3369407/1‑95 (MQ=255)


                   cGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGA.GCT.CGGGGGGG.G.CGCTTCTgccc                                                                               >  2:1875534/1‑101 (MQ=255)


                   cGGTGGCGTGGGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGA.GCT.CGGGGGGG.G.CGCTTCTgccc                                                                               <  2:2836786/101‑1 (MQ=255)


                     gTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGA.GCT.CGGGGGGG.G.CGCTTCTgcccgc                                                                             >  2:3614218/1‑101 (MQ=255)


                      ttGCGTGTGGCAGGGGGCGGAGGGGGCGGGGGCCGCGGCGAGGAGGGGGGTCAGGACGGGGGCAAGTGTCGCCGA.GCTGG.GGGGGG.G.CGCTTCTgcccgcc                                                                            <  2:4178882/101‑1 (MQ=255)


                       tGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGA.GCT.CGGGGGGG.G.CGCTTCTgcccgccc                                                                           >  1:2242259/1‑101 (MQ=255)


                       tGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGA.GCT.CGGGGGGG.G.CGCTTCTgcccgccc                                                                           >  1:180632/1‑101 (MQ=255)


                        gCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGA.GCT.CGGGGGGG.G.GGCCTCTGCCCGcccc                                                                          >  2:518372/1‑101 (MQ=255)


                          gtgtGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCATGACGGGGGCAAGTGTCGCCGA.GCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCgg                                                                        <  2:2738370/101‑1 (MQ=255)


                          gtgtGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGA.GCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCgg                                                                        >  1:3885661/1‑101 (MQ=255)


                          gtgtGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGA.GCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCgg                                                                        >  1:2653219/1‑101 (MQ=255)


                            gagCCAAGTGACGGCGCGGGCAGGGGCCGTGGCGAGGACGGGGGTCAGGAGGGGGGCAAGTGCCGCCGA.GCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAg                                                                      <  2:4714887/99‑1 (MQ=255)


                             tgCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGGCGGGGGTCAGGACGGGGGCAAGTGTCGCCGA.GCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGc                                                                     <  1:3049890/101‑1 (MQ=255)


                             ggCCAAGTGACGCAGCGGGCAGGGGCCGCGGGGAGGACGGGGGTCAGGACGGGGCCAAGTGTCGCCGA.GCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGc                                                                     <  2:1010416/100‑1 (MQ=255)


                              gCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGA.GCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGCt                                                                    <  1:3202009/101‑1 (MQ=255)


                              gCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGA.GCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGCt                                                                    <  2:2611827/101‑1 (MQ=255)


                               ccAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGA.GCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGCtt                                                                   >  1:861647/1‑101 (MQ=255)


                                cAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGA.GCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTg                                                                  >  1:4480678/1‑101 (MQ=255)


                                 aaGTGACGCAGCGGGGAGGGGCGGGGGCGGGGAGGGGGGTCAGGAGGGGGGCAAGTGTCGCCGA.GCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGc                                                                 <  2:132844/101‑1 (MQ=255)


                                  aGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGA.GCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCg                                                                >  1:4553238/1‑101 (MQ=255)


                                    tGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGA.GCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGtt                                                              >  2:878405/1‑101 (MQ=255)


                                    tGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGA.GCT.C.GGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGttt                                                             >  1:447073/1‑101 (MQ=255)


                                      aCGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGA.GCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGtttt                                                            >  2:4050074/1‑101 (MQ=255)


                                        gcggggggCAGGGGCCGGGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGA.GCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCa                                                          <  2:3493656/96‑1 (MQ=255)


                                        gCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGA.GCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCa                                                          <  2:2362169/101‑1 (MQ=255)


                                          cgCGGGCCGGGGCCGGGGGGGGGCGGGGGGTGCGGAGGGGGGAATGTGTCCCCGG.GCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGc                                                        <  1:4490447/100‑1 (MQ=255)


                                            cGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGA.GCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCt                                                      <  2:2211465/101‑1 (MQ=255)


                                             gggCGGGGGCGGGGGCGGGGAGGGGGGTCGGGACGGGGGCAAGTGTCGCCGA.GCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCtt                                                     <  1:4418578/101‑1 (MQ=255)


                                             gggCGGGGGCCGGGGGGGGGAGGGGGGTCAGGGCGGGGGCAAGTGTCGCCGA.GCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCtt                                                     <  2:3337748/101‑1 (MQ=255)


                                             gggCGGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGA.GCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCtt                                                     <  2:2862114/101‑1 (MQ=255)


                                             gggCAGGGGCCGCGGGGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGA.GCTGG.GGGGGG.G.CGCTTCTGCCCGCCCCGGAGCCTGCGTTTTCAGCCtt                                                     <  1:4000470/101‑1 (MQ=255)


                                              ggggggggCCGCGGGGGGGAGGGGGGTGAGGGGAGGGGTAATTGTCGCCGAGGGG.G.GGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTc                                                    <  1:4070250/97‑1 (MQ=255)


                                              ggCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGA.GCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTc                                                    <  2:2719660/101‑1 (MQ=255)


                                               gggggggCCGGGGGGAGGCGGGGGGTGGGGAGGGGGGCAGGTGTCGCCGA.GCG.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCg                                                   <  1:4530309/98‑1 (MQ=255)


                                               gCAGGGGCGGCGGGGAGGACGGGGGTCAGGAGGGGGGCAAGTGTCGCCGA.GCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCg                                                   <  2:2475290/101‑1 (MQ=255)


                                                ggggggCGGCGGCGAGGGGGGGGGTGGGGACGGGGGCAAGTGTCGCCGA.GCG.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGa                                                  <  1:254783/99‑1 (MQ=255)


                                                gaGGGGCCGCGGCGGGGAGGGGGGTCAGGAGGGGGGCAAGTGTCGCCGA.GCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGa                                                  <  1:2092073/100‑1 (MQ=255)


                                                 gggggCGGCGGGGGGGCGGGGGGTCAGGACGGGGGCAAGTGGCGCCGA.CCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGAc                                                 <  1:4671446/100‑1 (MQ=255)


                                                 aGGGGCCGCGGCGAGGACGGGGGTTAGGACGGGGGCAAGTGTCGCCGA.GCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGAc                                                 >  2:4466875/1‑101 (MQ=255)


                                                 aGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGA.GCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGAc                                                 >  1:2833168/1‑101 (MQ=255)


                                                 aGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGA.GCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGAc                                                 >  1:3497379/1‑101 (MQ=255)


                                                  ggggCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGA.GCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGAcg                                                >  1:1615091/1‑101 (MQ=255)


                                                   gggCCGGGGGGGGGAGGGGGGGGAGGCCGGGGGCAGGTGTCGCCGA.GCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGAcgg                                               <  1:3614437/101‑1 (MQ=255)


                                                      ccGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGA.GCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGAcggcgg                                            >  1:249379/1‑101 (MQ=255)


                                                       cgcgGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGA.GCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGAcggcggc                                           >  1:4249912/1‑101 (MQ=255)


                                                        gcggcgAGGACGGGGGTGAGGACGGGGGCAAGTGTCGCCGA.GCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGAcggcggcg                                          <  1:2334732/101‑1 (MQ=255)


                                                        gcggcgAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGA.GCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGAcggcggcg                                          <  2:587918/101‑1 (MQ=255)


                                                         ggggggggAAGGGGGGCACGGGGGGAGCAAGCGTCACCGA.GCT.G.GGGGGG.GCCGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGAcggcggcgg                                         <  1:1099309/95‑1 (MQ=255)


                                                           gggggggcgcgggTCAGGACGGGGGCAAGTGTCGCCGA.GCTGG.GGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGgcggcggcg                                       <  1:4600111/91‑1 (MQ=255)


                                                            cgAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGA.GCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGa                                      >  2:1378787/1‑101 (MQ=255)


                                                            cgAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGA.GCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGa                                      <  2:249354/101‑1 (MQ=255)


                                                            cgAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGA.GCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTGGAGCCTTCGACGGCGGCGGCGa                                      >  1:1564810/1‑101 (MQ=255)


                                                             gAGGGCGGGGGTCAGGACGGGGGCAAGTGTCGCCGA.GCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGaa                                     <  1:3746477/101‑1 (MQ=255)


                                                              aGGACGGGGGTCGGGACGGGGCCAAGCGTCGCCGC.GGG.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAg                                    <  1:3811747/101‑1 (MQ=255)


                                                              aGGACGGGGGTCAGGACGGGGCCAAGTGTCGCCGA.GCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAg                                    <  2:4552618/101‑1 (MQ=255)


                                                                gACGGGGGTCGGGACGGGGGCAAGTGTCGCCGA.GCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGgc                                  <  1:1882872/101‑1 (MQ=255)


                                                                gACGGGGGTCAGGACGGGGGCAAGTGTCGCCGA.GCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGgc                                  <  2:478554/101‑1 (MQ=255)


                                                                 aggggggTCAGGACGGGGGCAAGTGTCGCCGA.GCG.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGgcg                                 <  1:4259399/99‑1 (MQ=255)


                                                                 aCGGGGGTCAGGACGGGGGCAAGTGTCGCCGA.GCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGGAGgcg                                 >  2:1619982/1‑101 (MQ=255)


                                                                    ggggggaGGGCGGGGGCAGGTGTCGGCGA.GGG.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGgcggtg                              <  1:1052460/95‑1 (MQ=255)


                                                                    ggggTCAGGACGGGGGCAAGTGTCTCCGA.GCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGgcggtg                              <  2:2013140/101‑1 (MQ=255)


                                                                    ggggTCAGGACGGGGGCAAGTGTCGCCGA.GCTGG.GGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGgcggtg                              <  2:392381/101‑1 (MQ=255)


                                                                    ggggTCAGGACGGGGGCAAGTGTCGCCGA.GCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGgcggtg                              >  2:3757682/1‑101 (MQ=255)


                                                                         ctggACGGGGGGAGGTGTCGCCGA.GCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTa                         <  1:681073/99‑1 (MQ=255)


                                                                          aGGACGGGGGCAAGTGTCGCCGG.GCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTaa                        <  1:1516071/101‑1 (MQ=255)


                                                                          aGGACGGGGGCAAGTGTCGCCGA.GCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTaa                        <  2:2058362/101‑1 (MQ=255)


                                                                            ggggggggCAAGTGTCGCCGG.GCTGG.GGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAg                      <  1:3699318/98‑1 (MQ=255)


                                                                            gaggggggCAAGTGTCGCCGA.GCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAg                      <  1:2002114/98‑1 (MQ=255)


                                                                            gACGGGGGCAAGTGTCGCCGA.GCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAg                      >  2:3783799/1‑101 (MQ=255)


                                                                             aCGGGGGCAAGTGTCGCCGA.GCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGGGCGGTAAAGc                     >  1:1203952/1‑101 (MQ=255)


                                                                               gggggGAGGTGGCGCCGG.GCTGG.GGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGt                   <  1:546349/101‑1 (MQ=255)


                                                                               gggggCAAGTGTCGCCGA.GCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGt                   >  2:738527/1‑101 (MQ=255)


                                                                               gggggCAAGTGTCGCCGA.GCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGt                   >  2:2791046/1‑101 (MQ=255)


                                                                               ggggaCAAGTGTCGCCGA.GCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGGAAAGCgg                   >  2:4020207/1‑100 (MQ=255)


                                                                                ggggCAAGTGTCGCCGA.GCTGG.GGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCATg                  <  1:2711646/101‑1 (MQ=255)


                                                                                ggggCAAGTGTCGCCGA.GCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTg                  >  2:943802/1‑101 (MQ=255)


                                                                                ggggCAAGTGTCGCCGA.GCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTg                  >  2:4134180/1‑101 (MQ=255)


                                                                                  ggCAAGTGTCGCCGA.GCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCt                >  2:387974/1‑101 (MQ=255)


                                                                                     aaGTGTCGCCGA.GCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCTGaa             <  2:2242135/101‑2 (MQ=255)


                                                                                     aaGTGTCGCCGA.GCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCTGAg             >  1:488573/1‑101 (MQ=255)


                                                                                      aGTGTCGCCGA.GCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCTGAGc            >  1:4398635/1‑101 (MQ=255)


                                                                                       gtggCGCCGA.GCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCTGAGct           <  1:1540182/101‑1 (MQ=255)


                                                                                         gtCGCCGA.GCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCACCGAAGGCGGTGCGGTAAAGCGTGCTGAGctct         >  1:4106803/1‑101 (MQ=255)


                                                                                          gcGCCGA.GCG.C.GGGGGGNG.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCTGAGCTCTg        <  2:1819156/100‑1 (MQ=255)


                                                                                             ccGA.GCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCTGAGCTCTGCCa     >  2:2924166/1‑101 (MQ=255)


                                                                                              cGA.GCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTGGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCTGAGCTCTGCCAg    >  2:3055332/1‑101 (MQ=255)


                                                                                                a.GCT.CGGGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCTGAGCTCTGCCAgcg  >  1:1433496/1‑101 (MQ=255)


GGCGACAAGAGCCAAGCACCGGTTGCGTGTGCCAAGTGACGCAGCGGGCAGGGGCCGCGGCGAGGACGGGGGTCAGGACGGGGGCAAGTGTCGCCGA.GCT.C.GGGGGG.G.CGCTTCTGCCCGCCCCGGAGCTTGCGTTTTCAGCCTTCGACGGCGGCGGCGAAGGCGGTGCGGTAAAGCGTGCTGAGCTCTGCCAGCG  >  NC_014034/586817‑587012

Base quality scores:


ATCG

< 3 ≤


ATCG

< 12 ≤


ATCG

< 22 ≤


ATCG

< 32 ≤


ATCG

< 40 ≤


ATCG