BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,243,253	G→T	20.4%	V115V (GTC→GTA)	nfi ←	endonuclease V

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,243,253	0	G	T	20.4%	184.2 / 31.5	93	V115V (GTC→GTA)	nfi	endonuclease V
Reads supporting (aligned to +/- strand): ref base G (59/15); new base T (12/7); total (71/22)
Fisher's exact test for biased strand distribution p-value = 1.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

Only 100 of 106 total aligned reads displayed.

AGAGCCGTTTTTTCGCCACGCCAATGGTCGGCACATCCACCAATAAGCCAAAATGGCTGGCGACGCCAAGACGGCGAGGATGCGAGATCCCATGACCATCGACAAACACTAAATCCGGCTTTTGCGACAGCA  >  minE/2243192‑2243323

agagCCGTTTTTTCGCCACGCCAATGGTCGGCACATCCACCAATAAGCCAAAATGTCTGGCGACGCCaa                                                                 <  1:2285851/69‑1 (MQ=255)

         ttttCGCTACGCCAATGGTCGGCACATCCACCAATAAGCCAAAATGGCTGGCTACGCCAAGACGGCGAg                                                        >  1:1517756/1‑69 (MQ=255)

         ttttCGCTACGCCAATGGTCGGCACATCCACCAATAAGCCAAAATGGCTGGCTACGCCAAGACGGCGAg                                                        >  1:2183922/1‑69 (MQ=255)

         ttttCGCTACGCCAATGGTCGGCACATCCACCAATAAGCCAAAATGGCTGGCTACGCCAAGACGGCGAg                                                        >  1:2180075/1‑69 (MQ=255)

         ttttCGCTACGCCAATGGTCGGCACATCCACCAATAAGCCAAAATGGCTGGCTACGCCAAGACGGCGAg                                                        >  1:2368963/1‑69 (MQ=255)

         ttttCGCTACGCCAATGGTCGGCACATCCACCAATAAGCCAAAATGGCTGGCTACGCCAAGACGGCGAg                                                        >  1:2504255/1‑69 (MQ=255)

         ttttCGCTACGCCAATGGTCGGCACATCCACCAATAAGCCAAAATGGCTGGCTACGCCAAGACGGCGAg                                                        >  1:1807434/1‑69 (MQ=255)

         ttttCGCTACGCCAATGGTCGGCACATCCACCAATAAGCCAAAATGGCTGGCTACGCCAAGACGGCGAg                                                        >  1:3163085/1‑69 (MQ=255)

         ttttCGCTACGCCAATGGTCGGCACATCCACCAATAAGCCAAAATGGCTGGCTACGCCAAGACGGCGAg                                                        >  1:1329473/1‑69 (MQ=255)

         ttttCGCTACGCCAATGGTCGGCACATCCACCAATAAGCCAAAATGGCTGGCTACGCCAAGACGGCGAg                                                        >  1:720447/1‑69 (MQ=255)

         ttttCGCTACGCCAATGGTCGGCACATCCACCAATAAGCCAAAATGGCTGGCTACGCCAAGACGGCGAg                                                        >  1:378381/1‑69 (MQ=255)

         ttttCGCTACGCCAATGGTCGGCACATCCACCAATAAGCCAAAATGGCTGGCTACGCCAAGACGGCGAg                                                        >  1:684733/1‑69 (MQ=255)

         ttttCGCTACGCCAATGGTCGGCACATCCACCAATAAGCCAAAATGGCTGGCTACGCCAAGACGGCGAg                                                        >  1:39435/1‑69 (MQ=255)

         ttttCGCCACGCCAATGGTCGGCACATCCACCAATAAGCCAAAATGGCTGGCGACGCCAAGACGGCGAg                                                        >  1:2849295/1‑69 (MQ=255)

                    ccaATGGTCGGCACATCCACCAATAAGCCAAAATGGCTGGCGACGCCAAGACGGCGAGGATGCgaga                                               >  1:2823557/1‑67 (MQ=255)

                    ccaATGGTCGGCACATCCACCAATAAGCCAAAATGGCTGGCGACGCCAAGACGGCGAGGATGCGAGAt                                              >  1:90508/1‑68 (MQ=255)

                    ccaATGGTCGGCACATCCACCAATAAGCCAAAATGGCTGGCGACGCCAAGACGGCGAGGATGCGAGAt                                              >  1:2207627/1‑68 (MQ=255)

                    ccaATGGTCGGCACATCCACCAATAAGCCAAAATGGCTGGCGACGCCAAGACGGCGAGGATGCGAGAt                                              >  1:1980720/1‑68 (MQ=255)

                    ccaATGGTCGGCACATCCACCAATAAGCCAAAATGGCTGGCGACGCCAAGACGGCGAGGATGCGAGAt                                              >  1:594958/1‑68 (MQ=255)

                    ccaATGGTCGGCACATCCACCAATAAGCCAAAATGGCTGGCGACGCCAAGACGGCGAGGATGCGAGAt                                              >  1:537408/1‑68 (MQ=255)

                     caATGGTCGGCACATCCACCAATAAGCCAAAATGGCTGGCGACGCCAAGACGGCGAGGATGCGAGATcc                                            <  1:1324130/69‑1 (MQ=255)

                     caATGGTCGGCACATCCACCAATAAGCCAAAATGGCTGGCGACGCCAAGACGGCGAGGATGCGAGATcc                                            <  1:927123/69‑1 (MQ=255)

                     caATGGTCGGCACATCCACCAATAAGCCAAAATGGCTGGCGACGCCAAGACGGCGAGGATGCGAGATcc                                            <  1:239255/69‑1 (MQ=255)

                     caATGGGCGGCACATCCACCAATAAGCCAAAATGGCTGGCGACGCCAAGACGGCGAGGATGCGAGATcc                                            <  1:921436/69‑1 (MQ=255)

                         ggTCGGCACATCCACCAATAAGCCAAAATGGCTGGCGACGCCAAGACGGCGAg                                                        <  1:1742736/53‑1 (MQ=255)

                         ggTCGGCACATCCACCAATAAGCCAAAATGGCTGGCGACGCCAAGACGGCGAg                                                        <  1:11153/53‑1 (MQ=255)

                          gTCGGCACATCCACCAATAAGCCAAAATGGCTGGCGACGCCAAGACGGCGAGGATGCGAGATCCCATg                                        >  1:1380516/1‑68 (MQ=255)

                          gTCGGCACATCCACCAATAAGCCAAAATGGCTGGCGACGCCAAGACGGCGAGGATGCGAGATCCCATg                                        >  1:3166337/1‑68 (MQ=255)

                            cGGCACATCCACCAATAAGCCAAAATGGCTGGCGACGCCAAGACGGCGAGGATGCGAGATCCCATGAcc                                     <  1:1696251/69‑1 (MQ=255)

                            cGGCACATCCACCAATAAGCCAAAATGGCTGGCGACGCCAAGACGGCGAGGATGCGAGATCCCATGAcc                                     <  1:2971893/69‑1 (MQ=255)

                            cGGCACATCCACCAATAAGCCAAAATGGCTGGCGACGCCAAGACGGCGAGGATGCGAGATCCCATGAcc                                     <  1:2042841/69‑1 (MQ=255)

                             ggCACATCCACCAATAAGCCAAAATGGCTGGCGACGCCAAGACGGCGAGGATGCGAGATCCCATGAcc                                     <  1:1592650/68‑1 (MQ=255)

                             ggCACATCCACCAATAAGCCAAAATGGCTGGCGACGCCAAGACGGCGAGGATGCGAGATCCCATGACCa                                    >  1:2786239/1‑69 (MQ=255)

                             ggCACATCCACCAATAAGCCAAAATGGCTGGCGACGCCAAGACGGCGAGGATGCGAGATCCCATGACCa                                    >  1:2476185/1‑69 (MQ=255)

                             ggCACATCCACCAATAAGCCAAAATGGCTGGCGACGCCAAGACGGCGAGGATGCGAGATCCCATGACCa                                    >  1:2401627/1‑69 (MQ=255)

                             ggCACATCCACCAATAAGCCAAAATGGCTGGCGACGCCAAGACGGCGAGGATGCGAGATCCCATGACCa                                    >  1:770578/1‑69 (MQ=255)

                             ggCACATCCACCAATAAGCCAAAATGGCTGGCGACGCCAAGACGGCGAGGATGCGAGATCCCATGACCa                                    >  1:1929743/1‑69 (MQ=255)

                             ggCACATCCACCAATAAGCCAAAATGGCTGGCGACGCCAAGACGGCGAGGATGCGAGATCCCATGACCa                                    >  1:1918851/1‑69 (MQ=255)

                                   tCCACCAATAAGCCAAAATGGCTGGCGACGCCAAGACGGCGAGGATGCGAGATCCCATGACCATCGAc                               >  1:2487373/1‑68 (MQ=255)

                                   tCCACCAATAAGCCAAAATGGCTGGCGACGCCAAGACGGCGAGGATGCGAGATCCCATGACCATCGACa                              >  1:1078969/1‑69 (MQ=255)

                                   tCCACCAATAAGCCAAAATGGCTGGCGACGCCAAGACGGCGAGGATGCGAGATCCCATGACCATCGACa                              >  1:139136/1‑69 (MQ=255)

                                   tCCACCAATAAGCCAAAATGGCTGGCGACGCCAAGACGGCGAGGATGCGAGATCCCATGACCATCGACa                              >  1:1308496/1‑69 (MQ=255)

                                   tCCACCAATAAGCCAAAATGGCTGGCGACGCCAAGACGGCGAGGATGCGAGATCCCATGACCATCGACa                              >  1:3288353/1‑69 (MQ=255)

                                   tCCACCAATAAGCCAAAATGGCTGGCGACGCCAAGACGGCGAGGATGCGAGATCCCATGACCATCGACa                              >  1:1125487/1‑69 (MQ=255)

                                   tCCACCAATAAGCCAAAATGGCTGGCGACGCCAAGACGGCGAGGATGCGAGATCCCATGACCATCGACa                              >  1:505167/1‑69 (MQ=255)

                                   tCCACCAATAAGCCAAAATGGCTGGCGACGCCAAGACGGCGAGGATGCGAGATCCCATGACCATCGACa                              >  1:580586/1‑69 (MQ=255)

                                   tCCACCAATAAGCCAAAATGGCTGGCGACGCCAAGACGGCGAGGATGCGAGATCCCATGACCATCGACa                              >  1:582878/1‑69 (MQ=255)

                                   tCCACCAATAAGCCAAAATGGCTGGCGACGCCAAGACGGCGAGGATGCGAGATCCCATGACCATCGACa                              <  1:3118133/69‑1 (MQ=255)

                                   tCCACCAATAAGCCAAAATGGCTGGCGACGCCAAGACGGCGAGGATGCGAGATCCCATGACCATCGACa                              >  1:612978/1‑69 (MQ=255)

                                   tCCACCAATAAGCCAAAATGGCTGGCGACGCCAAGACGGCGAGGATGCGAGATCCCATGACCATCGACa                              >  1:663617/1‑69 (MQ=255)

                                   tCCACCAATAAGCCAAAATGGCTGGCGACGCCAAGACGGCGAGGATGCGAGATCCCATGACCATCGACa                              >  1:666399/1‑69 (MQ=255)

                                   tCCACCAATAAGCCAAAATGGCTGGCGACGCCAAGACGGCGAGGATGCGAGATCCCATGACCATCGACa                              >  1:715394/1‑69 (MQ=255)

                                   tCCACCAATAAGCCAAAATGGCTGGCGACGCCAAGACGGCGAGGATGCGAGATCCCATGACCATCGACa                              >  1:95569/1‑69 (MQ=255)

                                   tCCACCAATAAGCCAAAATGGCTGGCGACGCCAAGACGGCGAGGATGCGAGATCCCATGACCATCGACa                              >  1:1003195/1‑69 (MQ=255)

                                   tCCACCAATAAGCCAAAATGGCTGGCGACGCCAAGACGGCGAGGATGCGAGATCCCATGACCATCGACa                              >  1:2358211/1‑69 (MQ=255)

                                   tCCACCAATAAGCCAAAATGGCTGGCGACGCCAAGACGGCGAGGATGCGAGATCCCATGACCATCGACa                              >  1:2362780/1‑69 (MQ=255)

                                   tCCACCAATAAGCCAAAATGGCTGGCGACGCCAAGACGGCGAGGATGCGAGATCCCATGACCATCGACa                              >  1:827334/1‑69 (MQ=255)

                                   tCCACCAATAAGCCAAAATGGCTGGCGACGCCAAGACGGCGAGGATGCGAGATCCCATGACCATCGACa                              >  1:2745890/1‑69 (MQ=255)

                                   tCCACCAATAAGCCAAAATGGCTGGCGACGCCAAGACGGCGAGGATGCGAGATCCCATGACCATCGACa                              >  1:2735728/1‑69 (MQ=255)

                                   tCCACCAATAAGCCAAAATGGCTGGCGACGCCAAGACGGCGAGGATGCGAGATCCCATGACCATCGACa                              >  1:1851237/1‑69 (MQ=255)

                                     caccaATAAGCCAAAATGGCTGGCGACGCCAAGACGGCGAGGATGCGAGATCCCATGACCATCGACaaa                            >  1:1048473/1‑69 (MQ=255)

                                          aTAAGCCAAAATGGCTGGCTACGCCAAGACGGCGAGGATGCGAGATCCCATGACCATCGACAAACACTa                       <  1:3361888/69‑1 (MQ=255)

                                          aTAAGCCAAAATGGCTGGCTACGCCAAGACGGCGAGGATGCGAGATCCCATGACCATCGACAAACACTa                       <  1:3244259/69‑1 (MQ=255)

                                          aTAAGCCAAAATGGCTGGCTACGCCAAGACGGCGAGGATGCGAGATCCCATGACCATCGACAAACACTa                       <  1:2022591/69‑1 (MQ=255)

                                          aTAAGCCAAAATGGCTGGCTACGCCAAGACGGCGAGGATGCGAGATCCCATGACCATCGACAAACACTa                       <  1:645879/69‑1 (MQ=255)

                                           tAAGCCAAAATGGCTGGCTACGCCAAGACGGCGAGGATGCGAGATCCCATGACCATCGACAAACACTa                       <  1:812948/68‑1 (MQ=255)

                                           tAAGCCAAAATGGCTGGCTACGCCAAGACGGCGAGGATGCGAGATCCCATGACCATCGACAAACACTa                       <  1:719046/68‑1 (MQ=255)

                                            aaGCCAAAATGGCTGGCGACTCCAAGACGGCGAGGATGCGAGATCCCATGACCATCGACAAACACTaaa                     >  1:3380161/1‑69 (MQ=255)

                                            aaGCCAAAATGGCTGGCGACGCCAAGACGGCGAGGATGCGAGAt                                              >  1:2749753/1‑44 (MQ=255)

                                            aaGCCAAAATGGCTGGCGACGCCAAGACGGCGAGGATGCGAGATCCCATGACCATCGACAAACACTaaa                     >  1:3426960/1‑69 (MQ=255)

                                            aaGCCAAAATGGCTGGCGACGCCAAGACGGCGAGGATGCGAGATCCCATGACCATCGACAAACACTaaa                     >  1:1785852/1‑69 (MQ=255)

                                            aaGCCAAAATGGCTGGCGACGCCAAGACGGCGAGGATGCGAGATCCCATGACCATCGACAAACACTaaa                     >  1:2868714/1‑69 (MQ=255)

                                             aGCCAAAATGGCTGGCGACGCCAAGACGGCGAGGATGCGAGATCCCATGACCATCGa                                >  1:2513095/1‑57 (MQ=255)

                                             aGCCAAAATGGCTGGCGACGCCAAGACGGCGAGGATGCGAGATCCCATGACCATCGACAAACACTAAAt                    >  1:3355141/1‑69 (MQ=255)

                                             aGCCAAAATGGCTGGCGACGCCAAGACGGCGAGGATGCGAGATCCCATGACCATCGACAAACACTAAAt                    >  1:2268799/1‑69 (MQ=255)

                                             aGCCAAAATGGCTGGCGACGCCAAGACGGCGAGGATGCGAGATCCCATGACCATCGACAAACACTAAAt                    >  1:875840/1‑69 (MQ=255)

                                             aGCCAAAATGGCTGGCGACGCCAAGACGGCGAGGATGCGAGATCCCATGACCATCGACAAACACTAAAt                    >  1:3060549/1‑69 (MQ=255)

                                             aGCCAAAATGGCTGGCGACGCCAAGACGGCGAGGATGCGAGATCCCATGACCATCGACAAACACTAAAt                    >  1:450475/1‑69 (MQ=255)

                                             aGCCAAAATGGCTGGCGACGCCAAGACGGCGAGGATGCGAGATCCCATGACCATCGACAAACACTAAAt                    >  1:408900/1‑69 (MQ=255)

                                             aGCCAAAATGGCTGGCGACGCCAAGACGGCGAGGATGCGAGATCCCATGACCATCGACAAACAATAAAt                    >  1:1370937/1‑69 (MQ=255)

                                                       gctggcGACGCCAAGACGGCGAGGATGCGAGATCCCa                                          >  1:236864/1‑37 (MQ=255)

                                                       gctggcGACGCCAAGACGGCGAGGATGCGAGATCCCATGACCATCGACAAACACTAAATCCGGCtttt           >  1:1101245/1‑68 (MQ=255)

                                                       gctggcGACGCCAAGACGGCGAGGATGCGAGATCCCATGACCATCGACAAACACTAAATCCGGCtttt           >  1:1432618/1‑68 (MQ=255)

                                                       gctggcGACGCCAAGACGGCGAGGATGCGAGATCCCATGACCATCGACAAACACTAAATCCGGCtttt           >  1:1565444/1‑68 (MQ=255)

                                                       gctggcGACGCCAAGACGGCGAGGATGCGAGATCCCATGACCATCGACAAACACTAAATCCGGCtttt           >  1:1569654/1‑68 (MQ=255)

                                                       gctggcGACGCCAAGACGGCGAGGATGCGAGATCCCATGACCATCGACAAACACTAAATCCGGCttt            >  1:1909862/1‑67 (MQ=255)

                                                       gctggcGACGCCAAGACGGCGAGGATGCGAGATCCCATGACCATCGACAAACACTAAATCCGGCttt            >  1:1272752/1‑67 (MQ=255)

                                                         tggcGACGCCAAGACGGCGAGGATGCGAGATCCCATGACCATCGACAAACACTAAATCCGGCTTTTGc         <  1:385079/68‑1 (MQ=255)

                                                             gACGCCAAGACGGCGAGGATGCGAGATCCCATGACCATCGACAAACACTAAGTCCGGCTTTTGCGACAg    >  1:3369212/1‑69 (MQ=255)

                                                             gACGCCAAGACGGCGAGGATGCGAGATCCCATGACCATCGACAAACACTAAATCCGGCTTTTGCGACAg    >  1:2553227/1‑69 (MQ=255)

                                                             gACGCCAAGACGGCGAGGATGCGAGATCCCATGACCATCGACAAACACTAAATCCGGCTTTTGCGACAg    >  1:834751/1‑69 (MQ=255)

                                                             gACGCCAAGACGGCGAGGATGCGAGATCCCATGACCATCGACAAACACTAAATCCGGCTTTTGCGACAg    >  1:3038649/1‑69 (MQ=255)

                                                             gACGCCAAGACGGCGAGGATGCGAGATCCCATGACCATCGACAAACACTAAATCCGGCTTTTGCGACAg    <  1:2993493/69‑1 (MQ=255)

                                                             gACGCCAAGACGGCGAGGATGCGAGATCCCATGACCATCGACAAACACTAAATCCGGCTTTTGCGACAg    >  1:2652114/1‑69 (MQ=255)

                                                             gACGCCAAGACGGCGAGGATGCGAGATCCCATGACCATCGACAAACACTAAATCCGGCTTTTGCGACAg    <  1:2253008/69‑1 (MQ=255)

                                                             gACGCCAAGACGGCGAGGATGCGAGATCCCATGACCATCGACAAACACTAAATCCGGCTTTTGCGACAg    >  1:1875052/1‑69 (MQ=255)

                                                             gACGCCAAGACGGCGAGGATGCGAGATCCCATGACCATCGACAAACACTAAATCCGGCTTTTGCGACAg    <  1:1377269/69‑1 (MQ=255)

                                                             gACGCCAAGACGGCGAGGATGCGAGATCCCATGACCATCGACAAACACTAAATCCGGCTTTTGCGACAg    >  1:1315605/1‑69 (MQ=255)

                                                             gACGCCAAGACGGCGAGGATGCGAGATCCCATGACCATCGACAAACACTAAATCCGGCTTTTGCGACAg    >  1:1284661/1‑69 (MQ=255)

                                                             gACGCCAAGACGGCGAGGATGCGAGATCCCATGACCATCGACAAACACTAAATCCGGCTTTTGCGACAg    <  1:1018133/69‑1 (MQ=255)

AGAGCCGTTTTTTCGCCACGCCAATGGTCGGCACATCCACCAATAAGCCAAAATGGCTGGCGACGCCAAGACGGCGAGGATGCGAGATCCCATGACCATCGACAAACACTAAATCCGGCTTTTGCGACAGCA  >  minE/2243192‑2243323

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