BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,725,155	C→T	47.5%	T5T (ACG→ACA)	yheT ←	predicted hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,725,155	0	C	T	47.5%	‑0.6 / 31.9	40	T5T (ACG→ACA)	yheT	predicted hydrolase
Reads supporting (aligned to +/- strand): ref base C (14/7); new base T (12/7); total (26/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

AATTGCTAAAGCCGCGCATGGGGATGAATTCAGCACTGCTGCTGAATTCATTGGCATCGGTCGTCGTTATCTGCGCCATCAGTTGCTTTGGCCTTCCAGCAGCATCTGCTCAAG  >  minE/2725091‑2725204

aaTTGCTAAAGCCGCGCATGGGGATGAATTCAGCACTGCTGCTGAATTCATTGGCATCGGTCGTCGTTa                                               <  1:3143714/69‑1 (MQ=255)

aaTTGCTAAAGCCGCGCATGGGGATGAATTCAGCACTGCTGCTGAATTCATTGGCATCGGTCGTCGTTa                                               <  1:2889350/69‑1 (MQ=255)

aaTTGCTAAAGCCGCGCATGGGGATGAATTCAGCACTGCTGCTGAATTCATTGGCATCGGTCGTCGTTa                                               <  1:2473940/69‑1 (MQ=255)

 aTTGCTAAAGCCGCGCATGGGGATGAATTCAGCACTGCTTCTGAATTCATTGGCATCGGTCGTCGTTa                                               <  1:2682233/68‑1 (MQ=255)

  ttGCTAAAGCCGCGCATAGGGGTGAATTCAGCACTGTTGCTAAATTCATTGGCATCGGTCGTTGTTATc                                             >  1:3133944/1‑69 (MQ=255)

  ttGCTAAAGCCGCGCATAGGGGTGAATTCAGCACTGCTGCTAAATTCATTGGCATCGGTCGTTGTTAt                                              >  1:658477/1‑68 (MQ=255)

  ttGCTAAAGCCGCGCATAGGGGTGAATTCAGCACTGCTGCTAAATTCATTGGCATCGGTCGTTGTTAt                                              >  1:2701984/1‑68 (MQ=255)

  ttGCTAAAGCCGCGCATAGGGGTGAATTCAGCACTGCTGCTAAATTCATTGGCATCGGTCGTTGTTAt                                              >  1:258930/1‑68 (MQ=255)

  ttGCTAAAGCCGCGCATAGGGGTGAATTCAGCACTGCTGCTAAATTCATTGGCATCGGTCGTTGTTATc                                             >  1:2808831/1‑69 (MQ=255)

  ttGCTAAAGCCGCGCATAGGGGTGAATTCAGCACTGCTGCTAAATTCATTGGCATCGGTCGTTGTTATc                                             >  1:997100/1‑69 (MQ=255)

  ttGCTAAAGCCGCGCATAGGGGTGAATTCAGCACTGCTGCTAAATTCATTGGCATCGGTCGTTGTTATc                                             >  1:410223/1‑69 (MQ=255)

  ttGCTAAAGCCGCGCATAGGGGTGAATTCAGCACTGCTGCTAAATTCATTGGCATCGGTCGTTGTTATc                                             >  1:2375973/1‑69 (MQ=255)

  ttGCTAAAGCCGCGCATAGGGGTGAATTCAGCACTGCTGCTAAATTCATTGGCATCGGTCGTTGTTATc                                             >  1:2072966/1‑69 (MQ=255)

  ttGCTAAAGCCGCGCATAGGGGTGAATTCAGCACTGCTGCTAAATTCATTGGCATCGGTCGTTGTTATc                                             >  1:1113169/1‑69 (MQ=255)

  ttGCTAAAGCCGCGCATAGGGGTGAATTCAGCACTGCTGCTAAATTCATTGGCATCGGTCGTTGTTATc                                             >  1:1125271/1‑69 (MQ=255)

  ttGCTAAAGCCGCGCATAGGGGTGAATTCAGCACTGCTGCTAAATTCATTGGCATCGGTCGTTGTTATc                                             >  1:1322719/1‑69 (MQ=255)

       aaaGCCGCGCATGGGGATGAATTCAGCACTGCTGCTGAATTCATTGGCATCGgtcgtc                                                   <  1:1747542/58‑1 (MQ=255)

       aaaGCCGCGCATGGGGATGAATTCAGCACTGCTGCTGAATTCATTGGCATCGgtcgtc                                                   <  1:2930899/58‑1 (MQ=255)

       aaaGCCGCGCATGGGGATGAATTCAGCACTGCTGCTGAATTCATTGGCATCGgtcgtc                                                   <  1:698384/58‑1 (MQ=255)

                   aggggTGAATTCAGCACTGCTGCTAAATTCATTGGCATCGGTCGTTGTTATCTGCGCCATCAGTTGCtt                            <  1:3278466/68‑1 (MQ=255)

                   aggggTGAATTCAGCACTGCTGCTAAATTCATTGGCATCGGTCGTTGTTATCTGCGCCATCAGTTGCtt                            <  1:2802760/68‑1 (MQ=255)

                   aggggTGAATTCAGCACTGCTGCTAAATTCATTGGCATCGGTCGTTGTTATCTGCGCCATCAGTTGCtt                            <  1:793685/68‑1 (MQ=255)

                   aggggTGAATTCAGCACTGCTGCTAAATTCATTGGCATCGGTCGTTGTTATCTGCGCCATCAGTTGCtt                            <  1:2881956/68‑1 (MQ=255)

                   aggggTGAATTCAGCACTGCTGCTAAATTCATTGGCATCGGTCGTTGTTATCTGCGCCATCAGTTGCtt                            <  1:2883217/68‑1 (MQ=255)

                   aggggTGAATTCAGCACTGCTGCTAAATTCATTGGCATCGGTCGTTGTTATCTGCGCCATCAGTTGCtt                            <  1:1584636/68‑1 (MQ=255)

                   aggggTGAATTCAGCACTGCTGCTAAATTCATTGGCATCGGTCGTTGTTATCTGCGCCATCAGTTGCtt                            <  1:3120930/68‑1 (MQ=255)

                    gggATGAATTCAGCACTGCTGCTGAATTCATTGGCATCGGTCGTCGTTATCTGCGCCATCAGTTGCttt                           >  1:1575934/1‑69 (MQ=255)

                         gAATTCAGCACTGCTGCTGAATTCATTGGCATCGgtcgtc                                                   <  1:928271/40‑1 (MQ=255)

                            ttCAGCACTGCTGCTGAATTCATTGGCATCGGTCGTCGTTATCTGCGCCATCAGTTGCTTTGGCCTTcc                   <  1:524757/69‑1 (MQ=255)

                            ttCAGCACTGCTGCTGAATTCATTGGCATCGGTCGTCGTTATCTGCGCCATCAGTTGCTTTGGCCTTcc                   <  1:2051710/69‑1 (MQ=255)

                             tCAGCACTGCTGCTGAATTCATTGGCATCGGTCGTCGTTATCTGCGCCATCAGTTGCTTTGGCCTTcc                   <  1:454133/68‑1 (MQ=255)

                                   ctgctgctgAATTCATTGGCATCGGTCGTCGTTATCTGCGCCATCAGTTGCTTTGGCCTTCCagcagc             >  1:927394/1‑68 (MQ=255)

                                   ctgctgctgAATTCATTGGCATCGGTCGTCGTTATCTGCGCCATCAGTTGCTTTGGCCTTCCagcagc             >  1:2680902/1‑68 (MQ=255)

                                   ctgctgctgAATTCATTGGCATCGGTCGTCGTTATCTGCGCCATCAGTTGCTTTGGCCTTCCAgcagca            >  1:2135082/1‑69 (MQ=255)

                                   ctgctgctgAATTCATTGGCATCGGTCGTCGTTATCTGCGCCATCAGTTGCTTTGGCCTTCCAgcagca            >  1:3423214/1‑69 (MQ=255)

                                   ctgctgctgAATTCATTGGCATCGGTCGTCGTTATCTGCGCCATCAGTTGCTTTGGCCTTCCAgcagca            >  1:3252501/1‑69 (MQ=255)

                                   ctgctgctgAATTCATTGGCATCGGTCGTCGTTATCTGCGCCATCAGTTGCTTTGGCCTTCCAgcagca            >  1:1013976/1‑69 (MQ=255)

                                   ctgctgctgAATTCATTGGCATCGGTCGTCGTTATCTGCGCCATCAGTTGCTTTGGCCTTCCAgcagca            >  1:1642689/1‑69 (MQ=255)

                                   ctgctgctgAATTCATTGGAATCGGTCGTCGTTATCTGCGCCATCAGTTGCTTTGGCCTTCCAgcagca            >  1:1615846/1‑69 (MQ=255)

                                           gAATTCATTGGCATCGGTCGTCGTTATCTGCGCCATCAGTTGCTTTGGCCTTCCAGCAGCATCTGCTCa    >  1:339188/1‑69 (MQ=255)

                                             aTTCATTGGCATCGGTCGTCGTTATCTGCGCCATCAGTTGCTTTGGCCTTCCAGCAGCATCTGCTCaa   >  1:2487786/1‑68 (MQ=255)

                                             aTTCATTGGCATCGGTCGTCGTTATCTGCGCCATCAGTTGCTTTGGCCTTCCAGCAGCATCTGCTCAAg  >  1:1153107/1‑69 (MQ=255)

                                             aTTCATTGGCATCGGTCGTCGTTATCTGCGCCATCAGTTGCTTTGGCCTTCCAGCAGCATCTGCTCAAg  >  1:3206576/1‑69 (MQ=255)

                                             aTTCATTGGCATCGGTCGTCGTTATCTGCGCCATCAGTTGCTTTGGCCTTCCAGCAGCATCTGCTCAAg  >  1:2771612/1‑69 (MQ=255)

AATTGCTAAAGCCGCGCATGGGGATGAATTCAGCACTGCTGCTGAATTCATTGGCATCGGTCGTCGTTATCTGCGCCATCAGTTGCTTTGGCCTTCCAGCAGCATCTGCTCAAG  >  minE/2725091‑2725204

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