BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	125,180	G→A	33.2%	H262H (CAC→CAT)	gcd ←	glucose dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	125,180	0	G	A	33.2%	84.8 / 53.8	102	H262H (CAC→CAT)	gcd	glucose dehydrogenase
Reads supporting (aligned to +/- strand): ref base G (16/52); new base A (18/16); total (34/68)
Fisher's exact test for biased strand distribution p-value = 3.99e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

Only 100 of 102 total aligned reads displayed.

ACTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTT  >  minE/125121‑125233

aCTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGaaa                                              >  1:690617/1‑69 (MQ=255)

aCTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGaaa                                              >  1:602357/1‑69 (MQ=255)

aCTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGaaa                                              >  1:1713968/1‑69 (MQ=255)

aCTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGaaa                                              >  1:1994441/1‑69 (MQ=255)

aCTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGaaa                                              >  1:1441309/1‑69 (MQ=255)

aCTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGaaa                                              >  1:1923513/1‑69 (MQ=255)

aCTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGaaa                                              >  1:2828240/1‑69 (MQ=255)

aCTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGaaa                                              >  1:3279672/1‑69 (MQ=255)

aCTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGaa                                               >  1:2847457/1‑68 (MQ=255)

aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag                                              >  1:3083329/1‑68 (MQ=255)

aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag                                              >  1:588566/1‑68 (MQ=255)

aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag                                              >  1:2992560/1‑68 (MQ=255)

aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag                                              >  1:2896890/1‑68 (MQ=255)

aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag                                              >  1:1910709/1‑68 (MQ=255)

aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag                                              >  1:1946691/1‑68 (MQ=255)

aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag                                              >  1:2066344/1‑68 (MQ=255)

aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag                                              >  1:2624873/1‑68 (MQ=255)

aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag                                              >  1:1058307/1‑68 (MQ=255)

aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag                                              >  1:2463468/1‑68 (MQ=255)

aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag                                              >  1:1263375/1‑68 (MQ=255)

aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag                                              >  1:337114/1‑68 (MQ=255)

aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag                                              >  1:774835/1‑68 (MQ=255)

aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag                                              >  1:134501/1‑68 (MQ=255)

aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag                                              >  1:3350012/1‑68 (MQ=255)

aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaag                                              >  1:324447/1‑68 (MQ=255)

