BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,246,332	G→C	27.7%	G56G (GGC→GGG)	znuA ←	zinc transporter subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,246,332	0	G	C	27.7%	73.8 / 29.7	54	G56G (GGC→GGG)	znuA	zinc transporter subunit
Reads supporting (aligned to +/- strand): ref base G (6/33); new base C (2/13); total (8/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGACTAAGTCCGCGTTCTGTAAGCGTTTTACATCCGATGGGCGCAGTGAATAATCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGCAATGGCAGAAGCGATGAACCC  >  minE/1246265‑1246394

cGACTAAGTCCGCGTTCTGTAAGCGTTTTACATCCGATGGGCGCAGTGAAT‑ATCATGTTCTGAAGCGcc                                                              <  1:217089/69‑1 (MQ=255)

 gACTAAGTCCGCGTTCTGTAAGCGTTTTACATCCGATGGGCGCAGTGAATAATCATGTTCTGAAGCGcc                                                              <  1:538383/69‑1 (MQ=255)

 gACTAAGTCCGCGTTCTGTAAGCGTTTTACATCCGATGGGCGCAGTGAATAATCATGTTCTGAAGCGcc                                                              <  1:2942213/69‑1 (MQ=255)

        tCCGCGTTCTGTAAGCGTTTTACATCCGATGGGCGCAGTGAATAATCATGTTCTGAAGCGCCGTCAgg                                                        <  1:1014615/68‑1 (MQ=255)

        tCCGCGTTCTGTAAGCGTTTTACATCCGATGGGCGCAGTGAATAATCATGTTCTGAAGCGCCGTCAgg                                                        <  1:3132150/68‑1 (MQ=255)

        tCCGCGTTCTGTAAGCGTTTTACATCCGATGGGCGCAGTGAATAATCATGTTCTGAAGCGCCGTCAgg                                                        <  1:2856305/68‑1 (MQ=255)

        tCCGCGTTCTGTAAGCGTTTTACATCCGATGGGCGCAGTGAATAATCATGTTCTGAAGCGCCGTCAgg                                                        <  1:2636033/68‑1 (MQ=255)

        tCCGCGTTCTGTAAGCGTTTTACATCCGATGGGCGCAGTGAATAATCATGTTCTGAAGCGCCGTCAgg                                                        <  1:2173329/68‑1 (MQ=255)

        tCCGCGTTCTGTAAGCGTTTTACATCCGATGGGCGCAGTGAATAATCATGTTCTGAAGCGCCGTCAgg                                                        <  1:1460262/68‑1 (MQ=255)

             gTTCTGTAAGCGTTTTACATCCGATGGGCGCAGTGAATAATCATGTTCTGAAGCGCCGTCAGGAAGTaa                                                  <  1:2257861/69‑1 (MQ=255)

             gTTCTGTAAGCGTTTTACATCCGATGGGCGCAGTGAATAATCATGTTCTGAAGCGCCGTCAGGAAGTaa                                                  <  1:1004677/69‑1 (MQ=255)

             gTTCTGTAAGCGTTTTACATCCGATGGGCGCAGTGAATAATCATGTTCTGAAGCGCCGTCAGGAAGTaa                                                  <  1:1666979/69‑1 (MQ=255)

             gTTCTGTAAGCGTTTTACATCCGATGGGCGCAGTGAATAATCATGTTCTGAAGCGCCGTCAGGAAGTaa                                                  <  1:121197/69‑1 (MQ=255)

                        gTTTTACATCCGATGGGCGCAGTGAATAATCATGTTCTGAAGCGCCGTCAgg                                                        <  1:1545669/52‑1 (MQ=255)

                         ttttACATCCGATGGGCGCAGTGAATAATCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCtctgtttc                                      <  1:1880761/69‑1 (MQ=255)

                             aCATCCGATGGGCGCAGTGAATAATCATGTTTTGAAGCGCCGTCAGGAAGTAAAACCtctgtttctgt                                   <  1:2668689/68‑1 (MQ=255)

                             aCATCCGATGGGCGCAGTGAATAATCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCtctgtttctgt                                   <  1:2064587/68‑1 (MQ=255)

                                 ccGATGGGCGCAGTGAATAATCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTAccc                              <  1:1220211/69‑1 (MQ=255)

                                 ccGATGGGCGCAGTGAATAATCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTAccc                              <  1:2307036/69‑1 (MQ=255)

                                  cGATGGGCGCAGTGAATAATCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTAcccc                             >  1:1446423/1‑69 (MQ=255)

                                  cGATGGGCGCAGTGAATAATCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTAcccc                             >  1:2717327/1‑69 (MQ=255)

                                         cgcAGTGAATAATCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTACCCGATCAGc                       <  1:439987/68‑1 (MQ=255)

                                         cgcAGTGAATAATCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGc                       <  1:888729/68‑1 (MQ=255)

                                         cgcAGTGAATAATCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGc                       <  1:1325376/68‑1 (MQ=255)

                                         cgcAGTGAATAATCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGc                       <  1:328388/68‑1 (MQ=255)

                                         cgcAGTGAATAATCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGc                       <  1:1038728/68‑1 (MQ=255)

                                         cgcAGTGAATAATCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGc                       <  1:2558948/68‑1 (MQ=255)

                                         cgcAGTGAATAATCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGc                       <  1:2391271/68‑1 (MQ=255)

                                         cgcAGTGAATAATCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGc                       <  1:2212824/68‑1 (MQ=255)

                                         cgcAGTGAATAATCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGc                       <  1:2194651/68‑1 (MQ=255)

                                         cgcAGTGAATAATCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGc                       <  1:1160504/68‑1 (MQ=255)

                                               gAATAATCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGc                       <  1:627760/62‑1 (MQ=255)

                                               gAATAATCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGc                       <  1:19539/62‑1 (MQ=255)

                                                aataatCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGc                       <  1:3394776/61‑1 (MQ=255)

                                                  taatCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGCa                      >  1:188292/1‑60 (MQ=255)

                                                  taatCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGCa                      >  1:646836/1‑60 (MQ=255)

                                                  taatCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGCa                      >  1:1019831/1‑60 (MQ=255)

                                                  taatCAAGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGCa                      >  1:1022/1‑60 (MQ=255)

                                                    atCATGTTCTGAAGCCCCGTCAGGTAGTAAAACCTCCGTTTCTGTTACCCCATCAGCAATGGCAGAAGc           <  1:1301558/69‑1 (MQ=255)

                                                    atCATGTTCTGAAGCCCCGTCAGGTAGTAAAACCTCCGTTTCTGTTACCCCATCAGCAATGGCAGAAGc           <  1:2875961/69‑1 (MQ=255)

                                                    atCATGTTCTGAAGCCCCGTCAGGTAGTAAAACCTCCGTTTCTGTTACCCCATCAGCAATGGCAGAAGc           <  1:272318/69‑1 (MQ=255)

                                                    atCATGTTCTGAAGCCCCGTCAGGTAGTAAAACCTCCGTTTCTGTTACCCCATCAGCAATGGCAGAAGc           <  1:1254785/69‑1 (MQ=255)

                                                    atCATGTTCTGAAGCCCCGTCAGGTAGTAAAACCTCCGTTTCTGTTACCCCATCAGCAATGGCAGAAGc           <  1:2045452/69‑1 (MQ=255)

                                                    atCATGTTCTGAAGCCCCGTCAGGTAGTAAAACCTCCGTTTCTGTTACCCCATCAGCAATGGCAGAAGc           <  1:842386/69‑1 (MQ=255)

                                                    atCATGTTCTGAAGCCCCGTCAGGTAGTAAAACCTCCGTTTCTGTTACCCCATCAGCAATGGCAGAAGc           <  1:1153876/69‑1 (MQ=255)

                                                    atCATGTTCTGAAGCCCCGTCAGGTAGTAAAACCTCCGTTTCTGTTACCCCATCAGCAATGGCAGAAGc           <  1:2472727/69‑1 (MQ=255)

                                                    atCATGTTCTGAAGCCCCGTCAGGTAGTAAAACCTCCGTTTCTGTTACCCCATCAGCAATGGCAGAAGc           <  1:2441270/69‑1 (MQ=255)

                                                    atCATGTTCTGAAGCCCCGTCAGGTAGTAAAACCTCCGTTTCTGTTACCCCATCAGCAATGGCAGAAGc           <  1:691073/69‑1 (MQ=255)

                                                    atCATGTTCTGAAGCCCCGTCAGGTAGTAAAACCTCCGTTTCTGTTACCCCATCAGCAATGGCAGAAGc           <  1:721927/69‑1 (MQ=255)

                                                    atCATGTTCTGAAGCCCCGTCAGGTAGTAAAACCTCCGTTTCTGTTACCCCATCAGCAATGGCAGAAGc           <  1:838085/69‑1 (MQ=255)

                                                      cATGTTCTGAAGCCCCGTCAGGTAGTAAAACCTCCGTTTCTGTTACCCCATCAGCAATGGCAGAAGc           <  1:2931890/67‑1 (MQ=255)

                                                            cTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGc                       <  1:1491809/49‑1 (MQ=255)

                                                             tGAAGCCCCGTCAGGTAGTAAAACCTCCGTTTCTGTTACCCCATCAGCAATGGCAGAAGCGATGAAccc  >  1:2812376/1‑69 (MQ=255)

                                                             tGAAGCCCCGTCAGGTAGTAAAACCTACGTTTCTGTTACCCCATCAGCAATGGCAGAAGCGATGAAccc  >  1:1064795/1‑69 (MQ=255)

CGACTAAGTCCGCGTTCTGTAAGCGTTTTACATCCGATGGGCGCAGTGAATAATCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGCAATGGCAGAAGCGATGAACCC  >  minE/1246265‑1246394

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