BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,212,787	A→G	35.7%	V97V (GTT→GTC)	rng ←	ribonuclease G

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,212,787	0	A	G	35.7%	31.4 / 31.5	56	V97V (GTT→GTC)	rng	ribonuclease G
Reads supporting (aligned to +/- strand): ref base A (28/8); new base G (13/7); total (41/15)
Fisher's exact test for biased strand distribution p-value = 3.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

CAGGCGCGCACCTTTAGTGCCAAGCGGATCTTTCACCACCTGCACCATCAGATCTTGCCCCTGACGAACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCACCCGCCACACATTCGGTGTG  >  minE/2212721‑2212852

cAGGCGCGCACCTTTAGTGCCAAGCGGATCTTTCACCACCTGCACCATCAGATCTTGCCCCTGACGaa                                                                  <  1:2044767/68‑1 (MQ=255)

cAGGCGCGCACCTTTAGTGCCAAGCGGATCTTTCACCACCTGCACCATCAGATCTTGCCCCTGACGaa                                                                  <  1:2860066/68‑1 (MQ=255)

cAGGCGCGCACCTTTAGTGCCAAGCGGATCTTTCACCACCTGCACCATCAGATCTTGCCCCTGACGaa                                                                  <  1:1711616/68‑1 (MQ=255)

