BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	302,598	C→T	32.3%	S388S (TCC→TCT)	amtB →	ammonium transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	302,598	0	C	T	32.3%	61.0 / 34.1	71	S388S (TCC→TCT)	amtB	ammonium transporter
Reads supporting (aligned to +/- strand): ref base C (29/19); new base T (18/5); total (47/24)
Fisher's exact test for biased strand distribution p-value = 1.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

GTGTGACGATGGGCCATCAGTTGCTGGTACAGCTGGAAAGCATCGCCATTACGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCGGATCTGACGGTTGGTCTGCGTGTACCGGAAGAGC  >  minE/302534‑302665

gtgtGACGATGGGCCATCAGTTGCTGGTACAGCTGGAAAGCATCGCCATTACGATTGTCTGGTCTGgtg                                                                 <  1:3081424/69‑1 (MQ=255)

gtgtGACGATGGGCCATCAGTTGCTGGTACAGCTGGAAAGCATCGCCATTACGATTGTCTGGTCTGgtg                                                                 <  1:1339572/69‑1 (MQ=255)

gtgtGACGATGGGCCATCAGTTGCTGGTACAGCTGGAAAGCATCGCCATTACGATTGTCTGGTCTGgtg                                                                 <  1:1650577/69‑1 (MQ=255)

gtgtGACGATGGGCCATCAGTTGCTGGTACAGCTGGAAAGCATCGCCATTACGATTGTCTGGTCTGgtg                                                                 <  1:2358798/69‑1 (MQ=255)

