BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,415,848	C→G	23.1%	noncoding (31/77 nt)	proM ←	tRNA‑Pro

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,415,848	0	C	G	23.1%	115.9 / 26.1	52	noncoding (31/77 nt)	proM	tRNA‑Pro
Reads supporting (aligned to +/- strand): ref base C (27/13); new base G (6/6); total (33/19)
Fisher's exact test for biased strand distribution p-value = 3.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.31e-01

GAAGCGGGGTTCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTCAATTCTTTAAA  >  minE/2415785‑2415916

gAAGCGGGGTTCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTgcg                                                               >  1:3122587/1‑71 (MQ=255)

gAAGCGGGGTTAAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTgcg                                                               >  1:3793516/1‑71 (MQ=255)

   gCGGGGTTCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTgcg                                                               >  1:577056/1‑68 (MQ=255)

      gggTTCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAAcc                                                                      <  1:5097850/58‑1 (MQ=255)

       ttttCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCag                                                        >  1:3108423/3‑70 (MQ=12)

         ttCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCag                                                        <  1:2367851/69‑2 (MQ=12)

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

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

         ttCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTgc                                                                <  1:4029939/61‑1 (MQ=255)

           cAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTGCGCTa                                                            >  1:3907592/1‑63 (MQ=255)

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

            aaaaTTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCa                                                                     <  1:1771869/53‑1 (MQ=255)

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

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

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

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

               aTTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTGCGCTa                                                            >  1:1358257/1‑59 (MQ=255)

               aTTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTGCGCTa                                                            >  1:2281846/1‑59 (MQ=255)

               aTTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTGCGCTa                                                            >  1:4671372/1‑59 (MQ=255)

                ttGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAAcc                                                                      <  1:4283347/48‑1 (MQ=255)

                ttGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAAcc                                                                      <  1:1723109/48‑1 (MQ=255)

                  ggTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTAc                                                           <  1:3823222/57‑1 (MQ=12)

                      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)

                          atagagGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTAc                                                           <  1:3747175/47‑1 (MQ=9)

                             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:1273409/71‑1 (MQ=255)

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

                                            cGACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAAtttt                   <  1:2270917/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:2421618/1‑70 (MQ=255)

                                                           aaaCCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTCAATTCttt     >  1:5201225/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)

GAAGCGGGGTTCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTCAATTCTTTAAA  >  minE/2415785‑2415916

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