cAGGCGCGCACCTTTAGTGCCAAGCGGATCTTTCACCACCTGCACCATCAGATCTTGCCCCTGACGaa                                                                  <  1:2950738/68‑1 (MQ=255)

 aGGCGCGCACCTTTAGTGCCAAGCGGATCTTTCACCACCTGCACCATCAGATCTTGCCCCTGACGaa                                                                  <  1:1298960/67‑1 (MQ=255)

 aGGCGCGCACCTTTAGTGCCAAGCGGATCTTTCACCACCTGCACCATCAGATCTTGCCCCTCACGaa                                                                  <  1:711402/67‑1 (MQ=255)

    cgcgcACCTTTAGTGCGAAGCGGATCTTTCACCACCTGCACCATCAGATCTTGCCCCTGACGAACCAGt                                                             >  1:2825616/1‑69 (MQ=255)

    cgcgcACCTTTAGTGCCAAGCGGATCTTTCCCCACCTGCACCATCAGATCTTGCCCCTGACGAACCAGt                                                             >  1:3364946/1‑69 (MQ=255)

    cgcgcACCTTTAGTGCCAAGCGGATCTTTCACCACCTGCACCATCAGATCTTGCCCCTGCCGAACCAGt                                                             >  1:3361058/1‑69 (MQ=255)

    cgcgcACCTTTAGTGCCAAGCGGATCTTTCACCACCTGCACCATCAGATCTTGCCCCTGACGAACCAGt                                                             >  1:301415/1‑69 (MQ=255)

    cgcgcACCTTTAGTGCCAAGCGGATCTTTCACCACCTGCACCATCAGATCTTGCCCCTGACGAACCAGt                                                             >  1:514839/1‑69 (MQ=255)

    cgcgcACCTTTAGTGCCAAGCGGATCTTTCACCACCTGCACCATCAGATCTTGCCCCTGACGAACCAGt                                                             >  1:9101/1‑69 (MQ=255)

    cgcgcACCTTTAGTGCCAAGCGGATCTTTCACCACCTGCACCATCAGATCTTGCCCCTGACGAACCAGt                                                             >  1:1868342/1‑69 (MQ=255)

    cgcgcACCTTTAGTGCCAAGCGGATCTTTCACCACCTGCACCATCAGATCTTGCCCCTGACGAACCAGt                                                             >  1:1102759/1‑69 (MQ=255)

         aCCTTTAGTGCCAAGCGGATCTTTCACCACCTGCACCATCAGTTCTTGCCCCTGACGAACCAGTTCCga                                                        <  1:12038/69‑1 (MQ=255)

         aCCTTTAGTGCCAAGCGGATCTTTCACCACCTGCACCATCAGATCTTGCCCCTGACGAACCAGTTCCga                                                        <  1:3079237/69‑1 (MQ=255)

         aCCTTTAGTGCCAAGCGGATCTTTCACCACCTGCACCATCAGATCTTGCCCCTGACGAACCAGTTCCga                                                        <  1:2944983/69‑1 (MQ=255)

          ccTTTAGTGCCAAGCGGATCTTTCACCACCTGCACCATCAGATCTTGCCCCTGACGAACCAGTTCCga                                                        <  1:46089/68‑1 (MQ=255)

           cTTTAGTGCCAAGCGGATCTTTCACCACCTGCACCATCAGATCTTGCCCCTGACGAACCAGTTCCga                                                        <  1:2851933/67‑1 (MQ=255)

                         ggATCTTTCACCACCTGCACCATCAGATCTTGCCCCTGACGAACCAGTTCCGAGATGTCGCGCACCGTg                                        <  1:2403578/69‑1 (MQ=255)

                         ggATCTTTCACCACCTGCACCATCAGATCTTGCCCCTGACGAACCAGTTCCGAGATGTCGCGCACCGTg                                        <  1:1623074/69‑1 (MQ=255)

                         ggATCTTTCACCACCTGCACCATCAGATCTTGCCCCTGACGAACCAGTTCCGAGATGTCGCGCACCGTg                                        <  1:2261462/69‑1 (MQ=255)

                            tCTTTCACCACCTGCACCGTCAGATCTTGCCCCTGACGAACCAGTTCCGAGATGTCGCGCACCGTGAAt                                     >  1:3295329/1‑69 (MQ=255)

                            tCTTTCACCACCTGCACCATCAGATCTTTCCCCTGACGAACCAGTTCCGAGATGTCGCGCACCGTGaa                                      >  1:783522/1‑68 (MQ=255)

                            tCTTTCACCACCTGCACCATCAGATCTTGCCCCTGACGaa                                                                  >  1:2344721/1‑40 (MQ=255)

                            tCTTTCACCACCTGCACCATCAGATCTTGCCCCTGACGAACCAGTTCCGAGATGTCGCGCACCGTGaa                                      >  1:1739513/1‑68 (MQ=255)

                            tCTTTCACCACCTGCACCATCAGATCTTGCCCCTGACGAACCAGTTCCGAGATGTCGCGCACCGTGAAt                                     >  1:1117724/1‑69 (MQ=255)

                            tCTTTCACCACCTGCACCATCAGATCTTGCCCCTGACGAACCAGTTCCGAGATGTCGCGCACCGTGAAt                                     >  1:3413775/1‑69 (MQ=255)

                            tCTTTCACCACCTGCACCATCAGATCTTGCCCCTGACGAACCAGTTCCGAGATGTCGCGCACCGTGAAt                                     >  1:218771/1‑69 (MQ=255)

                            tCTTTCACCACCTGCACCATCAGATCTTGCCCCTAACGAACCAGTTCCGAGATGTCGCGCACCGTGAAt                                     >  1:220872/1‑69 (MQ=255)

                                      ccTGCACCATCAGATCTTGCCCCTGACGAACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGtt                           >  1:1113943/1‑69 (MQ=255)

                                      ccTGCACCATCAGATCTTGCCCCTGACGAACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGtt                           >  1:3159877/1‑69 (MQ=255)

                                              aTCAGATCTTGCCCCTGACGAACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCAccc                   >  1:751579/1‑69 (MQ=255)

                                              aTCAGATCTTGCCCCTGACGAACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCAccc                   >  1:3262307/1‑69 (MQ=255)

                                              aTCAGATCTTGCCCCTGACGAACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCAccc                   >  1:612994/1‑69 (MQ=255)

                                              aTCAGATCTTGCCCCTGACGAACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCAccc                   >  1:3180179/1‑69 (MQ=255)

                                              aTCAGATCTTGCCCCTGACGAACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCAccc                   >  1:2400552/1‑69 (MQ=255)

                                              aTCAGATCTTGCCCCTGACGAACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCAccc                   >  1:2756089/1‑69 (MQ=255)

                                              aTCAGATCTTGCCCCTGACGAACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCAcc                    >  1:1447839/1‑68 (MQ=255)

                                                    tCTTGCCCCTGACGGACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCAcc                    >  1:1415844/1‑62 (MQ=255)

                                                    tCTTGCCCCTGACGGACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCACCCGCcaca             >  1:992543/1‑69 (MQ=255)

                                                    tCTTGCCCCTGACGGACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCACCCGCcaca             >  1:111242/1‑69 (MQ=255)

                                                    tCTTGCCCCTGACGGACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCACCCGCcaca             >  1:866520/1‑69 (MQ=255)

                                                    tCTTGCCCCTGACGGACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCACCCGCcaca             >  1:790580/1‑69 (MQ=255)

                                                    tCTTGCCCCTGACGGACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCACCCGCcaca             >  1:220521/1‑69 (MQ=255)

                                                    tCTTGCCCCTGACGGACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCACCCGCcaca             >  1:1134658/1‑69 (MQ=255)

                                                    tCTTGCCCCTGACGGACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCACCCGCcaca             >  1:1827720/1‑69 (MQ=255)

                                                    tCTTGCCCCTGACGGACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCACCCGCcaca             >  1:1200644/1‑69 (MQ=255)

                                                    tCTTGCCCCTGACGGACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCACCCGCcaca             >  1:3361842/1‑69 (MQ=255)

                                                    tCTTGCCCCTGACGGACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCACCCGCcaca             >  1:2896735/1‑69 (MQ=255)

                                                    tCTTGCCCCTGACGGACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCACCCGCcaca             >  1:3157752/1‑69 (MQ=255)

                                                    tCTTGCCCCTGACGGACCAGTTCCGAGATGACGCGCACCGTGAATTGCt                                 >  1:1228422/1‑49 (MQ=39)

                                                     cTTGCCCCTGACGAACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCACCCGCcacac            >  1:464001/1‑69 (MQ=255)

                                                     cTTGCCCCTGACGAACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCACCCGCcacac            >  1:175622/1‑69 (MQ=255)

                                                           ccTGACGGACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCACCCGCCACACATTcag      <  1:3364406/69‑3 (MQ=255)

                                                           ccTGACGGACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCACCCGCCACACATTcag      <  1:563820/69‑3 (MQ=255)

                                                           ccTGACGGACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCACCCGCCACACATTcag      <  1:2498330/69‑3 (MQ=255)

                                                           ccTGACGGACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCACCCGCCACACATTcag      <  1:335975/69‑3 (MQ=255)

                                                           ccTGACGGACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCACCCGCCACACATTcag      <  1:1054350/69‑3 (MQ=255)

                                                           ccTGACGGACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCACCCGCCACACATTcag      <  1:974353/69‑3 (MQ=255)

                                                            cTGACGGACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCACCCGCCACACATTcag      <  1:2350614/68‑3 (MQ=255)

                                                               aCGAACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCACCCGCCACACATTCGGtgtg  >  1:784404/1‑69 (MQ=255)

                                                               aCGAACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCACCCGCCACACATTCGGtgtg  >  1:210377/1‑69 (MQ=255)

CAGGCGCGCACCTTTAGTGCCAAGCGGATCTTTCACCACCTGCACCATCAGATCTTGCCCCTGACGAACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCACCCGCCACACATTCGGTGTG  >  minE/2212721‑2212852

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