BRESEQ :: Evidence

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

Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_014034	222,045	0	G→C	7.8%	225.7 ‑0.8 ‑0.0	1519	P204P (CCC→CCG)	RCAP_rcc00185	hypothetical protein
Reads supporting (aligned to +/- strand): new base (44/69): ref base (576/766): total (650/869)
Fisher's exact test strand distribution p-value = 4.30e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 1947 total aligned reads displayed.


TGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCG  >  NC_014034/221945‑222145


  gggtccgtcggATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATGGGCCGGGTCCGGGGGATGAGGCGGGACCGCGGGGGCCGCCGGGTCCTTggg                                                                                                    >  2:4300111/7‑101 (MQ=2)


  ggatccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTggg                                                                                                    <  1:3852036‑M1/92‑1 (MQ=255)


   ggtcggtgggGTCGGTCGCGTGCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGa                                                                                                   <  2:394630/99‑1 (MQ=11)


   ggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGa                                                                                                   >  1:518570‑M1/9‑101 (MQ=255)


   gCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTcggg                                                                                                   <  2:821425/101‑2 (MQ=11)


    ggccgtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAt                                                                                                  <  2:3019469/97‑1 (MQ=11)


     gctcggtgattggtggggtccgtgggagcggtcgggtccgtggCGTCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtc                                                                                                 <  1:289508/94‑1 (MQ=11)


      ccggcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtcg                                                                                                <  2:2040325/97‑1 (MQ=11)


         tgggatcggtcgggtcggtgggatcggtcgggtccctgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGtccgtg                                                                                             <  2:2667609/99‑2 (MQ=11)


          gggatcggtcgggggggtgggatcggtcgggtccgtgGGATTGGTTGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcggtcg                                                                                            <  1:1341021/100‑1 (MQ=14)


          gGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcggtcg                                                                                            >  1:2366659‑M1/2‑101 (MQ=255)


           ggatgggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtcgg                                                                                           <  2:1982681/101‑1 (MQ=11)


           ggatcggtcggatccgtgggatcggtcggatccgtaGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGAtcggtcgg                                                                                           <  2:1585429/101‑1 (MQ=11)


            gatcggtcgggtccgtgggatcggtcgggtccgtggGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGgccgtcgga                                                                                          >  1:1576024/1‑100 (MQ=11)


             gtccgtcgggtccgtgggatcggtcgggtccgtcggatcgGGCGGGGCCGTGGGATCGGCCGGGCCCGTCGGATCGGTGGGGTCCGTGGGATCGGTCgggc                                                                                         >  1:2691363/5‑100 (MQ=11)


             atcggtcggatccgtgggatcggtcgggtccgtgggATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGGGGGATCGGTCGGAACCGTGGGAtcggtcggat                                                                                         >  2:1470578/1‑99 (MQ=11)


              tcggtcggatccgtgggatcggtcggatccgtgggaTCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATc                                                                                        <  1:3495103/101‑1 (MQ=11)


               cggtcgggtccgtggggtccgtcggatccgtgggatCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGGcc                                                                                       >  2:1657380/1‑101 (MQ=11)


               cgggcggggccgtgggctgggtcggggccgtgggctCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTcc                                                                                       <  1:3234614/101‑1 (MQ=11)


                ggttgggggcgtgggatcggtcggatccgtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCg                                                                                      <  2:2118793/92‑1 (MQ=11)


                ggtcgggtccgtgggatcggtcgggtccgtcggatcGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCg                                                                                      >  1:731980/1‑101 (MQ=11)


                  tcgggtccgtgggatcggtcgggtccgtgggatcggTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTg                                                                                    >  1:4550111‑M2/1‑100 (MQ=255)


                  tcggatccgtgggatcggtgggatcggtaggatcggTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGtc                                                                                    <  2:2720940/101‑2 (MQ=11)


                   cgggtccgtgggctcggtctggcgcgtcggctcggtCGGGTCCGTGTGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGtgg                                                                                   <  2:4148905/101‑3 (MQ=11)


                   cgggtccgtgggatcggtcgggtccgtgggatcggtCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgg                                                                                   >  2:493651‑M2/1‑99 (MQ=255)


                   cgggtccgtgggatcggtcgggtccgtgggatcggtCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGGCGGGTCCGTGGGGTTGGTCGGGTGCGtgg                                                                                   >  2:2493096/1‑99 (MQ=11)


                   cggatccgtggggtccgtcgggtccgtgggatcggtCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGGGGGGTcccggg                                                                                   >  2:2418678/1‑97 (MQ=11)


                    ggatccgtgggatcggtcgggtccgtgggatcggtcggATCCGTGGGGTCGGGCGGGTCCGTGGGATCGGGCGGGTCCGTGGGACCGGTCGGGTCCGTGgg                                                                                  >  2:1882450/4‑101 (MQ=14)


                     ggtccgggggatcggtcgggtccgtgggatcggtcgGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGa                                                                                 <  1:3426000/101‑1 (MQ=11)


                     gatccgtgggatcggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGGCGGGGCCGTGGGATCGGGCGGATCCGTGGGATGGGGCGGGTCCGTGGGa                                                                                 >  2:617608/3‑101 (MQ=14)


                      gtccgttgggtccttcggatccgtgggaccggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAt                                                                                <  2:3189754/101‑1 (MQ=11)


                      gtccgtgggatcggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggat                                                                                >  1:438064‑M2/1‑96 (MQ=255)


                      atccgtgggatcggtcgggtccgtgggatcggtcgggTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTgggat                                                                                >  1:824563‑M2/2‑96 (MQ=255)


                       tccgtgggatcggtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGAGGGATCGGTCGGGTCAGTGGGAtc                                                                               <  2:120643/101‑1 (MQ=14)


                         cgtgggatcggtcgggtccgtggggtcggtcggATCCGTGGGATCGGTCGGATCCGTGGGGTCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGtcgg                                                                             >  2:3787436/1‑101 (MQ=11)


                          gtgggatcggtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatcggt                                                                            <  2:726054‑M2/101‑10 (MQ=255)


                          gtgggatcggtcgggtccgtcggatcggtgggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGGCGGGTCCGGGGGGtcggg                                                                            >  2:3627051/1‑100 (MQ=11)


                           gggggttgggcggggccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtc                                                                           <  2:2939002/91‑1 (MQ=11)


                            gggatcggtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcg                                                                          <  2:3337218‑M2/101‑12 (MQ=255)


                             ggatcggtcggttccgtgggctctgtcggATCCGTGGGATCTGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcgg                                                                         <  2:1367194/101‑1 (MQ=11)


                             ggatcggtcgggtccgtcggatcggtgggATCCGTGGGATCGGTCGGGTCCGCGGGATCGGTCGGGTCCGTGAGATCGGCCGGGTCCGTGGGAgcggccgg                                                                         >  2:3686392/1‑101 (MQ=14)


                              ggtcgggcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcgga                                                                        <  1:1257604/99‑2 (MQ=11)


                                   gtggtgtccgtggggtccgtcggATCGGTGGGATCGGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGgtcgggtccgt                                                                   <  2:4004638/96‑1 (MQ=14)


                                      gggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcggatcggtggg                                                                <  1:2034988‑M2/101‑22 (MQ=255)


                                      gggtccgtcggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGGTAGGTCGGGCCCGTGGGATCGgtcgggtccgtcgg                                                                >  1:3971322/1‑101 (MQ=11)


                                      ggatccgtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGgtcgggtccgtcgg                                                                >  1:1830145/4‑101 (MQ=255)


                                      ggatccgtgggatcggtcggATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGgtcgggtccgtcgg                                                                >  1:3034834/4‑101 (MQ=11)


                                       gatccgtaggatgggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCgtggggtccgtcggg                                                               <  1:1003866/99‑1 (MQ=11)


                                        gtccgtgggatcggtcggGTCCGTGGGATCGGTCGGCTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggat                                                              >  2:1272040‑M2/1‑78 (MQ=255)


                                        atccgtgggatcggtcgggTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggat                                                              >  2:3105685‑M2/2‑78 (MQ=255)


                                          cggtgggatcggtcgggtCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcg                                                            >  1:2121087‑M2/3‑76 (MQ=255)


                                           cgtgggatcggtcggGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcag                                                           <  1:3661487‑M2/101‑27 (MQ=255)


                                               gggtccgtcggATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGGCGGGTCCGTCGGATCGGTggg                                                       >  1:2826741/7‑101 (MQ=11)


                                               ggggccgtcggATCGGGGGGGTCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTggg                                                       <  1:4449562/95‑1 (MQ=11)


                                                ggtccgtcggatCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcggatccgtgggatcggtcggg                                                      <  1:1084193‑M2/99‑32 (MQ=255)


                                                 atcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcgggt                                                     >  1:2451517‑M2/1‑69 (MQ=255)


                                                 atcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTgggatcggtcggatccgtgggatcggtcggac                                                     <  2:1751634‑M2/101‑33 (MQ=255)


                                                  gccgtgggATCTGTCCGATCCGTTGGGGCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGAtc                                                    <  2:553149/95‑1 (MQ=14)


                                                    ggtcggGTCCGTGGGATCGGTCTGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGAtcgg                                                  >  1:2961317/1‑101 (MQ=11)


                                                    ggtcggATCCGTGTGATCGGTCGGGTCCTTGGGATCGGGCGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcgggtccg                                                  >  1:983917‑M2/1‑66 (MQ=255)


                                                      tcggGTCCGGGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGAtcggtc                                                <  1:4290638/101‑1 (MQ=11)


                                                      tcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGCTCGGTCGGGTCCGTCGGATCGGTCGGAtccgtg                                                >  1:4051803/1‑100 (MQ=11)


                                                       cggGTCCGTAGGATCGGTCGGGTCCGTCTGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGAtcggtcg                                               <  1:375875/101‑1 (MQ=14)


                                                          ggccgtGGGGTGGGGCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcgggtccgtgggat                                            >  1:633025‑M2/3‑60 (MQ=255)


                                                           tCCGGGGGATCGGTCGGGTCCGTGGGATCGGACGGATCCGTGGGATCGTTCGGGTCCGCCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTc                                           <  1:1185119/101‑1 (MQ=14)


                                                           gccGTGGGGTCCGGCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATc                                           <  2:4489470/100‑1 (MQ=11)


                                                            gcgTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcggatccgtgggatcggtcgggtccgtgggatcg                                          >  1:4640564‑M2/2‑58 (MQ=255)


                                                             cGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTggggtcggtcgggtccgtgggatcggtcgggtccgtgggatcgg                                         <  2:4076970‑M2/101‑45 (MQ=255)


                                                              gtcggATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCgg                                        >  1:2207289/4‑100 (MQ=17)


                                                              gTGGGGTCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGt                                        <  2:2592702/101‑1 (MQ=14)


                                                               tGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtcggatcggtcgggtccgtgggatcggtc                                       >  1:2539936‑M2/1‑55 (MQ=255)


                                                               nggngtcggtcggGGCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTg                                       <  1:4103683/96‑1 (MQ=14)


                                                                 ggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgg                                     <  1:3323288‑M2/101‑49 (MQ=255)


                                                                  gATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcgga                                    >  1:1654178‑M2/1‑52 (MQ=255)


                                                                  gATCGGTCGGGTCCGTGGGATCCGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGa                                    <  2:3886071/101‑1 (MQ=17)


                                                                   gtcgtgcggGTCCGTGGGGTTGGTCGGAGCTGTGGGATCGGTCGGGTCGGTAGGCTCGGTCGGGTCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAt                                   <  2:2265894/95‑1 (MQ=18)


                                                                   aTCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAt                                   >  2:3509535/1‑101 (MQ=255)


                                                                    tcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtc                                  <  1:3137713‑M2/101‑52 (MQ=255)


                                                                    tcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatc                                  <  1:2756016‑M2/101‑52 (MQ=255)


                                                                     cggtcggGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATTGGTGGGATCGGTGGGTTCGGTGGGACCGCTCGGGGCCGGGGGAccg                                 >  2:4208307/1‑101 (MQ=11)


                                                                     cggtcggGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcgggtccgggggatcggtcgggtcc                                 <  1:1734414‑M2/101‑53 (MQ=255)


                                                                        ttggttCCGTGGGGTCGCTCGGTGAGGTGGGATCGGGCGGGTCCGTAGGATCGGTCCGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtc                              <  2:3437221/96‑1 (MQ=18)


                                                                        tcggGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcggatccgtg                              >  1:4175601‑M2/1‑46 (MQ=255)


                                                                         ggggtCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgg                             >  1:904153‑M2/2‑45 (MQ=255)


                                                                                gcgggACCGGTCGGACCCGGGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGt                      <  2:915994/99‑1 (MQ=17)


                                                                                   ggATCGGTGGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgg                   >  2:4396997‑M2/1‑35 (MQ=255)


                                                                                     aTCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTgggat                 <  1:2634231/101‑3 (MQ=14)


                                                                                       cggtcggATCCGTGGGCGCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTaggatcg               <  2:2805045/101‑5 (MQ=14)


                                                                                       cggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCGGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTaggatcg               >  1:3923348/1‑97 (MQ=14)


                                                                                        cgtgggATCCGTGGGATCGGTCGGGTCCGTGGGGTCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTaggatcgg              <  1:1463148/97‑6 (MQ=14)


                                                                                           ggggtCCGGGGGATGGGTCGGGTCGGTTGGATCGGTGGGATCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCtggtgg           <  2:2196585/97‑1 (MQ=17)


                                                                                             ggtCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGa         >  1:3580600/3‑101 (MQ=255)


                                                                                             ggtCCGGGGGAGCGGTCGGGTCAGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGa         <  2:2870216/99‑1 (MQ=255)


                                                                                              gtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtggggtccgtcggatccgtgggatcggtcggatccgtgggaa        <  2:1102170‑M2/100‑78 (MQ=255)


                                                                                              gtCCCTGGGGTCCGTCGGGTCGGTAGGATTGGTCGGGTCCGTGGGGTCGGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGCCtgctggcga        <  2:3259017/100‑4 (MQ=17)


                                                                                               tCCGTGGGATCGGTGGGGTCCGTCGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGAtc       <  1:583313/101‑1 (MQ=255)


                                                                                                 cgggggAGCGGGGGGGTCCGTAGGATCGGGCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGAtcgg     <  1:1979261/98‑1 (MQ=21)


                                                                                                  gtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggt    <  2:4590773‑M2/101‑82 (MQ=255)


                                                                                                   tgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtggggtccgtc   >  2:2346640‑M2/1‑19 (MQ=255)


                                                                                                   ggggatcggtcgggtccgtcggagcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtc   >  1:2289059‑M2/2‑19 (MQ=255)


TGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCG  >  NC_014034/221945‑222145

Base quality scores:


ATCG

< 3 ≤


ATCG

< 14 ≤


ATCG

< 22 ≤


ATCG

< 31 ≤


ATCG

< 40 ≤


ATCG