BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,253,560	A→G	69.4%	T293T (ACT→ACC)	aspS ←	aspartyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,253,560	0	A	G	69.4%	46.3 / 21.7	52	T293T (ACT→ACC)	aspS	aspartyl‑tRNA synthetase
Reads supporting (aligned to +/- strand): ref base A (10/6); new base G (19/17); total (29/23)
Fisher's exact test for biased strand distribution p-value = 5.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.99e-01

TACAGCAAACTCAACAGATTTCAGCAGATCAGCAACGTCAGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCGTTCTGCTTCCGCA  >  minE/1253521‑1253623

tACAGCAAACTCAACAGATTTCAGCAGATCAGCAACGTCAGTCAGTTCCATCGGGTTACGCAGATCCgg                                    <  1:163092/69‑1 (MQ=255)

tACAGCAAACTCAACAGATTTCAGCAGATCAGCAACGTCAGTCAGTTCCATCGGGTTACGCAGATCCgg                                    <  1:1676405/69‑1 (MQ=255)

 aCAGCAAACTCAACAGATTTCAGCAGATCAGCAACGTCAGTCAGTTCCATCGGGTTACGCAGATCCgg                                    >  1:973646/1‑68 (MQ=255)

 aCAGCAAACTCAACAGATTTCAGCAGATCAG‑AACGTCAGTCAGTTCCATCGGGTTACGCAGATCCGGtt                                  >  1:1007453/1‑69 (MQ=255)

       aaCTCAACAGATTTCAGCAGATCAGCAACGTCAGTCAGTTcc                                                        >  1:1380568/1‑42 (MQ=255)

       aaCTCAACAGATTTCAGCAGATCAGCAACGTCAGTCAGTTCCATCGGGTTACGCAGATCCGGTTTAta                              >  1:1422909/1‑67 (MQ=255)

       aaCTCAACAGATTTCAGCAGATCAGCAACGTCAGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATc                              >  1:945861/1‑68 (MQ=255)

       aaCTCAACAGATTTCAGCAGATCAGCAACGTCAGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATc                              >  1:418632/1‑68 (MQ=255)

        aCTCAACAGATTTCAGCAGATCAGCAACGTCAGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAg                            <  1:792072/69‑1 (MQ=255)

        aCTCAACAGATTTCAGCAGATCAGCAACGTCAGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAg                            <  1:726609/69‑1 (MQ=255)

        aCTCAACAGATTTCAGCAGATCAGCAACGTCAGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAg                            <  1:1917024/69‑1 (MQ=255)

                      agcagATCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg               <  1:2316198/68‑1 (MQ=255)

                      agcagATCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg               <  1:1025237/68‑1 (MQ=255)

                      agcagATCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg               <  1:2375964/68‑1 (MQ=255)

                      agcagATCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg               <  1:2828635/68‑1 (MQ=255)

                      agcagATCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg               <  1:2954045/68‑1 (MQ=255)

                      agcagATCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg               <  1:297186/68‑1 (MQ=255)

                      agcagATCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg               <  1:3182686/68‑1 (MQ=255)

                      agcagATCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg               <  1:3414039/68‑1 (MQ=255)

                      agcagATCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg               <  1:460090/68‑1 (MQ=255)

                      agcagATCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg               <  1:1575565/68‑1 (MQ=255)

                      agcagATCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg               <  1:1176772/68‑1 (MQ=255)

                      agcagATCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg               <  1:1851042/68‑1 (MQ=255)

                       gcagATCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg               <  1:3311873/67‑1 (MQ=255)

                       gcagATCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg               <  1:788737/67‑1 (MQ=255)

                       gcagATCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg               <  1:1957261/67‑1 (MQ=255)

                       gcagATCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAAACATAACGGCg               <  1:692211/67‑1 (MQ=255)

                       gcagATCGGCAACGTCGGTCAGTTCCATCAGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg               <  1:243877/67‑1 (MQ=255)

                           aTCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGc                >  1:1662891/1‑62 (MQ=255)

                           aTCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg               >  1:625734/1‑63 (MQ=255)

                           aTCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg               >  1:3396170/1‑63 (MQ=255)

                           aTCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg               >  1:1510916/1‑63 (MQ=255)

                           aTCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg               >  1:1461619/1‑63 (MQ=255)

                           aTCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg               >  1:1366377/1‑63 (MQ=255)

                           aTCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg               >  1:571514/1‑63 (MQ=255)

                           aTCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg               >  1:1081628/1‑63 (MQ=255)

                           aTCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg               >  1:129976/1‑63 (MQ=255)

                           aTCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg               >  1:1165223/1‑63 (MQ=255)

                           aTCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg               >  1:2101798/1‑63 (MQ=255)

                           aTCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg               >  1:3144496/1‑63 (MQ=255)

                           aTCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg               >  1:1698047/1‑63 (MQ=255)

                           aTCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg               >  1:2792467/1‑63 (MQ=255)

                           aTCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg               >  1:25941/1‑63 (MQ=255)

                           aTCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg               >  1:256278/1‑63 (MQ=255)

                           aTCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg               >  1:2525754/1‑63 (MQ=255)

                           aTCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg               >  1:2232308/1‑63 (MQ=255)

                           aTCGGCAACGTCGGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCg               >  1:2172140/1‑63 (MQ=255)

                             cAGCAACGTCAGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCGTTCTGCtt       >  1:2648270/1‑69 (MQ=255)

                             cAGCAACGTCAGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCGTTCTGCtt       >  1:2432146/1‑69 (MQ=255)

                             cAGCAACGTCAGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCGTTCTGCtt       >  1:2003440/1‑69 (MQ=255)

                             cAGCAACGTCAGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCGTTCTGCtt       >  1:2176278/1‑69 (MQ=255)

                                   cGTCAGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCGTTCTGCTTACGCa  <  1:2392341/68‑1 (MQ=255)

TACAGCAAACTCAACAGATTTCAGCAGATCAGCAACGTCAGTCAGTTCCATCGGGTTACGCAGATCCGGTTTATCAGAACCATAACGGCGTTCTGCTTCCGCA  >  minE/1253521‑1253623

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