BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_014034	2,566,625	C→G	3.5%	intergenic (+97/+97)	purD → / ← sdaA	phosphoribosylamine‑‑glycine ligase/L‑serine ammonia‑lyase

Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_014034	2,566,625	0	C→G	3.5%	23.9 ‑3.6 ‑0.0	459	intergenic (+97/+97)	purD/sdaA	phosphoribosylamine‑‑glycine ligase/L‑serine ammonia‑lyase
Reads supporting (aligned to +/- strand): new base (15/1): ref base (300/143): total (315/144)
Fisher's exact test strand distribution p-value = 2.75e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.72e-01

Only 100 of 486 total aligned reads displayed.


AGACGAAAC.GGGGGCGCTGCCCCC.CTC.CTG..CTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCAACCGCCTCAG  >  NC_014034/2566528‑2566725


agaCGAAACGGGGGGCGCTGCCCCC.CTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGccc                                                                                                      >  2:5400934/1‑101 (MQ=21)


agaCGAAACGGGGGGCGCTGCCCCC.CTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGccc                                                                                                      <  2:5783653/101‑1 (MQ=21)


 gaCGAAACGGGGGGCGCTGCCCCC.CTC.TTGGCCTTCGGTCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCg                                                                                                     <  1:6561591/101‑1 (MQ=21)


 gaCGAAACGGGGGGCGCTGCCCCC.CTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCg                                                                                                     <  1:5189968/101‑1 (MQ=21)


   cGAAACGGGGGGCGCTGCCCCC.CTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGaa                                                                                                   <  2:1336170/101‑1 (MQ=21)


   cGAAACGGGGGGCGCTGCCCCC.CTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGaa                                                                                                   >  2:3579234/1‑101 (MQ=21)


    gAAACGGGGGGCGCTGCCCCC.CTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGaaa                                                                                                  <  2:3432603/101‑1 (MQ=21)


    gAAACGGGGGGCGCTGCCCCC.CTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGaaa                                                                                                  <  2:8572887/101‑1 (MQ=21)


     aaacggggggCGCTGCCCCC.CTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTACTTTCATCTTGCCCGAAat                                                                                                 >  2:2036643/6‑101 (MQ=21)


        cgggggtCGCTGCCCCC.CTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATc                                                                                              >  1:4705207/3‑101 (MQ=25)


        cggggggCGCTGCCCCC.NTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATc                                                                                              <  2:1542981/99‑1 (MQ=25)


        cggggggCGCTGCCCCC.CTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATc                                                                                              <  1:7023015/99‑1 (MQ=25)


        cggggggCGCTGCCCCC.CTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTGGCCCGAAATATc                                                                                              <  1:8148520/99‑1 (MQ=21)


        cagggggCGCTG.CCCC.CTC.TTGGCCTTTGGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCgggg                                                                                        >  1:6752015/3‑101 (MQ=0)


         ggggggCGCTGCCCCC.CTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATcc                                                                                             >  1:3739418/2‑101 (MQ=25)


         ggggggCGCTGCCCCC.CTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATcc                                                                                             >  2:2802692/2‑101 (MQ=25)


         agggggCGCTG.CCCC.CTC.TTGGCCTTTGGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggggg                                                                                       >  2:2453432/2‑101 (MQ=0)


          gggggCGCTG.CCCC.CTC.TTGGCCTTTGGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCgggggg                                                                                      >  1:3994954/1‑100 (MQ=0)


           ggggCGCTGCCCCC.CTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCg                                                                                           <  1:1914339/101‑1 (MQ=25)


            gggCGCTGCCCCC.CTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCgg                                                                                          >  1:327226/1‑101 (MQ=25)


             ggCGCTGCCCCCTCTT.GGC..CTTTGGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggggggagc                                                                                   <  2:7410274/101‑5 (MQ=2)


              gcgcTGCCCCC.CTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCgggg                                                                                        <  1:4373604/101‑1 (MQ=37)


                gcTGCCCCC.CTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGa                                                                                      >  2:8807066/1‑101 (MQ=37)


                  tGCCCCC.CTCTCGGC.CTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGc                                                                                    <  2:870211/101‑1 (MQ=37)


                   gCCCCC.CTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGcc                                                                                   <  1:2260739/101‑1 (MQ=37)


                           ctcttggccTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTg                                                                             >  1:5517517/9‑101 (MQ=255)


                           ctcttggccTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTg                                                                             <  2:4951422/93‑1 (MQ=255)


                           ctcttggccTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTg                                                                             <  1:390175/93‑1 (MQ=255)


                           ctcttggccTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTg                                                                             >  1:2635413/9‑101 (MQ=255)


                            tcttggccTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGTTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgc                                                                            <  2:6982025/94‑1 (MQ=255)


                             cttggccTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgcg                                                                           >  1:4533035/7‑101 (MQ=255)


                             cttggccTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgcg                                                                           >  1:8611290/7‑101 (MQ=255)


                             cttggccTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgcg                                                                           <  2:230235/95‑1 (MQ=255)


                             cttggccTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgcg                                                                           >  1:312018/7‑101 (MQ=255)


                              ttggccTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgcgg                                                                          <  2:7283970/96‑1 (MQ=255)


                                      cGGGCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGc                                                                  >  1:8791151/1‑101 (MQ=255)


                                      cGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGc                                                                  >  2:4434243/1‑101 (MQ=255)


                                      cGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGc                                                                  >  2:8305152/1‑101 (MQ=255)


                                       atctaaTTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCa                                                                 <  2:5682487/97‑1 (MQ=255)


                                        gCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAg                                                                >  2:3703901/1‑101 (MQ=255)


                                         gcAATTCGCCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAgc                                                               <  2:6915705/100‑1 (MQ=255)


                                          cAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAgcg                                                              >  1:2492006/1‑101 (MQ=255)


                                           aaTTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCaagat                                                             >  1:5682694/1‑97 (MQ=255)


                                            aTTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGcc                                                            >  2:5851998/1‑101 (MQ=255)


                                            aTTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGcc                                                            >  1:1808494/1‑101 (MQ=255)


                                              tCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGcccc                                                          <  1:3627092/101‑1 (MQ=255)


                                              tCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGcccc                                                          >  1:2027374/1‑101 (MQ=255)


                                               cACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGccccc                                                         <  2:525401/101‑1 (MQ=255)


                                                        ggATAGTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGccc                                                <  2:6471164/101‑1 (MQ=255)


                                                          atatTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGa                                              >  1:7950746/1‑101 (MQ=255)


                                                            atTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGAGCCCGAcg                                            >  1:1027874/1‑101 (MQ=255)


                                                            atTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCACCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGAcg                                            <  2:3648492/101‑1 (MQ=255)


                                                               tCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGAcgcgc                                         >  1:5301417/1‑101 (MQ=255)


                                                                cGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGAcgcgcg                                        >  2:7776059/1‑101 (MQ=255)


                                                                  ggCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCgcggc                                      >  2:5137340/1‑101 (MQ=255)


                                                                   gCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCgcggcg                                     <  1:5786396/101‑1 (MQ=255)


                                                                    cAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGcggcgg                                    >  2:2607609/1‑101 (MQ=255)


                                                                     aaGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCggg                                   >  2:4737337/1‑101 (MQ=255)


                                                                     aaGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCggg                                   >  1:7369723/1‑101 (MQ=255)


                                                                     aaGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCggg                                   >  1:9059880/1‑101 (MQ=255)


                                                                      aGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCCGCGGGc                                  >  1:2817684/1‑101 (MQ=255)


                                                                       gATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCt                                 >  1:5032025/1‑101 (MQ=255)


                                                                       gATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCt                                 >  1:8413504/1‑101 (MQ=255)


                                                                        aTAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCtt                                <  1:3495009/101‑1 (MQ=255)


                                                                        aTAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCtt                                <  1:2337306/101‑1 (MQ=255)


                                                                        aTAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCtt                                <  1:1221413/101‑1 (MQ=255)


                                                                         tAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTg                               >  2:7913552/1‑101 (MQ=255)


                                                                          aaaGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgc                              >  2:5078052/1‑101 (MQ=255)


                                                                          aaaGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgc                              >  1:8987065/1‑101 (MQ=255)


                                                                           aagaagGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgcc                             >  1:1518723/1‑101 (MQ=255)


                                                                              aagGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgcccgc                          <  2:2108404/101‑1 (MQ=255)


                                                                                ggTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgcccgccc                        <  2:1808441/101‑1 (MQ=255)


                                                                                 gTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCa                       <  2:33302/101‑1 (MQ=255)


                                                                                 gTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCa                       >  2:8100206/1‑101 (MQ=255)


                                                                                 gTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCCAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCa                       <  1:6446572/101‑1 (MQ=255)


                                                                                   cTTTTCGTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAag                     >  1:7107288/1‑101 (MQ=255)


                                                                                   cTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAag                     >  2:3893908/1‑101 (MQ=255)


                                                                                       tCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAaggagg                 >  2:6039871/1‑101 (MQ=255)


                                                                                        ccTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCCAGGCggg                >  1:4247455/1‑101 (MQ=255)


                                                                                        ccTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAggg                >  2:4091151/1‑101 (MQ=255)


                                                                                        ccTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAggg                >  2:8002183/1‑101 (MQ=255)


                                                                                        ccTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAggg                <  1:5367222/101‑1 (MQ=255)


                                                                                        ccTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCCTGCCCCCCCAAGGAggg                >  2:8965958/1‑101 (MQ=255)


                                                                                         cTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAACGAAGgc               >  1:7513821/1‑101 (MQ=255)


                                                                                          tttCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGgcg              >  1:6645298/1‑101 (MQ=255)


                                                                                          tttCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGgcg              >  2:4625369/1‑101 (MQ=255)


                                                                                            tCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCa            >  1:6156198/1‑101 (MQ=255)


                                                                                             cATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCac           >  1:2949143/1‑100 (MQ=255)


                                                                                              aTCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCAAc          >  2:1029970/1‑101 (MQ=255)


                                                                                              aTCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCAAc          <  2:1213620/101‑1 (MQ=255)


                                                                                              aTCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCAAc          <  2:6359022/101‑1 (MQ=255)


                                                                                               tCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCAAcc         >  1:7735873/1‑101 (MQ=255)


                                                                                               tCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCAAcc         <  1:5925208/101‑1 (MQ=255)


                                                                                                cTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCAGCGGGCTTGCCCGCCCAAGGAGGGCGCAACCg        >  2:6764846/1‑101 (MQ=255)


                                                                                                  tGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCAACCGcc      >  2:661567/1‑101 (MQ=255)


                                                                                                  tGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCAACCGcc      >  1:1339130/1‑101 (MQ=255)


                                                                                                   gCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCAACCGCCt     >  2:5545077/1‑101 (MQ=255)


                                                                                                    cccGAAATATCCCGGGGGAGCCCGCTGGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCAACCGCCTc    <  2:29577/101‑1 (MQ=255)


                                                                                                     ccGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCAACCGCCTca   <  2:2895287/101‑1 (MQ=255)


                                                                                                     ccGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCAACCGCCTca   <  2:3182206/101‑1 (MQ=255)


AGACGAAAC.GGGGGCGCTGCCCCC.CTC.CTG..CTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCAACCGCCTCAG  >  NC_014034/2566528‑2566725

Base quality scores:


ATCG

< 3 ≤


ATCG

< 20 ≤


ATCG

< 27 ≤


ATCG

< 33 ≤


ATCG

< 41 ≤


ATCG