BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	2,415,862	A→G	14.3%	noncoding (17/77 nt)	proM ←	tRNA‑Pro

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	2,415,862	0	A	G	14.3%	140.8 / 13.2	56	noncoding (17/77 nt)	proM	tRNA‑Pro
Reads supporting (aligned to +/- strand): ref base A (28/20); new base G (5/3); total (33/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

TTCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTCAATTCTTTAAAGTCGTGGTGCGAG  >  Exported/2415794‑2415929

ttCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGc                                                                   >  1:5014009/1‑71 (MQ=12)

ttCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGc                                                                   >  1:468949/1‑71 (MQ=12)

  cAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTg                                                                 >  1:4132185/1‑71 (MQ=255)

      tttGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTgcgcc                                                             <  1:4180925/70‑2 (MQ=16)

      tttGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTgcgcc                                                             <  1:1308051/70‑2 (MQ=16)

      tttGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTgcg                                                               >  1:1922640/2‑69 (MQ=16)

      tttGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTgcg                                                               >  1:2772012/2‑69 (MQ=16)

             ggCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGcc                                                      <  1:1147009/71‑1 (MQ=255)

              gCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGcc                                                      >  1:1599202/1‑70 (MQ=255)

                gatagagGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGaa                                                   >  1:1185964/4‑71 (MQ=12)

                gagagaGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGa                                                    >  1:3291001/1‑70 (MQ=255)

                    gagGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTgcgcc                                                             <  1:4244743/57‑2 (MQ=12)

                           cGAACCTCCGACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCtttt                                        >  1:2900712/1‑71 (MQ=255)

                                  ccGACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCt                                                             >  1:561022/1‑43 (MQ=255)

                                  ccGACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTa                                                            >  1:2476629/1‑44 (MQ=255)

                                  ccGACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCg                                                     <  1:3706236/51‑1 (MQ=255)

                                  ccGACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAAttt                                 >  1:2233676/1‑71 (MQ=255)

                                   cGACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAAtttt                                <  1:2270917/71‑1 (MQ=255)

                                   cGACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAAtttt                                <  1:1273409/71‑1 (MQ=255)

                                   cGACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAAtttt                                <  1:4269871/71‑1 (MQ=255)

                                    gACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGaa                                                   >  1:5520604/1‑51 (MQ=255)

                                    gACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAAtttt                                >  1:291555/1‑70 (MQ=255)

                                         aCTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATAc                                               >  1:3173002/1‑50 (MQ=255)

                                          cTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGaaa                                                  <  1:2936474/46‑1 (MQ=255)

                                          cTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGaaa                                                  <  1:4971219/46‑1 (MQ=255)

                                          cTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTc                          >  1:1032025/1‑70 (MQ=255)

                                            ggTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATa                                                >  1:1049354/1‑46 (MQ=255)

                                            ggTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAAt                                   >  1:2381708/1‑59 (MQ=255)

                                             gTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTCaa                        <  1:3630416/69‑1 (MQ=255)

                                                ccAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTCAATTCtt                   >  1:3243/1‑71 (MQ=255)

                                                ccAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTCAATTCt                    <  1:2149250/70‑1 (MQ=255)

                                                ccAAACCAGTTGCGCTACCAAGCTGAGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTCAATTCt                    <  1:2672216/70‑1 (MQ=255)

                                                  aaaCCAGTTGCGCTACGAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTCAATTCttt                  >  1:2451534/1‑70 (MQ=255)

                                                  aaaCCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTCAATTCttt                  >  1:5201225/1‑70 (MQ=255)

                                                  aaaCCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTCAATTCttt                  >  1:2421618/1‑70 (MQ=255)

                                                  aaaCCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTCAATTCtt                   >  1:2426201/1‑69 (MQ=255)

                                                    aCCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGc                                            >  1:3526204/1‑42 (MQ=255)

                                                    aCCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTa                              >  1:270388/1‑56 (MQ=255)

                                                    aCCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTCAATTCTTTaaa               <  1:1735515/71‑1 (MQ=255)

                                                       aGTTTCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAAttt                                 >  1:5364138/1‑50 (MQ=255)

                                                       aGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAAttt                                 >  1:255122/1‑50 (MQ=255)

                                                       aGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAAttt                                 >  1:5595757/1‑50 (MQ=255)

                                                        gTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTCAATTCTTTAAAGTCg           <  1:5412913/71‑1 (MQ=255)

                                                        gTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTCAATTCTTTAAAGTCg           <  1:1000378/71‑1 (MQ=255)

                                                        gTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTCAATTCTTTAAAGTCg           <  1:1975606/71‑1 (MQ=255)

                                                        gTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTCAATTCTTTAAAGTCg           <  1:289974/71‑1 (MQ=255)

                                                           gcgcTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTCAATTCTTTAAAGTCgtgg        >  1:4462916/1‑71 (MQ=255)

                                                           gcgcTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTCAATTCTTTAAAGTCgtgg        >  1:5012254/1‑71 (MQ=255)

                                                           gcgcTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTCAATTCTTTAAAGTCgtgg        >  1:5018759/1‑71 (MQ=255)

                                                              cTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTCAATTCTTTAAAGTCgt          <  1:3102157/66‑1 (MQ=255)

                                                              cTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTCAATTCTTTAAAGTCGTGGTGc     >  1:11720/1‑71 (MQ=255)

                                                                aCCAAGCTGCGCTACTCGCCGAAATACTGCttttt                                       <  1:1743586/35‑1 (MQ=255)

                                                                aCCAAGCTGCGCTACTCGCCGAAATACTGCttttt                                       <  1:3913914/35‑1 (MQ=255)

                                                                 ccAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTCAATTCTTTAAAGTCGTGGTGCg    <  1:2909527/69‑1 (MQ=255)

                                                                 ccAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTCAATTCTTTAAAGTCGTGGTGCGAg  >  1:1376633/1‑71 (MQ=255)

                                                                  cAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTCAATTCTTTa                 <  1:3509937/55‑1 (MQ=255)

                                                                   aaGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTCAATTCTTTaaa               >  1:5012326/1‑56 (MQ=255)

                                                                   aaGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTCAATTCTTTaaa               >  1:5100500/1‑56 (MQ=255)

                                                                   aaGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTCAATTCTTTaaa               >  1:188979/1‑56 (MQ=255)

TTCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTCAATTCTTTAAAGTCGTGGTGCGAG  >  Exported/2415794‑2415929

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 27 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`