BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,318,190	A→G	33.3%	D141D (GAT→GAC)	yiiT ←	stress‑induced protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,318,190	0	A	G	33.3%	97.5 / 70.9	126	D141D (GAT→GAC)	yiiT	stress‑induced protein
Reads supporting (aligned to +/- strand): ref base A (45/39); new base G (24/18); total (69/57)
Fisher's exact test for biased strand distribution p-value = 8.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

Only 100 of 130 total aligned reads displayed.

TGACGAAACATTTACGTCGCCGACCGGGCCATATGACAGCGGAGCATTACTGGTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATCCCGCGATATGCCGGCATCAA  >  minE/2318124‑2318258

tGACGAAACATTTACGTCGCCGACCGGGCCATATGACAGCGGAGCATTACTGGTAAGCGGTTACTTATc                                                                    >  1:2027677/1‑69 (MQ=255)

tGACGAAACATTTACGTCGCCGACCGGGCCATATGACAGCGGAGCATTACTGGTAAGCGGTTACTTATc                                                                    >  1:2129365/1‑69 (MQ=255)

tGACGAAACATTTACGTCGCCGACCGGGCCATATGACAGCGGAGCATTACTGGTAAGCGGTTACTTATc                                                                    >  1:2982671/1‑69 (MQ=255)

 gACGAAACATTTACGTCGCCGACCGGGCCATATGACAGCGGAGCATTACTGGTAAGCGGTTACTTATc                                                                    >  1:2145109/1‑68 (MQ=255)

 gACGAAACATTTACGTCGCCGACCGGGCCATATGACAGCGGAGCATTACTGGTAAGCGGTTACTTATc                                                                    >  1:181053/1‑68 (MQ=255)

 gACGAAACATTTACGTCGCCGACCGGGCCATATGACAGCGGAGCATTACTGGTAAGCGGTTACTTATc                                                                    >  1:2934621/1‑68 (MQ=255)

 gACGAAACATTTACGTCGCCGACCGGGCCATATGACAGCGGAGCATTACTGGTAAGCGGTTACTTATc                                                                    >  1:2051251/1‑68 (MQ=255)

       aCATTTACGTCGCCGACCGGGCCATATGACAGCGGAGCATTACTGGTAAGCGGTTACTTATCGATAAAc                                                             >  1:1030190/1‑69 (MQ=255)

       aCATTTACGTCGCCGACCGGGCCATATGACAGCGGAGCATTACTGGTAAGCGGTTACTTATCGATAAAc                                                             >  1:2417281/1‑69 (MQ=255)

       aCATTTACGTCGCCGACCGGGCCATATGACAGCGGAGCATTACTGGTAAGCGGTTACTTATCGATAAAc                                                             >  1:2313933/1‑69 (MQ=255)

       aCATTTACGTCGCCGACCGGGCCATATGACAGCGGAGCATTACTGGTAAGCGGTTACTTATCGATAAAc                                                             >  1:92701/1‑69 (MQ=255)

            tACGTCGCCGACCGGGCCATATGACAGCGGAGCATTACTGGTAAGCGGTTACTTATCGATAAACGGCAc                                                        <  1:3268954/69‑1 (MQ=255)

            tACGTCGCCGACCGGGCCATATGACAGCGGAGCATTACTGGTAAGCGGTTACTTATCGATAAACGGCAc                                                        <  1:1043615/69‑1 (MQ=255)

            tACGTCGCCGACCGGGCCATATGACAGCGGAGCATTACTGGTAAGCGGTTACTTATCGATAAACGGCAc                                                        <  1:256436/69‑1 (MQ=255)

            tACGTCGCCGACCGGGCCATATGACAGCGGAGCATTACTGGTAAGCGGTTACTTATCGATAAACGGCAc                                                        <  1:2502041/69‑1 (MQ=255)

            tACGTCGCCGACCGGGCCATATGACAGCGGAGCATTACTGGTAAGCGGTTACTTATCGATAAACGGCAc                                                        <  1:2392794/69‑1 (MQ=255)

            tACGTCGCCGACCGGGCCATATGACAGCGGAGCATTACTGGTAAGCGGGTACTTATCGATAAACGGCAc                                                        <  1:1692538/69‑1 (MQ=255)

                                 tGACAGCGGAGCATTACTGGTAAGCGGTTACTTGTCGATAAACGGCACGATCAAGAGATCCGCACTCaa                                   >  1:2004207/1‑68 (MQ=255)

                                 tGACAGCGGAGCATTACTGGTAAGCGGTTACTTGTCGATAAACGGCACGATCAAGAGATCCGCACTCaa                                   >  1:2474422/1‑68 (MQ=255)

                                 tGACAGCGGAGCATTACTGGTAAGCGGTTACTTGTCGATAAACGGCACGATCAAGAGATCCGCACTCaa                                   >  1:1968154/1‑68 (MQ=255)

                                 tGACAGCGGAGCATTACTGGTAAGCGGTTACTTGTCGATAAACGGCACGATCAAGAGATCCGCACTCaa                                   >  1:2582260/1‑68 (MQ=255)

                                 tGACAGCGGAGCATTACTGGTAAGCGGTTACTTGTCGATAAACGGCACGATCAAGAGATCCGCACTCaa                                   >  1:2654892/1‑68 (MQ=255)

                                 tGACAGCGGAGCATTACTGGTAAGCGGTTACTTGTCGATAAACGGCACGATCAAGAGATCCGCACTCaa                                   >  1:2812690/1‑68 (MQ=255)

                                 tGACAGCGGAGCATTACTGGTAAGCGGTTACTTGTCGATAAACGGCACGATCAAGAGATCCGCACTCaa                                   >  1:3376532/1‑68 (MQ=255)

                                 tGACAGCGGAGCATTACTGGTAAGCGGTTACTTGTCGATAAACGGCACGATCAAGAGATCCGCACTCaa                                   >  1:3432316/1‑68 (MQ=255)

                                 tGACAGCGGAGCATTACTGGTAAGCGGTTACTTGTCGATAAACGGCACGATCAAGAGATCCGCACTCaa                                   >  1:38755/1‑68 (MQ=255)

                                 tGACAGCGGAGCATTACTGGTAAGCGGTTACTTGTCGATAAACGGCACGATCAAGAGATCCGCACTCaa                                   >  1:512373/1‑68 (MQ=255)

                                 tGACAGCGGAGCATTACTGGTAAGCGGTTACTTGTCGATAAACGGCACGATCAAGAGATCCGCACTCaa                                   >  1:564966/1‑68 (MQ=255)

                                 tGACAGCGGAGCATTACTGGTAAGCGGTTACTTGTCGATAAACGGCACGATCAAGAGATCCGCACTCaa                                   >  1:1938670/1‑68 (MQ=255)

                                 tGACAGCGGAGCATTACTGGTAAGCGGTTACTTGTCGATAAACGGCACGATCAAGAGATCCGCACTCaa                                   >  1:1281426/1‑68 (MQ=255)

                                 tGACAGCGGAGCATTACTGGTAAGCGGTTACTTGTCGATAAACGGCACGATCAAGAGATACGCACTCaa                                   >  1:3274868/1‑68 (MQ=255)

                                 tGACAGCGGAGCATTACTGGTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACt                                      >  1:331857/1‑66 (MQ=255)

                                 tGACAGCGGAGCATTACTGGTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCAt                                   >  1:2795155/1‑69 (MQ=255)

                                 tGACAGCGGAGCATTACTGGTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCAt                                   >  1:558360/1‑69 (MQ=255)

                                 tGACAGCGGAGCATTACTGGTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCAt                                   >  1:2392176/1‑69 (MQ=255)

                                 tGACAGCGGAGCATTACTGGTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCAt                                   >  1:2325765/1‑69 (MQ=255)

                                 tGACAGCGGAGCATTACTGGTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCAt                                   >  1:212330/1‑69 (MQ=255)

                                 tGACAGCGGAGCATTACTGGTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCAt                                   >  1:1843542/1‑69 (MQ=255)

                                 tGACAGCGGAGCATTACTGGTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCAt                                   >  1:2097293/1‑69 (MQ=255)

                                 tGACAGCGGAGCATTACTGGTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATACGCACTCAt                                   >  1:3366452/1‑69 (MQ=255)

                                   aCAGCGGAGCATTACTGGTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCt                                 >  1:1466921/1‑69 (MQ=255)

                                     aGCGGAGCATTACTGGTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTa                               <  1:3157704/69‑1 (MQ=255)

                                     aGCGGAGCATTACTGGTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTa                               <  1:702766/69‑1 (MQ=255)

                                     aGCGGAGCATTACTGGTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTa                               <  1:1432227/69‑1 (MQ=255)

                                         gAGCATTACTGGTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGa                           <  1:203813/69‑1 (MQ=255)

                                          aGCATTACTGGTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGat                          <  1:2991962/69‑1 (MQ=255)

                                          aGCATTACTGGTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGat                          <  1:589522/69‑1 (MQ=255)

                                          aGCATTACTGGTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGat                          <  1:979244/69‑1 (MQ=255)

                                          aGCATTACTGGTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGat                          <  1:1885328/69‑1 (MQ=255)

                                                aCTGGTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGTTCATccc                    >  1:1185307/1‑69 (MQ=255)

                                                aCTGGTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATccc                    >  1:1708931/1‑69 (MQ=255)

                                                aCTGGTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATccc                    >  1:2012536/1‑69 (MQ=255)

                                                aCTGGTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATccc                    >  1:563942/1‑69 (MQ=255)

                                                aCTGGTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATccc                    >  1:1135818/1‑69 (MQ=255)

                                                aCTGGTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATccc                    >  1:1045051/1‑69 (MQ=255)

                                                aCTGGTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATcc                     >  1:991363/1‑68 (MQ=255)

                                                aCTGGTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATcc                     >  1:1460590/1‑68 (MQ=255)

                                                   ggTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATCCcgc                  <  1:3025806/68‑1 (MQ=255)

                                                   ggTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATCCcgc                  <  1:1298341/68‑1 (MQ=255)

                                                   ggTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATCCcgc                  <  1:893472/68‑1 (MQ=255)

                                                   ggTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATCCcgc                  <  1:1896067/68‑1 (MQ=255)

                                                   ggTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATCCcgc                  <  1:241202/68‑1 (MQ=255)

                                                   ggTAAGCGGATACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATCCcgc                  <  1:99010/68‑1 (MQ=255)

                                                    gTAAGCGGTTACTTGTCGATAAACGGCACGATCAAGAGATCCGCACTCAACTTATTGATCATCCcgcg                 >  1:3099240/1‑68 (MQ=255)

                                                    gTAAGCGGTTACTTGTCGATAAACGGCACGATCAAGAGATCCGCACTCAACTTATTGATCATCCCGCGa                >  1:1813008/1‑69 (MQ=255)

                                                    gTAAGCGGTTACTTGTCGATAAACGGCACGATCAAGAGATCCGCACTCAACTTATTGATCATCCCGCGa                >  1:3143246/1‑69 (MQ=255)

                                                    gTAAGCGGTTACTTGTCGATAAACGGCACGATCAAGAGATCAGCACTCAACTTATTGATCATCCCGCGa                >  1:2516997/1‑69 (MQ=255)

                                                    gTAAGCGGATACTTGTCGATAAACGGCACGATCAAGAGATCCGCACTCAACTTATTGATCATCCCGCGa                >  1:1546877/1‑69 (MQ=255)

                                                      atgCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATCCCGCGa                <  1:789469/65‑1 (MQ=255)

                                                      aaGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATCCCGCGa                <  1:71199/67‑1 (MQ=255)

                                                      aaGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATCCCGCGa                <  1:643074/67‑1 (MQ=255)

                                                      aaGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATCCCGCGa                <  1:591391/67‑1 (MQ=255)

                                                      aaGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATCCCGCGa                <  1:2379561/67‑1 (MQ=255)

                                                      aaGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATCCCGCGa                <  1:512472/67‑1 (MQ=255)

                                                      aaGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATCCCGCGa                <  1:2052908/67‑1 (MQ=255)

                                                      aaGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATCCCGCGa                <  1:1373795/67‑1 (MQ=255)

                                                       aGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATCCCGCGa                <  1:2425916/66‑1 (MQ=255)

                                                           gTTACTTGTCGATAAACGGCACGATCAAGAGATCCGCACTCAACTTATTGTTCATCCCGCGATATGCCg         <  1:2254602/69‑1 (MQ=255)

                                                           gTTACTTGTCGATAAACGGCACGATCAAGAGATCCGCACTCAACTTATTGATCATCCCGCGATATGCCg         <  1:3233883/69‑1 (MQ=255)

                                                           gTTACTTGTCGATAAACGGCACGATCAAGAGATCCGCACTCAACTTATTGATCATCCCGCGATATGCCg         <  1:279447/69‑1 (MQ=255)

                                                           gTTACTTGTCGATAAACGGCACGATCAAGAGATCCGCACTCAACTTATTGATCATCCCGCGATATGCCg         <  1:1423617/69‑1 (MQ=255)

                                                           gTTACTTGTCGATAAACGGCACGATCAAGAGATCCGCACTCAACTTATTGATCATCCCGCGATATGCCg         <  1:1726714/69‑1 (MQ=255)

                                                           gTTACTTGTCGATAAACGGCACGATCAAGAGATCCGCACTCAACTTATTGATCATCCCGCGATATGCCg         <  1:1418964/69‑1 (MQ=255)

                                                           gTTACTTGTCGATAAACGGCACGATCAAGAGATCCGCACTCAACTTATTGATCATCCCGCGATATGCCg         <  1:2696547/69‑1 (MQ=255)

                                                           gTTACTTGTCGATAAACGGCACGATCAAGAGATCCGCACTCAACTTATTGATCATCCCGCGATATGCCg         <  1:2622665/69‑1 (MQ=255)

                                                           gTTACTTGTCGATAAACGGCACGATCAAGAGATCCGCACTCAACTTATTGATCATCCCGCGATATGCCg         <  1:459810/69‑1 (MQ=255)

                                                           gTTACTTGTCGATAAACGGCACGATCAAGAGATCCGCACTCAACTTATTGATCATCCCGCGATATGCCg         <  1:1873406/69‑1 (MQ=255)

                                                           gTTACTTGTCGATAAACGGCACGATCAAGAGATCCGCACTCAACTTATTGATCATCCCGCGATATGCCg         <  1:1029141/69‑1 (MQ=255)

                                                           gTTACTTGTCGATAAACGGCACGATCAAGAGATCCGCACTCAACTTATTGATCATCCCGCGATATGCCg         <  1:3128580/69‑1 (MQ=255)

                                                           gTTACTTGTCGATAAACGGCACGATCAAGAGATCCGCACTCAACTTATTGATCATCCCGCGATATGCCg         <  1:1310622/69‑1 (MQ=255)

                                                            ttACTTGTCGATAAACGGCACGATCAAGAGATCCGCACTCAACTTATTGATCATCCCGCGATATGCCg         <  1:2251773/68‑1 (MQ=255)

                                                              aCTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATCCCGCGATATGCCGGc       >  1:726889/1‑68 (MQ=255)

                                                              aCTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATCCCGCGATATGCCGGc       >  1:278587/1‑68 (MQ=255)

                                                              aCTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATCCCGCGATATGCCGGc       >  1:1705759/1‑68 (MQ=255)

                                                              aCTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATCCCGCGATATGCCGGc       >  1:2529629/1‑68 (MQ=255)

                                                                 tATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATCCCGCGa                <  1:580748/56‑1 (MQ=255)

                                                                 tATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATCCCGCGa                <  1:1235985/56‑1 (MQ=255)

                                                                 tATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATCCCGCGa                <  1:3059815/56‑1 (MQ=255)

                                                                  aTCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATCCCGCGATATGCCGGCATCaa  >  1:1799699/1‑69 (MQ=255)

                                                                  aTCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATCCCGCGATATGCCGGCATCaa  >  1:2574453/1‑69 (MQ=255)

TGACGAAACATTTACGTCGCCGACCGGGCCATATGACAGCGGAGCATTACTGGTAAGCGGTTACTTATCGATAAACGGCACGATGAGCAAATCCGCACTCATCTTATTGATCATCCCGCGATATGCCGGCATCAA  >  minE/2318124‑2318258

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