BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NZ_CP034943	2,981,253	C→T	20.0%	noncoding (1267/1550 nt)	EO946_RS15435 ←	16S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP034943	2,981,253	0	C	T	20.0%	‑2.4 / 8.4	25	noncoding (1267/1550 nt)	EO946_RS15435	16S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base C (2/18); new base T (3/2); total (5/20)
Fisher's exact test for biased strand distribution p-value = 3.77e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

TGTGTACAAGGCCCGGGAACGTATTCACCGCGGCATGCTGATCCGCGATTACTAGCGATTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTCGCTGCCCTTTGTTCTGTCCATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCATGATGATTTGACGTCATCCCCACCTTCCTCCGGTTTGTCACCGGCAGTCACCTTAGAGTGCCCAACTGAATG  >  NZ_CP034943/2981115‑2981380

tgggtACAAGGCCCGGGAACGTTTTCACCGCGGCATGCTGATCCGCGATTACTAGCGATTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTCg                                                                                                                                <  2:2499344/137‑1 (MQ=1)

  tggaCAAGGCCCGGGAACGTATTCACCGCGGCATGCTGATCCGCGATTACTAGCGATTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTCGCt                                                                                                                              <  2:1854081/137‑1 (MQ=1)

  tggaCAAGGCCCGGGAACGTATTCACCGCGGCATGCTGATCCGCGATTACTAGCGATTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTAAACCTCGCGGTTTCGCt                                                                                                                              <  1:2011558/137‑1 (MQ=1)

  tggaCAAGGCCCGGGAACGTATTCACCGCGGCATGCTGATCCGCGATTACGAGCGAGTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTCGCt                                                                                                                              <  2:2844849/137‑1 (MQ=0)

                            cggggCATGCTGATCCGCGATTACTAGCGAGTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTTGCTGCCCTTTGTTCTGTCCATTGTAGCAc                                                                                                    <  2:1604814/137‑1 (MQ=0)

                            cgaggCATGCTGATCCGCGATTACTAGCGATTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTCACCTCGCGGTTTCGCTGCCCTTTGTTCTGTCCATTGTAGCAc                                                                                                    <  1:1348945/137‑1 (MQ=1)

                                 caggctgctCCGCGATTACTAGCGATTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTGGTGGGATTGGCGTAACCTCGCGGTTTCGCTGCCCTCTGTTCTGTCCATTGTAGCACgtgtg                                                                                               <  1:2833393/137‑1 (MQ=1)

                                 caggCTGATCCGCGATTACTAGCGATTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTCGCTGCCCTTTGTTCTGTCCATTGTAGCACgtgtg                                                                                               <  1:1852228/137‑1 (MQ=1)

                                 cacgCTGATCCGCGATTACTAGCGATTCCAGCTTCACGCAGTCGAGTTGCAGCCTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTAAACCTCGCGTTTTCGCTGCCCTTTGTTCTGTCGATTGTAGCACgtgtg                                                                                               <  1:298053/137‑1 (MQ=1)

                                               atgacgaGCGATTCCAGCTTCACGCAGTCGCGTTGCAGACTGCGATCCGAACTGAGAACAGATGTGTGGGATTGGGTTAACCTCGCGGTTTCGCTGCCCTTTGTTCTGTCCCTTGTAGCACGTGTGTCGCCCAGGTCATa                                                                                 >  2:1087789/7‑140 (MQ=1)

                                               atgaCTAGCGATTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTCGCGGCCCTTTGTTCTGTCCATTGTAGCACGTGTGTAGCCCAGGTCATa                                                                                 <  2:671/137‑1 (MQ=1)

                                               ataaCTAGCGATTCCAGCTTCACTCATTAGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTTGCTGCCCTTTGTTCTGTCCATTGTAGCACGTGTGTAGCCCAGGTCATa                                                                                 >  2:219577/4‑140 (MQ=1)

                                               ataaCTAGCGATTCCAGCTTCACTAATTAGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTTGCTGCCCTTTGTTCTGTCCATTGTAGCACGTGTGTAGCCCAGGTCATa                                                                                 >  2:1591043/4‑140 (MQ=1)

                                               ataaCTAGCGATTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTTGCTGCCCTTTGTTCTGTCCATTGTAGCACGTGTGTAGCCCAGGTCATa                                                                                 >  2:219524/4‑140 (MQ=1)

                                                        gagtCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTCGCTGCCCTTTGTTCTGTCCATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCatg                                                                        <  2:2694089/137‑1 (MQ=1)

                                                        gagtCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTCGCTGCCCTTTGGTCTGTCCATTGGAGCACGTGTGTAGCCCCGGTCATAAGGGGCatg                                                                        <  2:767182/137‑1 (MQ=0)

                                                                          gcgcgTTGCAGACTGCGATCCGAACAGAGAACAGATTTGTGGGAGTGGCATACCCTCGCGGTTGCGCTGCCCTTTGTTCTGTGCCATGTAGCACGTGTGTAGCCCAGCTCATAAGGGGCATGATGATTTGACGTCATccg                                                      <  2:325874/136‑2 (MQ=1)

                                                                              gtagCAGACTGCGATCCGCACTGAGAACAGATTTGTGGGATGGGCTTAACCTCGCGGATTCGCTGCCCTTTGTTCTGTCCATTGTACCACGTGTGTAGCCCAGGTCATAAGGGGCATGATGATTTGACGTCATCCCCAcc                                                  <  2:381805/137‑1 (MQ=1)

                                                                                               gatcTGAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTTGCTGCCCTTTGTTCTGTCCATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCATGATGATTTGACGTCATCCCCACCTTCCTCCGGTTTGTCAc                                 <  2:1414216/137‑1 (MQ=1)

                                                                                               gaccTGAGAACAGATTTGTGGGATTGGCTTCACCTCGCGGTTTCGCTGCCCAGTGTTCGGTCCATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCATGATGATTTGACGTCATCCCCACCTTCCTCCGGTTTGGCAc                                 <  1:2425451/137‑1 (MQ=1)

                                                                                                     cgcaccgATCTGTGGGAGTGGCTCACCCACGCGGTGTCGCTGCCCTATGTGCTGGCCATTGTAGCCCCTGTGGACCCCAGGTCATAAGGGGCATGATGAGTGGACGTCATCCCCACCTTCCACCGGTGTGTCACCGGCAg                           <  2:1880676/134‑1 (MQ=0)

                                                                                                     agtaCAGATTTGTGGGATTGGCTTAACCTCGCGGTTTCGCTGCCCTTTGTTCTGTCCATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCATGATGATTTGACGTCATCCCCACCTTCCTCCGGTTTGTCACCGGCAg                           <  1:85813/137‑1 (MQ=1)

                                                                                                                    gagtGGCTTAACCTCGCGGTTTCGCTTCCCTTTGTTCTGTCCTTGGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCATGATGATTTGACGTCATCCCCACCTTCCTCCGGTTTGTCACCGGCAGTCACCTTAGAGTGcc            <  2:2516185/137‑1 (MQ=0)

                                                                                                                    gagtGGCTGCACCGCGCGGTTTCGCTGCCCTTTGTTCGGTCCATCGGAGCCCGTGTGTAGCCCAGGTCATAAGGGGCAAGATGACTTGACGTCATCCCCACCGTCCTCCGGTTTGTCACCGGCAGTCACCTTAGAGTGcc            <  2:955521/137‑1 (MQ=1)

                                                                                                                              aactagcgGTTTCGCAGAAATTTTTTATGAACATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCATGATGATTTGACGTCATCCCCACCTTCCTCCGGTTTGTCACCGGCAGTCACCTTAGAGTGCCCAACTGAATg  >  2:314594/6‑140 (MQ=1)

TGTGTACAAGGCCCGGGAACGTATTCACCGCGGCATGCTGATCCGCGATTACTAGCGATTCCAGCTTCACGCAGTCGAGTTGCAGACTGCGATCCGAACTGAGAACAGATTTGTGGGATTGGCTTAACCTCGCGGTTTCGCTGCCCTTTGTTCTGTCCATTGTAGCACGTGTGTAGCCCAGGTCATAAGGGGCATGATGATTTGACGTCATCCCCACCTTCCTCCGGTTTGTCACCGGCAGTCACCTTAGAGTGCCCAACTGAATG  >  NZ_CP034943/2981115‑2981380

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