BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	143,709	T→C	12.5%	Y53C (TAC→TGC)	pcnB ←	poly(A) polymerase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	143,709	0	T	C	12.5%	210.8 / 18.6	80	Y53C (TAC→TGC)	pcnB	poly(A) polymerase I
Reads supporting (aligned to +/- strand): ref base T (31/39); new base C (3/7); total (34/46)
Fisher's exact test for biased strand distribution p-value = 5.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

AACAGGTCGCGCACGCCGCCGCCAACCAGCCAGGCTTCGTATCCCGCTTTATTGAGCCTGTACATTACCTTCAGGGCATTTTCACTGATATCTTTGCGGGAAATAGCATGCTGCTCACGCGGGATCACCGT  >  Exported/143649‑143779

aaCAGGTCGCGCACGCCGCCGCCAACCAGCCAGGCTTCGTATCCCGCTTTATTGAGCCTGTACATTACCtt                                                              >  1:4627750/1‑71 (MQ=255)

 aCAGGTCGCGCACGCCGCCGCCAACCAGCCAGGCTTCGTATCCCGCTTTATTGAGCCTGTACATTACCTTc                                                             <  1:3051129/71‑1 (MQ=255)

  cAGGTCGCGCACGCCGCCGCCAACCAGCCAGGCTTCGTATCCCGCTTTATTGAGCCTGTACATTACCTTCa                                                            <  1:817628/71‑1 (MQ=255)

  cAGGTCGCGCACGCCGCCGCCAACCAGCCAGGCTTCGTATCCCGCTTTATTGAGCCTGTACATTACCTTCa                                                            <  1:3111643/71‑1 (MQ=255)

  cAGGTCGCGCACGCCGCCGCCAACCAGCCAGGCTTCGTATCCCGCTTTATTGAGCCTGTACATTACCTTCa                                                            <  1:1821573/71‑1 (MQ=255)

    ggTCGCGCACGCCGCCGCCAACCAGCCAGGCTTCGTATCCCGCTTTATTGAGCCTGTACATTACCTTc                                                             <  1:2645170/68‑1 (MQ=255)

       cgcgCACGCCGCCGCCAACCAGCCAGGCTTCGTATCCCGCTTTATTGAGCCTGTACATTACCTTCAGGGCa                                                       >  1:4699433/1‑71 (MQ=255)

           cACGCCGCCGCCAACCAGCCAGGCTTCGTATCCCGCTTTATTGAGCCTGTACATTACCTTCAGGGCAtttt                                                   <  1:977416/71‑1 (MQ=255)

           cACGCCGCCGCCAACCAGCCAGGCTTCGTATCCCGCTTTATTGAGCCTGTACATTACCTTCAGGGCAtttt                                                   <  1:3596701/71‑1 (MQ=255)

           cACGCCGCCGCCAACCAGCCAGGCTTCGTATCCCGCTTTATTGAGCCTGTACATTACCTTCAGGGCAtttt                                                   <  1:4877537/71‑1 (MQ=255)

           cACGCCGCCGCCAACCAGCCAGGCTTCGTATCCCGCTTTATTGAGCCTGTACATTACCTTCAGGGCAtttt                                                   <  1:1408063/71‑1 (MQ=255)

            aCGCCGCCGCCAACCAGCCAGGCTTCGTATCCCGCTTTATTGAGCCTGTACATTACCTTCAGGGCAtttt                                                   <  1:3160321/70‑1 (MQ=255)

                    gccAACCAGCCAGGCTTCGTATCCCGCTTTATTGAGCCTGTACATTACCTTCAGGGCa                                                       >  1:1021309/1‑58 (MQ=255)

                    gccAACCAGCCAGGCTTCGTATCCCGCTTTATTGAGCCTGTACATTACCTTCAGGGCATTTTCACTGatat                                          >  1:3358819/1‑71 (MQ=255)

                    gccAACCAGCCAGGCTTCGTATCCCGCTTTATTGAGCCTGTACATTACCTTCAGGGCATTTTCACTGatat                                          >  1:3833584/1‑71 (MQ=255)

                    gccAACCAGCCAGGCTTCGTATCCCGCTTTATTGAGCCTGTACATTACCTTCAGGGCATTTTCACTGatat                                          >  1:2822889/1‑71 (MQ=255)

                    gccAACCAGCCAGGCTTCGTATCCCGCTTTATTGAGCCTGTACATTACCTTCAGGGCATTTTCACTGata                                           >  1:2827582/1‑70 (MQ=255)

                    gccAACCAGCCAGGCTTCGTATCCCGCTTTATTGAGCCTGTACATTACCTTCAGGGCATTTTCACTGata                                           >  1:3418805/1‑70 (MQ=255)

                    gccAACCAGCCAGGCTTCGTATCCCGCTTTATTGAGCCTGTACATTACCTTCAGGGCATTTTCACTGata                                           >  1:1975105/1‑70 (MQ=255)

                     ccAACCAGCCAGGCTTCGTATCCCGCTTTATTGAGCCTGTACATTACCTTCAGGGCATTTTCACTGATATc                                         <  1:1727076/71‑1 (MQ=255)

                      cAACCAGCCAGGCTTCGTATCCCGCTTTATTGAGCCTGt                                                                        <  1:717576/39‑1 (MQ=255)

                      cAACCAGCCAGGCTTCGTATCCCGCTTTATTGAGCCTGt                                                                        <  1:3956005/39‑1 (MQ=255)

                      cAACCAGCCAGGCTTCGTATCCCGCTTTATTGAGCCTGt                                                                        <  1:4735224/39‑1 (MQ=255)

                      cAACCAGCCAGGCTTCGTATCCCGCTTTATTGAGCCTGt                                                                        <  1:339253/39‑1 (MQ=255)

                            gccagGCTTCGTATCCCGCTTTATTGAGCCTGCACATTACCTTCAGGGCAtttt                                                   <  1:773341/54‑1 (MQ=255)

                            gccagGCTTCGTATCCCGCTTTATTGAGCCTGCACATTACCTTCAGGGCAtttt                                                   <  1:1534198/54‑1 (MQ=255)

                            gccagGCTTCGTATCCCGCTTTATTGAGCCTGCACATTACCTTCAGGGCAtttt                                                   <  1:4378646/54‑1 (MQ=255)

                            gccagGCTTCGTATCCCGCTTTATTGAGCCTGCACATTACCTTCAGGGCAtttt                                                   <  1:2837300/54‑1 (MQ=255)

                            gccagGCTTCGTATCCCGCTTTATTGAGCCTGCACATTACCTTCAGGGCAtttt                                                   <  1:1262512/54‑1 (MQ=255)

                            gccagGCTTCGTATCCCGCTTTATTGAGCCTGCACATTACCTTCAGGGCAtttt                                                   <  1:2364069/54‑1 (MQ=255)

                            gccagGCTTCGTATCCCGCTTTATTGAGCCTGCACATTACCTTCAGGGCAtttt                                                   <  1:1449258/54‑1 (MQ=255)

                              cagGCTTCGTATCCCGCTTTATTGAGCCTGCACATTACCTTCAGGGCATTTTCa                                                 >  1:1495553/1‑54 (MQ=255)

                              cagGCTTCGTATCCCGCTTTATTGAGCCTGCACATTACCTTCAGGGCATTTTCa                                                 >  1:4025087/1‑54 (MQ=255)

                              cagGCTTCGTATCCCGCTTTATTGAGCCTGCACATTACCTTCAGGGCATTTTCa                                                 >  1:4046421/1‑54 (MQ=255)

                                ggCTTCGTATCCCGCTTTATTGAGCCTGTACATTACCTTCAGGGCATTTTCACTGATATCTTTGCggg                                 >  1:4836202/1‑68 (MQ=255)

                                ggCTTCGTATCCCGCTTTATTGAGCCTGTACATTACCTTCAGGGCATTTTCACTGATATCTTTGCGGGaaa                              >  1:3132561/1‑71 (MQ=255)

                                ggCTTCGTATCCCGCTTTATTGAGCCTGTACATTACCTTCAGGGCATTTTCACTGATATCTTTGCGGGaaa                              >  1:4407751/1‑71 (MQ=255)

                                    tCGTATCCCGCTTTATTGAGCCTGTACATTACCTTCAGGGCATTTTCACTGat                                            <  1:4295635/53‑1 (MQ=255)

                                    tCGTATCCCGCTTTATTGAGCCTGTACATTACCTTCAGGGCATTTTCACTGat                                            <  1:3028077/53‑1 (MQ=255)

                                    tCGTATCCCGCTTTATTGAGCCTGTACATTACCTTCAGGGCATTTTCACTGat                                            <  1:4234120/53‑1 (MQ=255)

                                    tCGTATCCCGCTTTATTGAGCCTGTACATTACCTTCAGGGCATTTTCACTGat                                            <  1:4219833/53‑1 (MQ=255)

                                    tCGTATCCCGCTTTATTGAGCCTGTACATTACCTTCAGGGCATTTTCACTGat                                            <  1:4857418/53‑1 (MQ=255)

                                    tCGTATCCCGCTTTATTGAGCCTGTACATTACCTTCAGGGCATTTTCACTGat                                            <  1:2029108/53‑1 (MQ=255)

                                    tCGTATCCCGCTTTATTGAGCCTGTACATTACCTTCAGGGCATTTTCACTGat                                            <  1:2275325/53‑1 (MQ=255)

                                    tCGTATCCCGCTTTATTGAGCCTGTACATTACCTTCAGGGCATTTTCACTGat                                            <  1:4464904/53‑1 (MQ=255)

                                     cGTATCCCGCTTTATTGAGCCTGTACATTACCTTCAGGGCATTTTCACTGATATCTTTGCGGGAAATAg                           <  1:899346/69‑1 (MQ=255)

                                     cGTATCCCGCTTTATTGAGCCTGTACATTACCTTCAGGGCATTTTCACTGATATCTTTGCGGGAAATAg                           <  1:1440713/69‑1 (MQ=255)

                                     cGTATCCCGCTTTATTGAGCCTGTACATTACCTTCAGGGCATTTTCACTGATATCTTTGCGGGAAATAg                           <  1:2354788/69‑1 (MQ=255)

                                     cGTATCCCGCTTTATTGAGCCTGTACATTACCTTCAGGGCATTTTCACTGATATCTTTGCGGGAAATAGCa                         >  1:1248868/1‑71 (MQ=255)

                                      gTATCCCGCTTTATTGAGCCTGTACATTACCTTCAGGGCATTTTCACTGatat                                          >  1:2509956/1‑53 (MQ=255)

                                      gTATCCCGCTTTATTGAGCCTGTACATTACCTTCAGGGCATTTTCACTGatat                                          >  1:3483688/1‑53 (MQ=255)

                                      gTATCCCGCTTTATTGAGCCTGTACATTACCTTCAGGGCATTTTCACTGatat                                          >  1:4096730/1‑53 (MQ=255)

                                      gTATCCCGCTTTATTGAGCCTGTACATTACCTTCAGGGCATTTTCACTGatat                                          >  1:1873598/1‑53 (MQ=255)

                                      gTATCCCGCTTTATTGAGCCTGTACATTACCTTCAGGGCATTTTCACTGatat                                          >  1:3176729/1‑53 (MQ=255)

                                      gTATCCCGCTTTATTGAGCCTGTACATTACCTTCAGGGCATTTTCACTGATATCTTTGCGGGAAATAGCa                         <  1:5064678/70‑1 (MQ=255)

                                         tCCCGCTTTATTGAGCCTGTACATTACCTTCAGGGCa                                                       >  1:3043010/1‑37 (MQ=255)

                                         tCCCGCTTTATTGAGCCTGTACATTACCTTCAGGGCAt                                                      >  1:1547585/1‑38 (MQ=255)

                                         tCCCGCTTTATTGAGCCTGTACATTACCTTCAGGGCAt                                                      >  1:2499805/1‑38 (MQ=255)

                                         tCCCGCTTTATTGAGCCTGTACATTACCTTCAGGGCAt                                                      >  1:829574/1‑38 (MQ=255)

                                         tCCCGCTTTATTGAGCCTGTACATTACCTTCAGGGCAt                                                      >  1:2830376/1‑38 (MQ=255)

                                         tCCCGCTTTATTGAGCCTGTACATTACCTTCAGGGCATTTTCACTGATATCTTTGCGGGAAATAg                           <  1:3528276/65‑1 (MQ=255)

                                         tCCCGCTTTATTGAGCCTGTACATTACCTTCAGGGCATTTTCACTGATATCTTTGCGGGAAATAg                           <  1:191599/65‑1 (MQ=255)

                                         tCCCGCTTTATTGAGCCTGTACATTACCTTCAGGGCATTTTCACTGATATCTTTGCGGGAAATAg                           <  1:4648923/65‑1 (MQ=255)

                                         tCCCGCTTTATTGAGCCTGTACATTACCTTCAGGGCATTTTCACTGATATCTTTGCGGGAAATAg                           <  1:3877081/65‑1 (MQ=255)

                                         tCCCGCTTTATTGAGCCTGTACATTACCTTCAGGGCATTTTCACTGATATCTTTGCGGGAAATAg                           <  1:3819931/65‑1 (MQ=255)

                                         acccGCTTTATTGAGCCTGTACATTACCTTCAGGGCAt                                                      >  1:3528345/2‑38 (MQ=255)

                                          cccGCTTTATTGAGCCTGTACATTACCTTCAGGGCATTTTCACTGATATCTTTGCGGGAAATAg                           <  1:1720391/64‑1 (MQ=255)

                                            cGCTTTATTGAGCCTGTACATTACCTTCAGGGCAttt                                                    >  1:1879902/1‑37 (MQ=255)

                                            cGCTTTATTGAGCCTGTACATTACCTTCAGGGCATTTTCACTGATATCTTTGCGGGAAATAGCATgctgct                  >  1:241917/1‑71 (MQ=255)

                                             gCTTTATTGAGCCTGTACATTACCTTCAGGGCATTTTCACTGATATCTTTGCGGGAAATAGCATgctgct                  >  1:75116/1‑70 (MQ=255)

                                             gCTTTATTGAGCCTGTACATTACCTTCAGGGCATTTTCACTGATATCTTTGCGGGAAATAGCATgctgct                  >  1:3899001/1‑70 (MQ=255)

                                             gCTTTATTGAGCCTGTACATTACCTTCAGGGCATTTTCACTGATATCTTTGCGGGAAATAGCATgctgct                  >  1:4337304/1‑70 (MQ=255)

                                             gCTTTATTGAGCCTGTACATTACCTTCAGGGCATTTTCACTGATATCTTTGCGGGAAATAGCATgctgct                  >  1:3736123/1‑70 (MQ=255)

                                              cTTTATTGAGCCTGTACATTACCTTCAGGGCATTTTCACTGATATCTTTGCGGGAAATAGCAtgc                      >  1:2306386/1‑65 (MQ=255)

                                                   ttGAGCCTGTACATTACCTTCAGGGCATTTTCACTGATATCTTTGCGGGAAATAGCATGCTGCTCAcgcg            <  1:5045726/70‑1 (MQ=255)

                                                     gAGCCTGTACATTACCTTCAGGGCATTTTCACTGATATCTTTGCGGGAAATAGCATGCTGCTCACGCggg          <  1:2522811/70‑1 (MQ=255)

                                                      aGCCTGTACATTACCTTCAGGGCATTTTCACTGATATCTTTGCGGGAAATAGCATGCTGCTCACGCGGGa         <  1:2734531/70‑1 (MQ=255)

                                                       gCCTGTACATTACCTTCAGGGCATTTTCACTGATATCTTTGCGGGAAATAGCATGCTGCTCACGCGGGATc       <  1:4954301/71‑1 (MQ=255)

                                                        ccTGTACATTACCTTCAGGGCATTTTCTCTGATATCTTTGCGGGAAATAGCATGCTGCTCACGCGGGATc       <  1:4728356/70‑1 (MQ=255)

                                                        ccTGTACATTACCTTCAGGGCATTTTCACTGATATCTTTGCGGGAAATAGCATGCTGCTCACGCGGGATCa      <  1:2710203/71‑1 (MQ=255)

                                                        ccTGTACATTACCTTCAGGGCATTCTCACTGATATCTTTGCGGGAAATAGCATGCTGCTCACGCGGGATCa      <  1:1721050/71‑1 (MQ=255)

                                                           gTACATTACCTTCAGGGCATTTTCACTGATATCTTTGCGGGAAATAGCATGCTGCTCACGCGGGATc       <  1:4830224/67‑1 (MQ=255)

                                                           gTACATTACCTTCAGGGCATTTTCACTGATATCTTTGCGGGAAATAGCATGCTGCTCACGCGGGATc       <  1:2981370/67‑1 (MQ=255)

                                                           gTACATTACCTTCAGGGCATTTTCACTGATATCTTTGCGGGAAATAGCATGCTGCTCACGCGGGATc       <  1:896965/67‑1 (MQ=255)

                                                            tACATTACCTTCAGGGCATTTTCACTGATATCTTTGCGGGAAATAg                           <  1:639716/46‑1 (MQ=255)

                                                            tACATTACCTTCAGGGCATTTTCACTGATATCTTTGCGGGAAATAGCATGCTGCTCACGCGGGATCACCGt  <  1:4058835/71‑1 (MQ=255)

                                                            tACATTACCTTCAGGGCATTTTCACTGATATCTTTGCGGGAAATAGCATGCTGCTCACGCGGGATCACCGt  <  1:61438/71‑1 (MQ=255)

AACAGGTCGCGCACGCCGCCGCCAACCAGCCAGGCTTCGTATCCCGCTTTATTGAGCCTGTACATTACCTTCAGGGCATTTTCACTGATATCTTTGCGGGAAATAGCATGCTGCTCACGCGGGATCACCGT  >  Exported/143649‑143779

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 14 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`