BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_014034	2,566,626	G→C	11.0%	intergenic (+98/+96)	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,626	0	G→C	11.0%	42.7 ‑3.6 ‑0.1	172	intergenic (+98/+96)	purD/sdaA	phosphoribosylamine‑‑glycine ligase/L‑serine ammonia‑lyase
Reads supporting (aligned to +/- strand): new base (18/1): ref base (108/45): total (126/46)
Fisher's exact test strand distribution p-value = 2.66e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.45e-01

Only 100 of 191 total aligned reads displayed.


GACGAAAC.GG.GGGCGCTGCCCCC.CTCCTG..CTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCAACCGCCTCAGC  >  NC_014034/2566529‑2566726


gaCGAAACGGG.GGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCg                                                                                                      <  1:4327812/101‑1 (MQ=21)


gaCGAAACGGG.GGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCg                                                                                                      <  1:1430547/101‑1 (MQ=21)


 aCGAAACGGG.GGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGa                                                                                                     >  1:3476249/1‑101 (MQ=21)


   gAAACGGG.GGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGaaa                                                                                                   <  1:4492940/101‑1 (MQ=21)


    aaaC.GGCGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAat                                                                                                  <  1:766337/101‑1 (MQ=25)


     aacggg.gggCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAata                                                                                                 <  2:2855141/97‑1 (MQ=25)


     aacggg.gggCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAata                                                                                                 >  2:2512064/5‑101 (MQ=25)


     aacggg.gggCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGGCTTTTCCTTTCATCTTTCCCGAAata                                                                                                 >  1:3728697/5‑101 (MQ=21)


      gcagg.gggCGCTG.CCCC.CTCTTGGCCTTTGGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggg                                                                                          >  1:1479095/4‑101 (MQ=0)


      gcagg.gggCGCTG.CCCC.CTCTTGGCCTTTGGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggg                                                                                          >  2:1564774/4‑101 (MQ=0)


      acggg.gggCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGACTTTTTCTTTCATCTTTCCCGAAaaat                                                                                                >  1:266360/4‑101 (MQ=18)


       cggg.gggCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATc                                                                                               <  1:2141213/99‑1 (MQ=25)


       cagg.gggCGCTG.CCCC.CTCTTGGCCTTTGGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCgggg                                                                                         >  1:1296958/3‑101 (MQ=0)


       cagg.gggCGCTG.CCCC.CTCTTGGCCTTTGGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCgggg                                                                                         >  2:3063307/3‑101 (MQ=0)


        ggg.gggCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATcc                                                                                              <  2:3092070/100‑1 (MQ=25)


        agg.gggCGCTGCCCCCTCTTGGC..CTTTGGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggggg                                                                                        >  1:3136529/2‑101 (MQ=0)


        agg.gggCGCTG.CCCC.CTCTTGGCCTTTGGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggggg                                                                                        >  1:4097678/2‑101 (MQ=0)


        agg.gggCGCTG.CCCC.CTCTTGGCCTTTGGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggggg                                                                                        >  2:2620134/2‑101 (MQ=0)


        agg.gggCGCTG.CCCC.CTCTTGGCCTTTGGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggggg                                                                                        >  2:2229699/2‑101 (MQ=0)


        agg.gggCGCTG.CCCC.CTCTTGGCCTTTGGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCCTGCCGCAAATATCCCggggg                                                                                        >  1:4292733/2‑101 (MQ=2)


         gg.gggCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCCCCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATccc                                                                                             <  2:1496807/101‑1 (MQ=25)


          g.gggCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCg                                                                                            <  1:4310501/101‑1 (MQ=25)


            gggCGCTG.CCCC.CTCTTGGCCTTTGGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggggggag                                                                                     >  2:548785/1‑98 (MQ=2)


             ggCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCggg                                                                                          <  1:1041361/101‑1 (MQ=37)


                   gCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGAGATTTCGGGCCAGATAAAAAAGGGCTTTTTCTTTTATCTTTTCCGGAAAATCCCCGGGGGGcc                                                                                    >  2:1755715/1‑101 (MQ=12)


                        ccctcttggccTTCGGCCAATTCCCCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCt                                                                                <  1:1867746/91‑1 (MQ=38)


                         cctcttggccTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTCGCCCGAAATATCCCGGGGGAGCCCGCtt                                                                               <  1:816745/92‑1 (MQ=37)


                           tcttggccTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgc                                                                             <  2:1920230/94‑1 (MQ=255)


                           tcttggccTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgc                                                                             <  2:3659272/94‑1 (MQ=255)


                             ttggccTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgcgg                                                                           >  1:1260221/6‑101 (MQ=255)


                             ttggccTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgcgg                                                                           >  2:45593/6‑101 (MQ=255)


                               ggccTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgcgggc                                                                         <  1:1825339/98‑1 (MQ=255)


                               ggccTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgcgggc                                                                         >  2:4424767/4‑101 (MQ=255)


                                   ttCGGCCAATTCCCCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCgggg                                                                     <  1:3789572/101‑1 (MQ=255)


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


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


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


                                      ggCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCa                                                                  >  2:2112354/1‑101 (MQ=255)


                                       gCCAATTCACCCCCCGGGAAATTTTGGGGAAGATAAAAAAGGGCTTTTTCTTTTATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTTCGGGCggggggcg                                                                 >  2:1659953/1‑98 (MQ=25)


                                        ccAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAgc                                                                <  1:2470017/101‑1 (MQ=255)


                                            ttCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGccc                                                            >  2:4514412/1‑101 (MQ=255)


                                               cccccccGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCg                                                         <  1:3378599/100‑1 (MQ=255)


                                                    ccGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAgcg                                                    <  2:1846416/101‑1 (MQ=255)


                                                     cGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAgcgc                                                   >  1:3907795/1‑101 (MQ=255)


                                                      gggATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGcc                                                  >  1:3653583/1‑101 (MQ=255)


                                                        gattttcggggAAAGAAAAAGAAGTTTTTTTCCTTTCATCTTCCCCAAAATACCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCg                                                <  1:895890/93‑1 (MQ=255)


                                                        gATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCg                                                >  2:2206216/1‑101 (MQ=255)


                                                        gATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCg                                                >  2:1849825/1‑101 (MQ=255)


                                                         atatTTCGGGCAAGATAAAGAAGGCCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGCGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCTGAAGCGCCCGa                                               >  2:2298931/1‑101 (MQ=255)


                                                          tatTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGAc                                              >  2:414178/1‑101 (MQ=255)


                                                           atTTCGGGCAAGATACAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGAcg                                             >  1:3933773/1‑101 (MQ=255)


                                                           atTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGAcg                                             >  2:1165890/1‑101 (MQ=255)


                                                             ttCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGAcgcg                                           >  2:2343191/1‑101 (MQ=255)


                                                                gggCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCgcgg                                        >  2:3731532/1‑101 (MQ=255)


                                                                gggCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCgcgg                                        >  1:4079646/1‑101 (MQ=255)


                                                                gggCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCgcgg                                        >  1:4645759/1‑101 (MQ=255)


                                                                  gCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCgcggcg                                      >  1:2564334/1‑101 (MQ=255)


                                                                  gCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCgcggcg                                      >  1:1529398/1‑101 (MQ=255)


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


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


                                                                    gatctAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCggg                                    <  2:3863356/97‑1 (MQ=255)


                                                                    aaGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCCGCCCCCCCCAACCCCCCCAccccccgccga                                    >  1:882646/1‑91 (MQ=25)


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


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


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


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


                                                                        tAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCCTg                                >  1:3503542/1‑101 (MQ=255)


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


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


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


                                                                         aaaGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGAGCCCCCGAAGCACCCGACGCGCGGCGGGCGTgc                               >  1:3595745/1‑101 (MQ=255)


                                                                         aaaGAAGGTCTGTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgc                               <  2:75792/101‑1 (MQ=255)


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


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


                                                                           agaagGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGGGGCGGGCTTgccc                             >  2:3031894/1‑101 (MQ=255)


                                                                            gaagGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGCAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgcccg                            >  2:4105022/1‑101 (MQ=255)


                                                                             aagGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGGAGCGCCCCCGAAGCGCCCGACGCGCGGCCGGGTTgccccc                           >  2:1600504/1‑99 (MQ=255)


                                                                             aagGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgcccgc                           >  2:2410286/1‑101 (MQ=255)


                                                                               ggTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgcccgccc                         >  2:725416/1‑101 (MQ=255)


                                                                                 tCTTTCCCTTTCATCTGGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCaa                       <  2:712653/101‑1 (MQ=255)


                                                                                  cTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAag                      <  2:3475991/101‑1 (MQ=255)


                                                                                   ttttCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAagg                     >  1:2637512/1‑101 (MQ=255)


                                                                                    tttCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAagga                    >  1:302387/1‑101 (MQ=255)


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


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


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


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


                                                                                          ttCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGgcgc              >  1:2876540/1‑101 (MQ=255)


                                                                                          ttCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGgcgc              >  1:1586826/1‑101 (MQ=255)


                                                                                            cagcTTGCCCGAAATACCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCaa            <  1:2052508/98‑1 (MQ=255)


                                                                                            cATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCaa            >  2:1147358/1‑101 (MQ=255)


                                                                                              tCTTGCCCAAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCAAcc          <  2:4083621/101‑1 (MQ=255)


                                                                                                 tccccGAAATATCCCGGGGGAGCCCGCTCGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCAACCGcc       <  1:1071456/99‑1 (MQ=255)


                                                                                                 tGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCAACCGcc       <  2:4227869/101‑1 (MQ=255)


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


                                                                                                   cgcGAAATATCCCGGGGGAGCCCGATTGCGGGCGGGGGCCGCCCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCAACCGCCTc     <  2:512725/99‑1 (MQ=255)


                                                                                                   cccGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCAACCGCCTc     >  1:873177/1‑101 (MQ=255)


                                                                                                    ccGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCAACCGCCTca    >  2:1054570/1‑101 (MQ=255)


                                                                                                     cGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCAACCGCCTcag   <  2:2564160/101‑1 (MQ=255)


                                                                                                      gAAATATCCCGGGGCAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCTCGGCGGGCTTTCCCGCCCAAGGAGGGCGCAACCGCCTcagc  >  2:4623669/1‑101 (MQ=255)


GACGAAAC.GG.GGGCGCTGCCCCC.CTCCTG..CTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCAACCGCCTCAGC  >  NC_014034/2566529‑2566726

Base quality scores:


ATCG

< 3 ≤


ATCG

< 18 ≤


ATCG

< 27 ≤


ATCG

< 33 ≤


ATCG

< 41 ≤


ATCG