BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,899,164:1	+C	11.9%	coding (1389/1413 nt)	PP_2552 →	DOPA decarboxylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,899,164	1	.	C	11.9%	143.8 / 14.6	42	coding (1389/1413 nt)	PP_2552	DOPA decarboxylase
Reads supporting (aligned to +/- strand): ref base . (19/18); new base C (2/3); total (21/21)
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.93e-01

AAGGGCTGGGCCGAGCGGCTGAATGCATCCGGCGCTGCTTATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGC  >  NC_002947/2899027‑2899306

aaGGGCTGGGCCGAGCGGCTGAATGCATCCGGCGCTGCTTATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTg                                                                                                                                             >  2:89674/1‑141 (MQ=255)

aaGGGCTGGGCCGAGCGGCTGAATGCATCCGGCGCTGCTTATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTg                                                                                                                                             <  1:89674/141‑1 (MQ=255)

aaGGGCTGGGCCGAGCGGCTGAATGCATCCGGCGCTGCTTATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGt                                                                                                                                     >  7:186208/1‑149 (MQ=255)

  gggCTGGGCCGAGCGGCTGAATGCATCCGGCGCTGCTTATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGTCCTGCA‑GACGTgc                                                                                                                                   <  4:112760/149‑2 (MQ=255)

        ggCCGAGCGGCTGAATGCATCCGGCGCTGCTTATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAgg                                                                                                                             <  8:221584/149‑1 (MQ=255)

          ccGAGCGGCTGAATGCATCCGGCGCTGCTTATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGc                                                                                                                           <  7:42873/149‑1 (MQ=255)

                  cTGAATGCATCCGGCGCTGCTTATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCt                                                                                                                   <  8:186208/149‑1 (MQ=255)

                     aaTGCATCCGGCGCTGCTTATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGCGGGCACGT‑CTGCAGGACGTGATCAAgg                                                                                                                             <  1:198431/136‑1 (MQ=255)

                     aaTGCATCCGGCGCTGCTTATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGCGGGCACGT‑CTGCAGGACGTGATCAAgg                                                                                                                             >  2:198431/1‑136 (MQ=255)

                            ccGGCGCTGCTTATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTAcc                                                                                                         >  1:335607/1‑149 (MQ=255)

                               gcgcTGCTTATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGt                                                                                                      <  2:335607/149‑1 (MQ=255)

                                 gctgctTATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTAt                                                                                                    <  7:100/149‑1 (MQ=255)

                                      ttATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCg                                                                                               >  4:79275/1‑149 (MQ=255)

                                      ttATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCg                                                                                               >  2:181766/1‑149 (MQ=255)

                                        aTGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCgg                                                                                              <  1:181766/148‑1 (MQ=255)

                                        aTGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCgg                                                                                              <  3:79275/148‑1 (MQ=255)

                                                    cTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGccc                                                                                      >  1:148652/1‑144 (MQ=255)

                                                    cTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGccc                                                                                      <  2:148652/144‑1 (MQ=255)

                                                                       gATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGt                                                                                                      <  3:326250/109‑1 (MQ=255)

                                                                       gATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGt                                                                                                      >  4:326250/1‑109 (MQ=255)

                                                                                           gTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCa                                          >  6:300023/1‑149 (MQ=255)

                                                                                           gTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCa                                          >  6:129602/1‑149 (MQ=255)

                                                                                           gTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCa                                          <  5:300023/149‑1 (MQ=255)

                                                                                                        cgaGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGaa                                 <  7:35476/145‑1 (MQ=255)

                                                                                                        cgaGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGaa                                 >  8:35476/1‑145 (MQ=255)

                                                                                                         gaGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGcccc                                                              >  4:27154/1‑115 (MQ=255)

                                                                                                                   aTGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACggc                  >  8:133141/1‑149 (MQ=255)

                                                                                                                      tgCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGccc                                                               >  4:253123/1‑101 (MQ=255)

                                                                                                                      tgCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGccc                                                               <  3:253123/101‑1 (MQ=255)

                                                                                                                       gCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACggctggc              >  8:286746/1‑149 (MQ=255)

                                                                                                                       gCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACggctggc              >  2:192729/1‑149 (MQ=255)

                                                                                                                         aGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCtt            <  5:129602/149‑1 (MQ=255)

                                                                                                                            ggCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCAt         <  7:133141/149‑1 (MQ=255)

                                                                                                                                gtgGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAg                               <  1:53850/123‑1 (MQ=255)

                                                                                                                                gtgGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAg                               >  2:53850/1‑123 (MQ=255)

                                                                                                                                   ggCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCATCATAGGCTGc  >  6:67083/1‑149 (MQ=255)

                                                                                                                                    tcACGTCCTGCA‑GACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTAc                           <  6:116643/122‑1 (MQ=255)

                                                                                                                                    tcACGTCCTGCA‑GACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTAc                           >  5:116643/2‑123 (MQ=255)

                                                                                                                                    tcACGTCCTGCA‑GACGTGATCAAGGGCTGAAGCCTGTACGGGACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTc           >  2:17471/2‑139 (MQ=39)

                                                                                                                                    tcACGTCCTGCA‑GACGTGATCAAGGGCTGAAGCCTGTACGGGACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTc           <  1:17471/138‑1 (MQ=39)

                                                                                                                                        gT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACgg                   >  8:180715/1‑127 (MQ=255)

                                                                                                                                        gT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACgg                   <  7:180715/127‑1 (MQ=255)

AAGGGCTGGGCCGAGCGGCTGAATGCATCCGGCGCTGCTTATGTGACGCCGGCTACACTGGACGGGCGGTGGATGGTGCGGGTTTCGATTGGTGCGCTGCCGACCGAGCGGGGGGATGTGCAGCGGCTGTGGGCACGT‑CTGCAGGACGTGATCAAGGGCTGAAGCCTGTACGGTACCTGTATGAGCGGGCTTGCCCGCGAAGGGGCTGCAAAGCAGCCCCTGCAATCTATCAGTGGCTCACCGCCTGAAAGCTACCTTTACGGCTGGCTTCATAGGCTGC  >  NC_002947/2899027‑2899306

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