BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,419,735	T→C	21.9%	K76K (AAA→AAG)	ccmH ←	heme lyase, CcmH subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,419,735	0	T	C	21.9%	134.4 / 29.7	100	K76K (AAA→AAG)	ccmH	heme lyase, CcmH subunit
Reads supporting (aligned to +/- strand): ref base T (16/62); new base C (11/11); total (27/73)
Fisher's exact test for biased strand distribution p-value = 1.25e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 129 total aligned reads displayed.

AACGGCGGATCGTAAGTGACGAAGTTGCCGTAACGCGCCACCATATAATCGACAATCTCTTTTTTACTTTTACCTTCCTGCATCAGTTCATACACTTTCTGACGCAGGTCGGTGGCAATCATCG  >  minE/1419673‑1419796

aaCGGCGGATCGTAAGTGACGAAGTTGCCGTAACGCGCCACCATATAATCGACAATCTCTTttttact                                                          >  1:2627199/1‑68 (MQ=255)

aaCGGCGGATCGTAAGTGACGAAGTTGCCGTAACGCGCCACCATATAATCGACAATCTCTTttttact                                                          >  1:3296521/1‑68 (MQ=255)

aaCGGCGGATCGTAAGTGACGAAGTTGCCGTAACGCGCCACCATATAATCGACAATCTCTTttttact                                                          >  1:2163936/1‑68 (MQ=255)

aaCGGCGGATCGTAAGTGACGAAGTTGCCGTAACGCGCCACCATATAATCGACAATCTCTTttttact                                                          >  1:2964377/1‑68 (MQ=255)

aaCGGCGGATCGTAAGTGACGAAGTTGCCGTAACGCGCCACCATATAATCGACAATCTCTTttttact                                                          >  1:760950/1‑68 (MQ=255)

