BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,156,666	T→A	23.5%	Y24F (TAT→TTT)	PP_1909 ←	septum formation protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,156,666	0	T	A	23.5%	26.5 / 6.4	17	Y24F (TAT→TTT)	PP_1909	septum formation protein
Reads supporting (aligned to +/- strand): ref base T (5/8); new base A (2/2); total (7/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 = 5.60e-01

GCCTCGGCCTTTTGCCTGGCCAGCCGGCGCACCAGCTCTACGGGCGGTTCGCCATCCAGGCGTTGCTCGTCTATATCGGGGCTTGCCCAGGTGAAGGGCAGGTGCAGGCGTGCGAGCAGTTCGCGGCGATAGGCAGAGCTGGAAGCCAGTAACAAAGGAAGCATGGTAGACTCCGGGTTCGGTTGAACCAATTCTATCACGCACGATGCCGCCTAATTTCCTTTGACAGGGGCGGGGGGC  >  NC_002947/2156537‑2156776

gCCTCGGCCTTTTGCCTGGCCAGCCGGCGCACCAGCTCTACGGGCGGTTCGCCATCCAGGCGTTGCTCGTCTATATCGGGGCTTGCCCAGGTGAAGGGCAGGTGCAGGCGTGCGAGCAGTTCGCGGCGATAGGCag                                                                                                          >  2:190832/1‑136 (MQ=255)

       ccTTTTGCCTGGCCAGCCGGCGCACCAGCTCTACGGGCGGTTCGCCATCCAGGCGTTGCTCGTCTATATCGGGGCTTGCCCAGGTGAAGGGCAGGTGCAGGCGTGCGAGCAGTTCGCGGCGATAGGCAGAGCTGGa                                                                                                   <  1:271/136‑1 (MQ=255)

       ccTTTTGCCTGGCCAGCCGGCGCACCAGCTCTACGGGCGGTTCGCCATCCAGGCGTTGCTCGTCTATATCGGGGCTTGCCCAGGTGAAGGGCAGGTGCAGGCGTGCGAGCAGTTCGCGGCGATAGGCAGAGCTGGa                                                                                                   <  1:98927/136‑1 (MQ=255)

                    cAGCCGGCGCACCAGCTCTACGGGCGGTTCGCCATCCAGGCGTTGCTCGTCTATATCGGGGCTTGCCCAGGTGAAGGGCAGGTGCAGGCGTGCGAGCAGTTCGCGGCGAAAGGCAGAGCTGGAAGCCAGTAACaaa                                                                                      >  1:187057/1‑136 (MQ=255)

                    cAGCCGGCGCACCAGCTCTACGGGCGGTTCGCCATCCAGGCGTTGCTCGTCTATATCGGGGCTTGCCCAGGTGAAGGGCAGGTGCAGGCGTGCGAGCAGTTCGCGGCGAAAGGCAGAGCTGGAAGCCAGTAACaaa                                                                                      >  1:188357/1‑136 (MQ=255)

                     aGCCGGCGCACCAGCTCTACGGGCGGTTCGCCATCCAGGCGTTGCTCGTCTATATCGGGGCTTGCCCAGGTGAAGGGCAGGTGCAGGCGTGCGAGCAGTTCGCGGCGATAGGCAGAGCTGGAAGCCAGTAACAAAg                                                                                     <  2:182866/136‑1 (MQ=255)

                                   ctctACGGGCGGTTCGCCATCCAGGCGTTGCTCGTCTATATCGGGGCTTGCCCAGGTGAAGGGCAGGTGCAGGCGTGCGAGCAGTTCGCGGCGATAGGCAGAGCTGGAAGCCAGTAACAAAGGAAGCATGGTAGAc                                                                       <  2:212585/136‑1 (MQ=255)

                                                 cGCCATCCAGGCGTTGCTCGTCTATATCGGGGCTTGCCCAGGTGAAGGGCAGGTGCAGGCGTGCGAGCAGTTCGCGGCGATAGGCAGAGCTGGAAGCCAGTAACAAAGGAAGCATGGTAGACTCCGGGTTCGGTTg                                                         <  1:72089/136‑1 (MQ=255)

                                                                 cTCGTCTATATCGGGGCTTGCCCAGGTGAAGGGCAGGTGCAGGCGTGCGAGCAGTTCGCGGCGATAGGCACAGCTGGAAGACAGTAACAAAGGAAGCATGGAAGACTCCGGGTTCGGTTGAACCAATTCTATcacc                                         >  2:47068/1‑135 (MQ=255)

                                                                                   tGCCCAGGTGAAGGGCAGGTGCAGGCGTGCGAGCAGTTCGCGGCGATAGGCAGAGCTGGAAGCCAGTAACAAAGGAAGCATGGTAGACTCCGGGTTCGGTTGAACCAATTCTATCACGCACGATGACGCGTAAttt                       >  2:129815/1‑136 (MQ=255)

                                                                                                    ggTGCAGGCGTGCGAGCAGTTCGCGGCGATAGGCAGAGCTGGAAGCCAGTAACAAAGGAAGCATGGTAGACTCCGGGTTCGGTTGAACCAATTCTATCACGCACGATGCCGCCTAATTTCCTTTGACAGGGGCggg      <  1:190832/136‑1 (MQ=255)

                                                                                                     gTGCAGGCGTGCGAGCAGTTCGCGGCGAAAGGCAGAGCTGGAAGCCAGTAACAAAGGAAGCATGGTAGACTCCGGGTTCGGTTGAACCAATTCTATCACGCACGATGCCGCCTAATTTCCTTTGACAGGGGCgggg     <  2:187057/136‑1 (MQ=255)

                                                                                                     gTGCAGGCGTGCGAGCAGTTCGCGGCGAAAGGCAGAGCTGGAAGCCAGTAACAAAGGAAGCATGGTAGACTCCGGGTTCGGTTGAACCAATTCTATCACGCACGATGCCGCCTAATTTCCTTTGACAGGGGCgggg     <  2:188357/136‑1 (MQ=255)

                                                                                                       gCAGGCGTGCGAGCAGTTCGCGGCGATAGGCAGAGCTGGAAGCCAGTAACAAAGGAAGCATGGTAGACTCCGGGTTCGGTTGAACCAATTCTATCACGCACGATGCCGCCTAATTTCCTTTGACAGGGGCgggggg   <  1:83529/136‑1 (MQ=255)

                                                                                                        cAGGCGTGCGAGCAGTTCGCGGCGATAGGCAGAGCTGGAAGCCAGTAACAAAGGAAGCATGGTAGACTCCGGGTTCGGTTGAACCAATTCTATCACGCACGATGCCGCCTAATTTCCTTTGACAGGGGCGGGGGGc  >  2:21353/1‑136 (MQ=255)

                                                                                                                            ggcgATAGGCAGAGCTGGAAGCCAGTAACAAAGGAAGCATGGTAGACTCCGGGTTCGGTTGAACCAATTCTATCACGCACGATGCCGCCTAATTTCCttt                  >  2:149663/1‑100 (MQ=255)

                                                                                                                            ggcgATAGGCAGAGCTGGAAGCCAGTAACAAAGGAAGCATGGTAGACTCCGGGTTCGGTTGAACCAATTCTATCACGCACGATGCCGCCTAATTTCCttt                  <  1:149663/100‑1 (MQ=255)

GCCTCGGCCTTTTGCCTGGCCAGCCGGCGCACCAGCTCTACGGGCGGTTCGCCATCCAGGCGTTGCTCGTCTATATCGGGGCTTGCCCAGGTGAAGGGCAGGTGCAGGCGTGCGAGCAGTTCGCGGCGATAGGCAGAGCTGGAAGCCAGTAACAAAGGAAGCATGGTAGACTCCGGGTTCGGTTGAACCAATTCTATCACGCACGATGCCGCCTAATTTCCTTTGACAGGGGCGGGGGGC  >  NC_002947/2156537‑2156776

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