gtgtGACGATGGGCCATCAGTTGCTGGTACAGCTGGAAAGCATCGCCATTACGATTGTCTGGTCTGgtg                                                                 <  1:1303046/69‑1 (MQ=255)

 tgtgACGATGGGCCATCAGTTGCTGGTACAGCTGGAAAGCATCGCCATTACGATCGTCTGGTCCGgtgt                                                                <  1:1449256/69‑1 (MQ=255)

 tgtgACGATGGGCCATCAGTTGCTGGTACAGCTGGAAAGCATCGCCATTACGATCGTCTGGTCCGgtgt                                                                <  1:1526388/69‑1 (MQ=255)

 tgtgACGATGGGCCATCAGTTGCTGGTACAGCTGGAAAGCATCGCCATTACGATCGTCTGGTCCGgtgt                                                                <  1:582391/69‑1 (MQ=255)

 tgtgACGATGGGCCATCAGTTGCTGGTACAGCTGGAAAGCATCGCCATTACGATCGTCTGGTCCGgtgt                                                                <  1:787659/69‑1 (MQ=255)

 tgtgACGATGGGCCATCAGTTGCTGGTACAGCTGGAAAGCATCGCCATTACGATCGTCTGGTCCGgtgt                                                                <  1:3246505/69‑1 (MQ=255)

 tgtgACGATGGGCCATCAGTTGCTGGTACAGCTGGAAAGCATCGCCATTACGATCGTCTGGTCCGgtgt                                                                <  1:1074616/69‑1 (MQ=255)

 tgtgACGATGGGCCATCAGTTGCTGGTACAGCTGGAAAGCATCGCCATTACGATCGTCTGGTCCGgtgt                                                                <  1:323760/69‑1 (MQ=255)

  gtgACGATGGGCCATCAGTTGCTGGTACAGCTGGAAAGCATCGCCATTACGATCGTCTGGTCCGgtgt                                                                <  1:1321071/68‑1 (MQ=255)

                      gCTGGTACAGCTGGAAAGCATCGCCATTACGATCGTCTGGTCCGgt                                                                  >  1:2447874/1‑46 (MQ=255)

                      gCTGGTACAGCTGGAAAGCATCGCCATTACGATCGTCTGGTCCGgt                                                                  >  1:2874436/1‑46 (MQ=255)

                      gCTGGTACAGCTGGAAAGCATCGCCATTACGATCGTCTGGTCCGgt                                                                  >  1:596862/1‑46 (MQ=255)

                      gCTGGTACAGCTGGAAAGCATCGCCATTACGATCGTCTGGTCCGgt                                                                  >  1:921907/1‑46 (MQ=255)

                             cAGCTGGAAAGCATCGCCATTACGATTGTCTGGTCTGgtg                                                                 >  1:754222/1‑40 (MQ=38)

                             cAGCTGGAAAGCATCGCCATTACGATTGTCTGGTCTGgtg                                                                 >  1:1055988/1‑40 (MQ=38)

                             cAGCTGGAAAGCATCGCCATTACGATTGTCTGGTCTGgtg                                                                 >  1:2311200/1‑40 (MQ=38)

                             cAGCTGGAAAGCATCGCCATTACGATTGTCTGGTCTGgtg                                                                 >  1:368318/1‑40 (MQ=38)

                             cAGCTGGAAAGCATCGCCATTACGATTGTCTGGTCTGgtg                                                                 >  1:2580461/1‑40 (MQ=38)

                             cAGCTGGAAAGCATCGCCATTACGATTGTCTGGTCTGgtg                                                                 >  1:46963/1‑40 (MQ=38)

                             cAGCTGGAAAGCATCGCCATTACGATTGTCTGGTCTGgtg                                                                 >  1:1929373/1‑40 (MQ=38)

                             cAGCTGGAAAGCATCGCCATTACGATTGTCTGGTCTGgtg                                                                 >  1:476084/1‑40 (MQ=38)

                             cAGCTGGAAAGCATCGCCATTACGATTGTCTGGTCTGgtg                                                                 >  1:2669002/1‑40 (MQ=38)

                             cAGCTGGAAAGCATCGCCATTACGATTGTCTGGTCTGgtg                                                                 >  1:557340/1‑40 (MQ=38)

                             cAGCTGGAAAGCATCGCCATTACGATTGTCTGGTCTGgtg                                                                 >  1:1682626/1‑40 (MQ=38)

                             cAGCTGGAAAGCATCGCCATTACGATTGTCTGGTCTGgtg                                                                 >  1:1671485/1‑40 (MQ=38)

                             cAGCTGGAAAGCATCGCCATTACGATTGTCTGGTCTGgtg                                                                 >  1:3169240/1‑40 (MQ=38)

                             cAGCTGGAAAGCATCGCCATTACGATTGTCTGGTCTGgtg                                                                 >  1:2790570/1‑40 (MQ=38)

                             cAGCTGGAAAGCATCGCCATTACGATTGTCTGGTCTGgtg                                                                 >  1:832942/1‑40 (MQ=38)

                             cAGCTGGAAAGCATCGCCATTACGATTGTCTGGTCTGgtg                                                                 >  1:2967083/1‑40 (MQ=38)

                             cAGCTGGAAAGCATCGCCATTACGATTGTCTGGTCTGgtg                                                                 >  1:1502386/1‑40 (MQ=38)

                             cAGCTGGAAAGCATCGCCATTACGATTGTCTGGTCTGgtg                                                                 >  1:3038886/1‑40 (MQ=38)

                             cAGCTGGAAAGCATAGCCATTACGATTGTctggtctggt                                                                  >  1:1403793/1‑39 (MQ=37)

                               gCTGGAAAGCATCGCCATTACGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGctg                                   <  1:2958934/68‑3 (MQ=255)

                               gCTGGAAAGCATCGCCATTACGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCgg                                   <  1:3431900/68‑1 (MQ=255)

                               gCTGGAAAGCATCGCCATTACGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCgg                                   <  1:1004208/68‑1 (MQ=255)

                               gCTGGAAAGCATCGCCATTACGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCgg                                   <  1:2122950/68‑1 (MQ=255)

                               gCTGGAAAGCATCGCCATTACGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCgg                                   <  1:209595/68‑1 (MQ=255)

                               gCTGGAAAGCATCGCCATTACGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCgg                                   <  1:792873/68‑1 (MQ=255)

                               gCTGGAAAGCATCGCCATTACGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCgg                                   <  1:1807789/68‑1 (MQ=255)

                               gCTGGAAAGCATCGCCATTACGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCgg                                   <  1:1375692/68‑1 (MQ=255)

                               gCTGGAAAGCATCGCCATTACGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCgg                                   <  1:284418/68‑1 (MQ=255)

                               gCTGGAAAGCATCGCCATTACGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCgg                                   <  1:3287970/68‑1 (MQ=255)

                               gCTGGAAAGCATCGCCATTACGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCgg                                   <  1:1391430/68‑1 (MQ=255)

                                          tCGCCATTACGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCGGATCTGACGGTTg                       >  1:1545005/1‑69 (MQ=255)

                                                  aCGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCGGTTCTGACGGTTGGTCt                   >  1:1313426/1‑65 (MQ=255)

                                                  aCGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCGGATCTGAc                            >  1:2035163/1‑56 (MQ=255)

                                                  aCGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCGGATCTGACGGTTGGTCTGc                 >  1:2924543/1‑67 (MQ=255)

                                                  aCGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCGGATCTGACGGTTGGTCTGCg                >  1:2221304/1‑68 (MQ=255)

                                                  aCGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCGGATCTGACGGTTGGTCTGCg                >  1:898591/1‑68 (MQ=255)

                                                  aCGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCGGATCTGACGGTTGGTCTGCg                >  1:329631/1‑68 (MQ=255)

                                                  aCGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCGGATCTGACGGTTGGTCTGCg                >  1:688879/1‑68 (MQ=255)

                                                  aCGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCGGATCTGACGGTTGGTCTGCg                >  1:648458/1‑68 (MQ=255)

                                                  aCGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCGGATCTGACGGTTGGTCTGCg                >  1:170013/1‑68 (MQ=255)

                                                  aCGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCGGATCTGACGGTTGGTCTGCg                >  1:1032201/1‑68 (MQ=255)

                                                  aCGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCGGATCTGACGGTTGGTCTGCg                >  1:1128305/1‑68 (MQ=255)

                                                  aCGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCGGATCTGACGGTTGGTCTGCg                >  1:1134065/1‑68 (MQ=255)

                                                  aCGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCGGATCTGACGGTTGGTCTGCg                >  1:1170509/1‑68 (MQ=255)

                                                  aCGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCGGATCTGACGGTTGGTCTGCg                >  1:1322657/1‑68 (MQ=255)

                                                  aCGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCGGATCTGACGGTTGGTCTGCg                >  1:1328112/1‑68 (MQ=255)

                                                  aCGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCGGATCTGACGGTTGGTCTGCg                >  1:1458965/1‑68 (MQ=255)

                                                  aCGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCGGATCTGACGGTTGGTCTGCg                >  1:1751425/1‑68 (MQ=255)

                                                  aCGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCGGATCTGACGGTTGGTCTGCg                >  1:1912690/1‑68 (MQ=255)

                                                  aCGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCGGATCTGACGGTTGGTCTGCg                >  1:1979283/1‑68 (MQ=255)

                                                  aCGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCGGATCTGACGGTTGGTCTGCg                >  1:2165020/1‑68 (MQ=255)

                                                  aCGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCGGATCTGACGGTTGGTCTGCg                >  1:2293510/1‑68 (MQ=255)

                                                  aCGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCGGATCTGACGGTTGGTCTGCg                >  1:2879638/1‑68 (MQ=255)

                                                  aCGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCGGATCTGACGGTTGGTCTGCg                >  1:3048086/1‑68 (MQ=255)

                                                  aCGATCGTCTGGTCCGCTGTTGTGGCATTTATCGGCTACAAATTGGCGGATCTGACGGTTGGTCTGCg                >  1:111724/1‑68 (MQ=255)

                                                               ccGGTGTTGTGGCATTTATCGGCTACAAATTGGCGGATCTGACGGTTGGTCTGCGTGTACCGGAagag   <  1:2064689/68‑1 (MQ=255)

                                                               ccGGTGTTGTGGCATTTATCGGCTACAAATTGGCGGATCTGACGGTTGGTCTGCGTGTACCGGAAGAGc  <  1:3233119/69‑1 (MQ=255)

                                                               ccGGTGTTGTGGCATTTATCGGCTACAAATTGGCGGATCTGACGGTTGGTCTGCGTGTACCGGAAGAGc  <  1:159561/69‑1 (MQ=255)

GTGTGACGATGGGCCATCAGTTGCTGGTACAGCTGGAAAGCATCGCCATTACGATCGTCTGGTCCGGTGTTGTGGCATTTATCGGCTACAAATTGGCGGATCTGACGGTTGGTCTGCGTGTACCGGAAGAGC  >  minE/302534‑302665

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