BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	2,514,497	G→T	5.5%	intergenic (+12/‑35)	valU → / → valY	tRNA‑Val/tRNA‑Val

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	2,514,497	0	G	T	5.5%	1007.9 / 2.3	326	intergenic (+12/‑35)	valU/valY	tRNA‑Val/tRNA‑Val
Reads supporting (aligned to +/- strand): ref base G (124/176); new base T (15/10); total (140/186)
Fisher's exact test for biased strand distribution p-value = 9.16e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 330 total aligned reads displayed.

AGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCC  >  CP009273/2514430‑2514562

 gagCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAga                                                                  >  1:20280860/1‑67 (MQ=255)

 gagCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAg                                                                   <  1:25215716/67‑1 (MQ=255)

   gCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGc                                                                  <  1:12879360/66‑1 (MQ=255)

     aCCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCTCCAGCTata                                                              >  1:11598116/1‑66 (MQ=255)

      ccTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTACTCGCCATCTAAAt                                                             >  1:16308392/1‑68 (MQ=255)

          ccTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTGTGGCCTGCTATAtttag                                                         >  1:9305376/1‑66 (MQ=14)

             tACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGt                                                      <  1:14033661/68‑1 (MQ=255)

                aaGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTaaa                                                   >  1:2093321/1‑68 (MQ=255)

                aaGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTGGTCGCGAGATAAATTTATTGTaat                                                   >  1:29240474/1‑67 (MQ=14)

                aaGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTACTCGCGTGATAAATTTATTGTtaa                                                   >  1:3458883/1‑68 (MQ=14)

                aaGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTACTCGCCATCTAAAATTCTTGttat                                                   >  1:7217887/1‑65 (MQ=17)

                  ggaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCGCGCCAGCTAAATTTGTTGTaagca                                                 >  1:17096098/1‑65 (MQ=18)

                  ggaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTACTCGCCAGCTAAATTTGGTGTaaaaa                                                 >  1:4707573/1‑68 (MQ=18)

                    aggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCTGCTATATTTCTGGTAAAAAtg                                               >  1:14760429/1‑68 (MQ=21)

                    aggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTaaaaa                                                 <  1:27592220/66‑1 (MQ=255)

                    aggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAAtg                                               >  1:29355468/1‑68 (MQ=255)

                     gggggTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCGAGCTATATTTGTGGTTAAAAtgt                                              >  1:24292836/1‑68 (MQ=18)

                      ggggTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAAtgtg                                             >  1:7900937/1‑68 (MQ=255)

                       gggTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAAtgt                                              <  1:25491748/66‑1 (MQ=255)

                       gggTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAAtgt                                              <  1:12580147/66‑1 (MQ=255)

                       gggTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAAtgt                                              <  1:7949094/66‑1 (MQ=255)

                       gggTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAAtgt                                              <  1:15295447/66‑1 (MQ=255)

                       gggTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAAtgt                                              <  1:17617201/66‑1 (MQ=255)

                       gggTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAAtgt                                              <  1:3872986/66‑1 (MQ=255)

                       gggTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAAtgt                                              <  1:9063810/66‑1 (MQ=255)

                        ggTCGGCGGTTCGATCCCGTCATCACCCACCACTTACTCGCCATCTAAATTTGGTGTAAAAATGTGaa                                           >  1:5928786/1‑68 (MQ=21)

                         gTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAg                                          >  1:8747873/1‑68 (MQ=255)

                          tCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAAATTCTTGTTAAAATGTGAAGt                                         >  1:14234335/1‑68 (MQ=255)

                           gggcggTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGt                                         <  1:10703737/66‑1 (MQ=255)

                           cggcggTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGt                                         <  1:20490060/67‑1 (MQ=255)

                           cggcggTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGt                                         <  1:22049350/67‑1 (MQ=255)

                           cggcggTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGt                                         <  1:4989385/67‑1 (MQ=255)

                           cggcggTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGt                                         <  1:26129024/67‑1 (MQ=255)

                           cggcggTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGt                                         <  1:27036151/67‑1 (MQ=255)

                           cggcggTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGt                                         <  1:19998942/67‑1 (MQ=255)

                           cggcggTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGt                                         <  1:2669571/67‑1 (MQ=255)

                           cggcggTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGt                                         <  1:3080994/67‑1 (MQ=255)

                           cggcggTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGt                                         <  1:9543544/67‑1 (MQ=255)

                           cggcggTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGt                                         <  1:3258368/67‑1 (MQ=255)

                           cggcggTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGt                                         <  1:11997505/67‑1 (MQ=255)

                           cggcggTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGt                                         <  1:13334341/67‑1 (MQ=255)

                           cggcggTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGt                                         <  1:1347443/67‑1 (MQ=255)

                             gcggTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTATTTTAAAAATGTGAAGTAca                                      >  1:21698303/1‑67 (MQ=37)

                              cggTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTTAAAATGTGAATTaact                                     >  1:6894847/1‑65 (MQ=37)

                                  tCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAg                                  <  1:25999832/67‑1 (MQ=255)

                                  tCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGt                                 <  1:7722730/68‑1 (MQ=255)

                                  tCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGt                                 <  1:23595298/68‑1 (MQ=255)

                                   cGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTg                                <  1:17315834/68‑1 (MQ=255)

                                   cGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTg                                <  1:23271954/68‑1 (MQ=255)

                                      tCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGt                             <  1:29719203/68‑1 (MQ=255)

                                      tCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGt                             <  1:284462/68‑1 (MQ=255)

                                      tCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGt                             <  1:6737481/68‑1 (MQ=255)

                                       cccGTCATCACCCACCACTTTCTCGCCAGCTATATTTCTTGTAAAAATGTGAAGTACCGAATTGGGtt                            >  1:8926822/1‑67 (MQ=39)

                                       cccGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTg                            >  1:20738908/1‑68 (MQ=255)

                                       cccGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAAGGTGAAGTACCGAAGTGGGTg                            >  1:19204913/1‑68 (MQ=255)

                                        ccGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGa                           <  1:27867589/68‑1 (MQ=255)

                                        ccGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGa                           <  1:7378939/68‑1 (MQ=255)

                                        ccGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGa                           <  1:2498902/68‑1 (MQ=255)

                                          gTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGAtt                         >  1:25586161/1‑68 (MQ=255)

                                            catcaCCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTa                        >  1:7002954/1‑67 (MQ=255)

                                            catcaCCCACCACTTTCTCGCCAGATAAATTTGTTTTAAAAATGTGAACTAACGAAGTTGGTGAtaa                        >  1:10225408/1‑65 (MQ=255)

                                             atcaCCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAg                       <  1:11002403/67‑1 (MQ=255)

                                             atcaCCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGc                      <  1:16589514/68‑1 (MQ=255)

                                                acccgccACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCa                   >  1:27903125/1‑68 (MQ=255)

                                                acccaccACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCTAAGTGGGTGATTAGCTCa                   <  1:19787836/68‑1 (MQ=255)

                                                acccaccACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCa                   >  1:18549318/1‑68 (MQ=255)

                                                acccaccACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCa                   >  1:1380403/1‑68 (MQ=255)

                                                acccaccACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCa                   >  1:5582726/1‑68 (MQ=255)

                                                acccaccACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCa                   >  1:11474536/1‑68 (MQ=255)

                                                 cccaccACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCa                   >  1:19930996/1‑67 (MQ=255)

                                                 cccaccACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCa                   >  1:28852047/1‑67 (MQ=255)

                                                 cccaccACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCa                   >  1:22916892/1‑67 (MQ=255)

                                                 cccaccACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAg                  <  1:26500584/68‑1 (MQ=255)

                                                 cccaccACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAg                  <  1:22352519/68‑1 (MQ=255)

                                                 cccaccACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAg                  <  1:18884188/68‑1 (MQ=255)

                                                  ccaccaCTTTCTGGCGTGATAAATTTCTTGAAAAACTTTTAAGTACCCATTTGGGTGATTAGCTCAGc                 <  1:9861674/68‑1 (MQ=255)

                                                  ccaccaCTTTCTCGCCATCTAAATTTCTTTTAAAAATGTGAAGTCCCGAAGGGGGTGATTAGCTCAGc                 <  1:17832138/68‑1 (MQ=255)

                                                  ccaccaCTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGc                 <  1:18015154/68‑1 (MQ=255)

                                                  ccaccaCTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGc                 <  1:8377517/68‑1 (MQ=255)

                                                     ccacTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTgg              >  1:28112480/1‑68 (MQ=255)

                                                     ccacTTACTCGCCACCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGATAAGCTgg              >  1:10486896/1‑68 (MQ=255)

                                                      cacTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTggg             >  1:3519142/1‑68 (MQ=255)

                                                       acTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGga            <  1:6105820/68‑1 (MQ=255)

                                                       acTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGga            >  1:23822248/1‑68 (MQ=255)

                                                       acTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGga            >  1:15656267/1‑68 (MQ=255)

                                                       acTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGga            >  1:19503048/1‑68 (MQ=255)

                                                         tttCTCGCCAGCTAAATTTCTTTTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGgaga          <  1:27518240/68‑1 (MQ=255)

                                                            ctcGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGc        >  1:7755791/1‑67 (MQ=255)

                                                            ctcGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGc        >  1:18520012/1‑67 (MQ=255)

                                                            ctcGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCa       <  1:12356492/68‑1 (MQ=255)

                                                              cGCCAGCTAAATTTCTTTTAAAAATGTCAAGAACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCAcc     <  1:22751644/68‑1 (MQ=37)

                                                               gCCATCTAAATCTCTTTTTAAAAAGTCAACTACCCAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCt    <  1:55148/68‑1 (MQ=18)

                                                               gCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCAcc     >  1:14421467/1‑67 (MQ=255)

                                                                ccAGCTAAATTTCTTGTAAAAATGTCAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTc   <  1:12834399/68‑1 (MQ=255)

                                                                 cAGCTAAATTTCTTGTAAAATTGTGAAGTACCGAAGTGGGGGATTAGCTCAGCTGGGAGAGCACCTc   <  1:17258638/67‑1 (MQ=38)

                                                                 cAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTcc  <  1:15062548/68‑1 (MQ=255)

                                                                 cAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTcc  <  1:11320226/68‑1 (MQ=255)

                                                                 cAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTcc  <  1:5083818/68‑1 (MQ=255)

                                                                 cAGCTAAATTTCTTGTAAAAATGTCAAGTACCGACGGGGGTGATTAGCTCAGCTGGGAGAGCACCTc   <  1:10047088/67‑1 (MQ=25)

                                                                 cAGCTAAATATCATGTAAAAATGTGAAGAACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTcc  <  1:12013588/68‑1 (MQ=25)

AGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCC  >  CP009273/2514430‑2514562

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