BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,243,208	C→T	37.0%	V130V (GTG→GTA)	nfi ←	endonuclease V

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,243,208	0	C	T	37.0%	31.1 / 30.8	54	V130V (GTG→GTA)	nfi	endonuclease V
Reads supporting (aligned to +/- strand): ref base C (17/17); new base T (13/7); total (30/24)
Fisher's exact test for biased strand distribution p-value = 3.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

AGTGGGGCCAGCGCGCCCGGTTCGCTGGAGAGCGGTTCGAATTTACCGCAGAGCCGTTTTTTCGCCACGCCAATGGTCGGCACATCCACCAATAAGCCAAAATGGCTGGCGACGCCAAGACGGCGAG  >  minE/2243143‑2243269

aGTGGGGCCAGCGCGCCCGGTTCGCTGGAGAGCGGTTCGAATTTACCGCAGAGCCGTTTTTTcgccac                                                             >  1:969535/1‑68 (MQ=255)

aGTGGGGCCAGCGCGCCCGGTTCGCTGGAGAGCGGTTCGAATTTACCGCAGAGCCGTTTTTTcgccac                                                             >  1:432495/1‑68 (MQ=255)

  tGGGGCCAGCGCGCCCGGTTCGCTGGAGAGCGGTTCGAATTTACCGCAGAGCCGTTTTTTcgccacgcc                                                          <  1:2384593/69‑1 (MQ=255)

    gggCCAGCGCGCCCGGTTCGCTGGAGAGCGGTTCGAATTTACCGCAGAGCCGTTTTTTCGCCACGCCaa                                                        <  1:2835287/69‑1 (MQ=255)

         aGCGCGCCCGGTTCGCTGGAGAGCGGTTCGAATTTACCGCAGAGCCGTTTTTTCGCCACGCCAATGGTc                                                   <  1:2366183/69‑1 (MQ=255)

           cgcgcCCGGTTCGCTGGAGAGCGGTTCGAATTTACCGCAGAGCCGTTTTTTCGCCACGCCAATGGTCgg                                                 <  1:467024/69‑1 (MQ=255)

           cgcgcCCGGTTCGCTGGAGAGCGGTTCGAATTTACCGCAGAGCCGTTTTTTCGCCACGCCAATGGTCgg                                                 <  1:1790329/69‑1 (MQ=255)

           cgcgcCCGGTTCGCTGGAGAGCGGTTCGAATTTACCGCAGAGCCGTTTTTTCGCCACGCCAATGGGCgg                                                 <  1:2779122/69‑1 (MQ=255)

              gcCCGGTTCGCTGGAGAGCGGTTCGAATTTACCGCAGAGCCGTTTTTTCGCCACGCCAATGGTCGGcac                                              <  1:2635152/69‑1 (MQ=255)

              gcCCGGTTCGCTGGAGAGCGGTTCGAATTTACCGCAGAGCCGTTTTTTCGCCACGCCAATGGTCGGcac                                              <  1:2340450/69‑1 (MQ=255)

              gcCCGGTTCCCTGGAGAGCGGTTCGAATTTACCGCAGAGCCGTTTTTTCGCCACGCCAATGGTCGGcac                                              <  1:3178680/69‑1 (MQ=255)

               cccGGTTCGCTGGAGAGCGGTTCGAATTTACCGCAGAGCCGTTTTTTCGCCACGCCAATGGTCGGcac                                              <  1:503125/68‑1 (MQ=255)

               cccGGTTCGCTGGAGAGCGGTTCGAATTTACCGCAGAGCCGTTTTTTCGCCACGCCAATGGTCGGcac                                              <  1:2015350/68‑1 (MQ=255)

               cccGGTTCGCTGGAGAGCGGTTCGAATTTACCGCAGAGCCGTTTTTTCGCCACGCCAATGGTCGGcac                                              <  1:2859660/68‑1 (MQ=255)

                      cGCTGGAGAGCGGTTCGAATTTACCGCAGAGCCGTTTTTTCGCCACGCCAATGGTCGGCACATccacc                                       >  1:2124378/1‑68 (MQ=255)

                      cGCTGGAGAGCGGTTCGAATTTACCGCAGAGCCGTTTTTTCGCCACGCCAATGGTCGGCACATccacc                                       >  1:2125619/1‑68 (MQ=255)

                      cGCTGGAGAGCGGTTCGAATTTACCGCAGAGCCGTTTTTTCGCCACGCCAATGGTCGGCACATccacc                                       >  1:3069927/1‑68 (MQ=255)

                         tGGATAGCGGTTCGAATTTACCGCAGAGCCGTTTTTTCGCTACGCCAATGGTCGGCACATCCACCAATa                                   <  1:1771896/69‑1 (MQ=255)

                         tGGAGAGCGGTTCGAATTTACCGCAGAGCCGTTTTTTCGCTACGCCAATGGTCGGCACATCCACCAATa                                   <  1:2496383/69‑1 (MQ=255)

                         tGGAGAGCGGTTCGAATTTACCGCAGAGCCGTTTTTTCGCTACGCCAATGGTCGGCACATCCACCAATa                                   <  1:141944/69‑1 (MQ=255)

                         tGGAGAGCGGTTCGAATTTACCGCAGAGCCGTTTTTTCGCTACGCCAATGGTCGGCACATCCACCAATa                                   <  1:3286401/69‑1 (MQ=255)

                         tGGAGAGCGGTTCGAATTTACCGCAGAGCCGTTTTTTCGCTACGCCAATGGTCGGCACATCCACCAATa                                   <  1:631803/69‑1 (MQ=255)

                              agCGGTTCGAATTTACCGCAGAGCCGTTTTTTcgcc                                                               >  1:777575/1‑36 (MQ=255)

                              agCGGTTCGAATTTACCGCAGAGCCGTTTTTTCGCCACGCCAATGGTCgg                                                 >  1:121536/1‑50 (MQ=255)

                              agCGGTTCGAATTTACCGCAGAGCCGTTTTTTCGCCACGCCAATGGTCgg                                                 >  1:432184/1‑50 (MQ=255)

                              agCGGTTCGAATTTACCGCAGAGCCGTTTTTTCGCCACGCCAATGGTCGGCACATCCACCAATAAGcc                               >  1:86983/1‑68 (MQ=255)

                              agCGGTTCGAATTTACCGCAGAGCCGTTTTTTCGCCACGCCAATGGTCGGCACATCCACCAATAAGcc                               >  1:736623/1‑68 (MQ=255)

                              agCGGTTCGAATTTACCGCAGAGCCGTTTTTTCGCCACGCCAATGGTCGGCACATCCACCAATAAGcc                               >  1:1179804/1‑68 (MQ=255)

                              agCGGTTCGAATTTACCGCAGAGCCGTTTTTTCGCCACGCCAATGGTCGGCACATCCACCAATAAGcc                               >  1:1069291/1‑68 (MQ=255)

                              agCGGTTCGAATTTACCGCAGAGCCGTTTTTTCGCCACGCCAATGGTCGGCACATCCACCAATAAGcc                               >  1:330693/1‑68 (MQ=255)

                              agCGGTTCGAATTTACCGCAGAGCCGTTTTTTCGCCACGCCAATGGTCGGCACATCCACCAATAAGcc                               >  1:1620291/1‑68 (MQ=255)

                              agCGGTTCGAATTTACCGCAGAGCCGTTTTTTCGCCACGCCAATGGTCGGCACATCCACCAATAAGcc                               >  1:3178358/1‑68 (MQ=255)

                               gCGGTTCGAATTTACCGCAGAGCCGTTTTTTCGCCACGCCAATGGTCGGCACATCCACCAATAAGCCaa                             <  1:1110886/69‑1 (MQ=255)

                               gCGGTTCGAATTTACCGCAGAGCCGTTTTTTCGCCACGCCAATGGTCGGCACATCCACCAATAAGCCaa                             <  1:359471/69‑1 (MQ=255)

                                cGGTTCGAATTTACCGCAGAGCCGTTTTTTCGCCACGCCAATGGt                                                    >  1:790268/1‑45 (MQ=255)

                                cGGTTCGAATTTACCGCAGAGCCGTTTTTTCGCCACGCCAATGGTCGGCACCTCCACCAATAAGCCaaa                            >  1:1496367/1‑69 (MQ=255)

                                cGGTTCGAATTTACCGCAGAGCCGTTTTTTCGCCACGCCAATGGTCGGCACATCCACCAATAAGCCaa                             >  1:965546/1‑68 (MQ=255)

                                         tttACCGCAGAGCCGTTTTTTCGCCACGCCAATGGTCGGCACATCCACCAATAAGCCAAAAtggctggc                   <  1:2398984/69‑1 (MQ=255)

                                         tttACCGCAGAGCCGTTTTTTCGCCACGCCAATGGTCGGCACATCCACCAATAAGCCAAAAtggctggc                   <  1:3006019/69‑1 (MQ=255)

                                          ttACCGCAGAGCCGTTTTTTCGCCACGCCAATGGTCGGCACATCCACCAATAAGCCAAAAtggctggc                   <  1:2722536/68‑1 (MQ=255)

                                                cAGAGCCGTTTTTTCGCTACGCCAATGGTCGGCACATCCACCAATa                                   <  1:387095/46‑1 (MQ=255)

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

                                                         tttttCGCTACGCCAATGGTCGGCACATCCACCAATa                                   <  1:2571964/37‑1 (MQ=255)

                                                          ttttCGCTACGCCAATGGTCGGCACATCCACCAATa                                   >  1:705485/1‑36 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

                                                          ttttCGCCACGCCAATGGTCGGCACATCCACCAATAAGCCAAAAtggctg                     >  1:634270/1‑50 (MQ=255)

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

AGTGGGGCCAGCGCGCCCGGTTCGCTGGAGAGCGGTTCGAATTTACCGCAGAGCCGTTTTTTCGCCACGCCAATGGTCGGCACATCCACCAATAAGCCAAAATGGCTGGCGACGCCAAGACGGCGAG  >  minE/2243143‑2243269

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