BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_014034	2,249,374	+G	2.3%	intergenic (‑170/‑135)	parE ← / → RCAP_rcc02080	DNA topoisomerase IV subunit B/HAD superfamily hydrolase

Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_014034	2,249,374	2	.→G	2.3%	4.8 ‑2.1 ‑0.3	173	intergenic (‑170/‑135)	parE/RCAP_rcc02080	DNA topoisomerase IV subunit B/HAD superfamily hydrolase
Reads supporting (aligned to +/- strand): new base (4/0): ref base (86/83): total (90/83)
Fisher's exact test strand distribution p-value = 1.22e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 4.54e-01

Only 100 of 194 total aligned reads displayed.


GCTCTG.CCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CG.CCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGG  >  NC_014034/2249275‑2249472


gCTCTGCCCCCGTCCCTTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCa                                                                                                         <  2:1105257/101‑1 (MQ=255)


 ctctGCCCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..g                                                                                                      <  1:106933/101‑1 (MQ=255)


  tctgcccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..Gc                                                                                                     >  1:4224026/6‑101 (MQ=255)


  tctgcccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGGGAA.TGGCGCA..Gc                                                                                                     >  1:811589/6‑101 (MQ=255)


   ctgcccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..Gcc                                                                                                    <  1:1853884/97‑1 (MQ=255)


     gcccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCag                                                                                                  >  2:1915529/3‑101 (MQ=255)


     gcccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCag                                                                                                  <  1:2604459/99‑1 (MQ=255)


       ccccGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..CCCagag                                                                                                >  2:1681540/1‑101 (MQ=255)


         ccgtcccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAggg                                                                                              <  2:2269323/101‑1 (MQ=255)


             ccgatctcGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGc                                                                                          <  2:425147/94‑1 (MQ=255)


             cccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGc                                                                                          >  2:3872364/1‑101 (MQ=255)


              ccgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCa                                                                                         <  1:560713/101‑1 (MQ=255)


               cgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.g                                                                                       >  2:4102191/1‑101 (MQ=255)


               cgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.g                                                                                       >  1:84139/1‑101 (MQ=255)


               cgtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.g                                                                                       <  2:24645/101‑1 (MQ=255)


               cgcaagcggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGAT..GCCTTTCATCTTGCCCGAAATATCCC.GGGGGGTGAA.TGGCGAA..GCCAGAggggg                                                                                            >  2:2990221/8‑101 (MQ=18)


                gtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.g.c                                                                                     >  2:1123088/1‑101 (MQ=255)


                gtGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.g.c                                                                                     <  1:1001720/101‑1 (MQ=255)


                    cGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.gatcgg                                                                                  >  1:425204/1‑96 (MQ=255)


                    cGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CG.ccc                                                                                >  1:1918090/1‑101 (MQ=255)


                    cGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CG.ccc                                                                                >  1:2850463/1‑101 (MQ=255)


                     gggACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CG.cccc                                                                               >  1:354411/1‑101 (MQ=255)


                       gACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CG.cccccc                                                                             <  2:3593020/101‑2 (MQ=255)


                        aCT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGTGGGGCA.G.CG.cccccct                                                                            <  2:4129505/101‑3 (MQ=255)


                          t.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CG.ccccccttc                                                                          >  1:1710117/1‑97 (MQ=255)


                           ccccccGGGTTTTTTGGGGGAAGAGGAA...GGGG..GGTTTTCATTTTGCCAAAAAAATCCCGGGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CG.cccccctt                                                                           <  1:1422384/100‑4 (MQ=21)


                            cccccGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CG.ccccccttcg                                                                         <  2:3959378/101‑6 (MQ=255)


                            cccccGGGATATTTGAGGCAAGAGGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CG.ccccccttcg                                                                         <  2:450745/101‑6 (MQ=255)


                             ccccGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTcgc                                                                         <  1:1547977/101‑1 (MQ=255)


                             ccccGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTcgc                                                                         <  2:1633285/101‑1 (MQ=255)


                              cccGGGATATTTGCAGCAAGGTGAAACCGGTT..GCCTTTCATCTGGCCCGAAATATCCC.GGGGGGTGAA.TGGCGAA..GCCAGAGGGGGGCA.G.CG.ccctcctgc                                                                          <  1:291216/101‑7 (MQ=14)


                              cccGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTcgcg                                                                        <  2:3190971/101‑1 (MQ=255)


                               ccGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCgcgg                                                                       >  1:82437/1‑101 (MQ=255)


                                cGGGATATTTGCGGCAAGATGAA...GGGGCATGATTTCATCTTGCCCGAAATATCCC..AGGGGTGAATTGGCCAAAGGCCAAGAGGGGGCA.G.CG.CCCCCtgccc                                                                         >  2:1996119/1‑97 (MQ=2)


                                cGGGATATTTGCGGCAAGATGAA...GGGGCATGATTTCATCTTGCCCGAAATATCCC..AGGGGTGAATTGGCCAAAGGCCAAGAGGGGGCA.G.CG.CCCCCtgccc                                                                         >  2:1663196/1‑97 (MQ=0)


                                cGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCgcggg                                                                      <  2:3162182/101‑1 (MQ=255)


                                cGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCgcggg                                                                      >  1:2366116/1‑101 (MQ=255)


                                 gggATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCG.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCgcgggc                                                                     <  2:2916758/101‑1 (MQ=255)


                                 gggATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCgcgggc                                                                     >  2:3783320/1‑101 (MQ=255)


                                 gggATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCgcgggc                                                                     >  1:1763632/1‑101 (MQ=255)


                                  gaatatTGGAAGAAAGATGAG...GGGG..TTTTTTCCTCGTGCAAAAAATCTCCG.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.GCCC.CCCCCTTCgcgggcg                                                                    <  1:762483/99‑1 (MQ=255)


                                      atTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCcg                                                                >  1:3698587/1‑101 (MQ=255)


                                      atTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTCCGCGGGCGGCcg                                                                >  1:991321/1‑101 (MQ=255)


                                        ttGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGGGGGGGGCA.G.CGCCCCCCTTCGCGGCCGGCcgcg                                                              >  2:1268224/1‑101 (MQ=37)


                                          gAGGCAAGATGAA...GGGG..GTCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGAc                                                            <  2:1972989/101‑1 (MQ=255)


                                            ggCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGAcgc                                                          >  2:370932/1‑101 (MQ=255)


                                               aaGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGAcgcgcg                                                       >  2:1837501/1‑101 (MQ=255)


                                               aaGAGGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGAcgcgcg                                                       <  1:2440432/101‑1 (MQ=255)


                                                aGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACgcgcgc                                                      >  2:1588134/1‑101 (MQ=255)


                                                  aTGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGcgccg                                                    >  2:3618452/1‑101 (MQ=255)


                                                    gAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCa                                                  <  2:3521903/101‑1 (MQ=255)


                                                    gAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCa                                                  <  1:1461113/101‑1 (MQ=255)


                                                      a...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAAc                                                <  1:1899413/101‑1 (MQ=255)


                                                         ggggg..ggTTTTCTTCTTGCGAAAAATTCCGG.GGGGGGTAAA.TGGCGCG..GCCGGGGGGGGGCG.GCCC.CCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAAc                                                <  1:3132242/100‑1 (MQ=255)


                                                          gggg..ggTTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACa                                               <  2:899651/101‑1 (MQ=255)


                                                           ggg..ggCTTTCATCTTGCCGAAAATATCCG.GGGGGGTGAA.TGGCTCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAg                                              <  1:382674/101‑1 (MQ=37)


                                                           ggg..ggCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCCACAg                                              >  2:1053412/1‑101 (MQ=255)


                                                                ggCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTc                                           >  2:2568943/1‑101 (MQ=255)


                                                                 gCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCg                                          <  1:3207496/101‑1 (MQ=255)


                                                                 gCTTTCATCTTGCCGAAAATATCCC.GGGGGGGGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCg                                          <  1:3527858/101‑1 (MQ=255)


                                                                   tttCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTg                                        >  1:1186225/1‑101 (MQ=255)


                                                                     tCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTTGTGcc                                      >  2:2107527/1‑101 (MQ=255)


                                                                     tCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGcc                                      >  1:2272169/1‑101 (MQ=255)


                                                                        tCTTGCCGAAAATATCCG.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCt                                   <  1:2631767/101‑1 (MQ=255)


                                                                         cTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCtt                                  <  1:3529238/101‑1 (MQ=255)


                                                                          ttGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCttt                                 <  1:1265657/101‑1 (MQ=255)


                                                                           tGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTg                                >  1:2210826/1‑101 (MQ=255)


                                                                           tGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTg                                >  1:1651336/1‑101 (MQ=255)


                                                                           tGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTg                                >  2:2885210/1‑101 (MQ=255)


                                                                           tGCCAAAAATATCCG.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTg                                <  2:2924812/101‑1 (MQ=38)


                                                                           ggcagAAAATACCCG.GGGGGGTGAA.GGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTg                                <  2:811511/97‑1 (MQ=37)


                                                                            gCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGt                               >  2:3357811/1‑101 (MQ=255)


                                                                             ccGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTc                              >  1:3153813/1‑101 (MQ=255)


                                                                              cGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCa                             <  2:1009263/101‑1 (MQ=255)


                                                                                aaaaTATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGGGCCCCTTTGTCAtt                           >  1:1053140/1‑101 (MQ=255)


                                                                                 aaaTATCCC.GGGGGGTGAA.TGGCGCA..GCCAGGGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTc                          <  1:3785762/101‑1 (MQ=255)


                                                                                  aaTGTCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCt                         <  1:2259244/101‑1 (MQ=255)


                                                                                         c.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAc                  <  2:2787310/101‑1 (MQ=255)


                                                                                           ggggggTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAcc                 <  1:208656/101‑1 (MQ=255)


                                                                                           ggggggTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCAcc                 >  2:4010255/1‑101 (MQ=255)


                                                                                             ggggTGAA.TGGCGCA..GCCAGAGGGGGCCACC.CC.CCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCa               <  2:514208/101‑1 (MQ=255)


                                                                                              gggTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCAt              <  2:982947/101‑1 (MQ=255)


                                                                                              gggTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCAt              >  1:3355792/1‑101 (MQ=255)


                                                                                              gggTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCAt              <  2:1706323/101‑1 (MQ=255)


                                                                                               ggTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATg             <  2:3299626/101‑1 (MQ=255)


                                                                                                gTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGt            <  1:3028008/101‑1 (MQ=255)


                                                                                                gTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGt            <  2:1777126/101‑1 (MQ=255)


                                                                                                 tGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTc           >  2:1436569/1‑101 (MQ=255)


                                                                                                 tGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTc           >  1:1124005/1‑101 (MQ=255)


                                                                                                  ttccgatctcA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTc           <  2:4117340/92‑1 (MQ=255)


                                                                                                  gAA.GGGCGCA..GCCAGGGGGGGGCA.GCCC.CCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCa          <  1:3683295/101‑1 (MQ=255)


                                                                                                    a.TGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATa        <  1:3246538/101‑1 (MQ=255)


                                                                                                      tGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATaa       <  2:289922/101‑1 (MQ=255)


                                                                                                      tGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATaa       <  1:675004/101‑1 (MQ=255)


                                                                                                      tGGCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATaa       <  2:3818028/101‑1 (MQ=255)


                                                                                                       ggCGCA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAg      <  2:1882305/101‑1 (MQ=255)


                                                                                                         cgcA..GCCAGAGGGGGGGA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGATCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCt    >  2:1289080/1‑101 (MQ=255)


                                                                                                         cgcA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCt    <  1:2475575/101‑1 (MQ=255)


                                                                                                          gcA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTg   <  1:1251423/101‑1 (MQ=255)


                                                                                                           cA..GCCAGAGGGGGGCA.G.CGCCCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCagatcggaa  >  1:4117601/1‑93 (MQ=255)


GCTCTG.CCCCGTCCCGTGCCGGGACT.CCCCCGGGATATTTGAGGCAAGATGAA...GGGG..GGCTTTCATCTTGCCGAAAATATCCC.GGGGGGTGAA.TGGCGCA..GCCAGAGGGGGGCA.G.CG.CCCCCTTCGCGGGCGGCCGCGACGCGCGCCGCAACAGGTCGTGCCCCTTTGTCATTCTGTCGCACCCATGTCATAAGCTGG  >  NC_014034/2249275‑2249472

Base quality scores:


ATCG

< 3 ≤


ATCG

< 15 ≤


ATCG

< 24 ≤


ATCG

< 33 ≤


ATCG

< 40 ≤


ATCG