BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,139,855	G→T	17.4%	V56V (GTC→GTA)	PP_3638 ←	acyl‑CoA dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,139,855	0	G	T	17.4%	70.1 / 3.0	23	V56V (GTC→GTA)	PP_3638	acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand): ref base G (7/12); new base T (2/2); total (9/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCCGCGGCAGCAGGCGCTCACGCAGGGGCTGGTTGGGGCTGTGCAGCAGGTATTCGATGAAGGCACGCTGGCCCCAGCAGACGAAAGCCGCGGCCAATGAGCGGCTGGCAATGGCGGCGATGGCCTCGACCGCGTCGGTCACCTGCCCGCCCGTACCGCCCAGCGCTGGGTCGACCCCCGTGCGCAGCAGCTGCGATTCTGCGATCTGCGCCAGTACCTGCTGCGGCTCGCACAGGCCCCGGTCGATGG  >  NC_002947/4139727‑4139976

aGCCGCGGCAGCAGGCGCTCACGCAGGGGCTGGTTGGGGCTGTGCAGCAGGTATTCGATGAAGGCACGCTGGCCCCAGCAGACGAAAGCCGCGGCCAATGAGCGGCTGGCAATGGCGGCGATGGCCTCTACCGCGt                                                                                                                    >  2:485610/1‑136 (MQ=255)

aGCCGCGGCAGCAGGCGCTCACGCAGGGGCTGGTTGGGGCTGTGCAGCAGGTATTCGATGAAGGCACGCTGGCCCCAGCAGACGAAAGCCGCGGCCAATGAGCGGCTGGCAATGGCGGCGATGGCCTCTACCGCGt                                                                                                                    >  2:113178/1‑136 (MQ=255)

            agGCGCTCACGCAGGGGCTGGTTGGGGCTGTGCAGCAGGTATTCGATGAAGGCACGCTGGCCCCAGCAGACGAAAGCCGCGGCCAATGAGCGGCTGGCAATGGCGGCGATGGCCTCGACCGCGTCGGTCACCTgcc                                                                                                        <  1:394803/136‑1 (MQ=255)

               cgcTCACGCAGGGGCTGGTTGGGGCTGTGCAGCAGGTATTCGATGAAGGCACGCTGGCCCCAGCAGACGAAAGCCGCGGCCAATGAGCGGCTGGCAATGGCGGCGATGGCCTCGACCGCGTCGGTCACCTgcccgc                                                                                                     >  2:185644/1‑136 (MQ=255)

                 cTCACGCAGGGGCTGGTTGGGGCTGTGCAGCAGGTATTCGATGAAGGCACGCTGGCCCCAGCAGACGAAAGCCGCGGCCAATGAGCGGCTGGCAATGGCGGCGATGGCCTCGACCGCGTCGGTCACCTgcccgccc                                                                                                   <  2:281876/136‑1 (MQ=255)

                      gCAGGGGCTGGTTGGGGCTGTGCAGCAGGTATTCGATGAAGGCACGCTGGCCCCAGCAGACGAAAGCCGCGGCCAATGAGCGGCTGGCAATGGCGGCGATGGCCTCGACCGCGTCGGTCACCTGCCCGCCCGTAcc                                                                                              <  1:190659/136‑1 (MQ=255)

                           ggCTGGTTGGGGCTGTGCAGCAGGTATTCGATGAAGGCACGCTGGCCCCAGCAGACGAAAGCCGCGGCCAATGAGCGGCTGGCAATGGCGGCGATGGCCTCGACCGCGTCGGTCACCTGCCCGCCCGTACCGCCCa                                                                                         <  2:490939/136‑1 (MQ=255)

                             cTGGTTGGGGCTGTGCAGCAGGTATTCGATGAAGGCACGCTGGCCCCAGCAGACGAAAGCCGCGGCCAATGAGCGGCTGGCAATGGCGGCGATGGCCTCGACCGCGTCGGTCACCTGCCCGCCCGTACCGCCCAgc                                                                                       >  2:111988/1‑136 (MQ=255)

                                   ggggCTGTGCAGCAGGTATTCGATGAAGGCACGCTGGCCCCAGTAGACGAAAGCCGCGGCCAATGAGCGGCTGGCAATGGCGGCGATGGCCTCGACCGCGTCGGTCACCTGCCCGCCCGTACCGCCCAGCGCTggg                                                                                 >  2:460498/1‑136 (MQ=255)

                                   ggggCTGTGCAGCAGGTATTCGATGAAGGCACGCTGGCCCCAGCAGACGAAAGCCGCGGCCAATGAGCGGCTGGCAATGGCGGCGATGGCCTCGACCGCGTCGGTCACCTGCCCGCCCGTACCGCCCAGCGCTggg                                                                                 >  2:419736/1‑136 (MQ=255)

                                                       cGATGAAGGCACGCTGGCCCCAGCAGACGAAAGCCGCGGCCAATGAGCGGCTGGCAATGGCGGCGATGGCCTCGACCGCGTCGGTCACCTGCCCGCCCGTACCGCCCAGCGCTGGGTCGACCCCCGTGCGcagcag                                                             >  1:273608/1‑136 (MQ=255)

                                                       cGATGAAGGCACGCTGGCCCCAGCAGACGAAAGCCGCGGCCAATGAGCGGCTGGCAATGGCGGCGATGGCCTCGACCGCGTCGGTCACCTGCCCGCCCGTACCGCCCAGCGCTGGGTCGACCCCCGTGCGcagcag                                                             <  2:389766/136‑1 (MQ=255)

                                                       cGATGAAGGCACGCTGGCCCCAGCAGACGAAAGCCGCGGCCAATGAGCGGCTGGCAATGGCGGCGATGGCCTCGACCGCGTCGGTCACCTGCCCGCCCGTACCGCCCAGCGCAGGGTCGACCCCCGTGCgcagcac                                                             >  1:25067/1‑135 (MQ=255)

                                                                 aCGCTGGCCCCAGCAGACGAAAGCCGCGGCCAATGAGCGGCTGGCAATGGCGGCGATGGCCTCGACCGCGTCGGTCACCTGCCCGCCCGTACCGCCCAGCGCTGGGTCGACCCCCGTGCGCAGCAGCTGCGATTCt                                                   <  2:389360/136‑1 (MQ=255)

                                                                    cTGGCCCCAGCAGACGAAAGCCGCGGCCAATGAGCGGCTGGCAATGGCGGCGATGGCCTCGACCGCGTCGGTCACCTGCCCGCCCGTACCGCCCAGCGCTGGGTCGACCCCCGTGCGCAGCAGCTGCGATTCTGCg                                                <  2:272266/136‑1 (MQ=255)

                                                                    cTGGCCCCAGCAGACGAAAGCCGCGGCCAATGAGCGGCTGGCAATGGCGGCGATGGCCTCGACCGCGTCGGTCACCTGCCCGCCCGTACCGCCCAGCGCTGGGTCGACCCCCGTGCGCAGCAGCTGCGATTCTGCg                                                <  2:254894/136‑1 (MQ=255)

                                                                           cagcagACGAAAGCCGCGGCCAATGAGCGGCTGGCAATGGCGGCGATGGCCTCGACCGCGTCGGTCACCTGCCCGCCCGTACCGCCCAGCGCTGGGTCGACCCCCGTGCGCAGCAGCTGCGATTCTGCGATCTgcg                                         <  1:37757/136‑1 (MQ=255)

                                                                             gcagACGAAAGCCGCGGCCAATGAGCGGCTGGCAATGGCGGCGATGGCCTCGACCGCGTCGGTCACCTGCCCGCCCGTACCGCCCAGCGCTGGGTCGACCCCCGTGCGCAGCAGCTGCGATTCTGCGATCTGCGcc                                       <  1:269908/136‑1 (MQ=255)

                                                                                           cgGCCAATGAGCGGCTGGCAATGGCGGCGATGGCCTCGACCGCGTCGGTCACCTGCCCGCCCGTACCGCCCAGCGCTGGGTCGACCCCCGTGCGCAGCAGCTGCGATTCTGCGATCTGCGCCAGTACCTGCTGCgg                         >  2:100438/1‑136 (MQ=255)

                                                                                                                tGGCGGCGATGGCCTCGACCGCGTCGGTCACCTGCCCGCCCGTACCGCCCAGCGCTGGGTCGACCCCCGTGCGCAGCAGCTGCGATTCTGCGATCTGCGCCAGTACCTGCTGCGGCTCGCACAGGCCCCGGTCGAt    <  1:326327/136‑1 (MQ=255)

                                                                                                                  gcggcgATGGCCTCTACCGCGTCGGTCACCTGCCCGCCCGTACCGCCCAGCGCTGGGTCGACCCCCGTGCGCAGCAGCTGCGATTCTGCGATCTGCGCCAGTACCTGCTGCGGCTCGCACAGGCCCCGGTCGATgg  <  1:485610/136‑1 (MQ=255)

                                                                                                                  gcggcgATGGCCTCTACCGCGTCGGTCACCTGCCCGCCCGTACCGCCCAGCGCTGGGTCGACCCCCGTGCGCAGCAGCTGCGATTCTGCGATCTGCGCCAGTACCTGCTGCGGCTCGCACAGGCCCCGGTCGATgg  <  1:113178/136‑1 (MQ=255)

                                                                                                                  gcggcgATGGCCTCGACCGCGTCGGTCACCTGCCCGCCCGTACCGCCCAGCGCTGGGTCGACCCCCGTGCGCAGCAGCTGCGATTCTGCGATCTGCGCCAGTACCTGCTGCGGCTCGCACAGGCCCCGGTCGATgg  <  2:390150/136‑1 (MQ=255)

AGCCGCGGCAGCAGGCGCTCACGCAGGGGCTGGTTGGGGCTGTGCAGCAGGTATTCGATGAAGGCACGCTGGCCCCAGCAGACGAAAGCCGCGGCCAATGAGCGGCTGGCAATGGCGGCGATGGCCTCGACCGCGTCGGTCACCTGCCCGCCCGTACCGCCCAGCGCTGGGTCGACCCCCGTGCGCAGCAGCTGCGATTCTGCGATCTGCGCCAGTACCTGCTGCGGCTCGCACAGGCCCCGGTCGATGG  >  NC_002947/4139727‑4139976

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