aaCGGCGGATCGTAAGTGACGAAGTTGCCGTAACGCGCCACCATATAATCGACAATCTCTTttttact                                                          >  1:1590455/1‑68 (MQ=255)

   ggcggATCGTAAGTGACGAAGTTGCCGTAACGCGCCACCATATAATCGACAATCTCTTttttactttta                                                      <  1:402540/69‑1 (MQ=255)

   ggcggATCGTAAGTGACGAAGTTGCCGTAACGCGCCACCATATAATCGACAATCTCTTttttactttta                                                      <  1:596101/69‑1 (MQ=255)

   ggcggATCGTAAGTGACGAAGTTGCCGTAACGCGCCACCATATAATCGACAATCTCTTttttactttta                                                      <  1:110667/69‑1 (MQ=255)

   ggcggATCGTAAGTGACGAAGTTGCCGTAACGCGCCACCATATAATCGACAATCTCTTttttactttta                                                      <  1:1249776/69‑1 (MQ=255)

   ggcggATCGTAAGTGACGAAGTTGCCGTAACGCGCCACCATATAATCGACAATCTCTTttttactttta                                                      <  1:1513488/69‑1 (MQ=255)

   ggcggATCGTAAGTGACGAAGTTGCCGTAACGCGCCACCATATAATCGACAATCTCTTttttactttt                                                       <  1:3284002/68‑1 (MQ=255)

   ggcggATCGTAAGTGACGAAGTTGCCGTAACGCGCCACCATATAATCGACAATCTCTTttttactttt                                                       <  1:2375206/68‑1 (MQ=255)

   ggcggATCGTAAGTGACGAAGTTGCCGTAACGCGCCACCATATAATCGACAATCTCTTttttactttt                                                       <  1:1824228/68‑1 (MQ=255)

   ggcggATCGTAAGTGACGAAGTTGCCGTAACGCGCCACCATATAATCGACAATCTCTTttttactttt                                                       <  1:701521/68‑1 (MQ=255)

   ggcggATCGTAAGTGACGAAGTTGCCGTAACGCGCCACCATATAATCGACAATCTCTTttttactttt                                                       <  1:1665122/68‑1 (MQ=255)

   ggcggATCGTAAGTGACGAAGTTGCCGTAACGCGCCACCATATAATCGACAATCTCTTttttactttt                                                       <  1:3118359/68‑1 (MQ=255)

   ggcggATCGTAAGTGACGAAGTTGCCGTAACGCGCCACCATATAATCGACAATCTCTTttttactttt                                                       <  1:3066105/68‑1 (MQ=255)

       gATCGTAAGTGACGAAGTTGCCGTAACGCGCCACCATATAATCGACAATCTCTTTTTTACTTTTACCtt                                                  <  1:2039723/69‑1 (MQ=255)

       gATCGTAAGTGACGAAGTTGCCGTAACGCGCCACCATATAATCGACAATCTCTTTTTTACTTTTACCtt                                                  <  1:2345075/69‑1 (MQ=255)

       gATCGTAAGTGACGAAGTTGCCGTAACGCGCCACCATATAATCGACAATCTCTTTTTTACTTTTACCtt                                                  <  1:2025110/69‑1 (MQ=255)

       gATCGTAAGTGACGAAGTTGCCGTAACGCGCCACCATATAATCGACAATCTCTTTTTTACTTTTACCtt                                                  <  1:3042074/69‑1 (MQ=255)

       gATCGTAAGTGACGAAGTTGCCGTAACGCGCCACCATATAATCGACAATCTCTTTTTTACTTTTACCtt                                                  <  1:2942534/69‑1 (MQ=255)

       gATCGTAAGTGACGAAGTTGCCGTAACGCGCCACCATATAATCGACAATCTCTTTTTTACTTTTACCtt                                                  <  1:1288622/69‑1 (MQ=255)

         tCGTAAGTGACGAAGTTGCCGTAACGCGCCACCATATAATCGACAATCTCTTTTTTACTTTTACCTTcc                                                <  1:885360/69‑1 (MQ=255)

         tCGTAAGTGACGAAGTTGCCGTAACGCGCCACCATATAATCGACAATCTCTTTTTTACTTTTACCTTcc                                                <  1:2069626/69‑1 (MQ=255)

         tCGTAAGTGACGAAGTTGCCGTAACGCGCCACCATATAATCGACAATCTCTTTTTTACTTTTACCTTcc                                                <  1:295077/69‑1 (MQ=255)

         tCGTAAGTGACGAAGTTGCCGTAACGCGCCACCATATAATCGACAATCTCTTTTTTACTTTTACCTTcc                                                <  1:1916014/69‑1 (MQ=255)

         tCGTAAGTGACGAAGTTGCCGTAACGCGCCACCATATAATCGACAATCTCTTTTTTACTTTTACCTTcc                                                <  1:1184271/69‑1 (MQ=255)

         tCGTAAGTGACGAAGTTGCCGTAACGCGCCACCATATAATCGACAATCTCTTTTTTACTTTTACCTTcc                                                <  1:1590404/69‑1 (MQ=255)

                    gAAGTTGCCGTAACGCGCCACCATATAATCGACAATCTCTTTTTTACTTTTACCTTCCTGCATCAGtt                                      >  1:177699/1‑68 (MQ=255)

                    gAAGTTGCCGTAACGCGCCACCATATAATCGACAATCTCTTTTTTACTTTTACCTTCCTGCATCAGtt                                      >  1:1962300/1‑68 (MQ=255)

                    gAAGTTGCCGTAACGCGCCACCATATAATCGACAATCTCTTTCTTACTTTTGCCTTCCTGCATCAGtt                                      >  1:2570464/1‑68 (MQ=255)

                    gAAGTTGCCGTAACGCGCCACCATATAATCGACAATCTCTTTCTTACTTTTGCCTTCCTGCATCAGtt                                      >  1:603443/1‑68 (MQ=255)

                    gAAGTTGCCGTAACGCGCCACCATATAATCGACAATCTCTTTCTTACTTTTGCCTTCCTGCATCAGtt                                      >  1:2335753/1‑68 (MQ=255)

                    gAAGTTGCCGTAACGCGCCACCATATAATCGACAATCTCTTTCTTACTTTTGCCTTCCTGCATCAGtt                                      >  1:1092605/1‑68 (MQ=255)

                    gAAGTTGCCGTAACGCGCCACCATATAATCGACAATCTCTTTCTTACTTTTGCCTTCCTGCATCAGtt                                      >  1:31468/1‑68 (MQ=255)

                    gAAGTTGCCGTAACGCGCCACCATATAATCGACAATCTCTTTCTTACTTTTGCCTTCCTGCATCAGtt                                      >  1:1547760/1‑68 (MQ=255)

                    gAAGTTGCCGTAACGCGCCACCATATAATCGACAATCTCTTTCTTACTTTTGCCTTCCTGCATCAGtt                                      >  1:718581/1‑68 (MQ=255)

                    gAAGTTGCCGTAACGCGCCACCATATAATCGACAATCTCTTTCTTACTTTTGCCTTCCTGCATCAGtt                                      >  1:744771/1‑68 (MQ=255)

                    gAAGTTGCCGTAACGCGCCACCATATAATCGACAATCTCTTTCTTACTTTTGCCTTCCTGCATCAGtt                                      >  1:159982/1‑68 (MQ=255)

                     aaGTTGCCGTAACGCGCCACCATATAATCGACAATCTCTTTTTTACTTTTACCTTCCTGCATCAGTTCa                                    <  1:3035406/69‑1 (MQ=255)

                     aaGTTGCCGTAACGCGCCACCATATAATCGACAATCTCTTTTTTACTTTTACCTTCCTGCATCAGTTCa                                    <  1:1670858/69‑1 (MQ=255)

                     aaGTTGCCGTAACGCGCCACCATATAATCGACAATCTCTTTTTTACTTTTACCTTCCTGCATCAGTTCa                                    <  1:33829/69‑1 (MQ=255)

                     aaGTTGCCGTAACGCGCCACCATATAATCGACAATCTCTTTTTTACTTTTACCTTCCTGCATCAGTTCa                                    <  1:2303139/69‑1 (MQ=255)

                     aaGTTGCCGTAACGCGCCACCATATAATCGACAATCTCTTTTTTACTTTTACCTTCCTGCATCAGTTCa                                    <  1:543335/69‑1 (MQ=255)

                     aaGTTGCCGTAACGCGCCACCATATAATCGACAATCTCTTTTTTACTTTTACCTTCCTGCATCAGTTCa                                    <  1:2383421/69‑1 (MQ=255)

                     aaGTTGCCGTAACGCGCCACCATATAATCGACAATCTCTTTTTTACTTTTACCTTCCTGCATCAGTTCa                                    <  1:2059535/69‑1 (MQ=255)

                      aGTTGCCGTAACGCGCCACCATATAATCGACAATCTCTTTTTTACTTTTACCTTCCTGCATCAGTTCa                                    <  1:1514580/68‑1 (MQ=255)

                       gTTGCCGTAACGCGCCACCATATAATCGACAATCTCTTTTTTACTTTTACCTTCCTGCATCAGTTCAt                                   >  1:1045828/1‑68 (MQ=255)

                       gTTGCCGTAACGCGCCACCATATAATCGACAATCTCTTTTTTACTTTTACCTTCCTGCATCAGTTCATa                                  <  1:1001195/69‑1 (MQ=255)

                       gTTGCCGTAACGCGCCACCATATAATCGACAATCTCTTTTTTACTTTTACCTTCCTGCATCAGTTCATa                                  <  1:2888361/69‑1 (MQ=255)

                       gTTGCCGTAACGCGCCACCATATAATCGACAATCTCTTTTTTACTTTTACCTTCCTGCATCAGTTCATa                                  <  1:612217/69‑1 (MQ=255)

                        ttGCCGTAACGCGCCACCATATAATCGACAATCTCTTTTTTACTTTTACCTTCCTGCATCAGTTCATac                                 >  1:2757449/1‑69 (MQ=255)

                        ttGCCGTAACGCGCCACCATATAATCGACAATCTCTTTTTTACTTTTACCTTCCTGCATCAGTTCATac                                 >  1:2326432/1‑69 (MQ=255)

                         tGCCGTAACGCGCCACCATATAATCGACAATCTCTTTTTTACTTTTACCTTcc                                                <  1:2417185/53‑1 (MQ=255)

                         tGCCGTAACGCGCCACCATATAATCGACAATCTCTTTTTTACTTTTACCTTCCTGCATCAGTTCATaca                                <  1:3027505/69‑1 (MQ=255)

                         tGCCGTAACGCGCCACCATATAATCGACAATCTCTTTTTTACTTTTACCTTCCTGCATCAGTTCATaca                                <  1:2953822/69‑1 (MQ=255)

                         tGCCGTAACGCGCCACCATATAATCGACAATCTCTTTTTTACTTTTACCTTCCTGCATCAGTTCATaca                                <  1:93517/69‑1 (MQ=255)

                         tGCCGTAACGCGCCACCATATAATCGACAATCTCTTTTTTACTTTTACCTTCCTGCATCAGTTCATaca                                <  1:2550407/69‑1 (MQ=255)

                         tGCCGTAACGCGCCACCATATAATCGACAATCTCTTTTTTACTTTTACCTTCCTGCATCAGTTCATaca                                <  1:592354/69‑1 (MQ=255)

                         tGCCGTAACGCGCCACCATATAATCGACAATCTCTTTTTTACTTTTACCTTCCTGCATCAGTTCATaca                                <  1:2364670/69‑1 (MQ=255)

                         tGCCGTAACGCGCCACCATATAATCGACAATCTCTTTTTTACTTTTACCTTCCTGCATCAGTTCATaca                                <  1:2291835/69‑1 (MQ=255)

                         tGCCGTAACGCGCCACCATATAATCGACAATCTCTTTTTTACTTTTACCTTCCTGCATCAGTTCATaca                                <  1:2187697/69‑1 (MQ=255)

                         tGCCGTAACGCGCCACCATATAATCGACAATCTCTTTTTTACTTTTACCTTCCTGCATCAGTTCATaca                                <  1:635797/69‑1 (MQ=255)

                         tGCCGTAACGCGCCACCATATAATCGACAATCTCTTTTTTACTTTTACCTTCCTGCATCAGTTCATaca                                <  1:2017462/69‑1 (MQ=255)

                         tGCCGTAACGCGCCACCATATAATCGACAATCTCTTTTTTACTTTTACCTTCCTGCATCAGTTCATaca                                <  1:714494/69‑1 (MQ=255)

                         tGCCGTAACGCGCCACCATATAATCGACAATCTCTTTTTTACTTTTACCTTCCTGCATCAGTTCATaca                                <  1:1625201/69‑1 (MQ=255)

                         tGCCGTAACGCGCCACCATATAATCGACAATCTCTTTTTTACTTTTACCTTCCTGCATCAGTTCATaca                                <  1:1010229/69‑1 (MQ=255)

                         tGCCGTAACGCGCCACCATATAATCGACAATCTCTTTTTTACTTTTACCTTCCTGCATCAGTTCATaca                                <  1:109562/69‑1 (MQ=255)

                          gCCGTAACGCGCCACCATATAATCGACAATCTCTTTTTTACTTTTACCt                                                   <  1:328509/49‑1 (MQ=255)

                          gCCGTAACGCGCCACCATATAATCGACAATCTCTTTTTTACTTTTACCt                                                   <  1:2617668/49‑1 (MQ=255)

                          gCCGTAACGCGCCACCATATAATCGACAATCTCTTTTTTACTTTTACCt                                                   <  1:3153804/49‑1 (MQ=255)

                          gCCGTAACGCGCCACCATATAATCGACAATCTCTTTTTTACTTTTACCt                                                   <  1:890625/49‑1 (MQ=255)

                          gCCGTAACGCGCCACCATATAATCGACAATCTCTTTTTTACTTTTACCt                                                   <  1:1410134/49‑1 (MQ=255)

                          gCCGTAACGCGCCACCATATAATCGACAATCTCTTTTTTACTTTTACCt                                                   <  1:1516007/49‑1 (MQ=255)

                          gCCGTAACGCGCCACCATATAATCGACAATCTCTTTTTTACTTTTACCt                                                   <  1:352796/49‑1 (MQ=255)

                            cGTAACGCGCCACCATATAATCGACAATCTCTTttttactttta                                                      <  1:2376443/44‑1 (MQ=255)

                               aaCGCGCCACCATATAATCGACAATCTCTTTTTTACTTTTACCTTCCTGCATCAGTTCATACACTTTCt                          <  1:1633622/69‑1 (MQ=255)

                               aaCGCGCCACCATATAATCGACAATCTCTTTTTTACTTTTACCTTCCTGCATCAGTTCATACACTTTCt                          <  1:248604/69‑1 (MQ=255)

                                         caTATAATCGACAATCTCTTTTTTACTTTTACCTTCCTGCATCAGTTCATACACTTTCTGACGCAGGTc                <  1:1289300/69‑1 (MQ=255)

                                             taATCGACAATCTCTTTTTTACTTTTACCTTCCTGCATCAGTTCATACACTTTCTGACGCa                    >  1:2284661/1‑61 (MQ=255)

                                             taATCGACAATCTCTTTTTTACTTTTACCTTCCTGCATCAGTTCATACACTTTCTGACGCAGGTCGGt             >  1:2052239/1‑68 (MQ=255)

                                             taATCGACAATCTCTTTTTTACTTTTACCTTCCTGCATCAGTTCATACACTTTCTGACGCAGGTCGGt             >  1:1009535/1‑68 (MQ=255)

                                             taATCGACAATCTCTTTTTTACTTTTACCTTCCTGCATCAGTTCATACACTTTCTGACGCAGGTCGGt             >  1:88230/1‑68 (MQ=255)

                                             taATCGACAATCTCTTTTTTACTTTTACCTTCCTGCATCAGTTCATACACTTTCTGACGCAGGTCGGt             >  1:797635/1‑68 (MQ=255)

                                             taATCGACAATCTCTTTTTTACTTTTACCTTCCTGCATCAGTTCATACACTTTCTGACGCAGGTCGGt             >  1:3287238/1‑68 (MQ=255)

                                             taATCGACAATCTCTTTTTTACTTTTACCTTCCTGCATCAGTTCATACACTTTCTGACGCAGGTCGGt             >  1:615546/1‑68 (MQ=255)

                                             taATCGACAATCTCTTTTTTACTTTTACCTTCCTGCATCAGTTCATACACTTTCTGACGCAGGTCGGt             >  1:2154848/1‑68 (MQ=255)

                                                       tctcTTTCTTACTTTTGCCTTCCTGCATCAGTTCATACACTTTCTGGCGCAGGTCGGTGGCAATCATCg  <  1:1798670/69‑1 (MQ=255)

                                                       tctcTTTCTTACTTTTGCCTTCCTGCATCAGTTCATACACTTTCTGGCGCAGGTCGGTGGCAATCATCg  <  1:171895/69‑1 (MQ=255)

                                                       tctcTTTCTTACTTTTGCCTTCCTGCATCAGTTCATACACTTTCTGGCGCAGGTCGGTGGCAATCATCg  <  1:1702642/69‑1 (MQ=255)

                                                       tctcTTTCTTACTTTTGCCTTCCTGCATCAGTTCATACACTTTCTGGCGCAGGTCGGTGGCAATCATCg  <  1:1475009/69‑1 (MQ=255)

                                                       tctcTTTCTTACTTTTGCCTTCCTGCATCAGTTCATACACTTTCTGGCGCAGGTCGGTGGCAATCATCg  <  1:2556423/69‑1 (MQ=255)

                                                       tctcTTTCTTACTTTTGCCTTCCTGCATCAGTTCATACACTTTCTGGCGCAGGTCGGTGGCAATCATCg  <  1:2566374/69‑1 (MQ=255)

                                                       tctcTTTCTTACTTTTGCCTTCCTGCATCAGTTCATACACTTTCTGGCGCAGGTCGGTGGCAATCATCg  <  1:1331726/69‑1 (MQ=255)

                                                       tctcTTTCTTACTTTTGCCTTCCTGCATCAGTTCATACACTTTCTGGCGCAGGTCGGTGGCAATCATCg  <  1:2875215/69‑1 (MQ=255)

                                                        ctctTTCTTACTTTTGCCTTCCTGCATCAGTTCATACACTTTCTGGCGCAGGTCGGTGGCAATCATCg  <  1:2092013/68‑1 (MQ=255)

                                                        ctctTTCTTACTTTTGCCTTCCTGCATCAGTTCATACACTTTCTGGCGCAGGTCGGTGGCAATCATCg  <  1:2550091/68‑1 (MQ=255)

                                                        ctctTTCTTACTTTTGCCTTCCTGCATCAGTTCATACACTTTCTGGCGCAGGTCGGTGGCAATCATCg  <  1:2716300/68‑1 (MQ=255)

AACGGCGGATCGTAAGTGACGAAGTTGCCGTAACGCGCCACCATATAATCGACAATCTCTTTTTTACTTTTACCTTCCTGCATCAGTTCATACACTTTCTGACGCAGGTCGGTGGCAATCATCG  >  minE/1419673‑1419796

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