BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_014034	222,018	C→G	3.0%	P213P (CCG→CCC)	RCAP_rcc00185 ←	hypothetical protein

Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_014034	222,018	0	C→G	3.0%	85.1 ‑4.1 ‑0.9	1794	P213P (CCG→CCC)	RCAP_rcc00185	hypothetical protein
Reads supporting (aligned to +/- strand): new base (38/21): ref base (832/898): total (874/920)
Fisher's exact test strand distribution p-value = 1.66e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.36e-01

Only 100 of 2293 total aligned reads displayed.


GCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCG  >  NC_014034/221918‑222118


 gggatcggtcggatccgtgggatcggtcggatccgtaGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGAtcggtcg                                                                                                     >  2:4060058‑M1/38‑101 (MQ=255)


  ggatcggtcggatccgtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTTCGTGGGGTCCGTCGGATCCGTGGGAtcggtcgg                                                                                                    <  1:4985288‑M1/65‑1 (MQ=255)


        gtcggatccgtgggatctgtccgatccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGt                                                                                              >  2:3608345‑M1/31‑101 (MQ=255)


         tcggatccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTg                                                                                             >  2:2867353‑M1/30‑101 (MQ=255)


          cggatccgtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGAGCGGTCGGGTCCGTGGGATCGGTGGGGTCCGTGGGATCGGTCGGGTCCGTgg                                                                                            <  1:599320‑M1/73‑1 (MQ=255)


              tccgtggggtcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtc                                                                                        >  1:465648‑M1/25‑101 (MQ=255)


                    ggagcggtggggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcgg                                                                                  >  1:1654170‑M1/19‑101 (MQ=255)


                     gatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGa                                                                                 <  2:4385580‑M1/84‑1 (MQ=255)


                        cggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATcc                                                                              <  2:1043393‑M1/87‑1 (MQ=255)


                          gtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGt                                                                            <  2:3370080‑M1/89‑1 (MQ=255)


                           tGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTg                                                                           >  1:1100130/1‑101 (MQ=255)


                           tGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTg                                                                           <  1:5498622/101‑1 (MQ=255)


                            gggATCGGTCTGGGCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGtcg                                                                          <  1:5616936/101‑3 (MQ=11)


                               gtcggttggGTCGGGCGGGTCCGGCGGGTCCGTCGGGGCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAt                                                                       <  2:1678628/100‑1 (MQ=11)


                               gtccgtcggATCCGTTGGGTCCGTTGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAt                                                                       <  2:2283806/97‑1 (MQ=11)


                                 cggtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcg                                                                     >  1:2695729/1‑101 (MQ=255)


                                 ccgggggatcgctcgggtccgtgggctcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcg                                                                     <  2:3959704/96‑1 (MQ=14)


                                  ggtcggGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGCtcgg                                                                    >  1:924846/1‑101 (MQ=1)


                                  cgtgGGGTCGGTCGGGTCCGTGGGGTCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcgg                                                                    >  1:2606726‑M1/5‑101 (MQ=255)


                                   gtcggTTCTGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggt                                                                   <  2:4511478/101‑1 (MQ=11)


                                   gtcggATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAccggg                                                                   >  2:2181156/1‑100 (MQ=11)


                                   gtcggATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggt                                                                   >  1:2231073/1‑101 (MQ=11)


                                   gtcggATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGGCCGTGGGATCGGTCGGGTCCGTGGGAtcggt                                                                   <  1:5196645/101‑1 (MQ=11)


                                    tgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtc                                                                  <  1:1487915‑M1/99‑1 (MQ=255)


                                    gggggtccgtcgggtccgtgggatcggtcgggtccgtcggatcgGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc                                                                  <  2:4724732/93‑1 (MQ=11)


                                    ggGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc                                                                  >  1:4108048‑M1/3‑101 (MQ=255)


                                     gggatcggtcgggtccgtggagtcggtcgggtccgtgGGATCGGTCGGTTCCGTGGGATCGGTCGGTTCGGTGGGATCGGTCGGGTCCGTGGGAtcggtcg                                                                 <  1:3254806/100‑1 (MQ=11)


                                     gGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcg                                                                 >  2:2997595‑M1/2‑101 (MQ=255)


                                     cggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtgg                                                                 >  2:1770851/1‑99 (MQ=255)


                                       gatcggttggatccgtggggtcggtcgggtccgtggGGTCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtcgga                                                               <  2:1579820/101‑2 (MQ=11)


                                        gtcggtcggctccgcgggatcggtcgtgtccgtgcgaTCGGCGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGt                                                              <  2:5242925/100‑1 (MQ=14)


                                         tcggtcgggtccgtgggatcggtcgggtccgtggggTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGGCGGGTCCGTGGGATCGGTCGGGTc                                                             >  2:3262942/1‑101 (MQ=11)


                                         tcggtcgggtccgtgggatcggtcgggtccgtgggaTCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcggtagatcg                                                             >  2:2655326/1‑95 (MQ=14)


                                          cggtgggatgcgtgggatcggtcgggtccgtgggatCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTcc                                                            <  1:3604223/101‑1 (MQ=11)


                                          cggtcgggtccgtgggatcggtcgggtccgtgggatCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGatcc                                                            <  1:2334677‑M2/101‑5 (MQ=255)


                                          cggtcggatccgtgggatcggtcgggtccgtgggatCGGCCGGGTCCGTGGGGTCGGCCGGGTCCGCGGGGGCGGTCGGGTCGGTGGGGTAGGTCggggca                                                            >  2:853018/1‑98 (MQ=2)


                                          ccggcggatccgtagggtcggtgggctccgtggggtctgtCGGCTCCGTCGGCTCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTcc                                                            <  2:4845814/97‑1 (MQ=11)


                                           ggtcgggtccgtcggatcggtcgggtccgtgggatcGGTCGGGTCCGTCGGATCGGTGGGGTCCGTGGGATCGGTCGGGTCCGGGGGATCGGTCGGGTCCg                                                           >  1:402503/1‑101 (MQ=11)


                                           ggtcggatccgtgggatcggtcgggtccgtgggatcGGTCGGATCCGTGGGGTCGGTCGGGTCCGTGGGATCCGTCGGGGCCGTGGGATCTGTCGGGTACg                                                           >  2:594113/1‑101 (MQ=14)


                                           gggcgggtccgtgggatcgggcgggtccgtggggtctgtCGGCTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTAGGGTCCGTGGGATCGGTCGGGTCCg                                                           <  2:3159829/98‑1 (MQ=11)


                                             tcgggtccgtggggtcggtcgggtccgtcggatcggTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGtg                                                         <  1:3739548/101‑2 (MQ=14)


                                             tcgggtccgggggatcggtcgggtcggtgggatcggTCGGGTCCGTGGGGTCCGGCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGtg                                                         <  1:1210254/101‑2 (MQ=11)


                                             tcgcgtccgtgggatcggtcgggtccgggggatcggTCGGGTCCGTTGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGtg                                                         <  2:5134755/101‑2 (MQ=11)


                                              gggatccgtgggatcggtcgggtccgtaggatccgtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGtcg                                                        <  1:2926367/97‑3 (MQ=14)


                                              ctggtccttgggatcggtcgggtccgttggatcggtcgGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGtgg                                                        <  2:3092143/99‑3 (MQ=11)


                                              cgggtccgtgggatcggtcgggtccgtgggatcggtCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGtgg                                                        >  1:3069372/1‑99 (MQ=11)


                                              cgggtccgtgggatcggtcgggtccgtgggatcggtCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGatccgtgg                                                        >  2:1487822‑M2/1‑93 (MQ=255)


                                              cgggtccgtgggatcggtcggatccgtgggatcggtCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGatccgtag                                                        <  2:2068086‑M2/101‑9 (MQ=255)


                                               gggtccgtgggatcggtcgggtccgtggggtccgtcGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGgg                                                       <  2:4618579/101‑1 (MQ=14)


                                               gggtccgtgggatcggtcgggtccgtcggatcggtcGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGgg                                                       >  1:5328801/1‑101 (MQ=14)


                                               gggtccgtgggatcggtcggatccgtgggatcggtcGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGatccgtagg                                                       <  1:2812687‑M2/101‑10 (MQ=255)


                                                ggtccgtgggatcggtcgggtccgtgggatcggtcgGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGatccgagatc                                                      <  1:1719951‑M2/101‑11 (MQ=255)


                                                ggtccgtgggatcgctcgtgtccgtggggtccgtcgGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGa                                                      <  2:2397799/101‑1 (MQ=11)


                                                gatccgtgggatcggtcgggtccgtgggatcggtcggGTCCGTCGGATCGGTCGGGCCCGTGGGATCGGTCGGGCCCGTGGGACCGCTCGGGCCCGTGGGa                                                      >  1:3850455/3‑101 (MQ=11)


                                                 gtccgtgggatcggtcgggtccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGatccgtgggat                                                     <  1:1368218‑M2/101‑12 (MQ=255)


                                                 atccgtggggtccgtcgggtccgtgggatcggtcggGGCCGTCGGATCGGTCGGGTCCGTGGGGTCGGTCGGGGCCGTGGGGTGGGTAGGAGCCGTGggct                                                     >  2:566329/2‑99 (MQ=2)


                                                 atccgtgggatcggtcggatccgtgggatcggtcggaTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGgtccgtgggac                                                     <  1:3480234‑M2/100‑12 (MQ=255)


                                                  tccgtgggatcggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGGCCGTGGGATCGGTCGGGTCCggggggtc                                                    >  2:1875006/1‑94 (MQ=11)


                                                  tccgtgggatcggtcgggtccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGatccgtgggatc                                                    >  1:3153473‑M2/1‑89 (MQ=255)


                                                   ccgtgggatcggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGatccgtgggatcg                                                   >  1:5094476‑M2/1‑88 (MQ=255)


                                                     gtgggatcggtcgggtccgtcggatcggtgggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGAtcggt                                                 <  2:4244898/101‑1 (MQ=11)


                                                      tgggatcggtcggatccgtgggatcggtcggATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGGGGGAtcgggc                                                >  1:3952069/1‑99 (MQ=11)


                                                        ggatcggtcgggtccgtggggtgcgtaggAGCCGGGGGATCGGTCGGATCCGTGGGGTCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcgg                                              <  1:715557/101‑1 (MQ=11)


                                                        ggatcggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCGGTGGGGtcggtagg                                              >  2:3854849/1‑101 (MQ=11)


                                                          gtcggtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggt                                            <  2:4136321/100‑1 (MQ=11)


                                                            cggtcggttccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtcc                                          <  1:2132663/101‑1 (MQ=255)


                                                            cggtcgggtccgtgggatcggtcggATCCGTG.........GGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGgtcgggtccgtcgggtcc                                 >  1:1565522/1‑101 (MQ=11)


                                                             ggtcggatccgtaggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCgtggggtccg                                         >  1:5708758/1‑101 (MQ=11)


                                                               tggggtcggtggggtccgtgggATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGTTCGgtcgggtccgtc                                       <  1:4790845/99‑1 (MQ=11)


                                                                cgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTAGGGTCCGTGGGATCGgtcggatccgtgg                                      >  2:431105/1‑99 (MQ=11)


                                                                  ggtccgtggggtccgtcggAGCCGTGGGATCGGTCGGATCCGTGGGAACGGTCGCATCCGTCGGAGCGGAGGGGTCCGTAGGCGCggt.gggtcagttcgct                                   >  2:805771/1‑96 (MQ=9)


                                                                     ccgtgggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCgtggggtccgtcgggtcc                                 <  1:100380/101‑1 (MQ=14)


                                                                      cgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGGCGGGTCCGTGGGATCGGtcggctacgtgggggccg                                >  2:536046/1‑101 (MQ=11)


                                                                       gtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGgtccgtgggatcggtcgggtccgtgggatcggt                               >  1:82260‑M2/1‑68 (MQ=255)


                                                                       gtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTcggggccgcggggtccgt                               >  2:3482133/1‑101 (MQ=11)


                                                                       gtaggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTggggtccgtcgggtccgt                               >  1:4224916/4‑101 (MQ=14)


                                                                          ggatcggtcggATCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGGCGGATCGGTCGGCTAGGGGGGCTCGGTCGGGTGCGTGGGGTGCGTcgg                            >  2:5555418/1‑101 (MQ=2)


                                                                          ggatcggtcggATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTggg                            <  1:3366926/101‑1 (MQ=11)


                                                                            atcggtcggATCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGatccgtgggatcggtcgggtccgtgggatcggtcggat                          >  2:2731475‑M2/1‑63 (MQ=255)


                                                                             tcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGatccgtgggatcggtcgggtccgtgggatcggtcgggtc                         <  2:1980688‑M2/101‑40 (MQ=255)


                                                                              ccgtgggATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGAtcg                        <  2:254742/96‑1 (MQ=14)


                                                                               ggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGgtccgtgggatcggtcgggtccgtgggatcggtcgggtccg                       <  1:4355686‑M2/101‑42 (MQ=255)


                                                                                gtcggGTCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGatccgtgggatcggtcgggtccgtgggatcggtcgggtccgt                      >  2:3739442‑M2/1‑59 (MQ=255)


                                                                                gtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGAtcggt                      <  2:5453311/101‑1 (MQ=14)


                                                                                 tcggGTCCGTGGGATCGGGCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGAtcggtc                     <  1:4566544/101‑1 (MQ=11)


                                                                                   gggtCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGAtcggtcgg                   >  1:5404887/4‑101 (MQ=14)


                                                                                    ggtccGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGgtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtggga                  >  1:4482020‑M2/3‑55 (MQ=255)


                                                                                     gtcCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGgtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggat                 <  1:195957‑M2/100‑48 (MQ=255)


                                                                                     aTCCGTGGGATCGGTCGTATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAt                 <  2:1273984/101‑1 (MQ=11)


                                                                                      tCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatc                >  2:70385‑M2/1‑53 (MQ=255)


                                                                                      gcggtcggATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATc                <  1:4302373/95‑1 (MQ=17)


                                                                                       ccGTGGGGTCGGTCGGGTCCGTGGGATCGGTCGGCTCCGTGGGATCGGTCGGatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtggggtcc               >  2:1155708‑M2/1‑52 (MQ=255)


                                                                                        cGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCg              >  2:493194/1‑101 (MQ=14)


                                                                                          tGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGatccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtc            >  1:247680‑M2/1‑49 (MQ=255)


                                                                                           gggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcg           <  1:4894484‑M2/101‑54 (MQ=255)


                                                                                            ggATGGGTGGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTcgg          <  1:2029625/101‑1 (MQ=14)


                                                                                            ggATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTggg          <  2:5038015/101‑1 (MQ=14)


                                                                                              ctcggtcggGTCCGTTGGAGCGGTGGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGGGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAt        <  1:130119/100‑1 (MQ=14)


                                                                                               tcggtcggGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtc       <  2:1341548‑M2/101‑58 (MQ=255)


                                                                                                    cggGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgg  <  1:1853199‑M2/101‑63 (MQ=255)


GCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCG  >  NC_014034/221918‑222118

Base quality scores:


ATCG

< 3 ≤


ATCG

< 15 ≤


ATCG

< 23 ≤


ATCG

< 31 ≤


ATCG

< 40 ≤


ATCG