BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP009273	155,132	C→G	R161P (CGC→CCC)	pcnB ←	poly(A) polymerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	155,132	0	C	G	95.9%	154.0 / ‑1.2	49	R161P (CGC→CCC)	pcnB	poly(A) polymerase
Reads supporting (aligned to +/- strand): ref base C (1/1); new base G (30/17); total (31/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.46e-01

CGCCAACGTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCC  >  CP009273/155067‑155196

cgccAACGTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGGCTGTTGATAGTGAAATCgggg                                                                >  1:84619/1‑68 (MQ=255)

 gccAACGTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGGCTGTTGATAGTGAAATCggggc                                                               <  1:1168370/68‑1 (MQ=255)

  ccAACGTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGGCTGTTGATAGTGAAATCggggcg                                                              <  1:3393275/68‑1 (MQ=255)

      cGTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGGCTGTTGATAGTGAAATCGGGGCGCTgg                                                          >  1:2350269/1‑68 (MQ=255)

       gTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGGCTGTTGATAGTGAAATCGGGGCGCgggg                                                         <  1:4252325/68‑1 (MQ=255)

          aTCACGGACGGTAAAATCCGCTACGCTGTAATACAGGCTGTTGATAGTGAAATCGGGGCGCTGGGCGt                                                      <  1:1838927/68‑1 (MQ=255)

          aTCACGGACGGTAAAATCCGCTACGCTGTAATACAGGCTGTTGATAGTGAAATCGGGGCGCTGGGCGt                                                      >  1:719417/1‑68 (MQ=255)

           tCACGGACGGTAAAATCCGCTACGCTGTAATACAGGCTGTTGATAGTGAAATCGGGGCGCTGGGCGtc                                                     >  1:2217725/1‑68 (MQ=255)

            cACGGACGGTAAAATCCGCTACGCTGTAATACAGGCTGTTGATAGTGAAATCGGGGCGCTGGGCGtct                                                    >  1:1105900/1‑68 (MQ=255)

            cACGGACGGTAAAATCCGCTACGCTGTAATACAGGCTGTTGATAGTGAAATCGGGGCGCTGGGCGtct                                                    >  1:2776475/1‑68 (MQ=255)

            cACGGACGGTAAAATCCGCTACGCTGTAATACAGGCTGTTGATAGTGAAATCGGGGCGCTGGGCGtct                                                    >  1:1946138/1‑68 (MQ=255)

            cACGGACGGTAAAATCCGCTACGCTGTAATACAGGCTGTTGATAGTGAAATCGGGGCGCTGGGCGtct                                                    >  1:4280042/1‑68 (MQ=255)

             acggacggTAAAATCCGCTACGCTGTAATACAGGCTGTTGATAGTGAAATCGGGGCGCTGGGCGtctt                                                   <  1:3238746/68‑1 (MQ=255)

             acggacggTAAAATCCGCTACGCTGTAATACAGGCTGTTGATAGTGAAATCGGGGCGCTGGGCGtctt                                                   >  1:536222/1‑68 (MQ=255)

             acggacggTAAAATCCGCTACGCTGTAATACAGGCTGTTGATAGTGAAATCGGGGCGCTGGGCGtctt                                                   >  1:1769095/1‑68 (MQ=255)

               ggacggTAAACTCCGCTACGCTGTAATACAGGCTGTTGATAGTGAAATCGGGGCGCTGGGCGtcttct                                                 <  1:819032/68‑1 (MQ=255)

                   ggTAAAATCCGCTACGCTGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGa                                             <  1:29569/68‑1 (MQ=255)

                       aaaTCCGCTACGCTGTAATACAGGCTGTTGATAGTGAAATCGGGGCGCTGGGCGTCTTCTTCGATGGa                                         <  1:3721659/68‑1 (MQ=255)

                          tCCGCTACGCTGTAATACAGGCTGTTGATAGTGAAATCGGGGCGCTGGGCGTCTTCTTCGATGGAGcc                                      <  1:3984869/68‑1 (MQ=255)

                          tCCGCTACGCTGTAATACAGGCTGTTGATAGTGAAATCGGGGCGCTGGGCGTCTTCTTCGATGGAGcc                                      <  1:3915764/68‑1 (MQ=255)

                          tCCGCTACGCTGTAATACAGGCTGTTGATAGTGAAATCGGGGCGCTGGGCGTCTTCTTCGATGGAGcc                                      <  1:3670299/68‑1 (MQ=255)

                           ccGCTACGCTGTAATACAGGCTGTTGATAGTGAAATCGGGGCGCTGGGCGTCTTCTTCGATGGAGCCg                                     >  1:807088/1‑68 (MQ=255)

                           ccGCTACGCTGTAATACAGGCTGTTGATAGTGAAATCGGGGCGCTGGGCGTCTTCTTCGATGGAGCCg                                     >  1:967448/1‑68 (MQ=255)

                           ccGCTACGCTGTAATACAGGCTGTTGATAGTAAAATCGGGGCGCTGGGCGTCTTCTTCGATGGAGCCg                                     >  1:2379935/1‑68 (MQ=255)

                                aCGCTGTAATACAGGCTGTTGATAGTGAAATCGGGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAAt                                <  1:3219935/68‑1 (MQ=255)

                                aCGCTGTAATACAGGCTGTTGATAGTGAAATCGGGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAAt                                <  1:586574/68‑1 (MQ=255)

                                aCGCTGCAATACAGGCTGTTGATAGTGAAATCGGGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAAt                                >  1:2229791/1‑68 (MQ=255)

                                   cTGTAATACAGGCTGTTGATAGTGAAATCGGGGCGCTGGGCGTCTTCTTCGATGGTGCCGAAAAtgtt                             <  1:2395235/68‑1 (MQ=255)

                                       aaTACAGGCTGTTGATAGTGAAATCGGGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTcg                         >  1:1515233/1‑68 (MQ=255)

                                        aTACAGGCTGTTGATAGTGAAATCGGGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTcgc                        <  1:4100415/68‑1 (MQ=255)

                                        aTACAGGCTGTTGATAGTGAAATCGGGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTcgc                        <  1:3123678/68‑1 (MQ=255)

                                        aTACAGGCTGTTGATAGTGAAATCGGGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTcgc                        <  1:2223698/68‑1 (MQ=255)

                                        aTACAGGCTGTTGATAGTGAAATCGGGGCGCTGGGCGGCTTCTTCGATGGAGCCGAAAATGTTGTcgc                        <  1:2202496/68‑1 (MQ=255)

                                          aCAGGCTGTTGATAGTGAAATCGGGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCgcgc                      >  1:2207790/1‑68 (MQ=255)

                                          aCAGGCTGTTGATAGTGAAATCGGGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCgcgc                      >  1:3782320/1‑68 (MQ=255)

                                          aCAGGCTGTTGATAGTGAAATCGGGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCgcgc                      >  1:3904095/1‑68 (MQ=255)

                                          aCAGGCTGTTGATAGTGAAATCGGGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCgcgc                      >  1:4005504/1‑68 (MQ=255)

                                          aCAGGCTGTTGATAGTGAAATCGGGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCgcgc                      >  1:3352527/1‑68 (MQ=255)

                                            aGGCTGTTGATAGTGCTATCGGGGCGCTGGGCGTCTTCTTCTATGGAGACGAAAATGTTCTCgcgccg                    >  1:4687276/1‑66 (MQ=255)

                                            aGGCTGTTGATAGTGAAATCGGGGCGCTGGTCGTCTTCTTCGATGGAGCCGAAAATGTTGTCTCgcag                    >  1:603171/1‑68 (MQ=255)

                                            aGGCTGTTGATAGTGAAATCGGGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCgcag                    >  1:4232486/1‑68 (MQ=255)

                                            aGGCTGTTGATAGTGAAATCGGGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCgcag                    >  1:2280457/1‑68 (MQ=255)

                                            aGGCTGTTGATAGTGAAATCGGGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCgcag                    >  1:3068924/1‑68 (MQ=255)

                                             ggCTGTTGATAGTGAAATCGGGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCgcagc                   <  1:2885585/68‑1 (MQ=255)

                                              gCTGTTGATAGTGAAATCGGGGCGCTGGGCGTCTTCTTCGATGGAGCCGAATATGTTGTCGCgcagca                  >  1:732342/1‑68 (MQ=255)

                                              gCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCgcagca                  >  1:4368423/1‑68 (MQ=255)

                                                  ttGTTAGTGAAATCGGGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACAt              <  1:309104/68‑1 (MQ=255)

                                                       aGTGAAATCGGGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGt         >  1:2751374/1‑68 (MQ=255)

                                                           aaaTCGGGGCGCTGGGCGTCTTCTTCGATGGAGCCGAATGTGGTGTCGCGCAGCAACAGGCGGTTTTg     >  1:837546/1‑68 (MQ=255)

                                                             aTCGGGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGcc   >  1:4206156/1‑68 (MQ=255)

                                                             aTCGGGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGcc   >  1:1276119/1‑68 (MQ=255)

                                                              tCGGGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGccc  >  1:1625333/1‑68 (MQ=255)

CGCCAACGTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCAACATGCCGTTTTGCCC  >  CP009273/155067‑155196

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: `‑`