BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,598,461	T→C	17.3%	T1106A (ACC→GCC)	PP_4071 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,598,461	0	T	C	17.3%	47.2 / 4.0	23	T1106A (ACC→GCC)	PP_4071	membrane protein
Reads supporting (aligned to +/- strand): ref base T (11/8); new base C (2/2); total (13/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCATTGAGCGCACCGCGGCTCAGCAGGCGAAACATCGACCAAGGCCCTCGATACTCCAGGCTGCTGCTATTGCCGTTCTGCCTGAGCAGGGTCAGGTTGCTGCGCACGTGTTGCCCCAAGGTATTGGGCCAGACAATACCGGTGATCTGGCTTGGGCCGTGGGTGTAGGAAATCAGCTGGCCATCCAGATCCAGCAGGCTGGTGCGTTGA‑TTTGCACTAAGCCCCAACGGCTCGATACTGAACTGCACGCTGAGGT  >  NC_002947/4598340‑4598595

ccATTGAGCGCACCGCGGCTCAGCAGGCGAAACATCGACCAAGGCCCTCGATACTCCAGGCTGCTGCTATTGCCGTTCTGCCTGAGCAGGGTCAGGTTGCTGCGCACGTGTTGCCCCAAGGTATTGGGCCAGACaa                                                                                                                           >  1:50024/1‑136 (MQ=255)

 cATTGAGCGCACCGCGGCTCAGCAGGCGAAACATCGACCAAGGCCCTCGATACTCCAGGCTGCTGCTATTGCCGTTCTGCCTGAGCAGGGTCAGGTTGCTGCGCACGTGTTGCCCCAAGGCATTGGGCCAGACAAt                                                                                                                          >  2:189838/1‑136 (MQ=255)

 cATTGAGCGCACCGCGGCTCAGCAGGCGAAACATCGACCAAGGCCCTCGATACTCCAGGCTGCTGCTATTGCCGTTCTGCCTGAGCAGGGTCAGGTTGCTGCGCACGTGTTGCCCCAAGGCATTGGGCCAGACAAt                                                                                                                          >  2:131254/1‑136 (MQ=255)

      aGCGCACCGCGGCTCAGCAGGCGAAACATCGACCAAGGCCCTCGATACTCCAGGCTGCTGCTATTGCCGTTCTGCCTGAGCAGGGTCAGGTTGCTGCGCACGTGTTGCCCCAAGGTATTg                                                                                                                                     <  2:259075/120‑1 (MQ=255)

      aGCGCACCGCGGCTCAGCAGGCGAAACATCGACCAAGGCCCTCGATACTCCAGGCTGCTGCTATTGCCGTTCTGCCTGAGCAGGGTCAGGTTGCTGCGCACGTGTTGCCCCAAGGTATTg                                                                                                                                     >  1:259075/1‑120 (MQ=255)

                          gCGAAACATCGACCAAGGCCCTCGATACTCCAGGCTGCTGCTATTGCCGTTCTGCCTGAGCAGGGTCAGGTTGCTGCGCACGTGTTGCCCCAAGGTATTGGGCCAGACAATACCGGTGATCTGGCTTGGGCCGTgg                                                                                                 >  2:270976/1‑136 (MQ=255)

                               aCATCGACCAAGGCCCTCGATACTCCAGGCTGCTGCTATTGCCGTTCTGCCTGAGCAGGGTCAGGTTGCTGCGCACGTGTTGCCCCAAGGTATTGGGCCAGACAATACCGGTGATCTGGCTTGGGCCGTGGgtgtt                                                                                            >  2:135770/1‑135 (MQ=255)

                                       cAAGGCCCTCGATACTCCAGGCTGCTGCTATTGCCGTTCTGCCTGAGCAGGGTCAGGTTGCTGCGCACGTGTTGCCCCAAGGTATTGGGCCAGACAATACCGGTGATCTGGCTTGGGCCGTGGGTGTAGGAAATCa                                                                                    >  2:149347/1‑136 (MQ=255)

                                                         aGGCTGCTGCTATTGCCGTTCTGCCTGAGCAGGGTCAGGTTGCTGCGCACGTGTTGCCCCAAGGTATTGGGCCAGACAATACCGGTGATCTGGCTTGGGCCGTGGGTGTAGGAAATCAGCTGGCCatccagatcca                                                                  >  2:56594/1‑136 (MQ=255)

                                                         aGGCTGCTGCAATTGCCCATCTGCCAGAGCACCGACAGGTTGCTGCGAACTTCTAGCCCCAAGGTATAGGGCAAGACAATACCGGTGATCTGGAATGCGCCGTGGGTGTATGAAATCAGCTGGCCatccagatcca                                                                  >  1:84639/1‑136 (MQ=255)

                                                                   tATTGCCGTTCTGCATGAGCAGGGTCAGGATGCTGCGCACCTGTTGCCCCAAGGTATTGGGCCAGACAATACCGCTGATCTGGCTTGGGCCGTGGGTGTAGCAAATCAGCTGGCCATCAAGATCAAGCAGGCTGGt                                                        >  1:146579/1‑136 (MQ=21)

                                                                        ccGTTCTGCCTGAGCAGGGTCAGGTTGCTGCGCACGTGTTGCCCCAAGGTATTGGGCCAGACAATACCGGTGATCTGGCTTGGGCCGTGGGTGTAGGAAATCAGCTGGCCATCCAGATCCAGCAGGCTGGTGCGtt                                                   <  1:66525/136‑1 (MQ=255)

                                                                                        gggTCAGGTTGCTGCGCACGTGTTGCCCCAAGGTATTGGGCCAGACAATACCGGTGATCTGGCTTGGGCCGTGGGTGTAGGAAATCAGCTGGCCATCCAGATCCAGCAGGCTGGTGCGTTGA‑TTTGCACTAAGccc                                  >  2:194300/1‑136 (MQ=255)

                                                                                             cggTTGCTGCGCACGTGTTGCCCCAAGGTATTGGGCCAGACAATACCGGTGATCTGGCTTGGGCCGTGGGTGTAGGAAATCAGCTGGCCATCCAGATCCAGCAGGCTGGTGCGTTGA‑TTTGCACTAAGCCCCAACg                             <  1:134940/135‑1 (MQ=255)

                                                                                             aGGTTGCTGCGCACGTGTTGCCCCAAGGTATTGGGCCAGACAATACCGGTGATCTGGGTTGGGCCGTGGGTGTAGGAAATCAGCTGGCCATCCAGATCCAGCAGGCTGGTGCGTTGA‑TTTGCACTAAGCCCCAACg                             <  1:245451/136‑1 (MQ=255)

                                                                                             aGGTTGCTGCGCACGTGTTGCCCCAAGGTATTGGGCCAGACAATACCGGTGATCTGGCTTGGGCCGTGGGTGTAGGAAATCAGCTGGCCATCCAGATCCAGCAGGCTGGTGCGTTGA‑TTTGCACTAAGCCCCAACg                             <  1:135770/136‑1 (MQ=255)

                                                                                             aGGTTGCTGCGCACGTGTTGCCCCAAGGTATTGGGCCAGACAATACCGGTGATCTGGCTTGGGCCGTGGGTGTAGGAAATCAGCTGGCCATCCAGATCCAGCAGGCTGGTGCATTGA‑TTTGCACTAAGCCCCAACg                             <  1:56594/136‑1 (MQ=255)

                                                                                                  gctgcGCACGTGTTGCCCCAAGGTTTTGGGCCAGACAATACCGGTGATCTGGCTTGGGCCGTGGGTGTAGGAAATCAGCTGGCCATCCAGATCCAGCAGGCTGGTGCGTTGA‑TTTGCACTAAGCCCCAACGGATCg                        <  1:39246/136‑1 (MQ=255)

                                                                                                  gctgcGCACGTGTTGCCCCAAGGTATTGGGCCAGACAATACCGGTGATCTGGCTTGGGCCGTGGGTGTAGGAAATCAGCTGGCCATCCAGATCCAGCAGGCTGGTGCGTTGA‑TTTGCACTAAGCCCCAACGGCTCg                        <  1:189405/136‑1 (MQ=255)

                                                                                                  gctgcGCACGTGTTGCCCCAAGGCATTGGGCCAGACAATACCGGTGATCTGGCTTGGGCCGTGGGTGTAGGAAATCAGCTGGCCATCCAGATCCAGCAGGCTGGTGCGTTGA‑TTTGCACTAAGCCCCAACGGCTCg                        <  1:131254/136‑1 (MQ=255)

                                                                                                  gctgcGCACGTGTTGCCCCAAGGCATTGGGCCAGACAATACCGGTGATCTGGCTTGGGCCGTGGGTGTAGGAAATCAGCTGGCCATCCAGATCCAGCAGGCTGGTGCGTTGA‑TTTGCACTAAGCCCCAACGGCTCg                        <  1:189838/136‑1 (MQ=255)

                                                                                                      cgcACGTGTTGCCCCAAGGTATTGGGCCAGACAATACCGGTGATCTGGCTTGGGCCGTGGGTGTAGGAAATCAGCTGGCCATCCAGATCCAGCAGGCTGGTGCGTTGATTTTGCACTACGCCCAACCGGcgggatg                     >  2:11138/1‑130 (MQ=21)

                                                                                                                        gTATTGGGCCAGACAATACCGGTGATCTGGCTTGGGCCGTGGGTGTAGGAAATCAGCTGGCCATCCAGATCCAGCAGGCTGGTGCGTTGA‑TTTGCACTAAGCCCCAACGGCTCGATACTGAACTGCACGCTGAGGt  >  2:189195/1‑136 (MQ=255)

CCATTGAGCGCACCGCGGCTCAGCAGGCGAAACATCGACCAAGGCCCTCGATACTCCAGGCTGCTGCTATTGCCGTTCTGCCTGAGCAGGGTCAGGTTGCTGCGCACGTGTTGCCCCAAGGTATTGGGCCAGACAATACCGGTGATCTGGCTTGGGCCGTGGGTGTAGGAAATCAGCTGGCCATCCAGATCCAGCAGGCTGGTGCGTTGA‑TTTGCACTAAGCCCCAACGGCTCGATACTGAACTGCACGCTGAGGT  >  NC_002947/4598340‑4598595

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