BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,787,739	T→C	33.9%	P231P (CCA→CCG)	lysC ←	aspartokinase III

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,787,739	0	T	C	33.9%	38.5 / 28.7	56	P231P (CCA→CCG)	lysC	aspartokinase III
Reads supporting (aligned to +/- strand): ref base T (9/28); new base C (13/6); total (22/34)
Fisher's exact test for biased strand distribution p-value = 3.19e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

AGTTGCCATCTCTGCCGCTTCGGCAAACGCGATTTCATCAATGCGTTTTGCTGCGGAAACTACGCGTGGATCGGTGGTGTAGATGCCCGGGACGTCGGTCCAGATATCAACACGAGATGCGTGTAAA  >  minE/2787673‑2787799

aGTTGCCATCTCTGCCGCTTCGGCAAACGCGATTTCATCAATGCGTTTTGCTGCGGAAACTACGCGTgg                                                            <  1:2424083/69‑1 (MQ=255)

aGTTGCCATCTCTGCCGCTTCGGCAAACGCGATTTCATCAATGCGTTTTGCTGCGGAAACTACGCGTgg                                                            <  1:1735098/69‑1 (MQ=255)

aGTTGCCATCTCTGCCGCTTCGGCAAACGCGATTTCATCAATGCGTTTTGCTGCGGAAACTACGCGTgg                                                            <  1:1195522/69‑1 (MQ=255)

aGTTGCCATCTCTGCCGCTTCGGCAAACGCGATTTCATCAATGCGTTTTGCTGCGGAAACTACGCGTgg                                                            <  1:1115992/69‑1 (MQ=255)

aGTTGCCATCTCTGCCGCTTCGGCAAACGCGATTTCATCAATGCGTTTTGCTGCGGAAACTACGCGTgg                                                            <  1:2849807/69‑1 (MQ=255)

aGTTGCCATCTCTGCCGCTTCGGCAAACGCGATTTCATCAATGCGTTTTGCTGCGGAAACTACGCGTgg                                                            <  1:3384765/69‑1 (MQ=255)

aGTTGCCATCTCTGCCGCTTCGGCAAACGCGATTTCATCAATGCGTTTTGCTGCGGAAACTACGCGTgg                                                            <  1:3050048/69‑1 (MQ=255)

aGTTGCCATCTCTGCCGCTTCGGCAAACGCGATTTCATCAATGCGTTTTGCTGCGGAAACTACGCGTgg                                                            <  1:125831/69‑1 (MQ=255)

