BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,314,380	C→T	29.2%	noncoding (272/1522 nt)	PP_16SG ←	16S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,314,380	0	C	T	29.2%	‑1.5 / 14.0	24	noncoding (272/1522 nt)	PP_16SG	16S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base C (5/12); new base T (4/3); total (9/15)
Fisher's exact test for biased strand distribution p-value = 3.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTGCTTTACAATCCGAAGACCTTCTTCACACACGCGGCATGGCTGGATCAGGCTTTCGCCCATTGTCCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGACTGATCATCCTCTCAGACCAGTTACGGATCGTCGCCTTGGTGAGCCATTACCCCACCAACTAGCTAATCCGACCTAGGCTCATCTGATAGCGCAAGGCCCGAAGGTCCCCTGCTTTCTCCCGTAGGACGTATGCGGTATTAGCGTTCCTTTC  >  NC_002947/5314236‑5314514

agagCTTTACAACCCTAGGGCCTTCATCACTCACGCGGCATGGCTGGATCAGGCTTGCGCCCATTGTCCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCTGATCATCCTCTCAGACCAgct                                                                                                                                   >  1:1607834/4‑148 (MQ=1)

                aaaaCCTTCATCACTCACGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAAGATTCCCTACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCCGATCACCCTCTCAGGTCGGCTACGCATCGTTGCCTTg                                                                                                                   <  2:3517362/147‑1 (MQ=0)

                aaaaCCTTCATCACTCACGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAAGATTCCCTACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCCGATCACCCTCTCAGGTCGGCTACGCATCGTTGCCTTg                                                                                                                   >  2:2503381/4‑150 (MQ=0)

                aaaaCCTTCATCACTCACGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAAGATTCCCTACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCCGATCACCCTCTCAGGTCGGCTACGCATCGTTGCCTTg                                                                                                                   >  2:2503354/4‑150 (MQ=0)

                    ccgtcttcACACCCGCGGCATGGCTGGATCAGGCTTTCGCCCATTGTCCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGACTGATCATCCTCTCAGACCAGTTACGGATCGTCGCCTTGGTGa                                                                                                               <  1:221046/147‑1 (MQ=1)

                       tcatcactcacGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAAGATTCCCTACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCCGATCACCCTCTCAGGTCGGCTACGCATCGTTGCCTTGGTGAGcc                                                                                                            >  1:353759/4‑150 (MQ=0)

                       tcatcactcacGCGGCATGGCTGGATCAGGCTTGCGCCCATTGTCCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCTGATCATCCTCTCAGACCAGCTATGGATCGTCGCCTTGGTGAGcc                                                                                                            >  1:976611/4‑150 (MQ=1)

                       tcatcactcacGCGGCATGGCTGGATCAGGCTTGCGCCCATTGTCCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCTGATCATCCTCTCAGACCAGCTATGGATCGTCGCCTTGGTGAGcc                                                                                                            <  1:577680/147‑1 (MQ=1)

                       tcatcactcacGCGGCATGGCTGGATCAGGCTTGCGCCCATTGTCCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCTGATCATCCTCTCAGACCAGCTATGGATCGTCGCCTTGGTGAGcc                                                                                                            <  1:577660/147‑1 (MQ=1)

                       tcatcactcacGCGGCATGGCTGGATCAGGCTTGCGCCCATTGTCCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCTGATCATCCTCTCAGACCAGCTATGGATCGTCGCCTTGGTGAGcc                                                                                                            <  1:1992990/147‑1 (MQ=1)

                       tcatcactcacGCGGCATGGCTGGATCAGGCTTGCGCCCATTGTCCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCTGATCAGCCTCTCAGACCAGCTATGGATCGTCGCCTTGGTGAGcc                                                                                                            >  2:4382511/4‑150 (MQ=1)

                            actcacGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAAGATTCCCTACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCCGATCACCCTCTCAGGTCGGCTACGCATCGTTGCCTTGGTGAGCCGTTAc                                                                                                       <  1:57466/147‑1 (MQ=0)

                            actcacGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAAGATTCCCTACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCCGATCACCCTCTCAGGTCGGCTACGCATCGTTGCCTTGGTGAGCCGTTAc                                                                                                       >  1:4615606/4‑150 (MQ=0)

                            actcacGCGGCGTTGCTCCGTCAGACTTTCGTCCATTGCGGAAGATTCCCTACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCCGATCACCCTATCAGGTCGGCTACGCATCGTTGCCTTGGTGAGCCGTTAc                                                                                                       <  2:1540004/147‑1 (MQ=0)

                            actcacGCGGCATGGCTGGATCAGGCTTGCGCCCATTGTCCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCTGATCATCCTCTCAGACCAGCTATGGATCGTCGCCTTGGTGAGCCTTTAc                                                                                                       >  1:5067203/4‑150 (MQ=1)

                            actcacGCGGCATGGCTGGATCAGGCTTGCGCCCATTGTCCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCTGATCATCCTCTCAGACCAGCTATGGATCGTCGCCTTGGTGAGCCTTTAc                                                                                                       <  1:1055266/147‑1 (MQ=1)

                            actcacGCGGCATGGCTGGATCAGGCTTGCGCCCATTGTCCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCTGATCATCCTCTCAGACCAGCTATGGATCGTCGCCTTGGTGAGCCTTTAc                                                                                                       >  2:2412892/4‑150 (MQ=1)

                            actcacGCGGCATGGCTGGATCAGGCTTGCGCCCATTGTCCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCTGATCATCCTCTCAGACCAGCTATGGATCGTCGCCTTGGTGAGCCTTTAc                                                                                                       <  1:4296356/147‑1 (MQ=1)

                            actcacGCGGCATGGCTGGATCAGGCTTGCGCCCATTGTCCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCTGATCATCCTCTCAGACCAGCTATGGATCGTCGCCTTGGTGAGCCTTTAc                                                                                                       <  1:1786854/147‑1 (MQ=1)

                            actcacGCGGCATGGCTGGATCAGGCTTGCGCCCATTGTCCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCTGATCATCCTCTCAGACCAGCTATGGATCGTCGCCTTGGTGAGCCTTTAc                                                                                                       <  2:398785/147‑1 (MQ=1)

                            actcacGCGGCATGGCTGGATCAGGCTTGCGCCCATTGTCCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCTGATCATCCTCTCAGACCAGCTATGGATCGTCGCCTTGGTGAGCCTTTAc                                                                                                       <  2:4304367/147‑1 (MQ=1)

                            actcacGCGGCATGGCTGGATCAGGCTTGCGCCCATTGTCCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAGTGTGGCTGATCATCCTCTCAGACCAGCTATGGATCGTCGCCTTGGTGAGCCTTTAc                                                                                                       <  2:4578191/147‑1 (MQ=1)

                                                                                           cggaGGAGTCTGGACCGTGTCTCAGGTCCAGTGTGACTGATCATCCTCTCAGACCAGTTACGGATCGTCGCCTTGGTGAGCCATTACCCCACCAACTAGCTAATCCGACCTAGGCTCATCTGATAGCGCAAGGCCCGAAGGTCCCCTGCt                                        <  2:4482366/147‑1 (MQ=1)

                                                                                                                                 gaacatcCTCTCAGACCAGTGACGGATCGTCGCCTTGGTGAGCCATTACCCCACCAACTAGCTAATCCGACCTAGGCTCATCTGATAGCGCAAGGCCCGAAGGTCCCCTGCTTTCTCCCGTAGGAAGTATGCGGTATTAGCGTTCCTTTc  <  2:4512170/147‑1 (MQ=0)

AGTGCTTTACAATCCGAAGACCTTCTTCACACACGCGGCATGGCTGGATCAGGCTTTCGCCCATTGTCCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGACTGATCATCCTCTCAGACCAGTTACGGATCGTCGCCTTGGTGAGCCATTACCCCACCAACTAGCTAATCCGACCTAGGCTCATCTGATAGCGCAAGGCCCGAAGGTCCCCTGCTTTCTCCCGTAGGACGTATGCGGTATTAGCGTTCCTTTC  >  NC_002947/5314236‑5314514

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