BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_014034	2,209,677	G→C	72.0%	intergenic (+9/‑262)	cobO → / → cobQ1	cob(I)yrinic acid a,c‑diamide adenosyltransferase/cobyric acid synthase CobQ

Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_014034	2,209,677	0	G→C	72.0%	115.6 ‑3.1 ‑0.1	149	intergenic (+9/‑262)	cobO/cobQ1	cob(I)yrinic acid a,c‑diamide adenosyltransferase/cobyric acid synthase CobQ
Reads supporting (aligned to +/- strand): new base (27/80): ref base (18/24): total (45/104)
Fisher's exact test strand distribution p-value = 4.70e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.93e-01

Only 100 of 156 total aligned reads displayed.


TGATCGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAAG..GC.C..G.GGCGCT.G.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGC  >  NC_014034/2209577‑2209765


tGATCGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGg..g                                                                                                >  1:3306499/1‑101 (MQ=255)


   ncGNNNCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG..GC.C..g.g                                                                                        <  1:1452416/100‑1 (MQ=255)


    cGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG..GC.C..g.gg                                                                                       <  2:3634568/101‑1 (MQ=255)


    cGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG..GC.C..g.gg                                                                                       >  2:1689248/1‑101 (MQ=255)


    cGAGGCGGCCGATATAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG..GC.C..g.gg                                                                                       >  2:3542076/1‑101 (MQ=255)


     gAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG..GC.C..g.ggg                                                                                      >  2:4211116/1‑100 (MQ=255)


      aGGCGGCCTATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG..GC.C..g.gggg                                                                                     <  2:4012140/101‑3 (MQ=255)


       gtcggcCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG..GC.C..G.Ggggc                                                                                    <  1:2168742/99‑1 (MQ=255)


       ggcggtCCATCTGGTGGCGGACATGGCTCGGGTGAAGCGCCCGATCCGGGAAGGCGTCGAGGGCCATCCGGGGGTCGAGTTCTGATCTC.AAG..GC.C..G.Ggggc                                                                                    <  2:3978187/101‑1 (MQ=255)


       ggcggcCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG..GC.C..G.Ggggc                                                                                    <  1:4275425/101‑1 (MQ=255)


         cggcCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG..GC.C..G.Ggggcgc                                                                                  <  1:2231869/101‑3 (MQ=255)


          ggcCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG..GC.C..g.gggggggt                                                                                 >  2:3831272/1‑95 (MQ=255)


           gcCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG..GC.C..g.ggggggctg                                                                                >  1:1727672/1‑94 (MQ=255)


           gcCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG..GC.C..G.Ggggcgctg                                                                                <  2:1220697/101‑5 (MQ=255)


             cGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG..GC.C..G.Ggggcgctgcc                                                                              <  2:3476819/101‑7 (MQ=255)


               aTCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G.GGCGCT.G.cccc                                                                            >  1:2729496/1‑101 (MQ=255)


                tCTAGTGGCCGAGATGGGTCAGGTGAAGCCCCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG..GC.CGGG.GGCGCT.G.ccccc                                                                           <  2:3209662/101‑1 (MQ=255)


                tCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG..GC.CGGG.GGCGCT.G.ccccc                                                                           <  2:515466/101‑1 (MQ=255)


                 cTAGTGACCGAGATTGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG..GC.CGGG.GGCGCT.G.CCCCCg                                                                          <  2:1299032/101‑1 (MQ=255)


                 cTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G.GGCGCT.G.CCCCCg                                                                          >  1:1443117/1‑101 (MQ=255)


                 cTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG..GC.CGGG.GGCGCT.G.CCCCCg                                                                          <  1:3521589/101‑1 (MQ=255)


                  tAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG..GC.CGGG.GGCGCT.G.CCCCCgg                                                                         <  1:1173166/101‑1 (MQ=255)


                     tGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G.GGCGCT.G.CCCCCGGAcc                                                                      >  2:1031567/1‑101 (MQ=255)


                      gACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G.GGGGCT.G.CCCCCCGAccc                                                                     >  1:1024402/1‑101 (MQ=255)


                      gACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CC.CG.G.GGCGCT.G.CCCCCAGAccc                                                                     >  2:576334/1‑101 (MQ=255)


                      gACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTC.AAG..GC.CGGG.GGCGCT.G.CCCCCGGAccc                                                                     <  1:1107137/101‑1 (MQ=255)


                        ccGAGATGGGTCAGGTGNAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G.GGCGCT.G.CCCCCGGAccccc                                                                   <  2:83088/101‑1 (MQ=38)


                        ccGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G.GGCGCT.G.CCCCCGGAccccc                                                                   <  1:3927381/101‑1 (MQ=38)


                        ccGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G.GGCGCT.G.CCCCCGGAccccc                                                                   <  1:1946486/101‑1 (MQ=38)


                          gagaTGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G.GGCGCT.G.CCCCCGGACCCCCgg                                                                 >  2:2848499/1‑101 (MQ=38)


                           agaTGGGTCAGGTGAAGCACCCGTTCCGGGCGGGCATCAACGGCCAGCCGGGGGTCGAGTTCTGATCTCGCGG..CCGG..G.GGCGCT.G.CCCCCGGACCCCCggg                                                                <  1:3972133/101‑1 (MQ=25)


                           agaTGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G.GGCGCT.G.CCCCCGGACCCCCggg                                                                >  2:693281/1‑101 (MQ=38)


                             aggggTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G.GGCGCT.G.CCCCCGGACCCCCGGGat                                                              <  1:3153411/99‑1 (MQ=38)


                             aTGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G.GGCGCT.G.CCCCCGGACCCCCGGGat                                                              <  2:1071005/101‑1 (MQ=38)


                             aTGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G.GGCGCT.G.CCCCCGGACCCCCGGGat                                                              >  2:1474867/1‑101 (MQ=38)


                             aTGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G.GGCGCT.G.CCCCCGGACCCCCGGGat                                                              <  2:1544796/101‑1 (MQ=38)


                             aTGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G.GGCGCT.G.CCCCCGGACCCCCGGGat                                                              <  2:2686132/101‑1 (MQ=38)


                                ggTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G.GGCGCT.G.CCCCCGGACCCCCGGGATAtt                                                           <  1:334735/101‑1 (MQ=38)


                                ggTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G.GGCGCT.G.CCCCCGGACCCCCGGGATAtt                                                           <  1:3148985/101‑1 (MQ=38)


                                 gTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G.GGCGCT.G.CCCCCGGACCCCCGGGATAttt                                                          >  1:673469/1‑101 (MQ=38)


                                  tCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CC.G..GTGGCGCT.G.CCCCCGGACCCCCGGGATATTtg                                                         <  1:3972310/101‑1 (MQ=38)


                                     ggTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G.GGCGCT.G.CCCCCGGACCCCCGGGATATTTGTgg                                                      >  2:1458145/1‑101 (MQ=37)


                                     ggTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G.GGCGCT.G.CCCCCGGACCCCCGGGATATTTGTgg                                                      <  2:1146295/101‑1 (MQ=37)


                                      gTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G.GGCGCT.G.CCCCCGGACCCCCGGGATATTTGTGGc                                                     <  2:1019952/101‑1 (MQ=37)


                                           gCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G.GGCGCT.G.CCCCCGGACCCCCGGGATATTTGTGGCAAGAt                                                >  1:1325037/1‑101 (MQ=25)


                                            cccccGTGCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G.GGCGCT.G.CCCCCGGACCCCCGGGATATTTGTGGCAAGATg                                               <  2:1559291/99‑1 (MQ=25)


                                            cACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G.GGCGCT.G.CCCCCGGACCCCCGGGATATTTGTGGCAAGATg                                               <  1:2035916/101‑1 (MQ=25)


                                            cACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G.GGCGCT.G.CCCCCGGACCCCCGGGATATTTGTGGCAAGATg                                               <  2:2087361/101‑1 (MQ=25)


                                            cACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G.GGCGCT.G.CCCCCGGACCCCCGGGATATTTGTGGCAAGATg                                               >  2:4350153/1‑101 (MQ=25)


                                            cACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G.GGCGCT.G.CCCCCGGACCCCCGGGATATTTGTGGCAAGATg                                               <  1:3567246/101‑1 (MQ=25)


                                            cACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CC.CG.G.GGGGCTGG.CCCCCCAGCCCCCGGGA.ATTTTTggggaaaag                                               >  1:241764/1‑94 (MQ=255)


                                             aCCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G.GGCGCT.G.CCCCCGGACCCCCGGGGTATTTTTGGCAAGATGa                                              >  1:376223/1‑101 (MQ=25)


                                             aCCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G.GGCGCT.G.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGa                                              >  1:2119024/1‑101 (MQ=25)


                                             aCCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G.GGCGCT.G.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGa                                              <  2:3365205/101‑1 (MQ=25)


                                             aCCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G.GGCGCT.G.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGa                                              <  2:44732/101‑1 (MQ=25)


                                              cccGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G.GGCGCC.GCCCCCCGGGCCCCCGGG.TATTTGTGGccaaaggga                                             >  2:1234133/1‑93 (MQ=255)


                                               ccGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGC.GG.G..G.GGCGCT.G.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGaaa                                            <  1:4077201/101‑1 (MQ=25)


                                               ccGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G.GGCGCT.G.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGaaa                                            <  1:1484392/101‑1 (MQ=25)


                                               ccGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G.GGCGCT.G.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGaaa                                            <  2:3506344/101‑1 (MQ=25)


                                               ccGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G.GGCGCT.G.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGaaa                                            <  1:951799/101‑1 (MQ=25)


                                                cGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G.GGCGCT.G.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAAc                                           <  2:2735672/101‑1 (MQ=25)


                                                   tCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G.GGCGCT.G.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACaga                                        <  2:2118887/101‑1 (MQ=25)


                                                    ccGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G.GGCGCT.G.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACagag                                       <  1:875770/101‑1 (MQ=25)


                                                       ggcaggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G.GGCGCT.G.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGt                                    >  1:773206/1‑101 (MQ=25)


                                                        gcaggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G.GGCGCT.G.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGtt                                   >  1:2201858/1‑101 (MQ=25)


                                                         caggcaTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G.GGCGCT.G.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTa                                  <  1:2993664/101‑1 (MQ=25)


                                                             caTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGGCCGG.G..G.GGGGCG.G.CCCCCGGCCCCCCGGG.TATTTTGGGGAAGATGAAAAAGAGCGTTAgggg                              >  1:2620111/1‑100 (MQ=255)


                                                              aTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G.GGGGGC.G.CCCCCCGGCCCCCCGGGGATTTTTGGCCAGATGAAACAGAGCGTTAggggt                             >  2:3275667/1‑99 (MQ=255)


                                                              aTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G.GGCGCT.G.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGAt                             >  2:506197/1‑101 (MQ=25)


                                                              aTCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G.GGCGCT.G.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGAt                             <  1:1009476/101‑1 (MQ=25)


                                                               tCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G.GGCGCT.G.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGAtt                            >  2:2478557/1‑101 (MQ=25)


                                                                cAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G.GGCGCT.G.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGAttt                           <  1:2142802/101‑1 (MQ=25)


                                                                cAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G.GGCGCT.G.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGAttt                           <  2:3061464/101‑1 (MQ=25)


                                                                cAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G.GGCGCT.G.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGAttt                           <  1:580958/101‑1 (MQ=25)


                                                                 aaGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G.GGCGCT.G.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTg                          <  1:2933886/101‑1 (MQ=25)


                                                                 aaGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G.GGCGCT.G.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTg                          <  1:3542324/101‑1 (MQ=25)


                                                                   gggCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G.GGCGCT.G.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGt                        <  2:861902/101‑1 (MQ=25)


                                                                   gccccACCCGGGGGTCG.GCT.TGATCTCAAGG..CCGG..G.GGCGCT.G.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTaa                      <  1:4389627/98‑1 (MQ=17)


                                                                    ggCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G.GGCGCT.G.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTa                       <  1:2597893/101‑1 (MQ=25)


                                                                    ggCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G.GGCGCT.G.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTa                       <  2:4321209/101‑1 (MQ=25)


                                                                     gCCAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G.GGCGCT.G.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTaa                      <  2:2053768/101‑1 (MQ=25)


                                                                      ccAGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G.GGCGCT.G.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAAc                     <  2:3150179/101‑1 (MQ=25)


                                                                      ccAGCCGGGGGTCGAGTTCTGATCTCAAGG..CC.CG.G.GGGGGA.T.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGGTAGGGATTTGGTAAc                     >  2:3589180/1‑101 (MQ=25)


                                                                        aGCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G.GGCGCT.G.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCa                   <  2:1023402/101‑1 (MQ=25)


                                                                         gCCGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G.GGCGCT.G.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCAt                  <  2:1606282/101‑1 (MQ=25)


                                                                          ccGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G.GGCGCT.G.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATg                 <  2:3339823/101‑1 (MQ=25)


                                                                           cGGGGGTCGAGTTCTGATCTCAAGG..CCGG..G.GGCGCT.G.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGa                <  1:2785271/101‑1 (MQ=25)


                                                                               ggTCGAGTTCTGATCTCAAGG..CCGG..G.GGCGCT.G.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCgg            <  2:2729306/101‑1 (MQ=25)


                                                                               ggTCGAGTTCTGATCTCAAGG..CCGG..G.GGCGCT.G.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCgg            <  1:1689359/101‑1 (MQ=25)


                                                                                  cGAGTTCTGATCTCAAGG..CCGG..G.GGCGCT.G.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCa         <  2:1835853/101‑1 (MQ=25)


                                                                                   gAGTTCTGATCTC.AAG..GC.CGGG.GGCGCT.G.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCaa        >  1:1316381/1‑101 (MQ=25)


                                                                                    aGTTCTGATCTCAAGG..CCGG..G.GGCGCT.G.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAg       <  1:2392618/101‑1 (MQ=25)


                                                                                     gTTCTGATCTC.AAG..GC.CGGG.GGGGCC.C.CCCCCGGGCCCCCGGGGGATTTGTGGCAAGATGCAACACAGCGCTAGGGGTTTTGGAAACCCGATCGGGCGAGc      >  2:1717024/1‑101 (MQ=255)


                                                                                     gTTCTGATCTC.AAG..GC.CGGG.GGCGCT.G.CCCCCCGACCCCCCGGGTATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAAACATGATCGGGCAAGc      >  1:4231082/1‑101 (MQ=21)


                                                                                      ttCTGATCTCAAGG..CCGG..G.GGCGCT.G.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCt     <  1:2848706/101‑1 (MQ=25)


                                                                                      ttCTGATCTCAAGG..CCGG..G.GGCGCT.G.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCt     <  2:593024/101‑1 (MQ=25)


                                                                                       tCTGATCTCAAGG..CCGG..G.GGCGCT.G.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTg    <  2:1335158/101‑1 (MQ=25)


                                                                                         tGATCTCAAGG..CC.GG.G.GGCGCT.G.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCCAGCTGGc  >  2:2062942/1‑101 (MQ=21)


                                                                                         tGATCTCAAGG..CC.GG.G.GGCGCT.G.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGc  >  2:4134759/1‑101 (MQ=25)


                                                                                         tGATCTCAAGG..CC.GG.G.GGCGCT.G.CCCCCGGAACCCCGGGGTAATTTTGGCAAGATGAAACAGAGCGTTAGGGGTTCGGGAACCATGAACGGGCAAGCTGGc  >  2:3685757/1‑101 (MQ=19)


TGATCGAGGCGGCCGATCTAGTGACCGAGATGGGTCAGGTGAAGCACCCGTTCCGGGCAGGCATCAAGGGCCAGCCGGGGGTCGAGTTCTGATCTCAAAG..GC.C..G.GGCGCT.G.CCCCCGGACCCCCGGGATATTTGTGGCAAGATGAAACAGAGCGTTAGGGATTTGGTAACCATGATCGGGCAAGCTGGC  >  NC_014034/2209577‑2209765

Base quality scores:


ATCG

< 3 ≤


ATCG

< 20 ≤


ATCG

< 29 ≤


ATCG

< 34 ≤


ATCG

< 41 ≤


ATCG