BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,422,085	A→G	62.4%	R124R (CGT→CGC)	ccmF ←	heme lyase, CcmF subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,422,085	0	A	G	62.4%	34.5 / 35.9	61	R124R (CGT→CGC)	ccmF	heme lyase, CcmF subunit
Reads supporting (aligned to +/- strand): ref base A (15/8); new base G (22/16); total (37/24)
Fisher's exact test for biased strand distribution p-value = 6.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGAGGTAAAGAGAATGAACAGCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACACGGGCCACAATATCCAGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAG  >  minE/1422025‑1422142

agagGTAAAGAGAATGAACAGCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACACGGGCcac                                                   <  1:2625728/69‑1 (MQ=255)

agagGTAAAGAGAATGAACAGCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACACGGGCcac                                                   <  1:2261351/69‑1 (MQ=255)

      aaaGAGAATGAACAGCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACACGGGCCACAata                                               >  1:274328/1‑67 (MQ=255)

      aaaGAGAATGAACAGCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACACGGGCCACAata                                               >  1:2724758/1‑67 (MQ=255)

      aaaGAGAATGAACAGCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACACGGGCCACAata                                               >  1:2126540/1‑67 (MQ=255)

       aaGAGAATGAACAGCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACACGGGCCACAATATcc                                            <  1:991411/69‑1 (MQ=255)

       aaGAGAATGAACAGCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACACGGGCCACAATATcc                                            <  1:1439656/69‑1 (MQ=255)

              tGAACAGCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACACGGGCCACAATATCCAGCGGaa                                     <  1:1097097/69‑1 (MQ=255)

              tGAACAGCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACACGGGCCACAATATCCAGCGGaa                                     <  1:967848/69‑1 (MQ=255)

                aaCAGCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACACGGGCCACAATATCCAGCGGAAt                                    >  1:2094663/1‑68 (MQ=255)

                aaCAGCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACACGGGCCACAATATCCAGCGGAATa                                   >  1:2498122/1‑69 (MQ=255)

                aaCAGCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACACGGGCCACAATATCCAGCGGAATa                                   >  1:2456322/1‑69 (MQ=255)

                aaCAGCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACACGGGCCACAATATCCAGCGGAATa                                   >  1:2793764/1‑69 (MQ=255)

                aaCAGCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACACGGGCCACAATATCCAGCGGAATa                                   >  1:2258919/1‑69 (MQ=255)

                aaCAGCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACACGGGCCACAATATCCAGCGGAATa                                   >  1:2910830/1‑69 (MQ=255)

                aaCAGCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACACGGGCCACAATATCCAGCGGAATa                                   >  1:2004696/1‑69 (MQ=255)

                aaCAGCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACACGGGCCACAATATCCAGCGGAATa                                   >  1:3115248/1‑69 (MQ=255)

                aaCAGCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACACGGGCCACAATATCCAGCGGAATa                                   >  1:3168609/1‑69 (MQ=255)

                    gCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCt                               <  1:2016442/69‑1 (MQ=255)

                    gCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCt                               <  1:3244380/69‑1 (MQ=255)

                    gCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCt                               <  1:1328176/69‑1 (MQ=255)

                    gCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCt                               <  1:330218/69‑1 (MQ=255)

                    gCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCt                               <  1:2845568/69‑1 (MQ=255)

                    gCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCt                               <  1:3394754/69‑1 (MQ=255)

                    gCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCt                               <  1:615369/69‑1 (MQ=255)

                    gCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCt                               <  1:2280282/69‑1 (MQ=255)

                    gCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCt                               <  1:2763876/69‑1 (MQ=255)

                    gCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCt                               <  1:1440724/69‑1 (MQ=255)

                    gCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCt                               <  1:274988/69‑1 (MQ=255)

                    gCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGAACGCGGGCAACAATATCCAGCGGAATACGCt                               <  1:1419468/69‑1 (MQ=255)

                     cAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCt                               <  1:197918/68‑1 (MQ=255)

                      aaaaaGCCGACACTGACCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCt                               <  1:849496/67‑1 (MQ=255)

                              gACACTGACCATCCCCATTATCGCCAGTACACGGTCCACAATATCCAGCGGaa                                     <  1:36728/53‑1 (MQ=255)

                               acacTGACCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCt                               <  1:2485649/58‑1 (MQ=255)

                                   tGACCATCCCCATTATCGCCAGTACACGGGCCACAATATCCAGCGGAATACGCTGACTAAAAGTCGCCa                <  1:440148/69‑1 (MQ=255)

                                     aCCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCTGACTAAAAATCGCTAc               >  1:1186313/1‑68 (MQ=255)

                                     aCCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCTGACTAAAAATCGCTAc               >  1:1434402/1‑68 (MQ=255)

                                     aCCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCTGACTAAAAATCGCTAc               >  1:1500279/1‑68 (MQ=255)

                                     aCCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCTGACTAAAAATCGCTAc               >  1:903492/1‑68 (MQ=255)

                                     aCCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCTGACTAAAAATCGCTAc               >  1:807859/1‑68 (MQ=255)

                                     aCCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCTGACTAAAAATCGCTAc               >  1:691688/1‑68 (MQ=255)

                                     aCCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCTGACTAAAAATCGCTAc               >  1:617359/1‑68 (MQ=255)

                                     aCCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCTGACTAAAAATCGCTAc               >  1:1717680/1‑68 (MQ=255)

                                     aCCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCTGACTAAAAATCGCTAc               >  1:377029/1‑68 (MQ=255)

                                     aCCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCTGACTAAAAATCGCTAc               >  1:1745356/1‑68 (MQ=255)

                                     aCCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCTGACTAAAAATCGCTAc               >  1:1792434/1‑68 (MQ=255)

                                     aCCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCTGACTAAAAATCGCTAc               >  1:1891646/1‑68 (MQ=255)

                                     aCCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCTGACTAAAAATCGCTAc               >  1:29065/1‑68 (MQ=255)

                                     aCCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCTGACTAAAAATCGCTAc               >  1:2898629/1‑68 (MQ=255)

                                     aCCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCTGACTAAAAATCGCTAc               >  1:2863371/1‑68 (MQ=255)

                                     aCCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCTGACTAAAAATCGCTAc               >  1:2240015/1‑68 (MQ=255)

                                     aCCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCTGACTAAAAATCGCTAc               >  1:1288206/1‑68 (MQ=255)

                                     aCCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCTGACTAAAAATCGCTAc               >  1:2338288/1‑68 (MQ=255)

                                     aCCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCTGACTAAAAATCGCTAc               >  1:2760495/1‑68 (MQ=255)

                                     aCCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCTGACTAAAAATCGCTAc               >  1:1165456/1‑68 (MQ=255)

                                     aCCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCTGACTAAAAATCGCTAc               >  1:1240805/1‑68 (MQ=255)

                                     aCCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCTGACTAAAAACCGCTAc               >  1:36191/1‑68 (MQ=255)

                                      ccATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCt                               <  1:1897217/51‑1 (MQ=255)

                                           cccATTATCGCCAGTACACGGGCCACAATATCCAGCGGAATACGCTGACTAAAAATCGCCACCGCAAAg        >  1:2183644/1‑69 (MQ=255)

                                                 aTCGCCAGTACACGGGCCACAATATCCAGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAg  >  1:2787557/1‑69 (MQ=255)

                                                 aTCGCCAGTACACGGGCCACAATATCCAGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAg  >  1:2777035/1‑69 (MQ=255)

AGAGGTAAAGAGAATGAACAGCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACACGGGCCACAATATCCAGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAG  >  minE/1422025‑1422142

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 14 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`