BRESEQ :: Evidence

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

Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_014034	222,054	0	C→G	0.9%	11.2 ‑2.7 ‑0.0	1467	P201P (CCG→CCC)	RCAP_rcc00185	hypothetical protein
Reads supporting (aligned to +/- strand): new base (3/12): ref base (650/800): total (654/813)
Fisher's exact test strand distribution p-value = 6.76e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.72e-01

Only 100 of 2306 total aligned reads displayed.


TCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGT  >  NC_014034/221954‑222154


tgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGGCCGGGGGGTCGGTCGGGGCCGTGGGGTCGGCCGGGTCCGTGGGAgcggtc                                                                                                      >  2:1656896/3‑101 (MQ=11)


 gGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGtcggtcg                                                                                                     <  1:2878596‑M1/100‑1 (MQ=255)


  gtatcgatggggtccgtgggaccggtcgcgtccgtgggACCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcgg                                                                                                    <  2:17605/99‑1 (MQ=11)


  gggttggtcggggccgtgggatccggcggggccgggggggcCGGCGGGTTCGTGGGGTCGGGCGGGTCCCTGGGATCGGTCGGGTCCGTGGGAtcggtcgg                                                                                                    <  1:3082231/96‑1 (MQ=2)


   gttctgtcggggccgtgggatcggtcgggtccgtgggatcgGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCggg                                                                                                   <  2:590840/96‑1 (MQ=11)


    gtcggtcgggtccgtgggctccgtcggacccgtgggaTCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGt                                                                                                  <  2:5095306/100‑1 (MQ=11)


    atcggtgggatccgtgggatcggtcgggtccgtgggATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGGGGGATCGGTCGGGt                                                                                                  >  1:2380409/1‑101 (MQ=11)


    atcggtcgggtccgtggcatcggtcgggtccgtgggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcggtcggat                                                                                                  <  2:5195244/101‑3 (MQ=14)


     ggagtgcggtccgtgggagcggtcggggccgtggggtccgtcgGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTc                                                                                                 <  1:2647263/94‑1 (MQ=14)


      cggtcgggtccgtgggatcggtcgggtccgcgggatCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCGGTCGGGTCCGTGGGATCGGTCGGGGcc                                                                                                >  1:4896000/1‑101 (MQ=11)


      cggtcggatccgtgggatcggtcgggtccgtgggatCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGCTCGGTCGGATcc                                                                                                >  1:533861/1‑101 (MQ=11)


       ggtcgggtccgtgggatcggtcgggtccgtgggatcGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGatccg                                                                                               <  2:2622972‑M2/101‑6 (MQ=255)


        gtcgggtacgtgggatcggtcgggtccgtagggtcgGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGt                                                                                              <  2:4281950/101‑1 (MQ=11)


         tcgggtccgtgggatcggtcggatccgtgggatcggTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGatccgtg                                                                                             <  2:3788433‑M2/101‑8 (MQ=255)


         tcgggtccgtgggatcggtcggatccgtgggatcggTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGatccgta                                                                                             <  1:2267761‑M2/101‑8 (MQ=255)


         cgggagccgtgggatcggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTGGGGTCCGTGGGATCGGCCGGGTCCGTGGGAGCGGTCGGGTccttg                                                                                             >  2:2031728/7‑98 (MQ=11)


          cgggtccgtgggatcggtcgggtccgtgggatcggtCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGatccgtgg                                                                                            <  1:2036144‑M2/101‑9 (MQ=255)


          cgggtccgtgggatcggtcgggtccgtcggatcggtCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGGTCGGTCGGGTCCGtgg                                                                                            >  2:1953420/1‑99 (MQ=14)


          cggatccgtgggatcggtcgggtccgtgggatcggtCGGATCCGTGGGATCGGTTGGGGCCGTGGGATCCGCCGGACCCGCGGGATCAGTCTGGTCCgcgg                                                                                            >  2:1562558/1‑98 (MQ=14)


           gggtccgttcgatcggtcgggtccgtggggtccgtcGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGgg                                                                                           <  2:1859149/101‑1 (MQ=11)


           gggtccgtgggatcggttgggtccgtggggtccgtcGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGgg                                                                                           <  1:4804519/101‑1 (MQ=11)


           gggtccgtgggatcggtcgggtccgtggggtccgtcGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGgg                                                                                           >  2:3992583/1‑101 (MQ=14)


           gggtccgtgggatcggtcgggtccgtgggatcggtcGGATCCGTGGGATCGGTCGGGTCCGTGGGCTCGGTCGGATCCGTGGGGGCCGTGGGGTcggccgg                                                                                           >  1:1964455/1‑95 (MQ=11)


            gatccgtgggatcggtcgggtccgttggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGGCGGATCCGTGGGATCGGTCGGATCCGTGGGa                                                                                          >  1:3943457/3‑101 (MQ=11)


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


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


              tccgtgggatcggtcgggtccgtcggatcggtgggATCCGTGGGATCGGTCGGGGCCGGGGGATCGGTCGGGTCCGGCGGGTCGGTCGGGTCCggggggtg                                                                                        >  2:4339911/1‑94 (MQ=11)


              tccgtgggatcggtcggatccgtgggatcggtcggGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGatccgtgggatc                                                                                        <  1:1579621‑M2/101‑13 (MQ=255)


               ccgttggctcggtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcg                                                                                       <  1:2531574/101‑1 (MQ=11)


               ccgtgggatcggtcgggtccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGgtccgtgggatcg                                                                                       >  1:1105026‑M2/1‑88 (MQ=255)


                cgtggggtcggtcgggtccgtggggtccgtcggGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcgg                                                                                      <  1:5005178/101‑1 (MQ=11)


                cgtgggatcggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGgtccgtgggatcgg                                                                                      <  2:1332337‑M2/101‑15 (MQ=255)


                 gtgggatcggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGCtcggt                                                                                     >  1:658664/1‑101 (MQ=14)


                 gtgggatcggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGGCCgagcgggccgt                                                                                     >  2:5660233/1‑91 (MQ=11)


                  tgggatcggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGatccgtgggatcggtc                                                                                    >  2:4187659‑M2/1‑85 (MQ=255)


                  cgggatcggtcggttccgtggtatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGAtccgtg                                                                                    <  2:2334551/100‑2 (MQ=11)


                   gggatcggtcgggtccgtcggatcggtgggATCCGTGGGATCGGTCGGGTCCGGGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcg                                                                                   >  2:2714051/1‑101 (MQ=14)


                    ggatcgttcgggtccgtggaatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCCTGGGATCGGTCGGGTCCGTGGGAtcggtcgg                                                                                  <  1:3263142/101‑1 (MQ=14)


                       tcggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGatccgtgggatcggtcgggtc                                                                               >  2:3393829‑M2/1‑80 (MQ=255)


                       gcgggcgggtccgtgggatcggtcggaTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGatccgtgggatcggtcggatc                                                                               >  2:2877620‑M2/2‑80 (MQ=255)


                        cgggcgggtccggggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcggatcc                                                                              <  1:1362309/101‑1 (MQ=11)


                          gtcggatccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGgtccgtgggatcggtcgggtccgg                                                                            <  2:3668044‑M2/101‑25 (MQ=255)


                            cgggtccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGGTCGGTCGGATCCGTGGGATCGGTCGGatccgtgggatcggtcggatccgtgg                                                                          <  1:5509599‑M2/101‑27 (MQ=255)


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


                              ggtccgtaggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATGGGTCGGGTCGGTGGGGTCCgtggggtccgtcggg                                                                        >  2:3096678/1‑101 (MQ=11)


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


                                 ccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGgtccgtgggatcggtcggatccgtgggatcg                                                                     >  2:2230037‑M2/1‑70 (MQ=255)


                                    tgggatcggtcggATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTg                                                                  >  1:4350263/1‑101 (MQ=11)


                                     gggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGatccgtgggatcggtcgggtccgtgggatcggtcg                                                                 >  1:1662845‑M2/1‑66 (MQ=255)


                                      gggtccgtcggATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGGCCGTGGGATCGGTCGGGCCCGTCGGATCGGTggg                                                                >  1:773891/7‑101 (MQ=11)


                                      ggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGGGGGATCGGTCGGGTCCGGGGGGACCGTcgg                                                                >  2:35237/1‑101 (MQ=11)


                                        gtccgtcggATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGAt                                                              <  1:5097013/97‑1 (MQ=11)


                                        gtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGac                                                              >  1:5756475/5‑100 (MQ=11)


                                         tcggtcggATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGtc                                                             <  1:3505051/101‑1 (MQ=11)


                                         tccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGAtc                                                             >  2:5293375/4‑101 (MQ=11)


                                          cggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGgtccgtgggatcggtcgggtccgtgggatcggtggggtcc                                                            <  2:4069869‑M2/101‑41 (MQ=255)


                                              cgggtagtgcggATGGGTCTGGGTCGTGGGATCGGTCGGGTCCGTGGGGTCGGTGGGATCTGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGtccgtcg                                                        <  1:2147636/91‑1 (MQ=11)


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


                                               ggATCCGTGGGGTCGGTCGGGTCCGTGGGAGCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGAtccgtggg                                                       <  1:2978759/101‑1 (MQ=11)


                                               ggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGAtccgtggg                                                       <  1:1302496/101‑1 (MQ=11)


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


                                                 gtcCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggat                                                     >  2:1162392‑M2/2‑54 (MQ=255)


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


                                                      tGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCGGTGGGGTCCGTGGGAtcggtgggggccggg                                                >  2:5586102/1‑94 (MQ=11)


                                                       cggATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTgg                                               >  2:3906945/2‑101 (MQ=255)


                                                          gtcggtcggGTCGGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGAt                                            <  1:4668293/100‑1 (MQ=17)


                                                          aTCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggt                                            <  2:4548309‑M2/101‑57 (MQ=255)


                                                          aTCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGatccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcggat                                            >  1:2754301‑M2/1‑45 (MQ=255)


                                                          aTCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGAGCGGTCGGATCCGTGGGAt                                            >  2:212497/1‑101 (MQ=18)


                                                              gtggggTCCGTGGGATCGGTCGGGTCCGTGGGATAGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggt                                        >  1:875042/4‑101 (MQ=255)


                                                               gcgggTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtg                                       >  2:4643506‑M2/2‑40 (MQ=255)


                                                                cggGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgg                                      <  2:878358‑M2/101‑63 (MQ=255)


                                                                cggGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGgtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgg                                      <  1:4015967‑M2/101‑63 (MQ=255)


                                                                cggGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGatccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcggatccgtag                                      >  2:4419066‑M2/1‑39 (MQ=255)


                                                                cggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtggggtccgtcggatccgtgg                                      <  1:1859259‑M2/101‑63 (MQ=255)


                                                                 gggTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcgg                                     >  1:1742272/1‑101 (MQ=17)


                                                                    tccgtggggtccgtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatc                                  >  1:4207384‑M2/1‑35 (MQ=255)


                                                                    tccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtggggtcggtcggatccgtgggatc                                  <  1:3739670‑M2/101‑67 (MQ=255)


                                                                    tccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatc                                  >  2:1453447‑M2/1‑35 (MQ=255)


                                                                     ccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcg                                 >  2:500238‑M2/1‑34 (MQ=255)


                                                                     ccgtgggatcggtcggatccgtgggatcggtcggatccgtaggatcggtcgggtccgtaggatcggtcgggtccgtcggatcggtcggatccgtgggatcg                                 <  1:5227592‑M2/101‑68 (MQ=255)


                                                                     ccGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTcc                                 >  1:5651628/1‑101 (MQ=255)


                                                                      cGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCAGTCGGGTCCGTGGGATCGGACGGGTACGTAGGAGCGGTCGGGGCCg                                >  2:4763767/1‑101 (MQ=2)


                                                                      cGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCg                                <  2:5392439/101‑1 (MQ=11)


                                                                      cGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGGTCGGTCGGGTCCg                                >  2:3020036/1‑101 (MQ=17)


                                                                        tgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtggggtccgtcggatccgtgggatcggtc                              <  1:3636911‑M2/101‑71 (MQ=255)


                                                                          ggATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGGCCGTggg                            >  2:1268504/1‑101 (MQ=17)


                                                                            aTCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTgggat                          >  1:4285224/1‑99 (MQ=17)


                                                                              cggtcggATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTgggatcg                        <  1:966945/101‑5 (MQ=14)


                                                                               ggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccg                       >  2:3175766‑M2/1‑24 (MQ=255)


                                                                               ggtcggGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCTGGTCCGTcgggtccg                       >  1:1480221/1‑97 (MQ=11)


                                                                               ggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTaggatcgg                       >  2:3278181/1‑96 (MQ=14)


                                                                               ggtagggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccg                       <  2:4903656‑M2/101‑78 (MQ=255)


                                                                                 tgggATCCGTGGGATCGGTCGGGTCGGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTaggatcggtc                     <  1:4559626/99‑8 (MQ=14)


                                                                                 gcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtg                     >  1:4109274‑M2/2‑22 (MQ=255)


                                                                                  tgggtCCGTGGGATCGGTGGGGTCCGTGGGATCCGTCGGGTCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCtggtgg                    <  2:1058379/97‑1 (MQ=17)


                                                                                   ggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgttgg                   <  2:5472603‑M2/101‑82 (MQ=255)


                                                                                    gATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTgggatcggtcggg                  >  1:1219196/1‑91 (MQ=11)


                                                                                            ggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtggggtccgtcggatccgtgggatcggtcggatccgtgggatcggtcgg          >  1:3389840‑M2/1‑11 (MQ=255)


                                                                                                cggtcggATCCGTGGGATCGGTCGGGTCCGTGGGACCGGTGGGGTCCGTGGGATCGGTCGGGCCGGCGGGGGAGGTGGGGACGGTGGGGTCGGCCGCACCg      >  2:2484437/1‑101 (MQ=2)


TCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGT  >  NC_014034/221954‑222154

Base quality scores:


ATCG

< 3 ≤


ATCG

< 15 ≤


ATCG

< 23 ≤


ATCG

< 31 ≤


ATCG

< 40 ≤


ATCG