BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,313,394	A→T	46.6%	noncoding (1258/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,313,394	0	A	T	46.6%	‑5.1 / 37.8	30	noncoding (1258/1522 nt)	PP_16SG	16S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base A (5/11); new base T (9/5); total (14/16)
Fisher's exact test for biased strand distribution p-value = 1.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCACTCCCATGGTGTGACGGGCGGTGTGTACAAGGCCCGGGAACGTATTCACCGCGACATTCTGATTCGCGATTACTAGCGATTCCGACTTCACGCAGTCGAGTTGCAGACTGCGATCCGGACTACGATCGGTTTTGTGAGATTAGCTCCACCTCGCG‑GCTTGGCAACCCTCTGTACCGACCATTGTAGCACGTGTGTAGCCCAGGCCGTAAGGGCCATGATGACTTGACGTCATCCCCACCTTCCTCCGGTTTGTCACCGGCAGTCTC  >  NC_002947/5313249‑5313518

ccaaCTCCCATGGTGTGACGGGCGGTGTGTACAAGGCCCGGGAACGTATTCACCGCGGCATGCTGATCCGCGATTACTAGCGATTCCAACTTCATGCACTCGAGTTGCAGAGTGCAATCCGAACTGAGATGGCTTTTG‑GATATTAGCTca                                                                                                                          >  2:6004546/4‑149 (MQ=1)

ccaaCTCCCATGGTGTGACGGGCGGTGTGTACAAGGCCCGGGAACGTATTCACCGCGGCATGCTGATCCGCGATTACTAGCGATTCCAACTTCATGCACTCGAGTTGCAGAGTGCAATCCGAACTGAGATGGCTTTTG‑GAGATTAGCTca                                                                                                                          >  1:4844414/4‑149 (MQ=1)

ccaaCTCCCATGGTGTGACGGGCGGTGTGTACAAGGCCCGGGAACGTATTCACCGCGGCATGCTGATCCGCGATTACTAGCGATTCCAACTTCATGCACTCGAGTTGCAGAGTGCAATCCGAACTGAGATGGCTTTTG‑GAGATTAGCTca                                                                                                                          <  1:6120392/147‑2 (MQ=1)

ccaaCTCCCATGGTGTGACGGGCGGTGTGTACAAGGCCCGGGAACGTATTCACCGCGGCATGCTGATCCGCGATTACTAGCGATTCCAACTTCATGCACTCGAGTTGCAGAGTGCAATCCGAACTGAGATGGCTTTTG‑GAGATTAGCTca                                                                                                                          >  1:6471878/4‑149 (MQ=1)

             tgagACGGGCGGTGTGTACAAGGCCCGGGAACGTATTCACCGCGACATTCTGATTCGCGATTACTAGCGATTCCGACTTCACGCAGTCGAGTTGCAGACTGCGATCCGGACTACGATCGGTTTTGTGAGATTAGCTCCACCTCGCG‑GCtt                                                                                                             <  1:8477859/147‑1 (MQ=1)

                                                       cggcATGCTGATCCGCGATTACTAGCGATTCCAGCTTCATGTAGGCGAGTTGCAGCCTACAATCCGAACTGAGAATGGTTTTATGGGATTTGCTCGACCTCGCG‑GTTTTGCTGCCCTTTGTACCATCCATTGTAGCACGTGTGTAGCCCa                                                                   <  2:2040187/147‑1 (MQ=0)

                                                           atgcTGATCCGCGATTACTAGCGATTCCAGCTTCATGTAGGCGAGTTGCAGCCTACAATCCGAACTGAGAATGGTTTTATGGGATTTGCTCGACCTCGCG‑GTTTTGCTGCCCTTTGTACCATCCATTGTAGCACGTGTGTAGCCCAggtc                                                               <  1:1518647/147‑3 (MQ=0)

                                                           atgcTGATCCGCGATTACTAGCGATTCCAACTTCATGCACTCGAGTTGCAGAGTGCAATCCGAACTGAGATGGCTTTTG‑GAGATTAGCTCACACTCGCGTGCTT‑GCTGCCCACTGTCACCACCATTGTAGCACGTGTGTAGCCCAGCCCg                                                              >  1:8378354/4‑150 (MQ=0)

                                                                 atccgcgATTACTAGCGATTCCAACTTCATGCACTCGAGTTGCAGAGTGCAATCCGAACTGAGATGGCTTTTG‑GAGATTAGCTCACACTCGCGTGCTT‑GCTGCCCACTGTCACCACCATTGTAGCACGTGTGTAGCCCAGCCCGTAAggg                                                        <  1:1969153/147‑1 (MQ=0)

                                                                 atccgcgATTACTAGCGATTCCAACTTCATGCACTCGAGTTGCAGAGTGCAATCCGAACTGAGATGGCTTTTG‑GAGATTAGCTCACACTCGCGTGCTT‑GCTGCCCACTGTCACCACCATTGTAGCACGTGTGTAGCCCAGCCCGTAAggg                                                        <  1:5638953/147‑1 (MQ=0)

                                                                 atccgcgATTACTAGCGATTCCAACTTCATGCACTCGAGTTGCAGAGTGCAATCCGAACTGAGATGGCTTTTG‑GAGATTAGCTCACACTCGCGTGCTT‑GCTGCCCACTGTCACCACCATTGTAGCACGTGTGTAGCCCAGCCCGTAAggg                                                        <  1:7371031/147‑1 (MQ=0)

                                                                                     ccaaCTTCATGCACTCGAGTTGCAGAGTGCAATCCGAACTGAGATGGCTTTTG‑GAGATTAGCTCACACTCGCGTGCTT‑GCTGCCCACTGTCACCACCATTGTAGCACGTGTGTAGCCCAGCCCGTAAGGGCCATGAGGACTTGACGTCAt                                    <  1:8356669/147‑1 (MQ=0)

                                                                                     ccaaCTTCATGCACTCGAGTTGCAGAGTGCAATCCGAACTGAGATGGCTTTTG‑GAGATTAGCTCACACTCGCGTGCTT‑GCTGCCCACTGTCACCACCATTGTAGCACGTGTGTAGCCCAGCCCGTAAGGGCCATGAGGACTTGACGTCAt                                    <  2:771242/147‑1 (MQ=0)

                                                                                             atgtaggcGAGTTGCAGCCTACAATCCGAACTGAGAATGGTTTTATGGGATTTGCTCGACCTCGCG‑GTTTTGCTGCCCTTTGTACCATCCATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCATGATGATTTGACGTCATCCCCAcc                             <  1:21772/143‑1 (MQ=0)

                                                                                             atgtaggcGAGTTGCAGCCTACAATCCGAACTGAGAATGGTTTTATGGGATTTGCTCGACCTCGCG‑GTTTTGCTGCCCTTTGTACCATCCATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCATGATGATTTGACGTCATCCCCAcc                             <  1:2016106/143‑1 (MQ=0)

                                                                                              tgtaggcGAGTTGCAGCCTACAATCCGAACTGAGAATGGTTTTATGGGATTTGCTCGACCTCGCG‑GTTTTGCTGCCCTTTGTACCATCCATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCATGATGATTTGACGTCATCCCCACCt                            >  2:1976148/7‑150 (MQ=0)

                                                                                                taggcGAGTTGCAGCCTACAATCCGAACTGAGAATGGTTTTATGGGATTTGCTCGACCTCGCG‑GTTTTGCTGCCCTTTGTACCATCCATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCATGATGATTTGACGTCATCCCCACCTtc                          <  2:8834504/146‑1 (MQ=0)

                                                                                                taggcGAGTTGCAGCCTACAATCCGAACTGAGAATGGTTTTATGGGATTTGCTCGACCTCGCG‑GTTTTGCTGCCCTTTGTACCATCCATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCATGATGATTTGACGTCATCCCCACCTtc                          >  2:6873893/5‑150 (MQ=0)

                                                                                                taggcGAGTTGCAGCCTACAATCCGAACTGAGAATGGTTTTATGGGATTTGCTCGACCTCGCG‑GTTTTGCTGCCCTTTGTACCATCCATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCATGATGATTTGACGTCATCCCCACCTtc                          >  2:6350083/5‑150 (MQ=0)

                                                                                                taggcGAGTTGCAGCCTACAATCCGAACTGAGAATGGTTTTATGGGATTTGCTCGACCTCGCG‑GTTTTGCTGCCCTTTGTACCATCCATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCATGATGATTTGACGTCATCCCCACCTtc                          >  2:6350059/5‑150 (MQ=0)

                                                                                                cactCGAGTTGCAGAGTGCAATCCGAACTGAGATGGCTTTTG‑GAGATTAGCTCACACTCGCGTGCTT‑GCTGCCCACTGTCACCACCATTGTAGCACGTGTGTAGCCCAGCCCGTAAGGGCCATGAGGACTTGACGTCATCCCCACCTtcc                         <  2:3005396/147‑1 (MQ=0)

                                                                                                aaggcGAGTTGCAGCCTACAATCCGAACTGAGAATGGTTTTATGGGATTTGCTCGACCTCGCG‑GTTTTGCTGCCCTTTGTACCATCCATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCATGATGATTTGACGTCATCCCCACCTtc                          >  2:6873868/5‑150 (MQ=0)

                                                                                                 aggcGAGTTGCAGCCTACAATCCGAACTGAGAATGGTTTTATGGGATTTGCTCGACCTCGCG‑GTTTTGCTGCCCTTTGTACCATCCATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCATGATGATTTGACGTCATCCCCACCTtcc                         >  2:2977507/4‑150 (MQ=0)

                                                                                                 aggcGAGTTGCAGCCTACAATCCGAACTGAGAATGGTTTTATGGGATTTGCTCGACCTCGCG‑GTTTTGCTGCCCTTTGTACCATCCATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCATGATGATTTGACGTCATCCCCACCTtcc                         >  2:2819442/4‑150 (MQ=0)

                                                                                                 aggcGAGTTGCAGCCTACAATCCGAACTGAGAATGGTTTTATGGGATTTGCTCGACCTCGCG‑GTTTTGCTGCCCTTTGTACCATCCATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCATGATGATTTGACGTCATCCCCACCTtcc                         >  2:7438776/4‑150 (MQ=0)

                                                                                                 aggcGAGTTGCAGCCTACAATCCGAACTGAGAATGGTTTTATGGGATTTGCTCGACCTCGCG‑GTTTTGCTGCCCTTTGTACCATCCATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCATGATGATTTGACGTCATCCCCACCTtcc                         >  2:7438887/4‑150 (MQ=0)

                                                                                                             gagtGCAATCCGAACTGAGATGGCTTTTG‑GAGATTAGCTCACACTCGCGTGCTT‑GCTGCCCACTGTCACCACCATTGTAGCACGTGTGTAGCCCAGCCCGTAAGGGCCATGAGGACTTGACGTCATCCCCACCTTCCTCCAGCTTATCAc            >  1:8029712/4‑150 (MQ=0)

                                                                                                                 gcaaTCCGAACTGAGATGGCTTTTG‑GAGATTAGCTCACACTCGCGTGCTT‑GCTGCCCACTGTCACCACCATTGTAGCACGTGTGTAGCCCAGCCCGTAAGGGCCATGAGGACTTGACGTCATCCCCACCTTCCTCCAGCTTATCACTGGc        <  2:2183364/147‑1 (MQ=0)

                                                                                                                 gcaaTCCGAACTGAGATGGCTTTTG‑GAGATTAGCTCACACTCGCGTGCTT‑GCTGCCCACTGTCACCACCATTGTAGCACGTGTGTAGCCCAGCCCGTAAGGGCCATGAGGACTTGACGTCATCCCCACCTTCCTCCAGCTTATCACTGGc        <  2:7737934/147‑1 (MQ=0)

                                                                                                                        ggccTACGATCGGTTTTGTGAGATTAGCTCCACCTCGCG‑GCTTGGCAACCCTCTGTACCGACCATTGTAGCACGTGTGTAGCCCAGGCCGTAAGGGCCATGATGACTTGACGTCATCCCCACCTTCCTCCGGTTTGTCACCGGCAGtctc  <  1:6101721/147‑1 (MQ=1)

CCCACTCCCATGGTGTGACGGGCGGTGTGTACAAGGCCCGGGAACGTATTCACCGCGACATTCTGATTCGCGATTACTAGCGATTCCGACTTCACGCAGTCGAGTTGCAGACTGCGATCCGGACTACGATCGGTTTTGTGAGATTAGCTCCACCTCGCG‑GCTTGGCAACCCTCTGTACCGACCATTGTAGCACGTGTGTAGCCCAGGCCGTAAGGGCCATGATGACTTGACGTCATCCCCACCTTCCTCCGGTTTGTCACCGGCAGTCTC  >  NC_002947/5313249‑5313518

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