BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_014034	221,967	G→A	21.7%	D230D (GAC→GAT)	RCAP_rcc00185 ←	hypothetical protein

Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_014034	221,967	0	G→A	21.7%	0.0 ‑0.1 ‑0.7	950	D230D (GAC→GAT)	RCAP_rcc00185	hypothetical protein
Reads supporting (aligned to +/- strand): new base (110/96): ref base (400/342): total (512/438)
Fisher's exact test strand distribution p-value = 9.37e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.82e-01

Only 100 of 2295 total aligned reads displayed.


CCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATC.CGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGAT  >  NC_014034/221867‑222067


    ggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccg                                                                                                   >  1:1942731‑M1/86‑101 (MQ=255)


    ggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggggccgtgggatcggtcgggtccg                                                                                                   <  1:5181468‑M1/16‑1 (MQ=255)


      tcggatccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtg                                                                                                 >  1:3532299‑M1/84‑101 (MQ=255)


            ccGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcg                                                                                           <  2:5626631/101‑1 (MQ=255)


             cgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtggggtccgtcgggtccgtgggatcgg                                                                                          <  1:3751902‑M1/25‑1 (MQ=255)


              gtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggacgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggt                                                                                         <  1:5972708‑M1/26‑1 (MQ=255)


               tgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtc                                                                                        <  2:4171885‑M1/27‑1 (MQ=255)


               ggggatcggtcgggtccgtggggtccgtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtc                                                                                        >  2:6744052‑M1/75‑101 (MQ=255)


                  ggtcggtcgggtccgtcggatcggtgggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtcggatcggtcgggtccgtgggatcggtcggg                                                                                     >  2:4286991‑M1/72‑101 (MQ=255)


                   ttCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggt                                                                                    <  2:5402349/101‑1 (MQ=255)


                   ttCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggt                                                                                    <  1:2653756/101‑1 (MQ=255)


                     cGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtcc                                                                                  <  1:4022613/101‑1 (MQ=255)


                      ggtcgggtccgtcggatcggtgggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtggggtcggtcgggtccg                                                                                 <  2:3260620‑M1/34‑1 (MQ=255)


                      ggtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtggggtcggtcgggtccg                                                                                 <  1:6119635‑M1/34‑1 (MQ=255)


                       gtgtTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtccgt                                                                                >  1:3070928/1‑101 (MQ=255)


                          gggtccgtgggatcggtcgggtccgtcgaatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtggg                                                                             <  1:5064372‑M1/38‑1 (MQ=255)


                          gggtccgtcggatcggtgggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtcggatcggtcgggtccgtgggatcggtcgggtccgtggg                                                                             >  1:33036‑M1/64‑101 (MQ=255)


                             tccgtcggatcggtgggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtcggatcggtcgggtccgtgggatcggtcgggtccgtgggatc                                                                          <  1:3203715‑M1/41‑1 (MQ=255)


                             gtggtCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCggtcgggtccgtgggatcggtcgggtccgtgggatc                                                                          >  1:5306276/1‑101 (MQ=255)


                              ccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcg                                                                         <  1:1497892‑M1/42‑1 (MQ=255)


                              ccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgggggatcggtcgggtcggtgggatcg                                                                         <  2:1307999‑M1/42‑1 (MQ=255)


                                  cAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGggtccgtgggatcggtcgggtccgtgggatcggtcg                                                                     <  2:873622/101‑1 (MQ=255)


                                        ggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTC.Cg                                                              <  2:1157638‑M1/52‑1 (MQ=255)


                                        cGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATC.Cg                                                              >  2:6148266/1‑101 (MQ=255)


                                          ccGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATC.CGTg                                                            >  1:1239740/1‑101 (MQ=255)


                                           cggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTC.CGTgg                                                           <  1:2403550‑M1/55‑1 (MQ=255)


                                           cggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTC.CGTgg                                                           <  1:542653‑M1/55‑1 (MQ=255)


                                            ggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATC.CGTggg                                                          >  1:5547251‑M1/46‑101 (MQ=255)


                                             ggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATC.CGTGGGa                                                         <  2:2779992‑M1/57‑1 (MQ=255)


                                             gATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTGGGATC.CGTGGGa                                                         >  1:7058460/1‑101 (MQ=255)


                                                ccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATC.CGTGGGAtcg                                                      <  1:776304‑M1/60‑1 (MQ=255)


                                                 cgtcggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTC.CGTGGGAtcgg                                                     >  1:134468‑M1/41‑101 (MQ=255)


                                                      cATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGACGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATC.CGTGGGATCGGTCggg                                                >  2:3452329/1‑101 (MQ=255)


                                                       atcggtcgggtccgggggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATC.CGTGGGATCGGTCGGGt                                               >  1:2755424‑M1/35‑101 (MQ=255)


                                                         gggtcgggtgcgtgggaccggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTC.CGTGGGATCGGTCGGGTcc                                             >  1:62297‑M1/33‑101 (MQ=255)


                                                         gcgtcgggtccgtgggatcggtcgggtccgtcGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATC.CGTGGGATCGGTCGGGTcc                                             >  2:5512437‑M1/33‑101 (MQ=255)


                                                         cggtgggattcgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATC.CGTGGGATCGGTCGGGTcc                                             >  1:5907037‑M1/33‑101 (MQ=255)


                                                           gtcgggtccgtgggatcggtcgggtccgtcGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATC.CGTGGGATCGGTCGGGTCCGt                                           >  2:3134365‑M1/31‑101 (MQ=255)


                                                             ggTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATC.CGTGGGATCGGTCGGGTCCGTgg                                         >  1:6654342/1‑101 (MQ=255)


                                                              ggatccgtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTC.CGTGGGATCGGTCGGGTCCGTggg                                        <  1:4164008‑M1/74‑1 (MQ=255)


                                                               gatccgtgggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTC.CGTGGGATCGGTCGGGTCCGTGGGa                                       <  1:16912‑M1/75‑1 (MQ=255)


                                                               gatccgtgggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATC.CGTGGGATCGGTCGGATCCGTGGGa                                       <  1:4445926‑M1/75‑1 (MQ=255)


                                                               gatccgtgggatcggtcggatccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATC.CGTGGGATCGGTCGGGTCCGTGGGa                                       <  2:6840131‑M1/75‑1 (MQ=255)


                                                                 tccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATC.CGTGGGATCGGTCGGATCCGTGGGAtc                                     >  2:4914565‑M1/25‑101 (MQ=255)


                                                                 tccgtgggatcggtcggcgccgtgCGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATC.CGTGGGATCGGTCGGGTCCGTGGGAtc                                     >  2:5871228‑M1/25‑101 (MQ=255)


                                                                 tccgtaggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATC.CGTGGGATCGGTCGGATCCGTGGGAtc                                     >  2:628955‑M1/25‑101 (MQ=255)


                                                                  ccgtgggagcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTC.CGTGGGATCGGTCGGGTCCGTGGGAtcg                                    >  1:4273387‑M1/24‑101 (MQ=255)


                                                                   cgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCCGTCGGATC.CGTGGGATCGGTCGGATCCGGGGGAtcgg                                   <  2:2219753‑M1/79‑1 (MQ=255)


                                                                   cgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGGTCGGTCGGGTC.CGTGGGATCGGTCGGGTCCGTGGgatccg                                   <  1:4806175‑M1/79‑3 (MQ=255)


                                                                   cgtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATC.CGTGGGATCGGTCGGATCCGTGGGAtcgg                                   <  2:5074496‑M1/79‑1 (MQ=255)


                                                                     tCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTC.CGTAGGATCGGTCGGGTCCGTCGGGtccctg                                 >  2:5752656/1‑97 (MQ=11)


                                                                      gggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATC.CGTGGGATCGGTCGGGTCCGTGGGAtcggtcg                                <  2:3532102‑M1/82‑1 (MQ=255)


                                                                      gggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTC.CGTGGGATCGGTCGGGTCCGTGGGAtcggtcg                                <  1:3623097‑M1/82‑1 (MQ=255)


                                                                        gatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATC.CGTGGGATCGGTCGGGTCCGTGGGatcggtcggg                              <  2:379526‑M1/84‑2 (MQ=255)


                                                                         atcggtcgggtccgtcGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATC.CGTGGGATCGGTCGGGTCCGTGGgatcggtcgggt                             >  2:143559‑M1/17‑99 (MQ=255)


                                                                          tcggtctggtccgtgGGATCGGTCGGGTCCGGGGGATCGGTCGGGTCCGTGGGATCGGTCGGATC.CGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATc                            >  2:725480‑M1/16‑101 (MQ=255)


                                                                              tcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATC.CGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTg                        >  1:907099‑M1/12‑101 (MQ=255)


                                                                               gggCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATC.CGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTgg                       >  1:2952963/1‑101 (MQ=255)


                                                                                gggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATC.CGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTggg                      <  1:1087103‑M1/92‑1 (MQ=255)


                                                                                gggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTC.CGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTggg                      >  1:2531023‑M1/10‑101 (MQ=255)


                                                                                 ggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTC.CGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGa                     <  2:7010829‑M1/93‑1 (MQ=255)


                                                                                  gtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATC.CGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAt                    >  1:4350406‑M1/8‑101 (MQ=255)


                                                                                   tcggtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTC.CGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtc                   >  1:100103‑M1/7‑101 (MQ=255)


                                                                                   tccgtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATC.CGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtc                   <  2:2669392/98‑1 (MQ=11)


                                                                                   gtcgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTC.CGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtc                   >  1:5722011‑M1/7‑101 (MQ=255)


                                                                                    gcgttgcataggtcgggtccgtgggatcggtcgggtccgtcggatcGGTCGGGTC.CGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcg                  <  2:2092905/91‑1 (MQ=11)


                                                                                    ccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGATC.CGTGGGATCGGTGGGGTCCGTGGGATCGGTAGAATCCGTgggggca                  >  1:3868550/6‑97 (MQ=11)


                                                                                      gtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTC.CGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcggt                >  2:6413239‑M1/4‑101 (MQ=255)


                                                                                      gtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTC.CGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGCGCGTGGGAtcggc                >  2:478678/1‑100 (MQ=11)


                                                                                       tgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTC.CGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcggtc               <  2:2410830‑M1/99‑1 (MQ=255)


                                                                                        gggatcggtcggatccgtgggatcggtcggatccgtaGGATCGGTCGGGTC.CGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGAtcggtcg              >  1:4663466/2‑101 (MQ=11)


                                                                                        gGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATC.CGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtcg              <  1:6538979‑M1/100‑1 (MQ=255)


                                                                                          gatcggtcgggtccgtgggatgggtcggatccgtagGATCGGTCGGGTC.CGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCggg            <  2:1069972/101‑1 (MQ=11)


                                                                                          gatcggtcgggtccgtgggatcggtcgggtccgtggGGTCCGTCGGATC.CGTGGGATCGGTCGGAGCCGTGGGATCGGTCGGATCCGTGGGATCGGTCggg            >  2:1858092/1‑101 (MQ=14)


                                                                                           atcggtcgggtccgtgggatcggtcgggtccgtgggGTCCGTCGGATC.CGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGGTCGGTCGGGt           >  1:6358576/1‑101 (MQ=14)


                                                                                           atcggtcggatccgtgggatcggtcgggtccgtgggATCGGTCGGATC.CGTGGGATCGGGCGGGGCCGTGGGATCGGTCGGGTCCGTGGGGTCGGCCGGGt           >  1:2846379/1‑101 (MQ=11)


                                                                                            tccggcgggtccgtctgatcggtcgggtccgtcggatcggtCGGGTC.CGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTc          <  2:526099/96‑1 (MQ=11)


                                                                                            gcggtcgtgtgcctgtgctgggtcgggtccgtgggatCGGTCGGGTC.CGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATc          <  1:681579/100‑1 (MQ=11)


                                                                                              ggtcgggtccgtgggatcggtcgggtccgtggggtcCGTCGGATC.CGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGTTCCg        >  1:2215017/1‑101 (MQ=11)


                                                                                              ggtcgggtccgtgggatcggtcgggtccgtgggatcGGTCGGGTC.CGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTGCGTCGGATCCg        >  2:3044973/1‑101 (MQ=11)


                                                                                               gtcgggtgcgtgggatcggtcgggtccgtgggatcgGTCGGGTC.CGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGt       <  2:664020/101‑1 (MQ=11)


                                                                                               gtcgggtccgtgggatcggtcgggtccgtgggatcgGTCGGGTC.CGTGGGATCGGTGGGGTCCGTGGGGTCGGTGGGGTCCGTGGGGTGGGTGGGGTCCgg       >  2:4019289/1‑100 (MQ=11)


                                                                                               gtcggatccgtgggatcggtcgggtccgtgggatcgGTCGGGTC.CGTGGGATCGGTCGGGTCCGGGGGATCGGTCGGGTCCGTGGGTTCGGGCGGGTcggg       >  2:3606146/1‑98 (MQ=11)


                                                                                               gtcggatccgtgggatcggtcgggtccgtgggatcgGTCGGATC.CGTGGGGTCGGTCGGATCCGTGGGCTCCGTGGGATCGGTTGGGtcgggcggcacggg       >  2:4505129/1‑95 (MQ=11)


                                                                                               ggcgggtccgtgtgatcggtcgggtccgtggggtccgtCGGATC.CGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGt       <  2:1956977/99‑1 (MQ=14)


                                                                                                tcgggtccgtgggatcggtcgggtccgtgggatcggTCGGGTC.CGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTg      >  2:5960561‑M2/1‑100 (MQ=255)


                                                                                                tcgggtccgtgggatcggtcgggtccgtgggatcggTCGGGTC.CGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTg      >  1:3440809‑M2/1‑100 (MQ=255)


                                                                                                tcgggtccgtgggatcggtcgggtccgtcggatcggTCGGGTC.CGTGGGATCGGTCGGGTCCGTGGGATCGGTGGGGTCCGGGGGAACGGTTGGGTCCGtg      >  1:2092985/1‑100 (MQ=11)


                                                                                                tcgggtccgtgggatcggtcgggtccgtcggatcggTCGGGTC.CGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTg      >  1:4915837‑M2/1‑100 (MQ=255)


                                                                                                tcgggtccgtgggatcggtcgggtccgtcggatcggTCGGGTC.CGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGGTCCGTGGGATCGGTCGGGTCCggg      >  2:3905618/1‑99 (MQ=11)


                                                                                                tcggatccgtgggatccgtcgggtccgtgggctcggTCGGGTC.CGTGGGATCGTTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGtg      <  2:287702/101‑2 (MQ=11)


                                                                                                gcgggtcgattggatcggtcggagccgttggatcggtCGGATC.CGTGGGATCGGTCGGATCCGTGAGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTa      <  2:4639023/100‑1 (MQ=11)


                                                                                                gcgggtccgtgggatcggtcgggtccggggggtcggtCGGGTC.CGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGtg      <  1:455291/100‑2 (MQ=11)


                                                                                                 cgggtccgtgggatcggtcgggtccgtgggatcggtCGGATC.CGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTgg     <  1:6913612‑M2/101‑3 (MQ=255)


                                                                                                  gggtccgtgggatcggtcgggtccgtggggtccgtcGGATC.CGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGGCCggggg    >  1:1812361/1‑97 (MQ=11)


                                                                                                  gggtccgtgggatcggtcgggtccgtcggatcggtcGGGTC.CGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGGCCGTGgg    >  1:3938337/1‑101 (MQ=14)


                                                                                                   ggtccgtgggatcggtcgggtccgtgggatcggtcgGGTC.CGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTggga   <  2:177006‑M2/101‑5 (MQ=255)


                                                                                                    gtccgtggggtctgtcgggtccgtgggatcggtcggGTC.CGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggat  >  2:2550068‑M2/1‑96 (MQ=255)


                                                                                                    gtccgtgggagcggtcgggtccgtggggtcggtcggATC.CGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggat  >  2:1590424‑M2/1‑96 (MQ=255)


                                                                                                    gtccgtgggaacggtcgggcgcgtggcgtcgg.cggcGCTCGTGGGATCGGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAt  <  2:4377718/101‑1 (MQ=2)


CCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATC.CGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGAT  >  NC_014034/221867‑222067

Base quality scores:


ATCG

< 3 ≤


ATCG

< 17 ≤


ATCG

< 24 ≤


ATCG

< 32 ≤


ATCG

< 40 ≤


ATCG