aGTTGCCATCTCTGCCGCTTCGGCAAACGCGATTTCATCAATGCGTTTTGCTGCGGAAACTACGCGTgg                                                            <  1:3206092/69‑1 (MQ=255)

 gTTGCCATCTCTGCCGCTTCGGCAAACGCGATTTCATCAATGAGTTTTGCTGCGGAAACTACGCGTgg                                                            <  1:1302831/68‑1 (MQ=255)

          tctGCCGCTTCGGCAAACGCGATTTCATCAATGCGTTTTGCTGCGGAAACTACGCGTGGATCggtgg                                                    >  1:1976767/1‑67 (MQ=255)

          tctGCCGCTTCGGCAAACGCGATTTCATCAATGCGTTTTGCTGCGGAAACTACGAGTGGATCggtgg                                                    >  1:1390898/1‑67 (MQ=255)

              ccGCTTCGGCAAACGCGATTTCATCAATGCGTTTTGCTGCGGAAACTACGCGTGGATCGGTGGTGTa                                                <  1:3123923/67‑1 (MQ=255)

              ccGCTTCGGCAAACGCGATTTCATCAATGCGTTTTGCTGCGGAAACTACGCGTGGATCGGTGGTGTa                                                <  1:3261413/67‑1 (MQ=255)

              ccGCTTCGGCAAACGCGATTTCATCAATGCGTTTTGCTGCGGAAACTACGCGTGGATCGGTGGTGTa                                                <  1:705533/67‑1 (MQ=255)

              ccGCTTCGGCAAACGCGATTTCATCAATGCGTTTTGCTGCGGAAACTACGCGTGGATCGGTGGTGTa                                                <  1:632596/67‑1 (MQ=255)

              ccGCTTCGGCAAACGCGATTTCATCAATGCGTTTTGCTGCGGAAACTACGCGTGGATCGGTGGTGTa                                                <  1:1171939/67‑1 (MQ=255)

               cGCTTCGGCAAACGCGATTTCATCAATGCGTTTTGCTGCGGAAACTACGCGTGGATCGGTGGTGTa                                                <  1:714317/66‑1 (MQ=255)

                      gCAAACGCGATTTCATCAATGCGTTTTGCTGCGGAAACTACGCGTGGATCGGTGGTGTAGATGCCCgg                                       <  1:2202198/68‑1 (MQ=255)

                      gCAAACGCGATTTCATCAATGCGTTTTGCTGCGGAAACTACGCGTGGATCGGTGGTGTAGATGCCCgg                                       <  1:2433420/68‑1 (MQ=255)

                           cgcgATTTCATCAATGCGTTTTGCTGCGGAAACTACGCGTGGATCGGTGGTGTAGATGCCCGGGACGt                                  >  1:862496/1‑68 (MQ=255)

                           cgcgATTTCATCAATGCGTTTTGCTGCGGAAACTACGCGTGGATCGGTGGTGTAGATGCCCGGGACGt                                  >  1:3244879/1‑68 (MQ=255)

                           cgcgATTTCATCAATGCGTTTTGCTGCGGAAACTACGCGTGGATCGGTGGTGTAGATGCCCGGGACGt                                  >  1:2256455/1‑68 (MQ=255)

                           cgcgATTTCATCAATGCGTTTTGCTGCGGAAACTACGCGTGGATCGGTGGTGTAGATGCCCGGGACGt                                  >  1:1499954/1‑68 (MQ=255)

                           cgcgATTTCATCAATGCGTTTTGCTGCGGAAACTACGCGTGGATCGGTGGTGTAGATGCCCGGGACGTc                                 >  1:737343/1‑69 (MQ=255)

                           cgcgATTTCATCAATGCGTTTTGCTGCGGAAACTACGCGTGGATCGGTGGTGTAGATGCCCGGGACGTc                                 >  1:709748/1‑69 (MQ=255)

                           cgagATTTCATCAATGCGTTTTGCTGCGGAAACTACGCGTGGATCGGTGGTGTAGATGCCCGGGACGTc                                 >  1:2140333/4‑69 (MQ=255)

                              gATTTCATCAATGCGTTTTGCTGCGGAAACTACGCGTGGATCGGTGGTGTa                                                <  1:2433927/51‑1 (MQ=255)

                                              tttGCTGCGGAAACTACGCGTGGATCGGTGGTGTAGATGCCCGGGACGTCGGTCCAGATATCAACACga              <  1:125481/69‑1 (MQ=255)

                                              tttGCTGCGGAAACTACGCGTGGATCGGTGGTGTAGATGCCCGGGACGTCGGTCCAGATATCAACACga              <  1:709702/69‑1 (MQ=255)

                                              tttGCTGCGGAAACTACGCGTGGATCGGTGGTGTAGATGCCCGGGACGTCGGTCCAGATATCAACACga              <  1:500483/69‑1 (MQ=255)

                                              tttGCTGCGGAAACTACGCGTGGATCGGTGGTGTAGATGCCCGGGACGTCGGTCCAGATATCAACACga              <  1:2827756/69‑1 (MQ=255)

                                              tttGCTGCGGAAACTACGCGTGGATCGGTGGTGTAGATGCCCGGGACGTCGGTCCAGATATCAACACga              <  1:1180705/69‑1 (MQ=255)

                                              tttGCTGCGGAAACTACGCGTGGATCGGTGGTGTAGATGCCCGGGACGTCGGTCCAGATATCAACACga              <  1:2497400/69‑1 (MQ=255)

                                              tttGCTGCGGAAACTACGCGTGGATCGGTGGTGTAGATGCCCGGGACGTCGGTCCAGATATCAACACga              <  1:1882288/69‑1 (MQ=255)

                                              tttGCTGCGGAAACTACGCGTGGATCGGTGGTGTAGATGCCCGGGACGTCGGTCCAGATATCAACACga              <  1:2025439/69‑1 (MQ=255)

                                               ttGCTGCGGAAACTACGCGTGGATCGGTGGTGTAGATGCCCGGGACGTCGGTCCAGATATCAACACga              <  1:2990760/68‑1 (MQ=255)

                                               ttGCTGCGGAAACCACGCGCTGATCGGTGGTGTAGATGCCCGGGACGTCGGTCCa                           <  1:1302049/55‑1 (MQ=255)

                                               ttGCTGCGGAAACCACGCGCGGATCGGTGGTGTAGATGCCCGGGACGTCGGTCCa                           <  1:1232681/55‑1 (MQ=255)

                                               ttGCTGCGGAAACCACGCGCGGATCGGTGGTGTAGATGCCCGGGACGTCGGTCCa                           <  1:3297702/55‑1 (MQ=255)

                                               ttGCTGCGGAAACCACGCGCGGATCGGTGGTGTAGATGCCCGGGACGTCGGTCCa                           <  1:657812/55‑1 (MQ=255)

                                               ttGCTGCGGAAACCACGCGCGGATCGGTGGTGTAGATGCCCGGGACGTCGGTCCa                           <  1:876569/55‑1 (MQ=255)

                                               ttGCTGCGGAAACCACGCGCGGATCGGTGGTGTAGATGCCCGGGACGTCGGTCCa                           <  1:2781600/55‑1 (MQ=255)

                                                          accaCGCGCGGATCGGTGGTGTAGATGCCCGGGACGTCGGTCTAGATATCAACACGAGATGCGTGTaaa  >  1:1673181/4‑69 (MQ=255)

                                                          accaCGCGCGGATCGGTGGTGTAGATGCCCGGGACGTCGGTCCAGATATCAACACGAGATg          >  1:3144988/4‑61 (MQ=255)

                                                          accaCGCGCGGATCGGTGGTGTAGATGCCCGGGACGTCGGTCCAGATATCAACACGAGATGCGTGTaaa  >  1:890548/4‑69 (MQ=255)

                                                          accaCGCGCGGATCGGTGGTGTAGATGCCCGGGACGTCGGTCCAGATATCAACACGAGATGCGTGTaaa  >  1:955178/4‑69 (MQ=255)

                                                          accaCGCGCGGATCGGTGGTGTAGATGCCCGGGACGTCGGTCCAGATATCAACACGAGATGCGTGTaaa  >  1:3404999/4‑69 (MQ=255)

                                                          accaCGCGCGGATCGGTGGTGTAGATGCCCGGGACGTCGGTCCAGATATCAACACGAGATGCGTGTaaa  >  1:2862476/4‑69 (MQ=255)

                                                          accaCGCGCGGATCGGTGGTGTAGATGCCCGGGACGTCGGTCCAGATATCAACACGAGATGCGTGTaaa  >  1:2488947/4‑69 (MQ=255)

                                                          accaCGCGCGGATCGGTGGTGTAGATGCCCGGGACGTCGGTCCAGATATCAACACGAGATGCGTGTaaa  >  1:22266/4‑69 (MQ=255)

                                                          accaCGCGCGGATCGGTGGTGTAGATGCCCGGGACGTCGGTCCAGATATCAACACGAGATGCGTGTaaa  >  1:2144845/4‑69 (MQ=255)

                                                          accaCGCGCGGATCGGTGGTGTAGATGCCCGGGACGTCGGTCCAGATATCAACACGAGATGCGTGTaaa  >  1:213422/4‑69 (MQ=255)

                                                          accaCGCGCGGATCGGTGGTGTAGATGCCCGGGACGTCGGTCCAGATATCAACACGAGATGCGTGTaaa  >  1:1554143/4‑69 (MQ=255)

                                                          accaCGCGCGGATCGGTGGTGTAGATGCCCGGGACGTCGGTCCAGATATCAACACGAGATGCGTGTaaa  >  1:1232911/4‑69 (MQ=255)

                                                          accaCGCGCGGATCGGTGGTGTAGATGCCCGGGACGTCGGTCCAGATATCAACACGAGATGCGTGTaa   >  1:1222993/4‑68 (MQ=255)

AGTTGCCATCTCTGCCGCTTCGGCAAACGCGATTTCATCAATGCGTTTTGCTGCGGAAACTACGCGTGGATCGGTGGTGTAGATGCCCGGGACGTCGGTCCAGATATCAACACGAGATGCGTGTAAA  >  minE/2787673‑2787799

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