BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	2,514,509	G→T	6.5%	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.5%	1804.6 / 17.9	613	intergenic (+24/‑23)	valU/valY	tRNA‑Val/tRNA‑Val
Reads supporting (aligned to +/- strand): ref base G (208/348); new base T (9/47); total (218/395)
Fisher's exact test for biased strand distribution p-value = 1.14e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 618 total aligned reads displayed.

TACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGG  >  CP009273/2514443‑2514575

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

 aCAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTa                                                                  <  1:12178307/68‑1 (MQ=255)

   aaGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTGTCGCGAGCTAAATTTCTTTTaaa                                                                >  1:41133765/1‑68 (MQ=18)

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

   aaGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCGAAAATTGTTGTaat                                                                >  1:3644459/1‑67 (MQ=18)

   aaGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCGCGCGAGCGAAATTTGTTGaaaa                                                                >  1:4667795/1‑68 (MQ=17)

   aaGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCGCGCCAGATAAATTTCTTGTtaa                                                                >  1:46326114/1‑68 (MQ=18)

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

       aggGGGGCGGCGGTTCGATCCCGGCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTaaaaa                                                              <  1:9185525/66‑1 (MQ=255)

          gggTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCTTCTAAATTTATTTTTAAAATGTGa                                                         >  1:50146212/1‑68 (MQ=18)

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

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

           ggTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCGATAGAAATTTTTGGTTAGAATGTGaa                                                        >  1:8183483/1‑68 (MQ=14)

             tCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTGGCCAGATATATGTAGTGTAAGAAGGTGAAGt                                                      >  1:14702384/1‑68 (MQ=17)

             tCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCTCCAGCTAAATTTCTTGTAAAAATTTGAAGt                                                      >  1:25957631/1‑68 (MQ=255)

              cggcggTTCGATGCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGt                                                      <  1:44153118/67‑1 (MQ=255)

              cggcggTTCGATCCCGTCATCACCCACCTCTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGt                                                      <  1:42754831/67‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

                  ggTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGa                                                 <  1:45819802/68‑1 (MQ=255)

                  ggTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGa                                                 <  1:31351254/68‑1 (MQ=255)

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

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

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

                          cccGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTTATAATGTGAAGTACCGAAGTGGGTg                                         >  1:25372077/1‑68 (MQ=255)

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

                          cccGTCATCACCCACCACTTTCTCGACAGCTAAATTTCT‑‑TAAAAATGTTTAGTACCGAAGTTGGTGAt                                       >  1:40285272/1‑68 (MQ=25)

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

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

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

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

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

                              tcatcaCCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTa                                     >  1:2207763/1‑68 (MQ=255)

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

                               catcaCCCACCACTTTCTCGCCAGATAAATTTATTGTAAAAATGTGAAGTACCGAAGTGGGTTATTAg                                    >  1:5479087/1‑68 (MQ=255)

                                 tcaCCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGc                                   >  1:16349336/1‑67 (MQ=255)

                                   acccaccACTTTCTCGCCTGCTAAATTTCTTGTAAAAATGTGAAGTCCCGAAGTGGGTGATTAGCTCa                                <  1:22674436/68‑1 (MQ=255)

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

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

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

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

                                    cccaccACTTTCTCGCCACCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCa                                >  1:43858373/1‑67 (MQ=255)

                                    cccaccACTTTCTCGACAGCTAAATATCTTGTAAAACTGTGAAGTACCGAAGTGGGTGATTAGCTCAg                               <  1:21048018/68‑1 (MQ=255)

                                    cacacgacTTTCTCGCCATCTAAATTTCTTGTAAAACTGTGAAGTACCGAAGTGGGTGATTAGCTCAg                               <  1:6655980/66‑1 (MQ=255)

                                     ccacgaCTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGc                              <  1:5289428/68‑1 (MQ=255)

                                     ccaccaCTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCa                                >  1:13447211/1‑66 (MQ=255)

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

                                     ccaccaCTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGGGGGTGATTAGCTCAGc                              <  1:16007889/68‑1 (MQ=255)

                                     ccaccaCTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGGGGGTGATTAGCTCAGc                              <  1:2055011/68‑1 (MQ=255)

                                     ccaccaCTTTATCGCCAGCTAAATTTCTTGTTAAACGGTGACGTACCGAAGTGGGTGATTAGCTCAGc                              <  1:17173610/68‑1 (MQ=255)

                                     caaccacTTTCTCGCGTGTTAAATATCTTTTAAATTTGTGTAGTACCGTAGTGGGTGATTAGCTCAGc                              <  1:47070444/66‑1 (MQ=255)

                                       accacTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTg                            >  1:21604561/1‑68 (MQ=255)

                                       accacTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTg                            >  1:45306775/1‑68 (MQ=255)

                                         cacTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGCAGTGGGTGATTAGCTCAGCTggg                          <  1:19988441/68‑1 (MQ=255)

                                         cacTTTCTCCCCAGCTAAATTTCTTGTAAAAATGTCAAGTACCGAAGTGGGTGATTAGCTCAGCTggg                          <  1:779738/68‑1 (MQ=255)

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

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

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

                                               ctcGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTggg                                           >  1:3992237/1‑45 (MQ=255)

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

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

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

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

                                                  gCCAGCTAAATATCTTGTAAATCTGAGAAGTCCCGTCGGGGGTGATTAGCTCAGCTGGGAGAGCACCt                 <  1:15542505/68‑1 (MQ=17)

                                                    tatcTAACTTTCTTTTAAAACTGTGAAGTACCGTAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTcc               <  1:1833087/65‑1 (MQ=21)

                                                    tatcTAAATTTCTTTTTAAAATGAGTAGAACCGGCGTGGGTGATTAGCTCAGCTGGGAGAGCACCTcc               <  1:2344748/65‑1 (MQ=17)

                                                    cAGCTAAATTTCTTTTAAAAATGTCAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTcc               <  1:41864841/68‑1 (MQ=255)

                                                    cAGCTAAATTTCTTGTAAAAATGTGAAGTACCGCAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTcc               <  1:30055893/68‑1 (MQ=255)

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

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

                                                    cAGCTAAATTTCTTGAAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTcc               <  1:12997720/68‑1 (MQ=255)

                                                    cAGCTAAATTTCATTAAAAAATGACAAGTACCCAAGGGGGTGATTAGCTCAGCTGGGAGAGCACCTcc               <  1:33525990/68‑1 (MQ=17)

                                                    cAGCTAAATATCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTcc               <  1:49084581/68‑1 (MQ=255)

                                                    aatcTAAATTTCTTGTAAAACCGTGAAGTACCTAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTc                <  1:47445360/64‑1 (MQ=25)

                                                      tcTAACTTTCATTTAAAAACGTGTACTACCTACGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCt             <  1:36272117/67‑1 (MQ=12)

                                                      tcTAAATTTCAATTAAAAATGAGACCTACCGAAGGGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCt             <  1:11674018/67‑1 (MQ=17)

                                                      tatAAATTTCTTTTAAAACTGTGAAGTCCCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCt             <  1:42131754/66‑1 (MQ=21)

                                                       cTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCt             <  1:38245403/67‑1 (MQ=255)

                                                        tAAATTTCTTGTAAAACTGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCt             <  1:5912146/66‑1 (MQ=255)

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

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

                                                        tAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCGCAGCGGGGAGAGCACCTccc              <  1:16623625/65‑1 (MQ=255)

                                                        tAAATATCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTccc              <  1:39445532/65‑1 (MQ=255)

                                                         aaaTTTCTTTTAAATATGTGAACTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTa           <  1:18135716/67‑1 (MQ=21)

                                                          aaTTTCTTGAAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACa         <  1:38594727/68‑1 (MQ=255)

                                                            tttCTAGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAag       <  1:17665133/68‑1 (MQ=255)

                                                             ttCTTGTAAAAATGTGAAGTACCGACGTGGGGGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAagg      >  1:29537253/1‑68 (MQ=255)

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

                                                              tCTTGTAAAAATGTCACCTACCCCAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAagga     <  1:8741595/68‑1 (MQ=17)

TACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGG  >  CP009273/2514443‑2514575

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