BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_014034	2,566,555	C→T	100%	intergenic (+27/+167)	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,555	0	C→T	100.0%	513.8	183	intergenic (+27/+167)	purD/sdaA	phosphoribosylamine‑‑glycine ligase/L‑serine ammonia‑lyase
Reads supporting (aligned to +/- strand): new base (48/134): ref base (0/0): total (49/134)

Only 100 of 206 total aligned reads displayed.


TATGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAAC.GGGG.G.CGCTGCCCCC.CTC.CT..G..CTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCC  >  NC_014034/2566457‑2566644


tATGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCAGCCGCGACATCGGCTGGGGCGCGCTCTGAGACGAAAAA.GGGG.G.GGCGGCCCCc.ccc.ct                                                                                              >  2:2423052/1‑101 (MQ=255)


      aTGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..ggcctt                                                                                      <  2:1516213/101‑6 (MQ=38)


       tGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..ggccttc                                                                                     <  2:3120128/101‑7 (MQ=38)


       tGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..ggccttc                                                                                     <  2:1666402/101‑7 (MQ=38)


        ggTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..ggccttcg                                                                                    <  1:1347379/101‑8 (MQ=38)


         gTCGACGCGATCGTCTGGCACGAGGGGTTCTGCCGCCGCGACAGCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..ggccttcgg                                                                                   <  2:724236/101‑9 (MQ=25)


          tCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..ggccttcggc                                                                                  <  2:1880317/101‑10 (MQ=37)


           cGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..ggccttccgcc                                                                                 >  1:1942831/1‑91 (MQ=38)


           cGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGcc                                                                                 <  1:3105290/101‑1 (MQ=37)


              tggtaTCGACGGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAAt                                                                              <  1:836120/97‑1 (MQ=25)


                cGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.G.GGGTGCCCCC.CCC.TT..GGCCTTCGGCCCATTc                                                                            >  1:2241320/1‑101 (MQ=255)


                 gATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCa                                                                           <  1:400313/101‑1 (MQ=37)


                 gATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCa                                                                           <  2:3008801/101‑1 (MQ=37)


                  aTCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCAc                                                                          <  2:4250848/101‑1 (MQ=37)


                   tCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCAcc                                                                         <  1:3710347/101‑1 (MQ=37)


                      aCTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCAccccc                                                                      <  1:2423007/101‑1 (MQ=25)


                       cTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCAcccccc                                                                     <  1:619811/101‑1 (MQ=25)


                       cTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCCCTGAGACGAAAAC.GGGG.G.GGCCGCCCCC.CCC.TT..GGCCTTCCGCCCAATCAcccccc                                                                     >  1:626548/1‑101 (MQ=255)


                        tGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCACTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCg                                                                    <  1:3484577/101‑1 (MQ=25)


                         ggCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCgg                                                                   <  1:1271329/101‑1 (MQ=25)


                         ggCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCgg                                                                   <  2:2107644/101‑1 (MQ=25)


                           cccGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGa                                                                 <  2:3361022/101‑1 (MQ=25)


                            ccGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCCCTGAGACG.AAACGGGGG.G.GGCTGCCCCC.CCC.TT..GGCCTTCGGCCAACTCACCCCCCggcgt                                                                >  1:2260402/1‑98 (MQ=255)


                             cGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGata                                                               >  2:402940/1‑101 (MQ=25)


                              gAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGatat                                                              <  1:3598415/101‑1 (MQ=25)


                              gAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGatat                                                              <  2:4448925/101‑1 (MQ=25)


                              gAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGatat                                                              <  2:4552337/101‑1 (MQ=25)


                              gAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGatat                                                              <  1:2293697/101‑1 (MQ=25)


                               aGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.G.GGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCCGGAAAtt                                                             >  1:4528074/1‑101 (MQ=21)


                                ggggTTCTGCCGCCGCGACATCGGCTGGCGCGCGCCCTGAGACG.AAACGGGGG.G.GGCTGCCCCC.CTC.TT..GGCCCTCCGCCACTTCAcccccccgcagaatt                                                            >  2:1761936/1‑92 (MQ=19)


                                 gggTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTc                                                           <  1:4680640/101‑1 (MQ=21)


                                 ccgatCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTc                                                           <  1:480242/97‑1 (MQ=21)


                                   gTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCgg                                                         <  2:3424315/101‑1 (MQ=21)


                                   gTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCgg                                                         <  2:1329466/101‑1 (MQ=21)


                                      cTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCa                                                      <  2:3220211/101‑1 (MQ=21)


                                        gccgccGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAg                                                    <  2:3584721/101‑1 (MQ=21)


                                        gccgccGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGCCaga                                                    >  2:3688659/1‑99 (MQ=21)


                                         ccgccgCGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGa                                                   <  1:2875269/101‑1 (MQ=21)


                                           gccgcGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTCTTT..GC.CTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATa                                                 >  2:3639782/1‑101 (MQ=21)


                                           gccgcGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATa                                                 <  2:3280760/101‑1 (MQ=21)


                                           gccgcGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATa                                                 <  1:734598/101‑1 (MQ=21)


                                           gccgcGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATa                                                 <  2:500633/101‑1 (MQ=21)


                                           gccgcGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATa                                                 <  1:2919991/101‑1 (MQ=21)


                                            ccgcGACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATaa                                                <  2:3892920/101‑1 (MQ=21)


                                             cgcgACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATaaa                                               <  2:2464653/101‑1 (MQ=21)


                                             cgcgACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATaaa                                               <  1:415380/101‑1 (MQ=21)


                                              gcgACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAaag                                              <  1:2561639/101‑1 (MQ=21)


                                               cgACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAaaga                                             <  2:1075959/101‑1 (MQ=21)


                                                gACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCCATTCACCCCCCGGGATATTTCGGGCAAGATAaagaa                                            >  1:2991139/1‑101 (MQ=21)


                                                gACATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAaagaa                                            <  1:4697722/101‑1 (MQ=21)


                                                 aCATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAaagaag                                           <  1:2399048/101‑1 (MQ=21)


                                                  cATCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAgg                                          <  2:4613648/101‑1 (MQ=21)


                                                   aTCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGt                                         <  1:1860053/101‑1 (MQ=21)


                                                    tCGGCTGGCGCGCGCTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTc                                        <  1:3907223/101‑1 (MQ=21)


                                                      ggCTGGCGCGCGCTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtt                                      <  1:1065005/101‑1 (MQ=21)


                                                       gCTGGCGCGCGCTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCttt                                     <  2:1134929/101‑1 (MQ=21)


                                                        cTGGCGCGCGCTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCtttt                                    <  1:872705/101‑1 (MQ=21)


                                                         tGGCGCGCGCTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGCTATTTCGGGCCAGAAAAAGAAGGGCTTTTc                                   >  1:2592339/1‑101 (MQ=18)


                                                         tGGCGCGCGCTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTc                                   <  1:3839415/101‑1 (MQ=21)


                                                         tGGCGCGCGCTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTc                                   <  2:2288983/101‑1 (MQ=21)


                                                          ggCGCGCGCTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTcc                                  >  1:4448274/1‑101 (MQ=21)


                                                          ggCGCGCGCTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTcc                                  <  2:3467003/101‑1 (MQ=21)


                                                          ggCGCGCGCTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTcc                                  <  1:198833/101‑1 (MQ=21)


                                                          ggCGCGCGCTATGAGACG.AGAC.GGGG.GCCGCTGCCCCC.CTC.TTGGG..CTTCGGCCAATTCTCCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTcc                                  <  2:875624/101‑1 (MQ=21)


                                                          ggCGCGCGCGCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTcc                                  <  1:4392873/101‑1 (MQ=21)


                                                            cgcgcgcTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCtt                                <  2:2227890/101‑1 (MQ=21)


                                                            cgcgcgcTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCtt                                <  2:487191/101‑1 (MQ=21)


                                                            cgcgcgcTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCtt                                <  1:2036391/101‑1 (MQ=21)


                                                             gcgcgcTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCttt                               >  1:4676021/1‑101 (MQ=21)


                                                             gcgcgcTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCttt                               <  2:2717771/101‑1 (MQ=21)


                                                             gcgcgcTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCttt                               <  2:1689403/101‑1 (MQ=21)


                                                              cgcgcTCTGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTc                              <  2:138080/101‑1 (MQ=21)


                                                                  ctctGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCt                          >  2:4182655/1‑101 (MQ=21)


                                                                    ctGAGACG.AACC.GGGGCG.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTg                        <  2:620015/101‑1 (MQ=21)


                                                                    ctGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTg                        <  2:2156545/101‑1 (MQ=21)


                                                                     tGAGACG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGGTATTTGCGGCAAGATAAAGAAAGGCTTTTCCTTTCATCtttc                       >  2:4047009/1‑99 (MQ=17)


                                                                      gagaCG.AAACGGGGG.G.CGCTGCCCCC.CTCTTT..GC.CTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGcc                      <  1:1485951/101‑1 (MQ=21)


                                                                      gagaCG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGcc                      <  2:3940742/101‑1 (MQ=21)


                                                                      gagaCG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCCGGGTATTTTTGGCAAAAAAAAAAAAGGGTTTTTCTTTTCTCttttc                      >  1:1691504/1‑98 (MQ=2)


                                                                         aCG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCCGCATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGa                   >  2:2699577/1‑101 (MQ=21)


                                                                          cG.AAACGGGGG.G.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGaa                  <  1:4319021/101‑1 (MQ=21)


                                                                           cgaaccggggg.g.CCCGGCCCCC.CTC.TTG.G..CTTTCGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGaaa                 <  1:3492803/94‑1 (MQ=21)


                                                                             aaacggggg.g.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAat                <  1:1396642/96‑1 (MQ=21)


                                                                             aaacggggg.g.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAat                >  1:4631004/6‑101 (MQ=21)


                                                                               acggggg.g.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAatat              <  1:699687/98‑1 (MQ=25)


                                                                                cagggg.g.CGCTGCCCCCTCTT.GG..C..CTTTGGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCgggg       >  2:2979150/3‑101 (MQ=0)


                                                                                cagggg.g.CGCTG.CCCC.CTC.TT..GGCCTTTGGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCgggg       >  2:506050/3‑101 (MQ=0)


                                                                                 ggggg.g.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCCCCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATcc            <  1:300878/100‑1 (MQ=25)


                                                                                 ggggg.g.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATcc            >  2:1003043/2‑101 (MQ=25)


                                                                                 ggggg.g.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATcc            <  2:3403742/100‑1 (MQ=25)


                                                                                 ggggg.g.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATcc            >  2:2639398/2‑101 (MQ=25)


                                                                                 agggg.g.CGCTG.CCCC.CTC.TT..GGCCTTTGGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggggg      >  1:2800801/2‑101 (MQ=0)


                                                                                 agggg.g.CGCTG.CCCC.CTC.TT..GGCCTTTGGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggggg      >  2:174876/2‑101 (MQ=0)


                                                                                   ggg.g.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCg          <  1:3343035/101‑1 (MQ=25)


                                                                                   ggg.g.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCg          <  1:4588220/101‑1 (MQ=25)


                                                                                    gg.g.CGCTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCgg         <  1:4511503/101‑1 (MQ=37)


                                                                                       g.cgcTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCgggg       >  2:4465774/1‑101 (MQ=37)


                                                                                         cgcTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCggggg      <  1:4182901/101‑1 (MQ=37)


                                                                                          gcTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGa     >  2:1913922/1‑101 (MQ=37)


                                                                                           cTGCCCCC.CTC.TT..GGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAg    <  2:1898783/101‑1 (MQ=37)


TATGCGATGGTCGATGCGATCGACTGGCCCGAGGGGTTCTGCCGCCGCGACATCGGCTGGCGCGCGCTCTGAGACG.AAAC.GGGG.G.CGCTGCCCCC.CTC.CT..G..CTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCC  >  NC_014034/2566457‑2566644

Base quality scores:


ATCG

< 3 ≤


ATCG

< 18 ≤


ATCG

< 27 ≤


ATCG

< 33 ≤


ATCG

< 40 ≤


ATCG