BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,813,620	A→G	36.9%	F519F (TTT→TTC)	alaS ←	alanyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,813,620	0	A	G	36.9%	23.1 / 26.0	46	F519F (TTT→TTC)	alaS	alanyl‑tRNA synthetase
Reads supporting (aligned to +/- strand): ref base A (8/21); new base G (15/2); total (23/23)
Fisher's exact test for biased strand distribution p-value = 1.48e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAACCCGCAGCAAGTTTACCGATGTGACCAATCGCCTGGCCGTATTTCTGCGTATCTTCCACCGCAAAGGAGAAGTTAGCGCCTTTCAGTTCGCCTTTATCGCCAACCTGACCGCCGG  >  minE/1813555‑1813672

gaACCCGCAGCAAGTTTACCGATGTGACCAATCGCCTGGCCGTATTTCTGCGTATCTTCCACCGCaa                                                     >  1:3361672/1‑67 (MQ=255)

   cccGCAGCAAGTTTACCGATGTGACCAATCGCCTGGCCGTATTTCTGCGTATCTTCCACCGCAAAgg                                                  >  1:633718/1‑67 (MQ=255)

   cccGCAGCAAGTTTACCGATGTGACCAATCGCCTGGCCGTATTTCTGCGTATCTTCCACCGCAAAgg                                                  >  1:1135011/1‑67 (MQ=255)

   cccGCAGCAAGTTTACCGATGTGACCAATCGCCTGGCCGTATTTCTGCGTATCTTCCACCGCAAAGgag                                                <  1:1216172/69‑1 (MQ=255)

   cccGCAGCAAGTTTACCGATGTGACCAATCGCCTGGCCGTATTTCTGCGTATCTTCCACCGCAAAGgag                                                <  1:916702/69‑1 (MQ=255)

   cccGCAGCAAGTTTACCGATGTGACCAATCGCCTGGCCGTATTTCTGCGTATCTTCCACCGCAAAGgag                                                <  1:730688/69‑1 (MQ=255)

   cccGCAGCAAGTTTACCGATGTGACCAATCGCCTGGCCGTATTTCTGCGTATCTTCCACCGCAAAGgag                                                <  1:1713663/69‑1 (MQ=255)

   cccGCAGCAAGTTTACCGATGTGACCAATCGCCTGGCCGTATTTCTGCGTATCTTCCACCGCAAAGgag                                                <  1:1742205/69‑1 (MQ=255)

   cccGCAGCAAGTTTACCGATGTGACCAATCGCCTGGCCGTATTTCTGCGTATCTTCCACCGCAAAGgag                                                <  1:1768646/69‑1 (MQ=255)

   cccGCAGCAAGTTTACCGATGTGACCAATCGCCTGGCCGTATTTCTGCGTATCTTCCACCGCAAAGgag                                                <  1:2900569/69‑1 (MQ=255)

   cccGCAGCAAGTTTACCGATGTGACCAATCGCCTGGCCGTATTTCTGCGTATCTTCCACCGCAAAGgag                                                <  1:2707834/69‑1 (MQ=255)

    ccGCAGCAAGTTTACCGATGTGACCAATCGCCTGGCCGTATTTCTGCGTATCTTCCACCGCAAAGgag                                                <  1:444823/68‑1 (MQ=255)

     cGCAGCAAGTTTACCGATGTGACCAATCGCCTGGCCGTATTTCTGCGTATCTTCCACCGCAAAGgag                                                <  1:1970796/67‑1 (MQ=255)

                             aaTCGCCTGGCCGTATTTCTGCGTATCTTCCACCGCAAAGGAGAAGTTAGCGCCTTTCAGTTCGCCttt                      <  1:2820880/69‑1 (MQ=255)

                             aaTCGCCTGGCCGTATTTCTGCGTATCTTCCACCGCAAAGGAGAAGTTAGCGCCTTTCAGTTCGCCttt                      <  1:3182514/69‑1 (MQ=255)

                             aaTCGCCTGGCCGTATTTCTGCGTATCTTCCACCGCAAAGGAGAAGTTAGCGCCTTTCAGTTCGCCttt                      <  1:2250267/69‑1 (MQ=255)

                             aaTCGCCTGGCCGTATTTCTGCGTATCTTCCACCGCAAAGGAGAAGTTAGCGCCTTTCAGTTCGCCttt                      <  1:2025791/69‑1 (MQ=255)

                             aaTCGCCTGGCCGTATTTCTGCGTATCTTCCACCGCAAAGGAGAAGTTAGCGCCTTTCAGTTCGCCttt                      <  1:1258274/69‑1 (MQ=255)

                             aaTCGCCTGGCCGTATTTCTGCGTATCTTCCACCGCAAAGGAGAAGTTAGCGCCTTTCAGTTCGCCttt                      <  1:183344/69‑1 (MQ=255)

                             aaTCGCCTGGCCGTATTTCTGCGTATCTTCCACCGCAAAGGAGAAGTTAGCGCCTTTCAGTTCGCCttt                      <  1:1824423/69‑1 (MQ=255)

                             aaTCGCATGGCCGTATTTCTGCGTATCTTCCACCGCAAAGGAGAAGTTAGCGCCTTTCAGTTCGCCttt                      <  1:848604/69‑1 (MQ=255)

                                        cGTATTTCTGCGTATCTTGCACCGCAAAGGAGAAGTTAGCTCGTTTCAGTTCGCCTTTATCGCCAACCt           <  1:635252/69‑1 (MQ=255)

                                        cGTATTTCTGCGTATCTTCCACCGCAAAGGAGAAGTTAGCGCCTTTCAGTTCGCCTTTATCGCCAACCt           <  1:2393649/69‑1 (MQ=255)

                                         gTATTTCTGCGTATCTTCCACCGCGAAGGAGAAGTTAGCGCCTTTCAGTTCGCCTTTATCACCAACCTg          <  1:2484706/69‑1 (MQ=255)

                                         gTATTTCTGCGTATCTTCCACCGCGAAGGAGAAGTTAGCGCCTTTCAGTTCGCCTTTATCACCAACCTg          <  1:3435713/69‑1 (MQ=255)

                                         gTATTTCTGCGTATCTTCCACCGCAAAGGAGAAGTTAGCGCCTTTCAGTTCGCCTTTATCGCCAACCt           <  1:1870591/68‑1 (MQ=255)

                                            tttCTGCGTATCTTCCACCGCAAAGGAGAAGTTAGCGCCTTTCAGTTCGCCTTTATCGCCACCCTGAcc       >  1:13416/1‑69 (MQ=255)

                                            tttCTGCGTATCTTCCACCGCAAAGGAGAAGTTAGCGCCTTTCAGTTCGCCTTTATCGCCAACCTGAcc       >  1:3107227/1‑69 (MQ=255)

                                            tttCTGCGTATCTTCCACCGCAAAGGAGAAGTTAGCGCCTTTCAGTTCGCCTTTATCGCCAACCTGAcc       >  1:3313758/1‑69 (MQ=255)

                                              tCTGCGTATCTTCCACCGCAAAGGAGAAGTTAGCGCCTTTCAGTTc                            >  1:684099/1‑46 (MQ=255)

                                              tCTGCGTATCTTCCACCGCAAAGGAGAAGTTAGCGCCTTTCAGTTCGc                          >  1:975394/1‑48 (MQ=255)

                                              tCTGCGTATCTTCCACCGCAAAGGAGAAGTTAGCGCCTTTCAGTTCGc                          >  1:1802254/1‑48 (MQ=255)

                                                 gCGTATCTTCCACCGCGAAGGAGAAGTTAGCGCCTTTc                                 >  1:621710/1‑38 (MQ=255)

                                                 gCGTATCTTCCACCGCGAAGGAGAAGTTAGCGCCTTTCAg                               >  1:884970/1‑40 (MQ=255)

                                                 gCGTATCTTCCACCGCGAAGGAGAAGTTAGCGCCTTTCAGTTCGCCTTTATCACCAACCTGAccgccg   >  1:3334548/1‑68 (MQ=255)

                                                 gCGTATCTTCCACCGCGAAGGAGAAGTTAGCGCCTTTCAGTTCGCCTTTATCACCAACCTGACCGCCgg  >  1:3420148/1‑69 (MQ=255)

                                                 gCGTATCTTCCACCGCGAAGGAGAAGTTAGCGCCTTTCAGTTCGCCTTTATCACCAACCTGACCGCCgg  >  1:3434467/1‑69 (MQ=255)

                                                 gCGTATCTTCCACCGCGAAGGAGAAGTTAGCGCCTTTCAGTTCGCCTTTATCACCAACCTGACCGCCgg  >  1:1630804/1‑69 (MQ=255)

                                                 gCGTATCTTCCACCGCGAAGGAGAAGTTAGCGCCTTTCAGTTCGCCTTTATCACCAACCTGACCGCCgg  >  1:3214685/1‑69 (MQ=255)

                                                 gCGTATCTTCCACCGCGAAGGAGAAGTTAGCGCCTTTCAGTTCGCCTTTATCACCAACCTGACCGCCgg  >  1:2829757/1‑69 (MQ=255)

                                                 gCGTATCTTCCACCGCGAAGGAGAAGTTAGCGCCTTTCAGTTCGCCTTTATCACCAACCTGACCGCCgg  >  1:271304/1‑69 (MQ=255)

                                                 gCGTATCTTCCACCGCGAAGGAGAAGTTAGCGCCTTTCAGTTCGCCTTTATCACCAACCTGACCGCCgg  >  1:1317083/1‑69 (MQ=255)

                                                 gCGTATCTTCCACCGCGAAGGAGAAGTTAGCGCCTTTCAGTTCGCCTTTATCACCAACCTGACCGCCgg  >  1:1895460/1‑69 (MQ=255)

                                                 gCGTATCTTCCACCGCGAAGGAGAAGTTAGCGCCTTTCAGTTCGCCTTTATCACCAACCTGACCGCCgg  >  1:884503/1‑69 (MQ=255)

                                                 gCGTATCTTCCACCGCGAAGGAGAAGTTAGCGCCTTTCAGTTCGCCTTTATCACCAACCTGACCGCCgg  >  1:2487088/1‑69 (MQ=255)

                                                 gCGTATCTTCCACCGCGAAGGAGAAGTTAGCGCCTTTCAGTTCGCCTTTATCACCAACCTGACCGCCgg  >  1:2465006/1‑69 (MQ=255)

                                                 gCGTATCTTCCACCGCGAAGGAGAAG‑TAGCGCCTTTCAGTTCGCCTTTATCACCAACCTGACCGCCgg  >  1:685938/1‑68 (MQ=255)

GAACCCGCAGCAAGTTTACCGATGTGACCAATCGCCTGGCCGTATTTCTGCGTATCTTCCACCGCAAAGGAGAAGTTAGCGCCTTTCAGTTCGCCTTTATCGCCAACCTGACCGCCGG  >  minE/1813555‑1813672

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