BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,419,744	A→G	21.9%	G73G (GGT→GGC)	ccmH ←	heme lyase, CcmH subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,419,744	0	A	G	21.9%	135.5 / 29.7	100	G73G (GGT→GGC)	ccmH	heme lyase, CcmH subunit
Reads supporting (aligned to +/- strand): ref base A (16/62); new base G (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 = 9.50e-01

Only 100 of 108 total aligned reads displayed.

GGCGGATCGTAAGTGACGAAGTTGCCGTAACGCGCCACCATATAATCGACAATCTCTTTTTTACTTTTACCTTCCTGCATCAGTTCATACACTTTCTGACGCAGGTCGGTGGCAATCATCGAGTTGGAATCGGC  >  minE/1419676‑1419809

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                  aaGTTGCCGTAACGCGCCACCATATAATCGACAATCTCTTTTTTACTTTTACCTTCCTGCATCAGTTCa                                                 <  1:682307/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:2888361/69‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                       gCCGTAACGCGCCACCATATAATCGACAATCTCTTTTTTACTTTTACCt                                                                <  1:2213507/49‑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:1497774/1‑68 (MQ=255)

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

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

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

                                          taATCGACAATCTCTTTTTTACTTTTACCTTCCTGCATCAGTTCATACACTTTCTGACGCAGGTCGGt                          >  1:3287238/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:2556423/69‑1 (MQ=255)

                                                    tctcTTTCTTACTTTTGCCTTCCTGCATCAGTTCATACACTTTCTGGCGCAGGTCGGTGGCAATCATCg               <  1:1475009/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:2550091/68‑1 (MQ=255)

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

                                                                 tttacCTTCCTGCATCAGTTCATACACTTTCTGACGCAGGTCGGTGGCAATCATCGAGTTGGAATCGGc  <  1:3389009/69‑1 (MQ=255)

GGCGGATCGTAAGTGACGAAGTTGCCGTAACGCGCCACCATATAATCGACAATCTCTTTTTTACTTTTACCTTCCTGCATCAGTTCATACACTTTCTGACGCAGGTCGGTGGCAATCATCGAGTTGGAATCGGC  >  minE/1419676‑1419809

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