aCTTCCGGCGAAGCAGTTTCTGCTTTGGCTTCATGATAAGAGACGCCACGGCAGGTTACATGCTGGaa                                               >  1:2860879/1‑68 (MQ=255)

            gCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGtt                                  >  1:3145991/1‑69 (MQ=255)

                 ttCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCa                              >  1:164464/1‑68 (MQ=255)

                 ttCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAg                             <  1:2585977/69‑1 (MQ=255)

                 ttCTGCTTTGGCTTCATGATAAGAGACACCAAGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAg                             <  1:2699784/69‑1 (MQ=255)

                 ttCTGCTTTGGCTTCATGATAAGAGA‑‑‑CACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTc                          >  1:822324/1‑69 (MQ=255)

                     gCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTc                          <  1:1452121/68‑1 (MQ=255)

                     gCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTc                          <  1:1362334/68‑1 (MQ=255)

                     gCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTc                          <  1:3144600/68‑1 (MQ=255)

                      cTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGa                                                <  1:119935/45‑1 (MQ=255)

                        ttGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAgg                       <  1:3393325/68‑1 (MQ=255)

                        ttGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAgg                       <  1:2360497/68‑1 (MQ=255)

                        ttGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAgg                       <  1:2743821/68‑1 (MQ=255)

                        ttGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAgg                       <  1:3251234/68‑1 (MQ=255)

                        ttGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAgg                       <  1:3048684/68‑1 (MQ=255)

                        ttGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAgg                       <  1:319048/68‑1 (MQ=255)

                        ttGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAgg                       <  1:2064846/68‑1 (MQ=255)

                        ttGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAgg                       <  1:212890/68‑1 (MQ=255)

                             ttCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:1852745/68‑1 (MQ=255)

                             ttCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:1800757/68‑1 (MQ=255)

                             ttCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:1797281/68‑1 (MQ=255)

                             ttCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:1753021/68‑1 (MQ=255)

                             ttCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:3168844/68‑1 (MQ=255)

                             ttCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:2788455/68‑1 (MQ=255)

                             ttCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:1734844/68‑1 (MQ=255)

                             ttCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:2586276/68‑1 (MQ=255)

                             ttCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:736960/68‑1 (MQ=255)

                             ttCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:1505843/68‑1 (MQ=255)

                             ttCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:2217404/68‑1 (MQ=255)

                             ttCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:563938/68‑1 (MQ=255)

                                   aTAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCa           <  1:1746104/69‑1 (MQ=255)

                                   aTAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCa           <  1:519770/69‑1 (MQ=255)

                                   aTAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCa           <  1:2970460/69‑1 (MQ=255)

                                    tAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:2239916/61‑1 (MQ=255)

                                      agagACGCCACGGCAGGTTACATGCTGGAAGGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCAt          <  1:1322401/67‑1 (MQ=255)

                                      agagACGCCACGGCAGGTTACATGCTGGAAGGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCAt          <  1:1526929/67‑1 (MQ=255)

                                      agagACGCCACGGCAGGTTACATGCTGGAAGGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCAt          <  1:3433508/67‑1 (MQ=255)

                                      agagACGCCACGGCAGGTTACATGCTGGAAGGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCAt          <  1:411243/67‑1 (MQ=255)

                                      agagACGCCACGGCAGGTTACATGCTGGAAGGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCAt          <  1:3176260/67‑1 (MQ=255)

                                      agagACGCCACGGCAGGTTACATGCTGGAAGGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCAt          <  1:424362/67‑1 (MQ=255)

                                      agagACGCCACGGCAGGTTACATGCTGGAAGGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCAt          <  1:2986548/67‑1 (MQ=255)

                                      agagACGCCACGGCAGGTTACATGCTGGAAGGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCAt          <  1:2266631/67‑1 (MQ=255)

                                      agagACGCCACGGCAGGTTACATGCTGGAAGGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCAt          <  1:2272983/67‑1 (MQ=255)

                                      agagACGCCACGGCAGGTTACATGCTGGAAGGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCAt          <  1:2431273/67‑1 (MQ=255)

                                      agagACGCCACGGCAGGTTACATGCTGGAAGGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCAt          <  1:719336/67‑1 (MQ=255)

                                      agagACGCCACGGCAGGTTACATGCTGGAAGGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCAt          <  1:1881451/67‑1 (MQ=255)

                                      agagACGCCACGGCAGGTTACATGCTGGAAGGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCAt          <  1:2704519/67‑1 (MQ=255)

                                      agagACGCCACGGCAGGTTACATGCTGGAAGGACTCGTTGCTTTTCAGCTCAGGATCGTAATGCCAt          <  1:1840574/67‑1 (MQ=255)

                                       gagaCGCCACGGCAGGTTACATGCTGGAAGGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCAt          <  1:186527/66‑1 (MQ=255)

                                         gaCACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTtctc     <  1:828197/69‑1 (MQ=255)

                                         gaCACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTtctc     <  1:669712/69‑1 (MQ=255)

                                          acgccacGGCAGGTTACATGCTGGAAGGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCAt          <  1:1550887/60‑1 (MQ=255)

                                          acacCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:2165924/55‑1 (MQ=255)

                                          acacCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCa           <  1:746231/62‑1 (MQ=255)

                                          acacCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTtc       <  1:2232742/66‑1 (MQ=255)

                                          acacCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTctct    <  1:2416226/69‑1 (MQ=255)

                                           caccacGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:1258666/54‑1 (MQ=255)

                                            accacGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCttt  >  1:3320879/1‑69 (MQ=255)

                                            accacGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCttt  >  1:2888365/1‑69 (MQ=255)

                                            accacGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCttt  >  1:3273713/1‑69 (MQ=255)

                                            accacGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCttt  >  1:1921597/1‑69 (MQ=255)

                                                cGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTtctc     <  1:1913353/62‑1 (MQ=255)

                                                       ttACGTGCTTGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:397986/42‑1 (MQ=255)

                                                       ttACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:257548/42‑1 (MQ=255)

                                                       ttACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:811431/42‑1 (MQ=255)

                                                       ttACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:324414/42‑1 (MQ=255)

                                                       ttACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:2961013/42‑1 (MQ=255)

                                                       ttACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:2464512/42‑1 (MQ=255)

                                                       ttACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:2146761/42‑1 (MQ=255)

                                                       ttACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:2143799/42‑1 (MQ=255)

                                                       ttACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:1916184/42‑1 (MQ=255)

                                                       ttACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:1644124/42‑1 (MQ=255)

                                                       ttACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:1641187/42‑1 (MQ=255)

                                                       ttACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGt                  <  1:1158550/42‑1 (MQ=255)

                                                       ttACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCa           <  1:1747331/49‑1 (MQ=255)

ACTTCCGGCGAAGCGGTTTCTGCTTTGGCTTCATGATAAGAGACACCACGGCAGGTTACGTGCTGGAAAGACTCGTTGGTTTTCAGCTCAGGATCGTAATGCCATTTCTCTTT  >  minE/125121‑125233

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