BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	649,288	A→G	48.4%	L656L (CTA→CTG)	pepN →	aminopeptidase N

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	649,288	0	A	G	48.4%	‑3.2 / 54.8	62	L656L (CTA→CTG)	pepN	aminopeptidase N
Reads supporting (aligned to +/- strand): ref base A (12/20); new base G (14/16); total (26/36)
Fisher's exact test for biased strand distribution p-value = 6.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCGGATGAGCTACTGGCTATTTACAACGCGAATTACCAGAGCGAGTACCGTGTTGAGCATGAAGATATTGCAAAACG  >  minE/649220‑649353

ccGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCGGATGAGCTa                                                                   >  1:3089956/1‑69 (MQ=255)

ccGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCGGATGAGCTa                                                                   <  1:1787745/69‑1 (MQ=255)

ccGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCGGATGAGCTa                                                                   >  1:1384253/1‑69 (MQ=255)

ccGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCGGATGAGCTa                                                                   >  1:1256431/1‑69 (MQ=255)

ccGATTGCTATTGCCGAAGTACGCGAAGCACTCAATCGTACTCTGGCGACTGAACTGGCGGATGAGCTa                                                                   >  1:1746206/1‑69 (MQ=255)

            gCCGAAGTACGCGAAGCACTCCCTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTAc                                                       >  1:2375520/1‑69 (MQ=255)

            gCCGAAGTACGCGAAGCACTCACTCGTACTCTGTCGACTGAACTGGCAGATGAGCTGCTGGCTATTTAc                                                       >  1:777253/1‑69 (MQ=255)

            gCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTAc                                                       >  1:1081225/1‑69 (MQ=255)

            gCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTAc                                                       >  1:2371486/1‑69 (MQ=255)

            gCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTAc                                                       >  1:1781457/1‑69 (MQ=255)

            gCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTAc                                                       >  1:511292/1‑69 (MQ=255)

            gCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTAc                                                       >  1:1731785/1‑69 (MQ=255)

            gCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTAc                                                       >  1:1687972/1‑69 (MQ=255)

            gCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTAc                                                       >  1:2697056/1‑69 (MQ=255)

            gCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTAc                                                       >  1:2700809/1‑69 (MQ=255)

            gCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTAc                                                       >  1:1939681/1‑69 (MQ=255)

            gCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTAc                                                       >  1:3314504/1‑69 (MQ=255)

            gCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTAc                                                       >  1:2929220/1‑69 (MQ=255)

            gCCGAAGTACGCGAAGCACTAACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTAc                                                       >  1:1170248/1‑69 (MQ=255)

               gAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCGGATGAGCTACTGGCTATTTACAAc                                                    <  1:1209284/69‑1 (MQ=255)

                  gTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCGGATGAGCTACTGGCTATTTACAAcgcg                                                 <  1:1241995/69‑1 (MQ=255)

                  gTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCGGATGAGCTACTGGCTATTTACAAcgcg                                                 <  1:72064/69‑1 (MQ=255)

                     cgcgAAGCACTCACTCGTACTCTGGCGACTGAACTGGCGGATGAGCTACTGGCTATTTACAACGCGAAt                                              >  1:2768018/1‑69 (MQ=255)

                     cgcgAAGCACTCACTCGTACTCTGGCGACTGAACTGGCGGATGAGCTACTGGCTATTTACAACGCGAAt                                              >  1:2025204/1‑69 (MQ=255)

                        gAAGCACTCACTCGTACTCTGGCGACTGAACTGGCGGATGAGCTACTGGCTATTTACAACGCGAATTAc                                           <  1:3148309/69‑1 (MQ=255)

                           gCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTACAACGCGAATTACCag                                        <  1:2122453/69‑1 (MQ=255)

                           gCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTACAACGCGAATTACCag                                        <  1:3156675/69‑1 (MQ=255)

                           gCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTACAACGCGAATTACCag                                        <  1:2879092/69‑1 (MQ=255)

                           gCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTACAACGCGAATTACCag                                        <  1:2822558/69‑1 (MQ=255)

                           gCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTACAACGCGAATTACCag                                        <  1:27993/69‑1 (MQ=255)

                           gCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTACAACGCGAATTACCag                                        <  1:3154873/69‑1 (MQ=255)

                           gCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTACAACGCGAATTACCag                                        <  1:411763/69‑1 (MQ=255)

                           gCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTACAACGCGAATTACCag                                        <  1:2684228/69‑1 (MQ=255)

                           gCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTACAACGCGAATTACCag                                        <  1:2572033/69‑1 (MQ=255)

                           gCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTACAACGCGAATTACCag                                        <  1:742896/69‑1 (MQ=255)

                           gCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTACAACGCGAATTACCag                                        <  1:529008/69‑1 (MQ=255)

                           gCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTACAACGCGAATTACCag                                        <  1:1000810/69‑1 (MQ=255)

                           gCACTCACTCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTACAACGCGAATTACCag                                        <  1:700612/69‑1 (MQ=255)

                            cactcactCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTACAACGCGAATTACCag                                        <  1:2602838/68‑1 (MQ=255)

                            cactcactCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTACAACGCGAATTACCag                                        <  1:1350446/68‑1 (MQ=255)

                            cactcactCGTACTCTGGCGACTGAACTGGCAGATGAGCTGCTGGCTATTTACAACGCGAATTACCag                                        <  1:2983182/68‑1 (MQ=255)

                                   tcGTACTCTGGCGACTGAACTGGCGGATGAGCTACTGGTTATTTACAACGCGAATTACCAGAGCGAGTa                                <  1:1477375/69‑1 (MQ=255)

                                   tcGTACTCTGGCGACTGAACTGGCGGATGAGCTACTGGCTATTTACAACGCGAATTACCAGAGCGAGTa                                <  1:2162719/69‑1 (MQ=255)

                                   tcGTACTCTGGCGACTGAACTGGCGGATGAGCTACTGGCTATTTACAACGCGAATTACCAGAGCGAGTa                                <  1:1304859/69‑1 (MQ=255)

                                   tcGTACTCTGGCGACTGAACTGGCGGATGAGCTACTGGCTATTTACAACGCGAATTACCAGAGCGAGTa                                <  1:1753305/69‑1 (MQ=255)

                                    cGTACTCTGGCGACTGAACTGGCGGATGAGCTACTGGCTATTTACAACGCGAATTACCAGAGCGAGTa                                >  1:1364016/1‑68 (MQ=255)

                                    cGTACTCTGGCGACTGAACTGGCGGATGAGCTACTGGCTATTTACAACGCGAATTACCAGAGCGAGTAc                               >  1:890897/1‑69 (MQ=255)

                                    cGTACTCTGGCGACTGAACTGGCGGATGAGCTACTGGCTATTTACAACGCGAATTACCAGAGCGAGTAc                               >  1:3324481/1‑69 (MQ=255)

                                    cGTACTCTGGCGACTGAACTGGCGGATGAGCTACTGGCTATTTACAACGCGAATTACCAGAGCGAGTAc                               >  1:954750/1‑69 (MQ=255)

                                    cGTACTCTGGCGACTGAACTGGCGGATGAGCTACTGGCTATTTACAACGCGAATTACCAGAGCGAGTAc                               >  1:2700229/1‑69 (MQ=255)

                                    cGTACTCTGGCGACTGAACTGGCGGATGAGCTACTGGCTATTTACAACGCGAATTACCAGAGCGAGTAc                               >  1:2444096/1‑69 (MQ=255)

                                    cGTACTCTGGCGACTGAACTGGCGGATGAGCTACTGGCTATTTACAACGCGAATTACCAGAGCGAGTAc                               >  1:1752492/1‑69 (MQ=255)

                                                      cTGGCGGATGAGCTACTGGCTATTTACAACGCGAATTACCAGAGCGAGTACCGTGTTGAGCATGAAGat             >  1:2892202/1‑69 (MQ=255)

                                                      cTGGCGGATGAGCTACTGGCTATTTACAACGCGAATTACCAGAGCGAGTACCGTGTTGAGCATGAAGat             >  1:2880047/1‑69 (MQ=255)

                                                      cTGGCGGATGAGCTACTGGCTATTTACAACGCGAATTACCAGAGCGAGTACCGTGTTGAGCATGAAGat             >  1:2843020/1‑69 (MQ=255)

                                                      cTGGCGGATGAGCTACTGGCTATTTACAACGCGAATTACCAGAGCGAGTACCGTGTTGAGCATGAAGa              <  1:1300499/68‑1 (MQ=255)

                                                                  cTACTGGCTATTTACAACGCGAATTACCAGAGCGAGTACCGTGTTGAGCATGAAGATATTGCAAAACg  <  1:1073718/68‑1 (MQ=255)

                                                                  cTACTGGCTATTTACAACGCGAATTACCAGAGCGAGTACCGTGTTGAGCATGAAGATATTGCAAAACg  <  1:3002906/68‑1 (MQ=255)

                                                                  cTACTGGCTATTTACAACGCGAATTACCAGAGCGAGTACCGTGTTGAGCATGAAGATATTGCAAAACg  <  1:917501/68‑1 (MQ=255)

                                                                  cTACTGGCTATTTACAACGCGAATTACCAGAGCGAGTACCGTGTTGAGCATGAAGATATTGCAAAACg  <  1:1269571/68‑1 (MQ=255)

                                                                  cTACTGGCTATTTACAACGCGAATTACCAGAGCGAGTACCGTGTTGAGCATGAAGATATTGCAAAACg  <  1:1878557/68‑1 (MQ=255)

                                                                  cTACTGGCTATTTACAACGCGAATTACCAGAGCGAGTACCGTGTTGAGCATGAAGATATTGCAAAACg  <  1:1902361/68‑1 (MQ=255)

                                                                  cTACTGGCTATTTACAACGCGAATTACCAGAGCGAGTACCGTGTTGAGCATGAAGATATTGCAAAACg  <  1:630019/68‑1 (MQ=255)

                                                                  cTACTGGCTATTTACAACGCGAATTACCAGAGCGAGTACCGTGTTGAGCATGAAGATATTGCAAAACg  <  1:2415569/68‑1 (MQ=255)

                                                                  cTACTGGCTATTTACAACGCGAATTACCAGAGCGAGTACCGTGTTGAGCATGAAGATATTGCAAAACg  <  1:520776/68‑1 (MQ=255)

                                                                  cTACTGGCTATTTACAACGCGAATTACCAGAGCGAGTACCGTGTTGAGCATGAAGATATTGCAAAACg  <  1:3290174/68‑1 (MQ=255)

                                                                  cTACTGGCTATTTACAACGCGAATTACCAGAGCGAGTACCGTGTTGAGCATGAAGATATTGCAAAACg  <  1:3219924/68‑1 (MQ=255)

CCGATTGCTATTGCCGAAGTACGCGAAGCACTCACTCGTACTCTGGCGACTGAACTGGCGGATGAGCTACTGGCTATTTACAACGCGAATTACCAGAGCGAGTACCGTGTTGAGCATGAAGATATTGCAAAACG  >  minE/649220‑649353

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