BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_014034	3,609,126	T→A	1.2%	intergenic (+65/‑106)	clpB → / → dinB	chaperone ClpB/DNA polymerase IV

Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_014034	3,609,126	0	T→A	1.2%	3.1 ‑1.0 ‑0.0	419	intergenic (+65/‑106)	clpB/dinB	chaperone ClpB/DNA polymerase IV
Reads supporting (aligned to +/- strand): new base (4/1): ref base (210/202): total (214/205)
Fisher's exact test strand distribution p-value = 3.73e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.97e-01

Only 100 of 441 total aligned reads displayed.


CGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGT.GGGGGCTCTG.CC.CCC..GCTTCGCTCCCCCCGGGATATTT.G.G.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGA.GGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGTTAGAATTG  >  NC_014034/3609029‑3609225


cgcgTCCAGACCGAACCNCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.G.GCAAGGTAAAAGCAtt                                                                                                      <  2:1977871/101‑1 (MQ=37)


cgcgTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.G.GCAAGGTAAAAGCAtt                                                                                                      <  2:4655672/101‑1 (MQ=38)


 gcgTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.G.GCAAGGTAAAAGCATTg                                                                                                     <  2:4329663/101‑1 (MQ=38)


   gTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.G.GCAAGGTAAAAGCATTGcc                                                                                                   <  2:5506765/101‑1 (MQ=38)


    tCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.G.GCAAGGTAAAAGCATTGccc                                                                                                  <  2:2189342/101‑1 (MQ=38)


      cAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCC.CCC..GCTTCGCTCCCCCCGGGATATTT.G.GAGCAAGGTAAAAGCATTGACCtt                                                                                                >  2:1693283/1‑101 (MQ=37)


       aGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.G.GCAAGGTAAAAGCATTGCCCttt                                                                                               <  1:4338362/101‑1 (MQ=38)


       aGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.G.GCAAGGTAAAAGCATTGCCCttt                                                                                               <  1:4134385/101‑1 (MQ=38)


       aGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.G.GCAAGGTAAAAGCATTGCCCttt                                                                                               <  2:212721/101‑1 (MQ=38)


        gACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.G.GCAAGGTAAAAGCATTGCCCTTTa                                                                                              <  1:3731959/101‑1 (MQ=38)


        gACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.G.GCAAGGTAAAAGCATTGCCCTTTa                                                                                              >  2:6580168/1‑101 (MQ=37)


        gACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.G.GCAAGGTAAAAGCATTGCCCTTTa                                                                                              >  1:5933999/1‑101 (MQ=37)


           cGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.G.GCAAGGTAAAAGCATTGCCCTTTATCt                                                                                           <  1:6720885/101‑1 (MQ=38)


            gAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.G.GCAAGGTAAAAGCATTGCCCTTTATCtt                                                                                          >  1:3092540/1‑101 (MQ=37)


            gAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.G.GCAAGGTAAAAGCATGGCCCTTTATCtt                                                                                          <  1:5124517/101‑1 (MQ=37)


             aaCCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.G.GCAAGGTAAAAGCATTGCCCTTTATCTTg                                                                                         >  2:126943/1‑101 (MQ=37)


             aaCCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.G.GCAAGGTAAAAGCATTGCCCTTTATCTTg                                                                                         >  2:4934716/1‑101 (MQ=37)


             aaCCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.G.GCAAGGTAAAAGCATTGCCCTTTATCTTg                                                                                         <  1:1236688/101‑1 (MQ=38)


                 gcgcgTGCCGAACTGATGCGGGAGTGGGGGGCTCTGCCC.CCC..GCTTCGCTCCCCCCGGGATATTT.G.GAGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTc                                                                                     <  2:2914273/101‑1 (MQ=37)


                  cgcgTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.G.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCa                                                                                    <  2:5898947/101‑1 (MQ=38)


                  cgcgTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.G.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCa                                                                                    <  2:3773739/101‑1 (MQ=38)


                   gcgTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.G.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAg                                                                                   <  1:4034858/101‑1 (MQ=38)


                     gTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.G.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGat                                                                                 <  2:6188528/101‑1 (MQ=38)


                     gTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.G.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGat                                                                                 <  2:6312089/101‑1 (MQ=38)


                     gTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.G.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGat                                                                                 >  2:5418272/1‑101 (MQ=37)


                      tGCTGAACTGATACGGGAGTGGGGGGCTCTGCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.G.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGata                                                                                <  1:2928869/101‑1 (MQ=38)


                      tGCTGAACTGATACGGGAGTGGGGGGCTCTGCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.G.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGata                                                                                <  1:4960106/101‑1 (MQ=38)


                      tGCTGAACTGATACGGGAGTGGGGGGCTCTGCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.G.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGata                                                                                <  1:4462566/101‑1 (MQ=38)


                       gCTGAACTGATACGGGAGTGGGGGGCTCTGCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.G.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGatat                                                                               <  1:5677309/101‑1 (MQ=37)


                         tGAACTGATACGGGAGTGGGGGGCTCTGCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.G.GCAAGGTAAAAGAATTGCCCTTTTTCTGTCCTCCGATATcc                                                                             >  2:6948141/1‑101 (MQ=25)


                          gAACTGATACGGGAGTGGGGGGCTCTGCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.G.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATccc                                                                            >  2:5651371/1‑101 (MQ=37)


                           aaCTGATACGGGAGTGGGGGGCTCTGCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.G.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCg                                                                           <  2:63870/101‑1 (MQ=37)


                            aCTGATACGGGAGTGGGGGGCTCTGCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.G.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCgg                                                                          <  2:1865653/101‑1 (MQ=37)


                              tGATACGGGAGTGGGGGGCTCTGCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.G.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCgggg                                                                        >  1:3190135/1‑101 (MQ=37)


                               gATACGGGAGTGGGGGGCTCTGCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.G.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCggggg                                                                       >  2:1541269/1‑101 (MQ=37)


                                   cGGGAGTGGGGGGCTCTGCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.G.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGaa                                                                   >  1:5584373/1‑101 (MQ=25)


                                    gggAGTGGGGGGCTCTGCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.G.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAAt                                                                  <  2:2372058/101‑1 (MQ=37)


                                    gggAGTGGGGGGCTCTG.CCGCCG..GCTTCGCGCCCCCCGGGGGATTTGG.G.GCAAGGTAAAAGCATTGCGCTTTAGCTTGCCTCAGAGATCCCgggggcgtgt                                                                  >  2:2254919/1‑96 (MQ=255)


                                      gagtggggggCTCTGCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.G.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAAttt                                                                <  1:6164238/96‑1 (MQ=37)


                                       agtggggggCTCTGCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.G.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAAtttt                                                               <  1:3431142/97‑1 (MQ=37)


                                        gtggggggCTCTGCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.G.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTg                                                              >  1:6441932/4‑101 (MQ=37)


                                         tggggggCTCTGCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.G.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTgg                                                             >  1:947206/3‑101 (MQ=37)


                                           gggggCTCTGCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.G.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCNGGGGGTGAATTTTGgcg                                                           >  2:2726294/1‑101 (MQ=37)


                                            gggtctgtgccc.ccc..GCTTCGCTCCCCCCGTGATATTTGG.G.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGTATTTTGgcgc                                                          <  2:681704/91‑1 (MQ=37)


                                              ggCTCTGCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.G.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAg                                                        <  2:33441/101‑1 (MQ=38)


                                              ggCTCTGCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.G.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAg                                                        <  2:3953384/101‑1 (MQ=38)


                                              ggCTCTG.CC.CCCCGGCTTCGCTCCCCCCGGGATATTT.GCG.GCAAGATGAAGGGGCATGATTTCATCTTGCCCGAAATATCCCAGGGGTGAAttggccaaag                                                         >  2:4173107/1‑93 (MQ=2)


                                                ctctGCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.G.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGcc                                                      <  2:2056284/101‑1 (MQ=37)


                                                 tctgccc.ccc..GCTTCGCTCCCCCCGGGATATTTGG.G.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCa                                                     <  2:5169106/96‑1 (MQ=255)


                                                 tctgccc.ccc..GCTTCGCTCCCCCCGGGATATTTGG.G.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCa                                                     >  1:2577277/6‑101 (MQ=255)


                                                  ctgccc.ccc..GCTTCGCTCCCCCCGGGATATTTGG.G.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCaa                                                    <  1:1417652/97‑1 (MQ=255)


                                                    gccc.ccc..GCTTCGCTCCCCCCGGGATATTTGG.G.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCaaaa                                                  <  2:4526553/99‑1 (MQ=255)


                                                     ccc.ccc..GCTTCGCTCCCCCCGGGATATTTGG.G.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAag                                                 <  2:3998698/100‑1 (MQ=255)


                                                     ccc.ccc..GCTTCGCTCCCCCCGGGATATTTGG.G.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAag                                                 >  1:4197969/2‑101 (MQ=255)


                                                     ccc.ccc..GCTTCGCTCCCCCCGGGATATTTGG.G.GCAAGGTAAAAGAATTGCCCATTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAag                                                 <  1:2800169/100‑1 (MQ=37)


                                                       c.ccc..GCTTCGCTCCCCCCGGGATATTT.GAG.GCGAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAaga.g                                              >  1:6272088/1‑101 (MQ=255)


                                                         ccc..GCTTCGCTCCCCCCGGGATATTTGG.G.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGA.gg                                             <  1:1880583/101‑1 (MQ=255)


                                                              gCTTCGCTCCCCCCGGGATATTTGG.G.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGA.ggggg                                          >  2:3435436/1‑100 (MQ=255)


                                                               cTTCGCTCCCCCCGGGATATTTGG.G.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGA.gggggc                                         >  1:3467647/1‑99 (MQ=255)


                                                               cTTCGCTCCCCCCGGGATATTTGG.G.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGA.gggggc                                         >  2:5373264/1‑99 (MQ=255)


                                                                 tCGCTCCCCCCGGGATATTTGG.G.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGA.gggggcag                                       >  2:1072397/1‑97 (MQ=255)


                                                                  cGCTCCCCCCGGGATATTTGG.G.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGA.gggggcagc                                      >  1:1969760/1‑96 (MQ=255)


                                                                  cGCTCCCCCCGGGATATTTGG.G.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGA.gggggcagc                                      <  1:6948514/101‑6 (MQ=255)


                                                                  cGCTCCCCCCGGGATATTTGG.G.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGA.gggggcagc                                      <  2:2577181/101‑6 (MQ=255)


                                                                   gCTCCCCCCGGGATATTTGG.G.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAgcg                                      <  1:2330827/101‑1 (MQ=25)


                                                                   gCTCCCCCCGGGATATTTGG.G.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAgcg                                      <  1:482084/101‑1 (MQ=25)


                                                                    cTCCCCCCGGGATATTTGG.G.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAgcgc                                     >  2:6045013/1‑101 (MQ=21)


                                                                     tCCCCCCGGGATATTTGG.G.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGcc                                    >  2:2987206/1‑101 (MQ=21)


                                                                      ccccccGGGATATTTGG.G.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGccc                                   >  2:6720381/1‑101 (MQ=21)


                                                                      ccccccGGGATATTTGG.G.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGccc                                   >  1:3132701/1‑101 (MQ=21)


                                                                          ccGGGATATTTGG.G.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTc                               <  1:1675362/101‑1 (MQ=21)


                                                                           cGGGATATTTGG.G.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTcc                              >  2:4790549/1‑101 (MQ=21)


                                                                           cGGGATATTTGG.G.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTcc                              >  1:3560701/1‑101 (MQ=21)


                                                                            gggATATTTGG.G.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTccc                             >  1:1520467/1‑101 (MQ=21)


                                                                             ggATATTTGG.G.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTcccc                            >  1:753979/1‑101 (MQ=21)


                                                                               atatTTGG.G.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCa                          >  1:1054010/1‑101 (MQ=21)


                                                                                 atttgg.g.gCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACg                        >  1:1127044/6‑101 (MQ=255)


                                                                                 atttgg.g.gCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACg                        >  1:2187288/6‑101 (MQ=255)


                                                                                  tttgg.g.gCAAGGTAACAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACgg                       <  1:3278816/97‑1 (MQ=18)


                                                                                  tttgg.g.gCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACgg                       >  1:6514744/5‑101 (MQ=255)


                                                                                  tttgg.g.gCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACgg                       <  1:544611/97‑1 (MQ=255)


                                                                                    tgg.g.gCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGcc                     >  2:2667270/3‑101 (MQ=255)


                                                                                    tgg.g.gCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGcc                     >  1:2708988/3‑101 (MQ=255)


                                                                                     gg.g.gCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCg                    >  1:5748314/2‑101 (MQ=255)


                                                                                        g.gCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCggg                  >  2:1410556/1‑101 (MQ=255)


                                                                                           cAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTc                >  2:880871/1‑101 (MQ=255)


                                                                                            aaGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTcc               >  2:3630366/1‑101 (MQ=255)


                                                                                             aGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCt              >  1:2468526/1‑101 (MQ=255)


                                                                                             aGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCt              >  1:1693383/1‑101 (MQ=255)


                                                                                               gTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGc            >  2:5737793/1‑101 (MQ=255)


                                                                                               gTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGc            <  1:1787832/101‑1 (MQ=255)


                                                                                                 aaaaGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGTGCCCCCTCCCCCACGGCCGGGTCCTGCGt          >  2:2821114/1‑101 (MQ=255)


                                                                                                 aaaaGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGt          >  1:3877919/1‑101 (MQ=255)


                                                                                                   aaGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGTTa        >  1:7117077/1‑101 (MQ=255)


                                                                                                   aaGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGTTa        <  2:3164462/101‑1 (MQ=255)


                                                                                                     gCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGTTAGa      <  1:595196/101‑1 (MQ=255)


                                                                                                       aTTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGTTAGAAt    <  2:6496053/101‑1 (MQ=255)


                                                                                                         tGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGTCCGGGTCCTGCGGTAGAATTg  >  1:1271265/1‑101 (MQ=18)


                                                                                                         tGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGTTAGAATTg  <  1:5377978/101‑1 (MQ=255)


                                                                                                         tGCCCTTTATCTTGCCTCAGATAGCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGTTAGAATTg  <  2:7033611/101‑1 (MQ=255)


CGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGT.GGGGGCTCTG.CC.CCC..GCTTCGCTCCCCCCGGGATATTT.G.G.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGA.GGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGTTAGAATTG  >  NC_014034/3609029‑3609225

Base quality scores:


ATCG

< 3 ≤


ATCG

< 23 ≤


ATCG

< 31 ≤


ATCG

< 35 ≤


ATCG

< 41 ≤


ATCG