BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	648,659	A→G	37.7%	I447V (ATT→GTT)	pepN →	aminopeptidase N

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	648,659	0	A	G	37.7%	22.0 / 25.8	45	I447V (ATT→GTT)	pepN	aminopeptidase N
Reads supporting (aligned to +/- strand): ref base A (4/24); new base G (10/7); total (14/31)
Fisher's exact test for biased strand distribution p-value = 2.88e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCGGAAACCGAGCAGTACACCCTGACCA  >  minE/648596‑648713

gAAGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATtgtg                                                   <  1:2445973/69‑1 (MQ=255)

  aGATGCGTCGAATGTCGATCTCTCCCATTTCCGCTGTTGGTACAGCCAGTCCGGTACACCGATTGTGa                                                  <  1:343352/68‑1 (MQ=255)

  aGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGa                                                  <  1:2622113/68‑1 (MQ=255)

  aGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGa                                                  <  1:232628/68‑1 (MQ=255)

  aGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGa                                                  <  1:2317664/68‑1 (MQ=255)

  aGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGa                                                  <  1:3042759/68‑1 (MQ=255)

  aGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGa                                                  <  1:2196633/68‑1 (MQ=255)

  aGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGa                                                  <  1:3339753/68‑1 (MQ=255)

   gATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGa                                                  <  1:197479/67‑1 (MQ=255)

          cGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGa                                                  <  1:2850539/60‑1 (MQ=255)

                          ccATTTCCGCCGTTGGTACAGCCAGTCCGGCACGCCGGTTGTGACCGTCAAAGACGACTACAATCCgg                          >  1:3319778/1‑68 (MQ=255)

                          ccATTTCCGCCGTTGGTACAGCCAGTCCGGCACGCCGGTTGTGACCGTCAAAGACGACTACAATCCGGa                         >  1:775401/1‑69 (MQ=255)

                          ccATTTCCGCCGTTGGTACAGCCAGTCCGGCACGCCGGTTGTGACCGTCAAAGACGACTACAATCCGGa                         >  1:2515633/1‑69 (MQ=255)

                          ccATTTCCGCCGTTGGTACAGCCAGTCCGGCACGCCGGTTGTGACCGTCAAAGACGACTACAATCCGGa                         >  1:2894755/1‑69 (MQ=255)

                          ccATTTCCGCCGTTGGTACAGCCAGTCCGGCACGCCGGTTGTGACCGTCAAAGACGACTACAATCCGGa                         >  1:2102778/1‑69 (MQ=255)

                          ccATTTCCGCCGTTGGTACAGCCAGTCCGGCACGCCGGTTGTGACCGTCAAAGACGACTACAATCCGGa                         >  1:396484/1‑69 (MQ=255)

                          ccATTTCCGCCGTTGGTACAGCCAGTCCGGCACGCCGGTTGTGACCGTCAAAGACGACTACAATCCGGa                         >  1:729724/1‑69 (MQ=255)

                          ccATTTCCGCCGTTGGTACAGCCAGTCCGGCACGCCGGTTGTGACCGTCAAAGACGACTACAATCCGGa                         >  1:1805311/1‑69 (MQ=255)

                          ccATTTCCGCCGTTGGTACAGCCAGTCCGGCACGCCGGTTGTGACCGTCAAAGACGACTACAATCCGGa                         >  1:1720884/1‑69 (MQ=255)

                          ccATTTCCGCCGTTGGTACAGCCAGTCCGGCACGCCGGTTGTGACCGTCAAAGACGACTACAATCCGGa                         >  1:1695358/1‑69 (MQ=255)

                           cATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCGGaa                        >  1:2400515/1‑69 (MQ=255)

                           cATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCGGaa                        >  1:1662426/1‑69 (MQ=255)

                           cATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCGGaa                        >  1:2692537/1‑69 (MQ=255)

                               tCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCGGAAACCg                    >  1:3227722/1‑69 (MQ=255)

                                  gccgTTGGTACAGCCAGTCCGGCACGCCGGTTGTGACCGTCAAAGACGACTACAATCCGGAAACTGAg                  <  1:552993/68‑1 (MQ=255)

                                  gccgTTGGTACAGCCAGTCCGGCACGCCGGTTGTGACCGTCAAAGACGACTACAATCCGGAAACTGAg                  <  1:288400/68‑1 (MQ=255)

                                  gccgTTGGTACAGCCAGTCCGGCACGCCGGTTGTGACCGTCAAAGACGACTACAATCCGGAAACTGAg                  <  1:1545525/68‑1 (MQ=255)

                                  gccgTTGGTACAGCCAGTCCGGCACGCCGGTTGTGACCGTCAAAGACGACTACAATCCGGAAACTGAg                  <  1:142647/68‑1 (MQ=255)

                                  gccgTTGGTACAGCCAGTCCGGCACGCCGGTTGTGACCGTCAAAGACGACTACAATCCGGAAACTGAg                  <  1:2251609/68‑1 (MQ=255)

                                  gccgTTGGTACAGCCAGTCCGGCACGCCGGTTGTGACCGTCAAAGACGACTACAATCCGGAAACTGAg                  <  1:1822937/68‑1 (MQ=255)

                                    cgTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCGGAAACCGAGCAg               <  1:3139501/69‑1 (MQ=255)

                                    cgTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCGGAAACCGAGCAg               <  1:2001885/69‑1 (MQ=255)

                                    cgTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCGGAAACCGAGCAg               <  1:454686/69‑1 (MQ=255)

                                    cgTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCGGAAACCGAGCAg               <  1:217851/69‑1 (MQ=255)

                                    cgTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCGGAAACCGAGCAg               <  1:2008610/69‑1 (MQ=255)

                                    cgTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCGGAAACCGAGCAg               <  1:2106231/69‑1 (MQ=255)

                                     gTTGGTACAGCCAGTCCGGCACGCCGGTTGTGACCGTCAAAGACGACTACAATCCGGAAACTGAg                  <  1:522255/65‑1 (MQ=255)

                                         gTACAGCCAGTCCGGTACAGCGATTGTGACCGTCAAAGACGACTACAATCCGGAAACCGAGCAg               <  1:2475079/64‑1 (MQ=255)

                                                cAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCGGAAACCGAGCAGTACACCCTGAcc   <  1:2223693/69‑1 (MQ=255)

                                                cAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCGGAAACCGAGCAGTACACCCTGAcc   <  1:2496548/69‑1 (MQ=255)

                                                 aGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCGGAAACCGAGCAGTACACCCTGACCa  <  1:3172020/69‑1 (MQ=255)

                                                 aGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCGGAAACCGAGCAGTACACCCTGACCa  <  1:3152137/69‑1 (MQ=255)

                                                 aGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCGGAAACCGAGCAGTACACCCTGACCa  <  1:2233516/69‑1 (MQ=255)

                                                 aGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCAGGAAACCGAGCAGTACACCCTGACCa  <  1:2151700/69‑1 (MQ=255)

                                                  gTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCGGAAACCGAGCAGTACACCCTGACCa  <  1:547169/68‑1 (MQ=255)

GAAGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCGGAAACCGAGCAGTACACCCTGACCA  >  minE/648596‑648713

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