BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	421,648	A→G	14.4%	G277G (GGT→GGC)	asnB ←	asparagine synthetase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	421,648	0	A	G	14.4%	93.1 / 5.7	55	G277G (GGT→GGC)	asnB	asparagine synthetase B
Reads supporting (aligned to +/- strand): ref base A (4/43); new base G (4/4); total (8/47)
Fisher's exact test for biased strand distribution p-value = 1.10e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAATTTCGTGATGCACCGTGCCCAGATGGTTTGCCACTTCCTGGGCTGCTTTCAGATCCGGTGAACCCGGCAGACCTACAGCAAAGGAGTGTAACTGCGGCCACCAGGCTTCAGAGCGTTCCTGATCTTC  >  minE/421584‑421713

gaaTTTCGTGATGCACCGTGCCCAGATTGTTTGCCACTTCCTGGGCTGCTTTCAGATCCGGTGAGccc                                                                >  1:1422488/1‑68 (MQ=255)

gaaTTTCGTGATGCACCGTGCCCAGATGGTTTGCCACTTCCTGGGCTGCTTTCAGATCCGGTGAGccc                                                                >  1:905888/1‑68 (MQ=255)

gaaTTTCGTGATGCACCGTGCCCAGATGGTTTGCCACTTCCTGGGCTGCTTTCAGATCCGGTGAGccc                                                                >  1:508627/1‑68 (MQ=255)

