BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	414,813	A→G	42.8%	C91C (TGT→TGC)	ybeX ←	predicteed ion transport

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	414,813	0	A	G	42.8%	12.1 / 40.8	56	C91C (TGT→TGC)	ybeX	predicteed ion transport
Reads supporting (aligned to +/- strand): ref base A (14/18); new base G (7/17); total (21/35)
Fisher's exact test for biased strand distribution p-value = 4.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCGTCTGGTTGCGTTTCAGGGTAATCATCTGGGAGCGGGGGATCATGATGTCGCG  >  minE/414772‑414878

aCCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCGTCTGGTTGCGTTTc                                        <  1:1979308/69‑1 (MQ=255)

   gggAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCGTCTGGTTGCGTTTCAg                                      >  1:2768422/1‑68 (MQ=255)

   gggAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCGTCTGGTTGCGTTTCAg                                      >  1:3173207/1‑68 (MQ=255)

          gTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCGTCTGGTTGCGTTTCAGGGTAatca                              <  1:947299/69‑1 (MQ=255)

          gTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCGTCTGGTTGCGTTTCAGGGTAatca                              <  1:2716412/69‑1 (MQ=255)

          gTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCGTCTGGTTGCGTTTCAGGGTAatca                              <  1:2990370/69‑1 (MQ=255)

          gTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCGTCTGGTTGCGTTTCAGGGTAatca                              <  1:3144230/69‑1 (MQ=255)

          gTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCGTCTGGTTGCGTTTCAGGGTAatca                              <  1:1934927/69‑1 (MQ=255)

          gTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCGTCTGGTTGCGTTTCAGGGTAatca                              <  1:804500/69‑1 (MQ=255)

          gTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCGTCTGGTTGCGTTTCAGGGTAatca                              <  1:808406/69‑1 (MQ=255)

          gTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCGTCTGGTTGCGTTTCAGGGTAatca                              <  1:1450022/69‑1 (MQ=255)

          gTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCGTCTGGTTGCGTCTCAGGGTAatca                              <  1:1740184/69‑1 (MQ=255)

             aGTGGGCGGACTCAATGATGACATCAAGGCATTCGTCCAGCGTCTGGTTGCGTTTc                                        >  1:33577/1‑56 (MQ=255)

             aGTGGGCGGACTCAATGATGACATCAAGGCATTCGTCCAGCGTCTGGTTGCGTTTc                                        >  1:3222083/1‑56 (MQ=255)

             aGTGGGCGGACTCAATGATGACATCAAGGCATTCGTCCAGCGTCTGGTTGCGTTTc                                        >  1:350601/1‑56 (MQ=255)

             aGTGGGCGGACTCAATGATGACATCAAGGCATTCGTCCAGCGTCTGGTTGCGTTTc                                        >  1:380268/1‑56 (MQ=255)

             aGTGGGCGGACTCAATGATGACATCAAGGCATTCGTCCAGCGTCTGGTTGCGTTTc                                        >  1:696659/1‑56 (MQ=255)

             aGTGGGCGGACTCAATGATGACATCAAGGCATTCGTCCAGCGTCTGGTTGCGTTTc                                        >  1:1216053/1‑56 (MQ=255)

             aGTGGGCGGACTCAATGATGACATCAAGGCATTCGTCCAGCGTCTGGTTGCGTTTc                                        >  1:1679787/1‑56 (MQ=255)

                             gatgaCATCAAGACATTCGTCCAGCGTCTGGTTGCGTTTCAGGGTAATCATCTGGGAGCGGGGGATCat           >  1:936845/1‑69 (MQ=255)

                                 aCATCAAGGCATTCGTCCAGCGTCTGGTTGCGTTTCAGGGTAATCATCTGGGAGCGGGGGATCatgatg       <  1:1545975/69‑1 (MQ=255)

                                 aCATCAAGGCATTCGTCCAGCGTCTGGTTGCGTTTCAGGGTAATCATCTGGGAGCGGGGGATCatgatg       <  1:57612/69‑1 (MQ=255)

                                 aCATCAAGGCATTCGTCCAGCGTCTGGTTGCGTTTCAGGGTAATCATCTGGGAGCGGGGGATCatgatg       <  1:576233/69‑1 (MQ=255)

                                 aCATCAAGGCATTCGTCCAGCGTCTGGTTGCGTTTCAGGGTAATCATCTGGGAGCGGGGGATCatgatg       <  1:638404/69‑1 (MQ=255)

                                 aCATCAAGGCATTCGTCCAGCGTCTGGTTGCGTTTCAGGGTAATCATCTGGGAGCGGGGGATCatgatg       <  1:3346306/69‑1 (MQ=255)

                                 aCATCAAGGCATTCGTCCAGCGTCTGGTTGCGTTTCAGGGTAATCATCTGGGAGCGGGGGATCatgatg       <  1:3336612/69‑1 (MQ=255)

                                 aCATCAAGGCATTCGTCCAGCGTCTGGTTGCGTTTCAGGGTAATCATCTGGGAGCGGGGGATCatgatg       <  1:1921027/69‑1 (MQ=255)

                                 aCATCAAGGCATTCGTCCAGCGTCTGGTTGCGTTTCAGGGTAATCATCTGGGAGCGGGGGATCatgatg       <  1:747274/69‑1 (MQ=255)

                                 aCATCAAGGCATTCGTCCAGCGTCTGGTTGCGTTTCAGGGTAATCATCTGGGAGCGGGGGATCatgatg       <  1:279487/69‑1 (MQ=255)

                                 aCATCAAGGCATTCGTCCAGCGTCTGGTTGCGTTTCAGGGTAATCATCTGGGAGCGGGGGATCatgatg       <  1:2734848/69‑1 (MQ=255)

                                 aCATCAAGGCATTCGTCCAGCGTCTGGTTGCGTTTCAGGGTAATCATCTGGGAGCGGGGGATCatgatg       <  1:2366335/69‑1 (MQ=255)

                                 aCATCAAGGCATTCGTCCAGCGTCTGGTTGCGTTTCAGGGTAATCATCTGGGAGCGGGGGATCatgatg       <  1:2658372/69‑1 (MQ=255)

                                 aCATCAAGGCATTCGTCCAGCGTCTGGTTGCGTTTCAGGGTAATCATCTGGGAGCGGGGGATCatgatg       <  1:2585225/69‑1 (MQ=255)

                                 aCATCAAGGCATTCGTCCAGCGTCTGGTTGCGTTTCAGGGTAATCATCTGGGAGCGGGGGATCatgatg       <  1:2438232/69‑1 (MQ=255)

                                  cATCAAGGCATTCGTCGAGCGGCTGGTTGCGTTTCAGGGTAATCATCTGGGAGCGGGGGATCatgatg       <  1:730648/68‑1 (MQ=255)

                                  cATCAAGGCATTCGTCCAGCGTCTGGTTGCGTTTCAGGGTAATCATCTGGGAGCGGGGGATCatgatg       <  1:3101780/68‑1 (MQ=255)

                                  cATCAAGGCATTCGTCCAGCGTCTGGTTGCGTTTCAGGGTAATCATCTGGGAGCGGGGGATCatgatg       <  1:3105861/68‑1 (MQ=255)

                                     cAAGACATTCGTCCAGCGTCTGGTTGCGTTTCAGGGTAATCATCTGGGAGCGGGGGAt              <  1:809249/58‑1 (MQ=255)

                                     cAAGACATTCGTCCAGCGTCTGGTTGCGTTTCAGGGTAATCATCTGGGAGCGGGGGAt              <  1:1515084/58‑1 (MQ=255)

                                     cAAGACATTCGTCCAGCGTCTGGTTGCGTTTCAGGGTAATCATCTGGGAGCGGGGGAt              <  1:2760109/58‑1 (MQ=255)

                                     cAAGACATTCGTCCAGCGTCTGGTTGCGTTTCAGGGTAATCATCTGGGAGCGGGGGAt              <  1:2755948/58‑1 (MQ=255)

                                      aaGACATTCGTCCAGCGTGTGGTTGCGTTTCAGGGTAATCATCTGGGAGCg                    >  1:1685000/1‑51 (MQ=255)

                                      aaGACATTCGTCCAGCGTCTGGTTGCGTTTCATGGTAATCATCTGGGAGCGGGGGATCATGATGTcgcg  >  1:1069713/1‑69 (MQ=255)

                                      aaGACATTCGTCCAGCGTCTGGTTGCGTTTCAGGGTAATCATCTGGGCGCGGGGGATCATGATGTcgcg  >  1:2096558/1‑69 (MQ=255)

                                      aaGACATTCGTCCAGCGTCTGGTTGCGTTTCAGGGTAATCATCTGGGAGCGGGGGATCATGATGTcgcg  <  1:440003/69‑1 (MQ=255)

                                      aaGACATTCGTCCAGCGTCTGGTTGCGTTTCAGGGTAATCATCTGGGAGCGGGGGATCATGATGTcgcg  <  1:1612219/69‑1 (MQ=255)

                                      aaGACATTCGTCCAGCGTCTGGTTGCGTTTCAGGGTAATCATCTGGGAGCGGGGGATCATGATGTcgcg  >  1:372105/1‑69 (MQ=255)

                                      aaGACATTCGTCCAGCGTCTGGTTGCGTTTCAGGGTAATCATCTGGGAGCGGGGGATCATGATGTcgcg  >  1:3353892/1‑69 (MQ=255)

                                      aaGACATTCGTCCAGCGTCTGGTTGCGTTTCAGGGTAATCATCTGGGAGCGGGGGATCATGATGTcgcg  >  1:1581251/1‑69 (MQ=255)

                                      aaGACATTCGTCCAGCGTCTGGTTGCGTTTCAGGGTAATCATCTGGGAGCGGGGGATCATGATGTcgcg  >  1:1702833/1‑69 (MQ=255)

                                      aaGACATTCGTCCAGCGTCTGGTTGCGTTTCAGGGTAATCATCTGGGAGCGGGGGATCATGATGTcgcg  <  1:1936653/69‑1 (MQ=255)

                                      aaGACATTCGTCCAGCGTCTGGTTGCGTTTCAGGGTAATCATCTGGGAGCGGGGGATCATGATGTcgcg  <  1:766872/69‑1 (MQ=255)

                                      aaGACATTCGTCCAGCGTCTGGTTGCGTTTCAGGGTAATCATCTGGGAGCGGGGGATCATGATGTcgcg  >  1:28053/1‑69 (MQ=255)

                                      aaGACATTCGTCCAGCGTCTGGTTGCGTTTCAGGGTAATCATCTGGGAGCGGGGGATCATGATGTcgcg  >  1:2308890/1‑69 (MQ=255)

                                      aaGACATTCGTCCAGCGTCTGGTTGCGTTTCAGGGTAATCATCTGGGAGCGGGGGATCATGATGTcgcg  >  1:2475295/1‑69 (MQ=255)

                                      aaGACATTCGTCCAGAGTCTGGTTGCGTTTCAGGGTAATCATCTGGGAGCGGGGGATCATGATGTcgcg  >  1:3047747/1‑69 (MQ=255)

ACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCGTCTGGTTGCGTTTCAGGGTAATCATCTGGGAGCGGGGGATCATGATGTCGCG  >  minE/414772‑414878

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