BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_014034	2,566,613	T→C	3.6%	intergenic (+85/+109)	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,613	0	T→C	3.6%	10.3 ‑3.8 ‑0.6	220	intergenic (+85/+109)	purD/sdaA	phosphoribosylamine‑‑glycine ligase/L‑serine ammonia‑lyase
Reads supporting (aligned to +/- strand): new base (8/0): ref base (119/93): total (127/93)
Fisher's exact test strand distribution p-value = 2.21e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.45e-01

Only 100 of 246 total aligned reads displayed.


GCGCGCGCTCTGAGACGAAAC.GGGGGCGCTGCCCCCC.TCCTG..CTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCG  >  NC_014034/2566516‑2566713


 cgcgcgcTCTGAGACGAAACGGGGGGCGCTGCCCCCC.TCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCtt                                                                                                     <  2:487191/101‑1 (MQ=21)


 cgcgcgcTCTGAGACGAAACGGGGGGCGCTGCCCCCC.TCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCtt                                                                                                     <  2:1685975/101‑1 (MQ=21)


   cgcgcTCTGAGACGAAACGGGGGGCGCTGCCCCCC.TCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTc                                                                                                   <  2:138080/101‑1 (MQ=21)


       ctctGAGACGAAACGGGGGGCGCTGCCCCCC.TCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCt                                                                                               >  2:4182655/1‑101 (MQ=21)


        tctGAGACGAAACGGGGGGCGCTGCCCCCC.TCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCtt                                                                                              <  2:2700968/101‑1 (MQ=21)


        tctGAGACGAAACGGGGGGCGCTGCCCCCC.TCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCtt                                                                                              >  2:505281/1‑101 (MQ=21)


          tGAGACGAAACGGGGGGCGCTGCCCCCC.TCTTGGCCTTCGGCCAATTCACCCCCCGGGGTATTTGCGGCAAGATAAAGAAAGGCTTTTCCTTTCATCtttc                                                                                            >  2:4047009/1‑99 (MQ=17)


           gagaCGAAACGGGGGGCGCTGCCCCCC.TCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGcc                                                                                           <  2:3940742/101‑1 (MQ=21)


           gagaCGAAACGGGGGGCGCTGCCCCCC.TCTTGGCCTTCGGCCAATTCACCCCCCCGGGTATTTTTGGCAAAAAAAAAAAAGGGTTTTTCTTTTCTCttttc                                                                                           >  1:1691504/1‑98 (MQ=2)


             gaCGAAACGGGGGGCGCTGCCCCCC.TCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCg                                                                                         <  2:883993/101‑1 (MQ=21)


               cGAAACGGGGGGCGCTGCCCCCC.TCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGaa                                                                                       <  1:4319021/101‑1 (MQ=21)


                gAAACGGGGGGCGCTGCCCCCC.TCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGaaa                                                                                      >  1:2914501/1‑101 (MQ=21)


                  aacggggggCGCTGCCCCCC.TCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAata                                                                                    <  1:4113473/97‑1 (MQ=25)


                   acgtggggCGCTGCCCCCC.TCTTGGCCTTCGGCCAATTCCCCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAatat                                                                                   <  1:2573218/97‑1 (MQ=25)


                   acggggggCGCTTCCCCCCTTTTTGC.CTTCGGCCAATTTACCCCCCGGGATATTTCGGGCAATATTAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAatat                                                                                   <  1:1372020/98‑1 (MQ=25)


                   acggggggCGCTGCCCCCC.TCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAatat                                                                                   <  1:2995379/98‑1 (MQ=25)


                   acggggggCGCTGCCCCCC.TCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAatat                                                                                   <  1:699687/98‑1 (MQ=25)


                     ggggggCGCTGCCCCCC.TCTTGGCCTTCGGCCAATTCCCCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATcc                                                                                 <  1:300878/100‑1 (MQ=25)


                     ggggggCGCTGCCCCCC.TCTTGGCCTTCGGCCAATTCACCCCCCGGGGGATTTTGGGGCAGAGAAAAAAAGGGTTTTTCTTTTTTTTTTCCCCCAATATcc                                                                                 >  1:984749/2‑101 (MQ=2)


                     ggggggCGCTGCCCCCC.TCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATcc                                                                                 >  2:1003043/2‑101 (MQ=25)


                     ggggggCGCTGCCCCCC.TCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATcc                                                                                 >  2:2639398/2‑101 (MQ=25)


                      gggggCGCTGCCCCCC.TCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATccc                                                                                >  2:4490350/1‑101 (MQ=25)


                      gggggCGCTGCCCCCC.TCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATccc                                                                                >  2:1746777/1‑101 (MQ=25)


                      gggggCGCTG.CCCCC.TCTTGGCCTTTGGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCgggggg                                                                          >  2:641046/1‑100 (MQ=0)


                       ggggCGCTGCCCCCC.TCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCg                                                                               <  1:4588220/101‑1 (MQ=25)


                       ggggCGCTGCCCCCC.TCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCg                                                                               <  1:1536663/101‑1 (MQ=25)


                       ggggCGCTGCCCCCC.TCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCg                                                                               <  2:2213454/101‑1 (MQ=25)


                        gggCGCTGCCCCCC.TCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCgg                                                                              <  1:4687976/101‑1 (MQ=37)


                          gcgcTGCCCCCC.TCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCgggg                                                                            <  1:4436796/101‑1 (MQ=37)


                           cgcTGCCCCCC.TCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCggggg                                                                           <  2:3278284/101‑1 (MQ=37)


                            gcTGCCCCCC.TCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGa                                                                          >  2:1913922/1‑101 (MQ=37)


                             cTGCCCCCC.TCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGGCAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAg                                                                         >  2:2414426/1‑101 (MQ=25)


                              tGCCCCCC.TCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGc                                                                        <  1:3489476/101‑1 (MQ=37)


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


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


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


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


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


                                         ttggccTTCGGCCAATTCCCCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgcgg                                                              <  1:3975753/96‑1 (MQ=255)


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


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


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


                                           ggccTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgcgggc                                                            >  1:2051894/4‑101 (MQ=255)


                                             ccTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgcgggcgg                                                          >  2:3570027/2‑101 (MQ=255)


                                             ccTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTgcgggcgg                                                          >  1:4260008/2‑101 (MQ=255)


                                              cTTCGGCCAATTCCCCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGcgggcggg                                                         <  2:2925804/101‑1 (MQ=255)


                                              cTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGGCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGcgggcggg                                                         >  2:428186/1‑101 (MQ=255)


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


                                                tCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCggggg                                                       <  2:2990951/101‑1 (MQ=255)


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


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


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


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


                                                      aaTTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGCGGGGGGCAgcgc                                                 >  2:2371746/1‑101 (MQ=255)


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


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


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


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


                                                               cccGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAgc                                        >  2:1996107/1‑101 (MQ=255)


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


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


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


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


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


                                                                            gggCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCgcgg                           <  2:3446511/101‑1 (MQ=255)


                                                                             ggCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCgcggc                          <  1:1996224/101‑1 (MQ=255)


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


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


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


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


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


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


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


                                                                                   gatctgaagGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCtt                    <  2:600043/96‑1 (MQ=255)


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


                                                                                   aTAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGcct                    >  1:1431052/1‑99 (MQ=255)


                                                                                    tAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGGTTg                   >  1:3222290/1‑101 (MQ=255)


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


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


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


                                                                                      aagaagGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTtcc                 >  2:1626063/1‑101 (MQ=255)


                                                                                       agaagGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgccc                >  1:2134527/1‑101 (MQ=255)


                                                                                        gaagGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCttagatc               >  1:600107/1‑96 (MQ=255)


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


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


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


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


                                                                                             ttttttCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCaa          <  1:1883910/99‑1 (MQ=255)


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


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


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


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


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


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


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


                                                                                                tttCCTTTCATCTTGCCCGAAATATCCAGGGGGAGTCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTAGCCCGCCCAaaga       >  2:3168933/1‑101 (MQ=37)


                                                                                                 ttCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCTaggag      >  2:1895315/1‑101 (MQ=255)


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


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


                                                                                                     tttCATCTTGCCCGAAATACCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCAAAGGAGGgcg  <  1:1742882/101‑1 (MQ=255)


GCGCGCGCTCTGAGACGAAAC.GGGGGCGCTGCCCCCC.TCCTG..CTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCG  >  NC_014034/2566516‑2566713

Base quality scores:


ATCG

< 3 ≤


ATCG

< 18 ≤


ATCG

< 26 ≤


ATCG

< 33 ≤


ATCG

< 41 ≤


ATCG