BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_014034	2,566,613	T→C	3.9%	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.9%	14.2 ‑4.0 ‑0.2	256	intergenic (+85/+109)	purD/sdaA	phosphoribosylamine‑‑glycine ligase/L‑serine ammonia‑lyase
Reads supporting (aligned to +/- strand): new base (10/0): ref base (160/86): total (170/86)
Fisher's exact test strand distribution p-value = 1.82e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.86e-01

Only 100 of 281 total aligned reads displayed.


GCGCGCGCTCTGAGACGAAAC.GGGGG.CGCTGCCCCCCTCCTG..CTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCG  >  NC_014034/2566516‑2566713


gcgcgcgcTCTGAGACGAAACGGGGGG.CGCTGCCCCCCTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCt                                                                                                      <  2:4093268/101‑1 (MQ=21)


gcgcaccgtctGAGACGAAAC.GGGGGCCGCTGCCCCCCTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCt                                                                                                      <  1:943921/93‑1 (MQ=21)


  gcgcgcTCTGAGACGAAACGGGGGG.CGCTGCCCCCCTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCttt                                                                                                    <  2:214335/101‑1 (MQ=21)


  gcgcgcTCTGAGACGAAACGGGGGG.CGCTGCCCCCCTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCttt                                                                                                    <  1:4619477/101‑1 (MQ=21)


   cgcgcTCTGAGACGAAACGGGGGG.CGCTGCCCCCCTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTc                                                                                                   <  2:4004732/101‑1 (MQ=21)


   cgcgcTCTGAGACGAAACGGGGGG.CGCTGCCCCCCTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTc                                                                                                   <  2:1014071/101‑1 (MQ=21)


     cgcTCTGAGACGAAACGGGGGG.CGCTGCCCCCCTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCAt                                                                                                 <  1:1003710/101‑1 (MQ=21)


             gaCGAAACGGGGGG.CGCTGCCCCCCTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCg                                                                                         <  2:2053853/101‑1 (MQ=21)


                gAAACGGGGGG.CGCTGCCCCCCTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGaaa                                                                                      <  1:1183627/101‑1 (MQ=21)


                 aaacgggggg.CGCTGCCCCCCTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAat                                                                                     >  1:802128/6‑101 (MQ=21)


                   gcaggggg.CGCTG.CCCCCTCTTGGCCTTTGGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggg                                                                             >  1:4762757/4‑101 (MQ=0)


                   gcaggggg.CGCTG.CCCCCTCTTGGCCTTTGGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCggg                                                                             >  2:5474240/4‑101 (MQ=0)


                    cgggggg.CGCTGCCCCCCTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATc                                                                                  <  2:3847742/99‑1 (MQ=25)


                    caggggg.CGCTG.CCCCCTCTTGGCCTTTGGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCgggg                                                                            >  1:728053/3‑101 (MQ=0)


                      ggggg.CGCTGCCCCCCTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATccc                                                                                <  1:2242135/101‑1 (MQ=25)


                      ggggg.CGCTG.CCCCCTCTTGGCCTTTGGCCAATTCA.CCCCTGGGATATTTCGGGCAAGAT...GAA.ATCATGCCCCTTCATCTTGCCGCAAATATCCCgggggg                                                                          >  2:3245842/1‑100 (MQ=0)


                        ggg.CGCTGCCCCCCTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCgg                                                                              <  1:1746994/101‑1 (MQ=37)


                          g.cgcTGCCCCCCTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCgggg                                                                            >  1:5236898/1‑101 (MQ=37)


                          g.cgcTGCCCCCCTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCgggg                                                                            <  2:4236780/101‑1 (MQ=37)


                             gcTGCCCCCCTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGa                                                                          >  2:4478594/1‑101 (MQ=37)


                             gcTGCCCCCCTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGa                                                                          >  1:5471362/1‑101 (MQ=37)


                              cTGCCCCCCTCTTGGCCTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAg                                                                         >  1:553744/1‑101 (MQ=37)


                                  ccctcttggg..CTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCt                                                                   >  2:5275561/10‑101 (MQ=255)


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


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


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


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


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


                                                 cGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGc                                                      <  2:1296680/101‑1 (MQ=255)


                                                  ggCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGcg                                                     >  1:3883801/1‑100 (MQ=255)


                                                  ggCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGcc                                                     >  2:3973502/1‑100 (MQ=255)


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


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


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


                                                         tCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGACCGCTTGCGGGCGGGGGCAGCGcccc                                              >  1:2170162/1‑101 (MQ=255)


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


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


                                                                ccGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAgcg                                       >  1:4391008/1‑101 (MQ=255)


                                                                ccGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAgcg                                       >  1:3105209/1‑101 (MQ=255)


                                                                ccGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAgcg                                       >  1:1950834/1‑101 (MQ=255)


                                                                  gggATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGcc                                     <  1:169593/101‑1 (MQ=255)


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


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


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


                                                                    gATATTTCGGGCGAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCGGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCg                                   <  1:5366197/101‑1 (MQ=255)


                                                                    gATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCg                                   <  2:5036603/101‑1 (MQ=255)


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


                                                                     atatTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGa                                  <  2:3680679/101‑1 (MQ=255)


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


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


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


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


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


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


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


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


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


                                                                             ggCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGC.GGGGCAGGGGCCCCCGAGCCCCCGACGCccggcg                         >  1:3620813/1‑101 (MQ=37)


                                                                              gCAAGATAAAGAAGGTCTTTTCCTTTCAGCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCgcggcg                         <  1:336027/101‑1 (MQ=255)


                                                                               cAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGcggggg                        >  1:3247773/1‑101 (MQ=255)


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


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


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


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


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


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


                                                                                    tAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGCCGCGCGGCGGGCttt                   >  2:533794/1‑100 (MQ=255)


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


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


                                                                                    tAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGACTTg                   >  1:87118/1‑101 (MQ=255)


                                                                                     aaaGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCTCCCGACGCGCGGCGGGCTTgc                  >  1:3541890/1‑101 (MQ=255)


                                                                                      aagaagGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGACGCGCCCGACGCGGGGCGGGCTTccc                 >  1:4545012/1‑101 (MQ=255)


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


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


                                                                                       agaagGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCtgtccc                >  2:2023071/1‑96 (MQ=255)


                                                                                       agaagGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgccc                <  2:5168874/101‑1 (MQ=255)


                                                                                       agaagGTCTTTTCCTTTCATCTTGCCCGAAATACCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgccc                <  1:4984098/101‑1 (MQ=255)


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


                                                                                         aagGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGGTTgcccgc              >  1:2628142/1‑101 (MQ=255)


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


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


                                                                                          agGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTgcccgcc             >  1:4711779/1‑101 (MQ=255)


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


                                                                                             tCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGGCCaa          >  1:5736815/1‑101 (MQ=255)


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


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


                                                                                             tCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCaa          <  2:2816096/101‑1 (MQ=255)


                                                                                              cTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCaac         >  1:2882303/1‑100 (MQ=255)


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


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


                                                                                                 ttCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGGTTGCCCGCCCAaggag      >  2:5068179/1‑101 (MQ=255)


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


                                                                                                  tCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAaggagg     <  2:4183561/101‑1 (MQ=255)


                                                                                                  tCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCCCTTGCGGGCGGGGGCACCGCCCCCGCGGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAagaggg     >  1:1837167/1‑97 (MQ=37)


                                                                                                   ccTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCCAGGAggg    >  2:3151995/1‑101 (MQ=255)


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


                                                                                                   ccTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCCCCCAAGGAggg    >  2:4593039/1‑101 (MQ=255)


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


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


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


GCGCGCGCTCTGAGACGAAAC.GGGGG.CGCTGCCCCCCTCCTG..CTTCGGCCAATTCACCCCCCGGGATATTTCGGGCAAGATAAAGAAGGTCTTTTCCTTTCATCTTGCCCGAAATATCCCGGGGGAGCCCGCTTGCGGGCGGGGGCAGCGCCCCCGAAGCGCCCGACGCGCGGCGGGCTTGCCCGCCCAAGGAGGGCG  >  NC_014034/2566516‑2566713

Base quality scores:


ATCG

< 3 ≤


ATCG

< 20 ≤


ATCG

< 27 ≤


ATCG

< 33 ≤


ATCG

< 41 ≤


ATCG