BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	2,514,509	G→T	6.9%	intergenic (+24/‑23)	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,509	0	G	T	6.9%	1948.5 / 14.3	668	intergenic (+24/‑23)	valU/valY	tRNA‑Val/tRNA‑Val
Reads supporting (aligned to +/- strand): ref base G (241/362); new base T (15/50); total (256/412)
Fisher's exact test for biased strand distribution p-value = 7.26e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 678 total aligned reads displayed.

TACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGG  >  CP009273/2514443‑2514576

tACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGt                                                                    >  1:52492570/1‑68 (MQ=255)

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

tACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGt                                                                    >  1:38799637/1‑68 (MQ=255)

tACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTACTGGCCAGATAAATTtgttgt                                                                    >  1:36585168/1‑68 (MQ=17)

  cAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTACTCGCGAGCTAAATTTCTTGTaa                                                                  >  1:17012954/1‑68 (MQ=255)

   aaGGAGGGGGTCGGCGGTTCGCTCCCGTCATCACCCCCCACTTACGCGCCAGCGAAAATTCGGGTaaa                                                                 >  1:1626985/1‑68 (MQ=14)

   aaGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTGGCCATCGAAATTtgttgttat                                                                 >  1:27315541/1‑65 (MQ=17)

   aaGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCTCGAGATAAATTTATTGTtaa                                                                 >  1:36316391/1‑68 (MQ=17)

   aaGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTAGTCGCCAGCTATATGTGTTGTaat                                                                 >  1:3839150/1‑67 (MQ=14)

   aaGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTACTCGCGTGATAAATTTATTTTaaa                                                                 >  1:14695738/1‑68 (MQ=14)

   aaGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTACTCGCCTTATAAAATTGTTGTaat                                                                 >  1:55732108/1‑67 (MQ=14)

   aaGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTACTCGCCATCTAAATTTGTTGaaat                                                                 >  1:19041312/1‑67 (MQ=17)

   aaGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTACTCGCCAGCTAAATTTATTGTaaa                                                                 >  1:26051680/1‑68 (MQ=255)

   aaGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTACGCGCCAGCTAAATTTCTTGTaaa                                                                 >  1:22607869/1‑68 (MQ=255)

     ggaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCTCCAGCTATAATTGTTGTtaaaa                                                               >  1:48338579/1‑68 (MQ=17)

     ggaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTaaaaa                                                               >  1:50195928/1‑68 (MQ=255)

     ggaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTaaaaa                                                               >  1:20904559/1‑68 (MQ=255)

      gaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTGTGGCCAGAGATATTTCTTGTAAAAAt                                                              >  1:49158857/1‑68 (MQ=17)

       aggGGGGCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAAt                                                              <  1:11513417/67‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

               ggcggTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTAc                                                     <  1:39432881/68‑1 (MQ=255)

                  gggtCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGa                                                  <  1:13217875/65‑1 (MQ=255)

                  ggTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCg                                                   <  1:39812188/67‑1 (MQ=255)

                  ggTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGa                                                  >  1:5815665/1‑68 (MQ=255)

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

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

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

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

                          cccGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATTTGAAGTAAAGAAGtgtgtg                                          >  1:32913683/1‑68 (MQ=38)

                          cccGTCATCACCCACCACTTTCTCGCCAGCGATATTTCTTTTAAAAATTTGAATTACCGAAGTGGGt                                           >  1:31820620/1‑67 (MQ=37)

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

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

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

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

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

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

                               catcaCCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTa                                      <  1:18384254/67‑1 (MQ=255)

                               catcaCCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAg                                     <  1:53822816/68‑1 (MQ=255)

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

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

                                  cacccaccACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTc                                  <  1:45041480/68‑1 (MQ=255)

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

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

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

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

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

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

                                    cccaccACTTGCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGGGGGTGATTAGCTCAg                                <  1:31123906/68‑1 (MQ=255)

                                     ccaccaCTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAACTGGGTGATTAGCTCa                                 >  1:8928545/1‑66 (MQ=255)

                                     acaccacTTTCTCGCCATCTAAATATCTCTAAAAAATGTCAAGACCCGCAGTGGGTGATTAGCTCAGc                               <  1:12270194/67‑1 (MQ=255)

                                      cagcacTTTCTCGCCAGCTAACTTTCTTGTAAAAATGTGACGGACCGAAGTGGGTGATTAGCTCAGCt                              <  1:26455351/65‑1 (MQ=255)

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

                                         cacTTTCTCGCCAGCTAAATTTCTGGTAAAAATGAGACCTACCGAACTGGGGGATTAGCTCAGCTggg                           <  1:20314981/68‑1 (MQ=255)

                                         cacTTTCAAGCCTGCTAAATTTCATGAAAAAATGAGAACGACCTAAGGGGGGGATTAGCTCAGCTggg                           <  1:15730325/68‑1 (MQ=16)

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

                                           cTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGgag                         >  1:643485/1‑68 (MQ=255)

                                             ttCTCGCCATCTAAATATCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGagag                       <  1:54854967/68‑1 (MQ=255)

                                             ttCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGgaga                        >  1:30453771/1‑67 (MQ=255)

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

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

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

                                               ctcGACAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGc                      >  1:370212/1‑67 (MQ=255)

                                                 cGCCAGCTAAATTTCTTTTAAAATTGTGAAGTACCGTAGTGGGTGATTAGCTCAGCTGGGAGAGCAcc                   <  1:21673558/68‑1 (MQ=37)

                                                 cGCCAGCTAAATTTCTTGTAAAAAAGTGAACTACCGCAGTGGGTGATTAGCTCAGCTGGGAGAGCAcc                   <  1:37880402/68‑1 (MQ=37)

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

                                                  gCCAGCTAAATTTCTCTTAAAACCGAGAACTAACGACGTGGGTGATTAGCTCAGCTGGGAGAGCAcc                   <  1:54354034/67‑1 (MQ=18)

                                                  gCCAGCTAAATATCTTGTAAAACTGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCt                  <  1:33316910/68‑1 (MQ=255)

                                                   ccATCTAAATTTCTTGTAAAAATGTGAAGTCCCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCt                  <  1:12841939/67‑1 (MQ=255)

                                                    catcTAAATTTCTTGTAAAAATGTGAACTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTcc                <  1:41172899/65‑1 (MQ=255)

                                                    catcTAAATTTCATTTAAAACTGTGAAGACCCGACGGGGGTGATTAGCTCAGCTGGGAGAGCACCTcc                <  1:46738271/65‑1 (MQ=17)

                                                    catcTAAATATCTCGTAAAAACGACAAGTCCCGAAGGGGGTGATTAGCTCAGCTGGGAGAGCACCTcc                <  1:24959061/65‑1 (MQ=17)

                                                    cAGCTAACTTTCTCTAAAAACTGTCAACTCCCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTcc                <  1:35508606/68‑1 (MQ=17)

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

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

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

                                                    attcTAAATTTCTTGTTAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTcc                <  1:31327924/65‑1 (MQ=255)

                                                     atcTAAATTTCTTGTAAAAAAGTGAAGTACCGCCGTGGGTGATTAGCTCAGCTGGGAGAGCACCTccc               <  1:17654555/66‑1 (MQ=25)

                                                      gCTAAATATCTTGTAAAACTGTGAAGTCCCCAAGGGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCt              <  1:50867490/68‑1 (MQ=18)

                                                        tAAATTTCTTTTAAAATTGTGAAGTAACTTAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCt              <  1:49889467/66‑1 (MQ=18)

                                                        tAAATTTCTTTTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTccc               <  1:25609411/65‑1 (MQ=255)

                                                        tAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTccc               <  1:46168954/65‑1 (MQ=255)

                                                        tAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTccc               <  1:10632426/65‑1 (MQ=255)

                                                        tAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCt              <  1:53340769/66‑1 (MQ=255)

                                                        tAAATTTCTTGTAAAAATGTGAACTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTa            <  1:40731547/68‑1 (MQ=255)

                                                        tAAATTTCATGTAAAACTGTGAAGAAACGCAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTcc                <  1:28158771/64‑1 (MQ=18)

                                                        tAAATATCTTGTAAAAATGTGAAGTACCGACGGGGGTGATTAGCTCAGCTGGGAGAGCACCTccc               <  1:35648091/65‑1 (MQ=21)

                                                            tttCTTGTAAAAATGTGAAGTACCGCAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAag        <  1:33773288/68‑1 (MQ=255)

                                                             ttCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAagg       >  1:3253799/1‑68 (MQ=255)

                                                             ttCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAagg       >  1:56320971/1‑68 (MQ=255)

                                                             ttCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAagg       <  1:29459439/68‑1 (MQ=255)

                                                                ttGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAaggagg    >  1:42070958/1‑68 (MQ=255)

                                                                 tGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAggg   >  1:13786584/1‑68 (MQ=255)

TACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGG  >  CP009273/2514443‑2514576

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