BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,344,055	G→A	67.8%	S68S (AGC→AGT)	wzc ←	protein‑tyrosine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,344,055	0	G	A	67.8%	99.9 / 64.4	121	S68S (AGC→AGT)	wzc	protein‑tyrosine kinase
Reads supporting (aligned to +/- strand): ref base G (33/6); new base A (59/23); total (92/29)
Fisher's exact test for biased strand distribution p-value = 1.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.75e-01

Only 100 of 121 total aligned reads displayed.

TCCGATGCAGGCGGTTTGTTGGCTAATGCCGAACCAATGTCCTGCACTAACGAATTGCCGCTGTTTTGCTCGATTTGTACCAGTGCGTCGGCACTATAAATCGGCGTGGCGAAGAAGG  >  minE/1343996‑1344113

tCCGATGCAGGCGGTTTGTTGGCTAATGCCGAACCAATGTCCTGCACTAACGAATTGCCGCTGTTTTGc                                                   >  1:2291846/1‑69 (MQ=255)

  cGATGCTGGCGGTTTGTTGGCTAATGCCGAACCAATGTCCTGCACTAACGAATTGCCACTGtttt                                                     <  1:1562057/65‑1 (MQ=255)

  cGATGCTGGCGGTTTGTTGGCTAATGCCGAACCAATGTCCTGCACTAACGAATTGCCACTGtttt                                                     <  1:2369729/65‑1 (MQ=255)

  cGATGCTGGCGGTTTGTTGGCTAATGCCGAACCAATGTCCTGCACTAACGAATTGCCACTGtttt                                                     <  1:2209087/65‑1 (MQ=255)

  cGATGCTGGCGGTTTGTTGGCTAATGCCGAACCAATGTCCTGCACTAACGAATTGCCACTGtttt                                                     <  1:2201438/65‑1 (MQ=255)

  cGATGCTGGCGGTTTGTTGGCTAATGCCGAACCAATGTCCTGCACTAACGAATTGCCACTGtttt                                                     <  1:1973031/65‑1 (MQ=255)

  cGATGCTGGCGGTTTGTTGGCTAATGCCGAACCAATGTCCTGCACTAACGAATTGCCACTGtttt                                                     <  1:1949419/65‑1 (MQ=255)

  cGATGCTGGCGGTTTGTTGGCTAATGCCGAACCAATGTCCTGCACTAACGAATTGCCACTGtttt                                                     <  1:1855238/65‑1 (MQ=255)

  cGATGCTGGCGGTTTGTTGGCTAATGCCGAACCAATGTCCTGCACTAACGAATTGCCACTGtttt                                                     <  1:1789741/65‑1 (MQ=255)

  cGATGCTGGCGGTTTGTTGGCTAATGCCGAACCAATGTCCTGCACTAACGAATTGCCACTGtttt                                                     <  1:3265552/65‑1 (MQ=255)

  cGATGCTGGCGGTTTGTTGGCTAATGCCGAACCAATGTCCTGCACTAACGAATTGCCACTGtttt                                                     <  1:3311854/65‑1 (MQ=255)

  cGATGCTGGCGGTTTGTTGGCTAATGCCGAACCAATGTCCTGCACTAACGAATTGCCACTGtttt                                                     <  1:591038/65‑1 (MQ=255)

  cGATGCTGGCGGTTTGTTGGCTAATGCCGAACCAATGTCCTGCACTAACGAATTGCCACTGtttt                                                     <  1:1276873/65‑1 (MQ=255)

  cGATGCTGGCGGTTTGTTGGCTAATGCCGAACCAATGTCCTGCACTAACGAATTGCCACTGtttt                                                     <  1:68219/65‑1 (MQ=255)

  cGATGCTGGCGGTTTGTTGGCTAATGCCGAACCAATGTCCTGCACTAACGAATTGCCACTGtttt                                                     <  1:1025886/65‑1 (MQ=255)

  cGATGCTGGCGGTTTGTTGGCTAATGCCGAACCAATGTCCTGCACTAACGAATTGCCACTGtttt                                                     <  1:95429/65‑1 (MQ=255)

  cGATGCTGGCGGTTTGTTGGCTAATGCCGAACCAATGTCCTGCACTAACGAATTGCCACTGtttt                                                     <  1:1155792/65‑1 (MQ=255)

  cGATGCAGGCGGTTTGTTGGCTAATGCCGAACCAATGTCCTGCACTAACGAATTGCCGCTGTTTTGCTc                                                 <  1:1490606/69‑1 (MQ=255)

   gATGCTGGCGGTTTGTTGGCTAATGCCGAACCAATGTCCTGCACTAACGAATTGCCACTGtttt                                                     <  1:651933/64‑1 (MQ=255)

              tttGTTGGCTAATGCCGAACCAATGTCCTGCACTAACGAATTGCCGCTGTTTTGCTCGATTTGTACCAg                                     <  1:1858279/69‑1 (MQ=255)

              tttGTTGGCTAATGCCGAACCAATGTCCTGCACTAACGAATTGCCGCTGTTTTGCTCGATTTGTACCAg                                     <  1:3155345/69‑1 (MQ=255)

               ttgttgGCTAATGCCGAACCAATGTCCTGCACTAACGAATTGCCGCTGTTTTGCTCGATTTGTACCAg                                     <  1:1133080/68‑1 (MQ=255)

                 gttgGCTAATGCCGAACCAATGTCCTGCACTAACGAATTGCCACTGtttt                                                     <  1:3191437/50‑1 (MQ=255)

                      cTAATGGCGAACCAATGTCCTGCACTAACGAATTGCCACTGTTTTGCTCGATTTGTACCa                                      >  1:1412061/1‑60 (MQ=255)

                      cTAATGCCGAACCAATGTCCTGTACTAACGAATTGCCACTGTTTTGCTCGATTTGTACCa                                      >  1:712410/1‑60 (MQ=255)

                      cTAATGCCGAACCAATGTCCTGCACTAACGAATTGCCACTGTTTTGCTCGATTTGTACCa                                      >  1:2475711/1‑60 (MQ=255)

                      cTAATGCCGAACCAATGTCCTGCACTAACGAATTGCCACTGTTTTGCTCGATTTGTACCa                                      >  1:2507953/1‑60 (MQ=255)

                      cTAATGCCGAACCAATGTCCTGCACTAACGAATTGCCACTGTTTTGCTCGATTTGTACCa                                      >  1:845302/1‑60 (MQ=255)

                      cTAATGCCGAACCAATGTCCTGCACTAACGAATTGCCACTGTTTTGCTCGATTTGTACCa                                      >  1:960003/1‑60 (MQ=255)

                      cTAATGCCGAACCAATGTCCTGCACTAACGAATTGCCACTGTTTTGCTCGATTTGTACCa                                      >  1:2396726/1‑60 (MQ=255)

                      cTAATGCCGAACCAATGTCCTGCACTAACGAATTGCCACTGTTTTGCTCGATTTGTACCa                                      >  1:2831106/1‑60 (MQ=255)

                      cTAATGCCGAACCAATGTCCTGCACTAACGAATTGCCACTGTTTTGCTCGATTTGTACCa                                      >  1:284706/1‑60 (MQ=255)

                      cTAATGCCGAACCAATGTCCTGCACTAACGAATTGCCACTGTTTTGCTCGATTTGTACCa                                      >  1:835664/1‑60 (MQ=255)

                      cTAATGCCGAACCAATGTCCTGCACTAACGAATTGCCACTGTTTTGCTCGATTTGTACCa                                      >  1:300239/1‑60 (MQ=255)

                      cTAATGCCGAACCAATGTCCTGCACTAACGAATTGCCACTGTTTTGCTCGATTTGTACCa                                      >  1:801094/1‑60 (MQ=255)

                      cTAATGCCGAACCAATGTCCTGCACTAACGAATTGCCACTGTTTTGCTCGATTTGTACCa                                      >  1:784797/1‑60 (MQ=255)

                      cTAATGCCGAACCAATGTCCTGCACTAACGAATTGCCACTGTTTTGCTCGATTTGTACCa                                      >  1:1021552/1‑60 (MQ=255)

                      cTAATGCCGAACCAATGTCCTGCACTAACGAATTGCCACTGTTTTGCTCGATTTGTACCa                                      >  1:354648/1‑60 (MQ=255)

                      cTAATGCCGAACCAATGTCCTGCACTAACGAATTGCCACTGTTTTGCTCGATTTGTACCa                                      >  1:55028/1‑60 (MQ=255)

                      cTAATGCCGAACCAATGTCCTGCACTAACGAATTGCCACTGTTTTGCTCGATTTGTACCa                                      >  1:582758/1‑60 (MQ=255)

                      cTAATGCCGAACCAATGTCCTGCACTAACGAATTGCCACTGTTTTGCTCGATTTGTACCa                                      >  1:2373639/1‑60 (MQ=255)

                      cTAATGCCGAACCAATGTCCTGCACTAACGAATTGCCACTGTTTTGCTCGATTTGTACCa                                      >  1:227882/1‑60 (MQ=255)

                      cTAATGCCGAACCAATGTCCTGCACTAACGAATTGCCACTGTTTTGCTCGATTTGTACCa                                      >  1:1625876/1‑60 (MQ=255)

                      cTAATGCCGAACCAATGTCCTGCACTAACGAATTGCCACTGTTTTGCTCGATTTGTACCa                                      >  1:1062308/1‑60 (MQ=255)

                      cTAATGCCGAACCAATGTCCTGCACTAACGAATTGCCACTGTTTTGCTCGATTTGTACCa                                      >  1:1342392/1‑60 (MQ=255)

                      cTAATGCCGAACCAATGTCCTGCACTAACGAATTGCCACTGTTTTGCTCGATTTGTACCa                                      >  1:1318295/1‑60 (MQ=255)

                      cTAATGCCGAACCAATGTCCTGCACTAACGAATTGCCACTGTTTTGCTAGATTTGTACCa                                      >  1:2565892/1‑60 (MQ=255)

                         aTGCCGAACCAATGTCCTGCTCTAACGAATTGCCGCTGTTTTGCTCGATTTg                                           >  1:488411/1‑52 (MQ=255)

                         aTGCCGAACCAATGTCCTGCACTAACGAATTGCCGCTGTTTTGCTCGAttt                                            >  1:1358738/1‑51 (MQ=255)

                         aTGCCGAACCAATGTCCTGCACTAACGAATTGCCGCTGTTTTGCTCGATTTg                                           >  1:1186546/1‑52 (MQ=255)

                         aTGCCGAACCAATGTCCTGCACTAACGAATTGCCGCTGTTTTGCTCGATTTg                                           >  1:1034029/1‑52 (MQ=255)

                         aTGCCGAACCAATGTCCTGCACTAACGAATTGCCGCTGTTTTGCTCGATTTg                                           >  1:1827074/1‑52 (MQ=255)

                         aTGCCGAACCAATGTCCTGCACTAACGAATTGCCGCTGTTTTGCTCGATTTg                                           >  1:477372/1‑52 (MQ=255)

                         aTGCCGAACCAATGTCCTGCACTAACGAATTGCCGCTGTTTTGCTCGATTTg                                           >  1:2942402/1‑52 (MQ=255)

                         aTGCCGAACCAATGTCCTGCACTAACGAATTGCCGCTGTTTTGCTCGATTTg                                           >  1:924094/1‑52 (MQ=255)

                         aTGCCGAACCAATGTCCTGCACTAACGAATTGCCGCTGTTTTGCTCGATTTg                                           >  1:504777/1‑52 (MQ=255)

                         aTGCCGAACCAATGTCCTGCACTAACGAATTGCCGCTGTTTTGCTCGATTTg                                           >  1:2808166/1‑52 (MQ=255)

                         aTGCCGAACCAATGTCCTGCACTAACGAATTGCCGCTGTTTTGCTCGATTTg                                           >  1:2788843/1‑52 (MQ=255)

                             cGAACCAATGTCCTGCACTAACGAATTGCCACTGtttt                                                     >  1:264506/1‑38 (MQ=255)

                             cGAACCAATGTCCTGCACTAACGAATTGCCACTGTTTTGCTCGATTTGTACCAGTGCGTCGGCa                           >  1:557208/1‑64 (MQ=255)

                             cGAACCAATGTCCTGCACTAACGAATTGCCACTGTTTTGCTCGATTTGTACCAGTGCGTCGGCa                           >  1:101888/1‑64 (MQ=255)

                             cGAACCAATGTCCTGCACTAACGAATTGCCACTGTTTTGCTCGATTTGTACCAGTGCGTCGGCa                           >  1:1124162/1‑64 (MQ=255)

                             cGAACCAATGTCCTGCACTAACGAATTGCCACTGTTTTGCTCGATTTGTACCAGTGCGTCGGCa                           >  1:3375414/1‑64 (MQ=255)

                             cGAACCAATGTCCTGCACTAACGAATTGCCACTGTTTTGCTCGATTTGTACCAGTGCGTCGGCa                           >  1:3296236/1‑64 (MQ=255)

                             cGAACCAATGTCCTGCACTAACGAATTGCCACTGTTTTGCTCGATTTGTACCAGTGCGTCGGCa                           >  1:689189/1‑64 (MQ=255)

                             cGAACCAATGTCCTGCACTAACGAATTGCCACTGTTTTGCTCGATTTGTACCAGTGCGTCGGCa                           >  1:1627602/1‑64 (MQ=255)

                             cGAACCAATGTCCTGCACTAACGAATTGCCACTGTTTTGCTCGATTTGTACCAGTGCGTCGGCa                           >  1:1771724/1‑64 (MQ=255)

                             cGAACCAATGTCCTGCACTAACGAATTGCCACTGTTTTGCTCGATTTGTACCAGTGCGTCGGCa                           >  1:2626544/1‑64 (MQ=255)

                             cGAACCAATGTCCTGCACTAACGAATTGCCACTGTTTTGCTCGATTTGTACCAGTGCGTCGGCa                           >  1:2528704/1‑64 (MQ=255)

                             cGAACCAATGTCCTGCACTAACGAATTGCCACTGTTTTGCTCGATTTGTACCAGTGCGTCGGCa                           >  1:314782/1‑64 (MQ=255)

                             cGAACCAATGTCCTGCACTAACGAATTGCCACTGTTTTGCTCGATTTGTACCAGTGCGTCGGCa                           >  1:2675702/1‑64 (MQ=255)

                             cGAACCAATGTCCTGCACTAACGAATTGCCACTGTTTTGCTCGATTTGTACCAGTGCGTCGGCa                           >  1:1989288/1‑64 (MQ=255)

                             cGAACCAATGTCCTGCACTAACGAATTGCCACTGTTTTGCTCGATTTGTACCAGTGCGTCGGCa                           >  1:2477004/1‑64 (MQ=255)

                             cGAACCAATGTCCTGCACTAACGAATTGCCACTGTTTTGCTCGATTTGTACCAGTGCGTCGGCa                           >  1:1985086/1‑64 (MQ=255)

                             cGAACCAATGTCCTGCACTAACGAATTGCCACTGTTTTGCTCGATTTGTACCAGTGCGTCGGCa                           >  1:2981184/1‑64 (MQ=255)

                              gAACCAATGTCCTGCACTAACGAATTGCCGCTGTTTTGCTCGATTTGTACCAGTGCGTCGGCACtata                      >  1:1476125/1‑68 (MQ=255)

                              gAACCAATGTCCTGCACTAACGAATTGCCGCTGTTTTGCTCGATTTGTACCAGTGCGTCGGCACTATaa                     >  1:755548/1‑69 (MQ=255)

                              gAACCAATGTCCTGCACTAACGAATTGCCGCTGTTTTGCTCGATTTGTACCAGTGCGTCGGCACTATaa                     >  1:3245073/1‑69 (MQ=255)

                               aaCCAATGTCCTGCACTAACGAATTGCCACTGttttt                                                    >  1:2869879/1‑36 (MQ=38)

                               aaCCAATGTCCTGCACTAACGAATTGCCACTGTTTTTCTCGATTTGTACCAGTGc                                  >  1:479397/1‑55 (MQ=255)

                               aaCCAATGTCCTGCACTAACGAATTGCCACTGTTTTTCTCGATTTGTACCAGTGc                                  >  1:1601160/1‑55 (MQ=255)

                               aaCCAATGTCCTGCACTAACGAATTGCCACTGTTTTTCTCGATTTGTACCAGTGc                                  >  1:3108824/1‑55 (MQ=255)

                               aaCCAATGTCCTGCACTAACGAATTGCCACTGTTTTTCTCGATTTGTACCAGTGc                                  >  1:1339906/1‑55 (MQ=255)

                               aaCCAATGTCCTGCACTAACGAATTGCCACTGTTTTTCTCGATTTGTACCAGTGc                                  >  1:1229592/1‑55 (MQ=255)

                               aaCCAATGTCCTGCACTAACGAATTGCCACTGTTTTTCTCGATTTGTACCAGTGc                                  >  1:2174532/1‑55 (MQ=255)

                               aaCCAATGTCCTGCACTAACGAATTGCCACTGTTTTTCTCGATTTGTACCAGTGc                                  >  1:1264665/1‑55 (MQ=255)

                                       tgcTGCACTAACGAATTGCCGCTGTTTTGCTCGATTTGTACCAGTGCGTCGGCAc                          >  1:371472/3‑55 (MQ=255)

                                       tCCTGCACTAACGAATTGCCGCTGTTTTGCTCGATTTGTACCAGTGCGTCGGCACTATAAATCGGCGTg            >  1:864465/1‑69 (MQ=255)

                                       tCCTGCACTAACGAATTGCCGCTGTTTTGCTCGATTTGTACCAGTGCGTCGGCACTATAAATCGGCGTg            >  1:317722/1‑69 (MQ=255)

                                                aaCGAATTGCCGCTGTTTTGCTCGATTTGTACCAGTGCGTCGGCACTATAAATCGGCGTGGCGaagaag   >  1:21309/1‑69 (MQ=255)

                                                aaCGAATTGCCGCTGTTTTGCTCGATTTGTACCAGTGCGTCGGCACTATAAATCGGCGTGGCGaagaag   >  1:2763479/1‑69 (MQ=255)

                                                aaCGAATTGCCGCTGTTTTGCTCGATTTGTACCAGTGCGTCGGCACTATAAATCGGCGTGGCGaagaag   >  1:1363959/1‑69 (MQ=255)

                                                aaCGAATTGCCGCTGTTTTGCTCGATTTGTACCAGTGCGTCGGCACTATAAATCGGCGTGGCGaagaag   >  1:3086095/1‑69 (MQ=255)

                                                aaCGAATTGCCGCTGTTTTGCTCGATTTGTACCAGTGCGTCGGCACTATAAAACGGCGTGGCGaagaag   >  1:3227587/1‑69 (MQ=255)

                                                 aCGAATTGCCGCTGTTTTGCTCGATTTGTACCAGTGCGTCGGCACTATAAATCGGCGTGGCGAAGAAgg  >  1:2641571/1‑69 (MQ=255)

                                                 aCGAATTGCCGCTGTTTTGCTCGATTTGTACCAGTGCGTCGGCACTATAAATCGGCGTGGCGAAGAAgg  >  1:2618867/1‑69 (MQ=255)

                                                 aCGAATTGCCGCTGTTTTGCTCGATTTGTACCAGTGCGTCGGCACTATAAATCGGCGTGGCGAAGAAgg  >  1:622616/1‑69 (MQ=255)

                                                 aCGAATTGCCGCTGTTTTGCTCGATTTGTACCAGTGCGTCGGCACTATAAATCGGCGTGGCGAAGAAgg  >  1:383765/1‑69 (MQ=255)

                                                 aCGAATTGCCGCTGTTTTGCTCGATTTGTACCAGTGCGTCGGCACTATAAATCGGCGTGGCGAAGAAgg  >  1:930631/1‑69 (MQ=255)

                                                       tGCCGCTGTTTTGCTCGATTTGTACCAGTGCGTCGGCACTATAAATCGGCGTGGCgaaga     <  1:2386929/60‑1 (MQ=255)

TCCGATGCAGGCGGTTTGTTGGCTAATGCCGAACCAATGTCCTGCACTAACGAATTGCCGCTGTTTTGCTCGATTTGTACCAGTGCGTCGGCACTATAAATCGGCGTGGCGAAGAAGG  >  minE/1343996‑1344113

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