gaaTTTCGTGATGCACCGTGCCCAGATGGTTTGCCACTTCCTGGGCTGCTTTCAGATCCGGTGAGccc                                                                >  1:140142/1‑68 (MQ=255)

 aaTTTCGTGATGCACCGTGCCCAGATGGTTTGCCACTTCCTGGGCTGCTTTCAGATCCGGTGAAccc                                                                >  1:1005685/1‑67 (MQ=255)

 aaTTTCGTGATGCACCGTGCCCAGATGGTTTGCCACTTCCTGGGCTGCTTTCAGATCCGGTGAAccc                                                                >  1:2336217/1‑67 (MQ=255)

 aaTTTCGTGATGCACCGTGCCCAGATGGTTTGCCACTTCCTGGGCTGCTTTCAGATCCGGTGAAccc                                                                >  1:1746395/1‑67 (MQ=255)

    ttCGTGATGCACCGTGCCCAGATGGTTTGCCACTTCCTGGGCTGCTTTCAGATCCGGTGAACCCGGCAg                                                           <  1:2396803/69‑1 (MQ=255)

    ttCGTGATGCACCGTGCCCAGATGGTTTGCCACTTCCTGGGCTGCTTTCAGATCCGGTGAACCCGGCAg                                                           <  1:402607/69‑1 (MQ=255)

    ttCGTGATGCACCGTGCCCAGATGGTTTGCCACTTCCTGGGCTGCTTTCAGATCCGGTGAACCCGGCAg                                                           <  1:385906/69‑1 (MQ=255)

    ttCGTGATGCACCGTGCCCAGATGGTTTGCCACTTCCTGGGCTGCTTTCAGATCCGGTGAACCCGGCAg                                                           <  1:332762/69‑1 (MQ=255)

    ttCGTGATGCACCGTGCCCAGATGGTTTGCCACTTCCTGGGCTGCTTTCAGATCCGGTGAACCCGGCAg                                                           <  1:3111461/69‑1 (MQ=255)

    ttCGTGATGCACCGTGCCCAGATGGTTTGCCACTTCCTGGGCTGCTTTCAGATCCGGTGAACCCGGCAg                                                           <  1:2376660/69‑1 (MQ=255)

    ttCGTGATGCACCGTGCCCAGATGGTTTGCCACTTCCTGGGCTGCTTTCAGATCCGGTGAACCCGGCAg                                                           <  1:2303659/69‑1 (MQ=255)

    ttCGTGATGCACCGTGCCCAGATGGTTTGCCACTTCCTGGGCTGCTTTCAGATCCGGTGAACCCGGCAg                                                           <  1:188147/69‑1 (MQ=255)

    ttCGTGATGCACCGTGCCCAGATGGTTTGCCACTTCCTGGGCTGCTTTCAGATCCGGTGAACCCGGCAg                                                           <  1:1823585/69‑1 (MQ=255)

    ttCGTGATGCACCGTGCCCAGATGGTTTGCCACTTCCTGGGCTGCTTTCAGATCCGGTGAACCCGGCAg                                                           <  1:1672842/69‑1 (MQ=255)

    ttCGTGATGCACCGTGCCCAGATGGTTTGCCACTTCCTGGGCTGCTTTCAGATCCGGTGAACCCGGCAg                                                           <  1:1753283/69‑1 (MQ=255)

    ttCGTGATGAACCGTGCCCAGATGGTTTGCCACTTCCTGGGCTGCTTTCAGATCCGGTGAACCCGGCAg                                                           <  1:2819645/69‑1 (MQ=255)

     tCGTGATGCACCGTGCCCAGATGGTTTGCCACTTCCTGGGCTGCTTTCAGATCCGGTGAACCCGGCAg                                                           <  1:3356465/68‑1 (MQ=255)

       gTGATGCACCGTGCCCAGATGGTTTGCCACTTCCTGGGCTGCTTTCAGATCCGGTGAACCCGGCAGAc                                                         <  1:1817024/68‑1 (MQ=255)

       gTGATGCACCGTGCCCAGATGGTTTGCCACTTCCTGGGCTGCTTTCAGATCCGGTGAACCCGGCAGAc                                                         <  1:159001/68‑1 (MQ=255)

          aTGCACCGTGCCCAGATGGTTTGCCACTTCCTGGGCTGCTTTCAGATCCGGTGAACCCGGCAGACCTAc                                                     <  1:847729/69‑1 (MQ=255)

          aTGCACCGTGCCCAGATGGTTTGCCACTTCCTGGGCTGCTTTCAGATCCGGTGAACCCGGCAGACCTAc                                                     <  1:2460627/69‑1 (MQ=255)

          aTGCACCGTGCCCAGATGGTTTGCCACTTCCTGGGCTGCTTTCAGATCCGGTGAACCCGGCAGACCTAc                                                     <  1:360532/69‑1 (MQ=255)

          aTGCACCGTGCCCAGATGGTTTGCCACTTCCTGGGCTGCTTTCAGATCCGGTGAACCCGGCAGACCTAc                                                     <  1:2549448/69‑1 (MQ=255)

                cGTGCCCAGATGGTTTGCCACTTCCTGGGCTGCTTTCAGATCCGGTGAACCCGGCAGACCTAc                                                     <  1:797981/63‑1 (MQ=255)

                  tGCCCAGATGGTTTGCCACTTCCTGGGCTGCTTTCAGATCCGGTGAACCCGGCAGACCTACAGCAAAgg                                             <  1:2412334/69‑1 (MQ=255)

                  tGCCCAGATGGTTTGCCACTTCCTGGGCTGCTTTCAGATCCGGTGAACCCGGCAGACCTACAGCAAAgg                                             <  1:2501868/69‑1 (MQ=255)

                       aGATGGTTTGCCACTTCCTGGGCTGCTTTCAGATCCGGTGAACCCGGCAGACCTACAGCAAAGGAgtgt                                        <  1:2312765/69‑1 (MQ=255)

                       aGATGGTTTGCCACTTCCTGGGCTGCTTTCAGATCCGGTGAACCCGGCAGACCTACAGCAAAGGAgtgt                                        <  1:73853/69‑1 (MQ=255)

                       aGATGGTTTGCCACTTCCTGGGCTGCTTTCAGATCCGGTGAACCCGGCAGACCTACAGCAAAGGAgtgt                                        <  1:1679359/69‑1 (MQ=255)

                       aGATGGTTTGCCACTTCCTGGGCTGCTTTCAGATCCGGTGAACCCGGCAGACCTACAGCAAAGGAgtgt                                        <  1:2091092/69‑1 (MQ=255)

                       aGATGGTTTGCCACTTCCTGGGCTGCTTTCAGAACCGGTGAACCCGGCAGACCTACAGCAAAGGAgtgt                                        <  1:2445699/69‑1 (MQ=255)

                       aGATGGTTTACCACTTCCTGGGCTGCTTTCAGATCCGGTGAACCCGGCAGACCTACAGCAAAGGAgtgt                                        <  1:2498418/69‑1 (MQ=255)

                       aGATGGTTTACCACTTCCTGGGCTGCTTTCAGATCCGGTGAACCCGGCAGACCTACAGCAAAGGAgtgt                                        <  1:2674576/69‑1 (MQ=255)

                         aTGGTTTGCCACTTCCTGGGCTGCTTTCAGATCCGGTGAACCCGGCAGACCTACAGCAAAGGAgtgt                                        <  1:63631/67‑1 (MQ=255)

                                 ccACTTCCTGGGCTGCTTTCAGATCCGGTGAACCCGGCAGACCTACAGCAAAGGAGTGTAACTGCgga                               <  1:24914/68‑2 (MQ=255)

                                         tGGGCTGCTTTCAGATCCGGTGAACCCGGCAGACCTAc                                                     <  1:1970263/38‑1 (MQ=255)

                                                     aGATCCGGTGAGCCCGGCAGACCTACAGCAAAGGAGTGCAACTGCGGCCACCAGGCTTCAGAGCGTTc           <  1:1115636/68‑1 (MQ=255)

                                                     aGATCCGGTGAGCCCGGCAGACCTACAGCAAAGGAGTGCAACTGCGGCCACCAGGCTTCAGAGCGTTc           <  1:3402324/68‑1 (MQ=255)

                                                     aGATCCGGTGAGCCCGGCAGACCTACAGCAAAGGAGTGCAACTGCGGCCACCAGGCTTCAGAGCGTTc           <  1:1081859/68‑1 (MQ=255)

                                                     aGATCCGGTGAGCCCGGCAGACCTACAGCAAAGGAGTGCAACTGCGGCCACCAGGCTTCAGAGCGTTc           <  1:2331812/68‑1 (MQ=255)

                                                       aTCCGGTGAACCCGGCAGACCTACAGCAAAGGAGTGTAACTGCGGCCACCAGGCTTCAGAGCGTTCCTg        <  1:2527158/69‑1 (MQ=255)

                                                       aTCCGGTGAACCCGGCAGACCTACAGCAAAGGAGTGTAACTGCGGCCACCAGGCTTCAGAGCGTTCCTg        <  1:1337824/69‑1 (MQ=255)

                                                       aTCCGGTGAACCCGGCAGACCTACAGCAAAGGAGTGTAACTGCGGCCACCAGGCTTCAGAGCGTTCCTg        <  1:1308524/69‑1 (MQ=255)

                                                       aTCCGGTGAACCCGGCAGACCTACAGCAAAGGAGTGTAACTGCGGCCACCAGGCTTCAGAGCGTTCCTg        <  1:405778/69‑1 (MQ=255)

                                                       aTCCGGTGAACCCGGCAGACCTACAGCAAAGGAGTGTAACTGCGGCCACCAGGCTTCAGAGCGTTCCTg        <  1:495449/69‑1 (MQ=255)

                                                       aTCCGGTGAACCCGGCAGACCTACAGCAAAGGAGTGTAACTGCGGCCACCAGGCTTCAGAGCGTTCCTg        <  1:2483535/69‑1 (MQ=255)

                                                       aTCCGGTGAACCCGGCAGACCTACAGCAAAGGAGTGTAACTGCGGCCACCAGGCTTCAGAGCGTTCCTg        <  1:599487/69‑1 (MQ=255)

                                                       aTCCGGTGAACCCGGCAGACCTACAGCAAAGGAGTGTAACTGCGGCCACCAGGCTTCAGAGCGTTCCTg        <  1:1618323/69‑1 (MQ=255)

                                                         ccGGTGAACCCGGCAGACCTACAGCAAAGGAGTGTAACTGCGGCCACCAGGCTTCAGAGCGTTCCTGAt      >  1:2116902/1‑69 (MQ=255)

                                                             tGAACCCGGCAGACCTACAGCAAAGGAGTGTAACTGCGGCCACCAGGCTTCAGAGCGTTCCTGATCTTc  <  1:1403846/69‑1 (MQ=255)

                                                             tGAACCCGGCAGACCTACAGCAAAGGAGTGTAACTGCGGCCACCAGGCTTCAGAGCGGTCCTGATCTTc  <  1:2135243/69‑1 (MQ=255)

                                                             tGAACCAGGCAGACCTACAGCAAAGGAGTGTAACTGCGGCCACCAGGCTTCAGAGCGTTCCTGATCTTc  <  1:1820124/69‑1 (MQ=255)

GAATTTCGTGATGCACCGTGCCCAGATGGTTTGCCACTTCCTGGGCTGCTTTCAGATCCGGTGAACCCGGCAGACCTACAGCAAAGGAGTGTAACTGCGGCCACCAGGCTTCAGAGCGTTCCTGATCTTC  >  minE/421584‑421713

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