BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	2,514,509	G→T	5.3%	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	5.3%	1201.1 / 3.6	392	intergenic (+24/‑23)	valU/valY	tRNA‑Val/tRNA‑Val
Reads supporting (aligned to +/- strand): ref base G (139/220); new base T (7/25); total (147/245)
Fisher's exact test for biased strand distribution p-value = 8.45e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 395 total aligned reads displayed.

TACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGG  >  CP009273/2514443‑2514571

   aaGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCTACCACTTTCTCGCCAGCTAAATTTGTTGTaaa                                                            >  1:10451726/1‑68 (MQ=255)

   aaGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTGGCCTGATAAATTTGTTGTaaa                                                            >  1:16342209/1‑68 (MQ=17)

   aaGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCGTTCTAAATGTCTTGTaaa                                                            >  1:28108028/1‑68 (MQ=17)

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

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

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

   aaGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTACTTGTaaa                                                            >  1:21985241/1‑68 (MQ=255)

   aaGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCGCTCCAGCGAAAATTGTTGTaaa                                                            >  1:23181418/1‑68 (MQ=17)

   aaGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCGCGCCATCTAAATGTATTGTaaa                                                            >  1:21350339/1‑68 (MQ=17)

   aaGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTAGTCGCCAGATAAATTTCTT‑Taaaa                                                           >  1:1498239/1‑68 (MQ=17)

   aaGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTACTCGCCAGGTAAATTTGTTGTaaa                                                            >  1:24929698/1‑68 (MQ=18)

   aaGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTACTCGCCAGATAAATTTTTTGTaaa                                                            >  1:25754807/1‑68 (MQ=18)

   aaGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTACGCTCCAGATAAATTTCTTGTaaa                                                            >  1:25286468/1‑68 (MQ=17)

    aggaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTaaaa                                                           <  1:28588479/68‑1 (MQ=255)

     ggaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTGGCGCGTGCTAAATTTGTTGTaaaaa                                                          >  1:620061/1‑68 (MQ=17)

     ggaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTGGTaaaaa                                                          >  1:936242/1‑68 (MQ=255)

       aggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTGTTTTAAAAAtt                                                        >  1:19610478/1‑67 (MQ=21)

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

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

            gTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGaa                                                    <  1:30170058/67‑1 (MQ=255)

             tCGGCGGTTCGATCCCGTCATCACCCACCACTTTGTCGCCAGCTAAATTTAGTGAAATAAGGTGAAGt                                                  >  1:5413269/1‑68 (MQ=18)

             tCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGt                                                  >  1:5820998/1‑68 (MQ=255)

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

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

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

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

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

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

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

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

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

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

                     tCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTATTGTAAAAATGTGAAGTTACGAAGt                                          >  1:18481170/1‑68 (MQ=37)

                         tCCCGTCATCACCCACCACTTTCTCGCCAGCTAACTTTATTGTAAAAATGTGAAGTAAAGAAGTgtgg                                      >  1:23032669/1‑65 (MQ=37)

                         tCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTggg                                       <  1:7620517/67‑1 (MQ=255)

                         tCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGt                                      >  1:11499717/1‑68 (MQ=255)

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

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

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

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

                         tCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTAACGAAGTGGGt                                      >  1:19010126/1‑68 (MQ=255)

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

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

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

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

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

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

                                atcaCCAACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGc                               <  1:10609953/68‑1 (MQ=255)

                                   acccaccACTTTCTCGCCAGCTAAATTTCTTGTAAAAATTTGAAGTACCGAAGTGTGTGATAAGCtaa                            >  1:31240401/1‑66 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

                                             ttCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGagag                  <  1:1568682/68‑1 (MQ=255)

                                              tctcGCCAGCTAAATCGCTCGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCGCAGCTGGGAGAGc                 <  1:3602185/68‑1 (MQ=37)

                                                 cGCCAGCTAACTTTCTTGAAAAAACGAGAAGTACCCCCGTGGGTGATTAGCTCAGCTGGGAGAGCAcc              <  1:17870138/68‑1 (MQ=18)

                                                 cGCCAGATAAATTTCTTTTAAAAATGTGACGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCAcc              <  1:4813900/68‑1 (MQ=37)

                                                  gCCATCGAAATTTCTTTTAAATCTGTGAAGTACCCAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCt             <  1:2157501/68‑1 (MQ=21)

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

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

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

                                                  gCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCt             <  1:18618829/68‑1 (MQ=255)

                                                  gCCAGCAAAATTTCTTCTAAAAATGTGAACTACCGAAGTGGGTGATTATATCAGATGGGAGACCAca              >  1:26214515/1‑66 (MQ=255)

                                                   ccAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCt             >  1:10099606/1‑67 (MQ=255)

                                                   ccAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCt             >  1:18902871/1‑67 (MQ=255)

                                                    catcTAACTATCTTGTAAAAGTGAGAAGTCCCGCAGGGGGTGATTAGCGCAGCTGGGAGAGCACCTcc           <  1:19490281/65‑1 (MQ=12)

                                                    caccTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTcc           <  1:18025648/65‑1 (MQ=255)

                                                    cAGCTAAATTTCTTGTAAAACAGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTcc           <  1:19240169/68‑1 (MQ=255)

                                                    cAGCTAAATTTCTTGTAAAAATGTGAAGTCCCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTcc           <  1:6951415/68‑1 (MQ=255)

                                                    cAGCTAAATTTCTTGTAAAAATGTGAAGTCCCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTcc           <  1:367768/68‑1 (MQ=255)

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

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

                                                    cAGCTAAATCTCTTTTAAAACTGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTcc           <  1:9647551/68‑1 (MQ=25)

                                                    cAGCTAAATATCTTGTAAAAATGTGAAGTACCGCAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTc            <  1:12204203/67‑1 (MQ=255)

                                                    cAGCGAACTTTCTTTTAAAATTGTGAAGTACCGACGTGGGTGATTAGCTCAGCTGGGAGAGCACCTcc           <  1:15599015/68‑1 (MQ=21)

                                                    aatcTAAATATCATTAAAAAATGTGAACTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTcc           <  1:30509037/65‑1 (MQ=18)

                                                     atcTAACTTTCTTGTAAAACTGAGAAGTCACGAAGGGGGTGATTAGCTCAGCTGGGAGAGCACCTccc          <  1:6642779/66‑1 (MQ=17)

                                                      tcTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCt         <  1:24617941/67‑1 (MQ=255)

                                                      tcTAAATTTCTTGTAAAAATGTGAACAACCGCAGGGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCt         <  1:9167257/67‑1 (MQ=21)

                                                      gCTAAATTTCTTTTAAAATTGTGTAGTCCCTTAGGGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCt         <  1:17014903/68‑1 (MQ=17)

                                                        tAAATTTCTTTTAAAAAGGTGAACTACCCACGGGGGTGATTAGCTCAGCTGGGAGAGCACCTccc          <  1:21810688/65‑1 (MQ=17)

                                                        tAAATTTCTTGTAAAAATGTGAAGTACCGTAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTccc          <  1:4327738/65‑1 (MQ=255)

                                                        tAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTa       <  1:10396750/68‑1 (MQ=255)

                                                        tAAATTTCTTGTAAAAAGGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTccc          <  1:21709035/65‑1 (MQ=255)

                                                          aaTTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACa     <  1:17464100/68‑1 (MQ=255)

                                                          aaTTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACa     <  1:377214/68‑1 (MQ=255)

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

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

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

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

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

TACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTTCTCGCCAGCTAAATTTCTTGTAAAAATGTGAAGTACCGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGG  >  CP009273/2514443‑2514571

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