BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,763,636	A→G	36.1%	noncoding (208/1542 nt)	rrsG ←	16S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,763,636	0	A	G	36.1%	17.8 / 18.9	36	noncoding (208/1542 nt)	rrsG	16S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base A (18/5); new base G (2/11); total (20/16)
Fisher's exact test for biased strand distribution p-value = 4.09e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTAT  >  minE/1763575‑1763695

gAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGATGCCCGAGGGTcccc                                                      <  1:1328425/69‑1 (MQ=255)

gAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTcccc                                                      <  1:113188/69‑1 (MQ=255)

gAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTcccc                                                      <  1:1218165/69‑1 (MQ=255)

gAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTcccc                                                      <  1:1595820/69‑1 (MQ=255)

gAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTcccc                                                      <  1:1700426/69‑1 (MQ=255)

   ccGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCctc                                                   >  1:1709555/1‑69 (MQ=255)

   ccGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCctc                                                   >  1:2990410/1‑69 (MQ=255)

   ccGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCct                                                    <  1:289842/68‑1 (MQ=255)

   ccGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCctc                                                   >  1:1045819/1‑69 (MQ=40)

   ccGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCctc                                                   >  1:946530/1‑69 (MQ=40)

   ccGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCctc                                                   >  1:914026/1‑69 (MQ=40)

   ccGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCctc                                                   >  1:3399364/1‑69 (MQ=40)

   ccGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCctc                                                   >  1:2409185/1‑69 (MQ=40)

   ccGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCct                                                    >  1:714772/1‑68 (MQ=40)

            caccaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCtt                                          <  1:1102816/69‑1 (MQ=255)

            caccaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCtt                                          <  1:1802454/69‑1 (MQ=255)

            caccaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCtt                                          <  1:3368617/69‑1 (MQ=255)

            caccaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCtt                                          <  1:2910124/69‑1 (MQ=255)

            caccaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCtt                                          <  1:1806425/69‑1 (MQ=255)

                        aaTCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGc                              >  1:245631/1‑69 (MQ=35)

                        aaTCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGc                              >  1:269376/1‑69 (MQ=35)

                        aaTCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGc                              >  1:41985/1‑69 (MQ=35)

                        aaTCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGc                              >  1:342302/1‑69 (MQ=35)

                        aaTCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGc                              >  1:2961717/1‑69 (MQ=35)

                        aaTCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGc                              >  1:3025312/1‑69 (MQ=35)

                        aaTCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGc                              >  1:3189950/1‑69 (MQ=35)

                        aaTCCCATCTGGGCACATCCGATGGAAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGc                              >  1:2043438/1‑69 (MQ=25)

                            ccATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCACTCTTTGGTCTTGCGACATTATGCGGt                           >  1:257980/1‑68 (MQ=25)

                                     cacaTCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACAt                                   >  1:273300/1‑51 (MQ=33)

                                     cacaTCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACAt                                   >  1:814151/1‑51 (MQ=33)

                                     cacaTCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACAt                                   >  1:101855/1‑51 (MQ=33)

                                                  aaGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGtt    <  1:333714/69‑1 (MQ=35)

                                                  aaGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGtt    <  1:1935880/69‑1 (MQ=35)

                                                  aaGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGtt    <  1:2872874/69‑1 (MQ=35)

                                                    gagGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTAt  <  1:2698081/69‑1 (MQ=35)

                                                    gagGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTAt  <  1:2046487/69‑1 (MQ=35)

GAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTAT  >  minE/1763575‑1763695

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