BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,329,744	C→A	11.4%	R38R (CGG→AGG)	yidH →	inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,329,744	0	C	A	11.4%	108.9 / 2.9	35	R38R (CGG→AGG)	yidH	inner membrane protein
Reads supporting (aligned to +/- strand): ref base C (16/14); new base A (2/2); total (18/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.56e-01

CTTTGGTGTGCGCTGATGGGCGTGCCTGAACCACGGCCTGACTGGTCGCGGCGGCTGTTGGGGCAGGGTGAGGCGCCTGATCCGCGCTTCACGCTGGCCAACGAACGTACCTTCCTGGCCTGGGTCCGGACCGCGCTGGCGCTGCTGGGTGGCGCAATCGCCATCGAAACTTTCGCCGGCCACGTCCTCGAGACACCCCTGCGGCTGTGGCTGGTCGCTGGCCTGATGCTGCTGAGTGCATTGCTGAGCGCTGGCGCCTGC  >  NC_002947/2329618‑2329878

cTTTGGTGTGCGCTGATGGGCGTGCCTGAACCACGGCCTGACTGGTCGCGGCGGCTGTTGGGGCAGGGTGAGGCGCCTGATCCGCGCTTCACGCTGGCCAACGAACGTACCTTCCTGGCCTGGGTCCGGACCgcgc                                                                                                                               <  1:71144/136‑1 (MQ=255)

     gtgtGCGCTGATGGGCGTGCCTGAACCACGGCCTGACTGGTCGCGGCGGCTGTTGGGGCAGGGTGAGGCGCCTGATCCGCGCTTCACGCTGGCCAACGAACGTACCTTCCTGGCCTGGGTCCGGACCgcgctggcg                                                                                                                          <  2:432733/136‑1 (MQ=255)

             tGATGGGCGTGCCTGAACCACGGCCTGACTGGTCGCGGCGGCTGTTGGGGCAGGGTGAGGCGCCTGATCCGCGCTTCACGCTGGCCAACGAACGTACCTTCCTGGCCTGGGTCCGGACCGCGCTGGCGCTGCTggg                                                                                                                  >  1:198493/1‑136 (MQ=255)

             tGATGGGCGTGCCTGAACCACGGCCTGACTGGTCGCGGCGGCTGTTGGGGCAGGGTGAGGCGCCTGATCCGCGCTTCACGCTGGCCAACGAACGTACCTTCCTGGCCTGGGTCCGGACCGCGCTGGCGCTGCTggg                                                                                                                  >  1:585251/1‑136 (MQ=255)

             tGATGGGCGTGCCTGAACCACGGCCTGACTGGTCGCGGCGGCTGTTGGGGCAGGGTGAGGCGCCTGATCCGCGCTTCACGCTGGCCAACGAACGTACCTTCCTGGCCTGGGTCCGGACCGCGCTGGCGCTGCTggg                                                                                                                  >  1:502843/1‑136 (MQ=255)

                       gCCTGAACCACGGCCTGACTGGTCGCGGCGGCTGTTGGGGCAGGGTGAGGCGCCTGATCCGCGCTTCACGCTGGCCAACGAACGTACCTTCCTGGCCTGGGTCCGGACCGCGCTGGCGCTGCTGGGTGGCGGAATc                                                                                                        >  2:163892/1‑136 (MQ=255)

                                    ccTGACTGGTCGCGGCGGCTGTTGGGGCAGGGTGAGGCGCCTGATCCGCGCTTCACGCTGGCCAACGAACGTACCTTCCTGGCCTGGGTCCGGACCGCGCTGGCGCTGCTGGGTGGCGCAATCGCCATCGAAACtt                                                                                           >  2:656387/1‑136 (MQ=255)

                                     cTGACTGGTCGCGGCGGCTGTTGGGGCAGGGTGAGGCGCCTGATCCGCGCTTCACGCTGGCCAACGAACGTACCTTCCTGGCCTGGGTCAGGACCGCGCTGGCGCTGCTGGGTGGCGCAATCGCCATCGAAACttt                                                                                          >  1:673726/1‑136 (MQ=255)

                                     cTGACTGGTCGCGGCGGCTGTTGGGGCAGGGTGAGGCGCCTGATCCGCGCTTCACGCTGGCCAACGAACGTACCTTCCTGGCCTGGGTCAGGACCGCGCTGGCGCTGCTGGGTGGCGCAATCGCCATCGAAACttt                                                                                          >  1:612536/1‑136 (MQ=255)

                                          gggTCGCGGCGGCTGGTGGGGCAGGGTGAGGCGACTGATCCGCGCTTCACGCTTGCCAACGAACGTACCTTCCTGGCCTGGGTCCGGACCGCGCTGGCGCTGCTGGGTGGCGCAATCGCCATCGAAACTTTCGCCg                                                                                     <  1:196229/135‑1 (MQ=255)

                                           ggTCGCGGCGGCTGTTGGGGCAGGGTGAGGCGCCTGATCCGCGCTTCACGCTGGCCAACGAACGTACCTTCCTGGCCTGGGTCCGGACCGCGCTGGCGCTGCTGGGTGGCGCAATCGCCATCGAAACTTTCGCCgg                                                                                    >  1:613060/1‑136 (MQ=255)

                                               gcggcgGCTGTTGGGGCAGGGTGAGGCGCCTGATCCGCGCTTCACGCTGGCCAACGAACGTACCTTCCTGGCCTGGGTCCGGACCGCGCTGGCGCTGCTGGGTGGCGCAATCGCCATCGAAACTTTCGCCGGCCAc                                                                                >  1:432590/1‑136 (MQ=255)

                                                cggcggCTGTTGGGGCAGGGTGAGGCGCCTGATCCGCGCTTCACGCTGGCCAACGAACGTACCTTCCTGGCCTTGGTCCGGAACGCGCTGGCGCTGCTGGGTGGCGCAATCGCCATCGAAACATTAGCCGGCCAcc                                                                               >  2:672228/1‑135 (MQ=255)

                                                         ttGGGGCAGGGTGAGGCGCCTGATCCGCGCTTCACGCTGGCCAACGAACGTACCTTCCTGGCCTGGGTCCGGACCGCGCTGGCGCTGCTGGGTGGCGCAATCGCCATCGAAACTTTCGCCGGCCACGTCCTCgaga                                                                      <  1:218391/136‑1 (MQ=255)

                                                             ggCAGGGTGAGGCGCCTGATCCGCGCTGCACGCTGGCCAACGAACGTACCTTCCTGGCCTGGGTCCGGACCGCGCTGGCGCTGCTGGGTGGCGCAATCGCCATCGAAACTTTCGCCGGCCACGTCCTCGAGACAcc                                                                  <  2:161931/136‑1 (MQ=255)

                                                              gCAGGGTGAGGCGCCTGATCCGCGCTTCACGCTGGCCAACGAACGTACCTTCCTGGCCTGGGTCCGGACCGCGCTGGCGCTGCTGGGTGGCGCAATCGCCATCGAAACTTTCGCCGGCCACGTCCTCGAGACAccc                                                                 <  1:65738/136‑1 (MQ=255)

                                                                     gAGGCGCCTGATCCGCGCTTCACGCTGGCCAACGAACGTACCTTCCTGGCCTGGGTCCGGACCGCGCTGGCGCTGCTGGGTGGCGCAATCGCCATCGAAACTTTCGCCGGc                                                                                   <  2:652025/111‑1 (MQ=255)

                                                                     gAGGCGCCTGATCCGCGCTTCACGCTGGCCAACGAACGTACCTTCCTGGCCTGGGTCCGGACCGCGCTGGCGCTGCTGGGTGGCGCAATCGCCATCGAAACTTTCGCCGGc                                                                                   >  1:652025/1‑111 (MQ=255)

                                                                         cgcCTGATCCGCGCTTCACGCTGGCCAACGAACGTACCTTCCTGGCCTGGGTCAGGACCGCGCTGGCGCTGCTGGGTGGCGCAATCGCCATCGAAACTTTCGCCGGCCACGTCCTCGAGACACCCCTGCGGCtgtg                                                      <  2:612536/136‑1 (MQ=255)

                                                                         cgcCTGATCCGCGCTTCACGCTGGCCAACGAACGTACCTTCCTGGCCTGGGTCAGGACCGCGCTGGCGCTGCTGGGTGGCGCAATCGCCATCGAAACTTTCGCCGGCCACGTCCTCGAGACACCCCTGCGGCtgtg                                                      <  2:673726/136‑1 (MQ=255)

                                                                              gATCCGCGCTTCACGCTGGCCAACGAACGTACCTTCCTGGCCTGGGTCCGGACCGCGCTGGCGCTGCTGGGTGGCGCAATCGCCATCGAAACTTTCGCCGGCCCCGTCCTCGAGACACCCCTGCGGCTGTGGCTgg                                                 >  2:165939/1‑136 (MQ=255)

                                                                               aTCCGCGCTTCACTCTGGCCAACGAACGTACCTTCCTGGCCTGGTTCCGGACCGCGCTGGCGCTGCTGGGTGGCGCAATCGCCATCGAAACTTTCGCCGGCCACGTCCTCGAGACACCCCTCCGGCTGTGGCTGGt                                                <  2:665572/136‑1 (MQ=255)

                                                                               aTCCGCGCTTCACGCTGGCCAACGAACGTAGCTTCCTGGCCTGGGTCGGGTCCGCGCTGGCGGTGCTGGGTGGCGGAATCGCCATCGATACTTTCGCCGGGCACGACCTCGAGTGACGCCTGGGGCTGTGGCTGGt                                                <  2:663312/136‑1 (MQ=255)

                                                                                       ttCACGCTGGCCAACGAACGTACCTTCCTGGCCTGGGTCCGGACCGCGCTGGCGCTGCTGGGTGGCGCAATCGCCATCGAAACTTTCGCCGGCCACGTCCTCGAGACACCCCTGCGGCTGTGGCTGGTCGCTGGcc                                        >  1:674646/1‑136 (MQ=255)

                                                                                       ttCACGCTGGCCAACGAACGTACCTTCCTGGCCTGGGTCCGGACCGCGCTGGCGCTGCTGGGTGGCGCAATCGCCATCGAAACTTTCGCCGGCCACGTCCTCGAGACACCCCTGCGGCTGTGGCTGGTCGCTGGcc                                        >  1:141221/1‑136 (MQ=255)

                                                                                            gCTGGCCAACGAACGTACCTTCCTGGCCTGGGTCCGGACCGCGCTGGCGCTGCTGGGTGGCGCAATCGCCATCGAAACTTTCGCCGGCCACGTCCTCGAGACACCCCTGCGGCTGTGGCTGGTCGCTGGCCTGAtg                                   <  1:173354/136‑1 (MQ=255)

                                                                                              tGGCCAACGAACGTACCTTCCTGGCCTGGGTCCGGACCGCGCTGGCGCTGCTGGGTGGCGGAATCGCCATCGAAACTTTCGCCGGCCACGTCCTCGAGACACCCCTGCGGCTGTGGCTGGTCGCTGGCCTGAtgct                                 <  1:163892/136‑1 (MQ=255)

                                                                                              tGGCCAACGAACGTACCTTCCTGGCCTGGGTCCGGACCGCGCTGGCGCTGCTGGGTGGCGCAATCGCCATCGAAACTTTCGCCGGCCACGTCCTCTAGACACCCCTGCGGCTGTGGCTGGTCGCTGGCCTGAtgct                                 <  1:428305/136‑1 (MQ=255)

                                                                                              tGGCCAACGAACGTACCTTCCTGGCCTGGGTCCGGACCGCGCTGGCGCTGCTGGGTGGCGCAATCGCCATCGAAACTTTCGCCGGCCACGTCCTCGAGACACCCCTGCGGCTGTGGCTGGTCGCTGGCCTGAtgct                                 <  2:405844/136‑1 (MQ=255)

                                                                                              tGGCCAACGAACGTACCTTCCTGGCCTGGGTCCGGACCGCGCTGGCGCTGCTGGGTGGCGCAATCGCCATCGAAACTTTCGCCGGCCACGTCCTCGAGACACCCCTGCGGCTGTGGCTGGTCGCTGGCCTGAtgct                                 <  1:172729/136‑1 (MQ=255)

                                                                                                gCCAACGAACGTACCTTCCTGGCCTGGGTCCGGACCGCGCTGGCGCTGCTGGGTGGCGCAATCGCCATCGAAACTTTCGCCGGCCACGTCCTCGAGACACCCCTGCGGCTGTGGCTGGTCGCTGGCCTGAtgctgc                               <  2:445239/136‑1 (MQ=255)

                                                                                                             ccTTCCTGGCCTGGGTCCGGACCGCGCTGGCGCTGCTGGGTGGCGCAATCGCCATCGAAACTTTCGCCGGCCACGTCCTCGAGACACCCCTGCGGCTGTGGCTGGTCGCTGGCCTGATGCTGCTGAGTGCATTGCt                  >  1:549836/1‑136 (MQ=255)

                                                                                                             ccTTCCTGGCCTGGGTCCGGACCGCGCTGGCGCTGCTGGGTGGCGCAATCGCCATCGAAACTTTCGCCGGCCACGTCCTCGAGACACCCCTGCGGCTGTGGCTGGTCGCTGGCCTGATGCTGCTGAGTGCATTGCt                  >  1:157198/1‑136 (MQ=255)

                                                                                                             ccTTCCTGGCCTGGGTCCGGACCGCGCAGGCGCTGCTGGGAGGCGCAATCGCCAACGAAACTATCGCCGGCCACGTCCTCGAGACCCCCCTGCGGCTGTGGCTGGTCGCTGGCCTGATGCTACAGAGTGCATTTCt                  >  1:106597/1‑136 (MQ=255)

                                                                                                                     gcctggGTCCGGACCGCGCTGGCGCTGCTGGGTGGCGCAATCGCCATCGAAACTTTCGCCGGCCACGTCCTCGAGACACCCCTGCGGCTGTGGCTGGTCGCTGGCCTGATGCTGCTGAGTGCATTGCTGAGCGCTg          >  1:129605/1‑136 (MQ=255)

                                                                                                                             ccGGACCGCGCTGGCGCTGCTGGGTGGCGCAATCGCCATCGAAACTTTCGCCGGCCACGTCCTCGAGACACCCCTGCGGCTGTGGCTGGTCGCTGGCCTGATGCTGCTGAGTGCATTGCTGAGCGCTGGCgcctgc  >  2:401112/1‑136 (MQ=255)

CTTTGGTGTGCGCTGATGGGCGTGCCTGAACCACGGCCTGACTGGTCGCGGCGGCTGTTGGGGCAGGGTGAGGCGCCTGATCCGCGCTTCACGCTGGCCAACGAACGTACCTTCCTGGCCTGGGTCCGGACCGCGCTGGCGCTGCTGGGTGGCGCAATCGCCATCGAAACTTTCGCCGGCCACGTCCTCGAGACACCCCTGCGGCTGTGGCTGGTCGCTGGCCTGATGCTGCTGAGTGCATTGCTGAGCGCTGGCGCCTGC  >  NC_002947/2329618‑2329878

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