BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_014034	2,566,626	G→C	3.5%	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	3.5%	10.7 ‑3.6 ‑0.4	231	intergenic (+98/+96)	purD/sdaA	phosphoribosylamine‑‑glycine ligase/L‑serine ammonia‑lyase
Reads supporting (aligned to +/- strand): new base (8/0): ref base (137/86): total (145/86)
Fisher's exact test strand distribution p-value = 2.71e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 3.65e-01

Only 100 of 256 total aligned reads displayed.


GACGAAAC.GGGGGCGCTGCCCCC.CTCCTG..CTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCAACCGCCTCAGC  >  NC_014034/2566529‑2566726


  cGAAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGaa                                                                                                    <  1:4319021/101‑1 (MQ=21)


   gAAACGGGGGGCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGaaa                                                                                                   >  2:4211402/1‑101 (MQ=21)


    aaacggggggCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAat                                                                                                  <  1:209702/96‑1 (MQ=21)


     aacggggggCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAata                                                                                                 <  1:4113473/97‑1 (MQ=25)


     aacggggggCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGGCTTTTCCTTTCATCTTGCCCGAAata                                                                                                 >  1:552176/5‑101 (MQ=21)


      gcagggggCGCTG.CCCC.CTCTTGGCCTTTGGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggg                                                                                          >  1:2279811/4‑101 (MQ=0)


      acgtggggCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCCCCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAatat                                                                                                <  1:2573218/97‑1 (MQ=25)


       cagggggCGCTGCCCCCTCTTGGC..CTTTGGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCgggg                                                                                         >  2:2979150/3‑101 (MQ=0)


       cagggggCGCTG.CCCC.CTCTTGGCCTTTGGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCgggg                                                                                         >  2:506050/3‑101 (MQ=0)


        ggggggCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATcc                                                                                              >  2:1003043/2‑101 (MQ=25)


        agggggCGCTG.CCCC.CTCTTGGCCTTTGGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggggg                                                                                        >  1:4604638/2‑101 (MQ=0)


        agggggCGCTG.CCCC.CTCTTGGCCTTTGGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggggg                                                                                        >  1:2925338/2‑101 (MQ=0)


         gggggCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATccc                                                                                             <  1:942584/101‑1 (MQ=25)


         gggggCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATccc                                                                                             <  2:3215671/101‑1 (MQ=25)


         gggggCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATccc                                                                                             >  2:1746777/1‑101 (MQ=25)


          ggggCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCg                                                                                            <  2:2213454/101‑1 (MQ=25)


          ggggCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCg                                                                                            <  1:3343035/101‑1 (MQ=25)


           gggCGCTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCgg                                                                                           <  1:4687976/101‑1 (MQ=37)


             gcgcTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCgggg                                                                                         <  1:4436796/101‑1 (MQ=37)


             gcgcTGCCCCC.CTCTTGGCCTGCGGCCTAGTCTCCCCCCGGGATATTTCGGGCAATATAAAGAAGGTCTTTTCCTTTCATCTTGCACGAAATATCCCgggg                                                                                         <  2:1636853/101‑1 (MQ=21)


              cgcTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCggggg                                                                                        <  1:4182901/101‑1 (MQ=37)


               gcTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGa                                                                                       >  2:1913922/1‑101 (MQ=37)


                cTGCCCCC.CTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGGCAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAg                                                                                      >  2:2414426/1‑101 (MQ=25)


                       ccctcttggccTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCt                                                                                <  1:1003795/91‑1 (MQ=255)


                        cctcttggccTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCtt                                                                               <  1:2133041/92‑1 (MQ=255)


                        cctcttggccTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCtt                                                                               <  1:1833509/92‑1 (MQ=255)


                          tcttggccTTCGGCCAATTCTCCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgc                                                                             <  1:2797394/94‑1 (MQ=255)


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


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


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


                             tggccTTCGGCCAATTCACCCCCCGGGGTATTTTGGGCAAGATAAAGAAGGGCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCCCTTgtggg                                                                          >  1:4664144/5‑101 (MQ=25)


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


                               gccTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgcgggcg                                                                        <  1:2303930/99‑1 (MQ=255)


                                ccTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgcgggcgg                                                                       <  2:1942714/100‑1 (MQ=255)


                                   tCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTTATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCggggg                                                                    >  1:1105181/1‑101 (MQ=255)


                                    cGGCCAACTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGc                                                                   <  2:4520084/101‑1 (MQ=255)


                                     ggCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCa                                                                  <  2:863892/101‑1 (MQ=255)


                                      gCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTTCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGcgg                                                                 >  1:3666957/1‑99 (MQ=255)


                                      gCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGGAg                                                                 >  2:2839463/1‑101 (MQ=255)


                                         aaTTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAgcgc                                                              >  2:3497631/1‑101 (MQ=255)


                                         aaTTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCCCTTGCGGGCGGGGGGAgcgc                                                              >  2:1851767/1‑101 (MQ=255)


                                             cACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGGAGCGccccc                                                          >  2:3336327/1‑101 (MQ=255)


                                             cACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGccccc                                                          >  2:268428/1‑101 (MQ=255)


                                              aCCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCg                                                         <  1:1585892/101‑1 (MQ=255)


                                               ccccccGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGa                                                        <  2:2260269/101‑1 (MQ=255)


                                                    cGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAgcgc                                                   <  2:740043/101‑1 (MQ=255)


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


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


                                                      ggATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGccc                                                 <  1:699301/101‑1 (MQ=255)


                                                      ggATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGccc                                                 <  2:1578690/101‑1 (MQ=255)


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


                                                         tatTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCCAAGCGCCCGAc                                              >  2:4477629/1‑101 (MQ=255)


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


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


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


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


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


                                                                  cAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGcggcgg                                     <  1:383857/101‑1 (MQ=255)


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


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


                                                                    aGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGc                                   <  2:3585484/101‑1 (MQ=255)


                                                                    aGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGc                                   >  2:898488/1‑101 (MQ=255)


                                                                      aTAAAGAAGGTGTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCtt                                 <  1:2371879/101‑1 (MQ=255)


                                                                      aTAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCtt                                 >  2:634062/1‑101 (MQ=255)


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


                                                                      aTAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCtt                                 >  1:1576416/1‑101 (MQ=255)


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


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


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


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


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


                                                                        aaaGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGAGCGGCGGGCTTgc                               >  2:255543/1‑101 (MQ=255)


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


                                                                           gaagGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgcccg                            <  2:626484/101‑1 (MQ=255)


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


                                                                            aagGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGATTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTtcccgc                           >  2:410070/1‑101 (MQ=255)


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


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


                                                                                tCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCaa                       >  1:1503155/1‑101 (MQ=255)


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


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


                                                                                  ttttCCTTTCATCTTGCCCGAAATATCCCGGGGGACCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAcgg                     >  2:416715/1‑101 (MQ=255)


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


                                                                                    ttCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAaggag                   >  2:3050333/1‑101 (MQ=255)


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


                                                                                       cTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGgc                >  1:2837340/1‑101 (MQ=255)


                                                                                        tttCATCTTGCCCGAAATATCCCGGGGGCGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAACTCCCGACGCGCGGCGGGGTGGCCCGCCCAAGGAGGgcg               >  1:913522/1‑101 (MQ=38)


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


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


                                                                                          tCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCa             <  1:1485959/101‑1 (MQ=255)


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


                                                                                          tCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGAGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCACAAGGAGGGCGCa             >  2:1484237/1‑101 (MQ=255)


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


                                                                                            aTCTTGCCCGAAATATCCGGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCAAc           <  2:4448024/101‑1 (MQ=255)


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


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


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


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


                                                                                                     gAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGGGGGATTGCCAGCCCAAGGAGGGCGCAACCGCCAcagg  >  1:202164/1‑100 (MQ=38)


                                                                                                     gAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCAACCGCCTcagc  >  1:3512229/1‑101 (MQ=255)


GACGAAAC.GGGGGCGCTGCCCCC.CTCCTG..CTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCGCAACCGCCTCAGC  >  NC_014034/2566529‑2566726

Base quality scores:


ATCG

< 3 ≤


ATCG

< 18 ≤


ATCG

< 25 ≤


ATCG

< 33 ≤


ATCG

< 41 ≤


ATCG