BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,382,371	G→T	10.0%	intergenic (+7/‑29)	glyX → / → glyY	tRNA‑Gly/tRNA‑Gly

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	4,382,371	0	G	T	10.0%	751.4 / 7.6	296	intergenic (+7/‑29)	glyX/glyY	tRNA‑Gly/tRNA‑Gly
Reads supporting (aligned to +/- strand): ref base G (74/175); new base T (13/23); total (87/209)
Fisher's exact test for biased strand distribution p-value = 4.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 297 total aligned reads displayed.

TAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCA  >  CP009273/4382308‑4382436

tAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTTTG‑AAAg                                                               >  1:6760293/1‑68 (MQ=14)

tAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAg                                                               >  1:4781338/1‑68 (MQ=31)

tAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAg                                                               >  1:18652538/1‑68 (MQ=31)

 agagCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTTTG‑AAAGt                                                              >  1:25850588/1‑68 (MQ=18)

 agagCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGt                                                              <  1:15532565/68‑1 (MQ=33)

 agagCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGt                                                              <  1:2882043/68‑1 (MQ=33)

 agagCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGt                                                              <  1:22875959/68‑1 (MQ=33)

 agagCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGt                                                              <  1:23570449/68‑1 (MQ=33)

 agagCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGt                                                              <  1:2431722/68‑1 (MQ=33)

  gagCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTg                                                             >  1:14781821/1‑68 (MQ=35)

    gCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGc                                                            <  1:15158084/67‑1 (MQ=37)

    gCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCt                                                           <  1:3476053/68‑1 (MQ=39)

    gCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCt                                                           >  1:21220213/1‑68 (MQ=39)

     cACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTTTG‑CAAGt                                                              >  1:26090764/1‑64 (MQ=255)

     cACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTT‑AAAGTGATg                                                          >  1:21016225/1‑68 (MQ=255)

     cACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGt                                                              >  1:21965917/1‑64 (MQ=33)

      ccGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGt                                                         >  1:10349469/2‑68 (MQ=255)

      aCGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCt                                                           <  1:2694688/66‑1 (MQ=39)

      aCGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCt                                                           <  1:24458340/66‑1 (MQ=39)

      aCGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCt                                                           <  1:23619923/66‑1 (MQ=39)

      aCGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGt                                                         <  1:28296457/68‑1 (MQ=40)

      aCGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGt                                                         <  1:2154929/68‑1 (MQ=40)

      aCGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGt                                                         <  1:27467326/68‑1 (MQ=40)

      aCGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGt                                                         <  1:11253961/68‑1 (MQ=40)

      aCGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGt                                                         >  1:27727308/1‑68 (MQ=40)

      aCGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGt                                                         <  1:9310526/68‑1 (MQ=40)

      aCGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGt                                                         <  1:14033302/68‑1 (MQ=40)

      aCGACCTTGCCAAGGGCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGGGCTGt                                                         <  1:28705199/68‑1 (MQ=255)

         aCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑ATAGTGATGTaat                                                      >  1:27649666/1‑67 (MQ=14)

         aCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAg                                                      >  1:10882151/1‑68 (MQ=255)

         aCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGATTTAAg                                                      >  1:10699198/1‑68 (MQ=255)

          ccTTGCCACGGTCGGGGTCGCGAGTTCGAGGCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTCAgg                                                     <  1:1732191/68‑1 (MQ=14)

          ccTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAgg                                                     <  1:17971214/68‑1 (MQ=255)

          ccTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAgg                                                     <  1:11383300/68‑1 (MQ=255)

          ccTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAgg                                                     <  1:25461480/68‑1 (MQ=255)

          ccTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAgg                                                     <  1:2807337/68‑1 (MQ=255)

          ccTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAgg                                                     <  1:15724002/68‑1 (MQ=255)

           cTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGc                                                    >  1:16948482/1‑68 (MQ=255)

           cTTGCCAAGGTCGGGGTCGCGAGTTCGAGTATCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGc                                                    <  1:22730173/68‑1 (MQ=255)

               ccAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTT‑AAAGTGCTGTAAGGGACAg                                                >  1:17489604/1‑68 (MQ=255)

                cAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTT‑AAAGTGATTTAAGGAACAg                                                >  1:10484374/1‑67 (MQ=14)

                cAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAg                                                >  1:936559/1‑67 (MQ=255)

                cAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAg                                                >  1:1188427/1‑67 (MQ=255)

                cAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGa                                               <  1:9780717/68‑1 (MQ=255)

                cAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGa                                               <  1:16271047/68‑1 (MQ=255)

                cAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGa                                               <  1:12734525/68‑1 (MQ=255)

                cAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGa                                               <  1:19038112/68‑1 (MQ=255)

                cAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGAACAg                                                >  1:18568613/1‑67 (MQ=255)

                cAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGATTTAAGGGACAg                                                >  1:6804032/1‑67 (MQ=14)

                cAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATATT‑AAAGAGCTATAATGAACAg                                                >  1:13533854/1‑67 (MQ=11)

                cAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTACCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGa                                               <  1:11523595/68‑1 (MQ=255)

                  aGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGacc                                             <  1:2642713/68‑1 (MQ=255)

                  aGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGacc                                             <  1:26610754/68‑1 (MQ=255)

                  aGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGacc                                             <  1:442166/68‑1 (MQ=255)

                    gTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATATT‑AAAGTTCTGTAATgaccagacca                                            >  1:3556002/1‑67 (MQ=14)

                     tCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGaccac                                           <  1:14357329/67‑1 (MQ=255)

                     tCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGaccac                                           <  1:11709765/67‑1 (MQ=255)

                     tCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGaccac                                           <  1:20057726/67‑1 (MQ=255)

                         ggTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCaa                                       <  1:7849594/67‑1 (MQ=255)

                          gTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAg                                      <  1:28421652/67‑1 (MQ=255)

                          gTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAg                                      <  1:7714336/67‑1 (MQ=255)

                          gTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAg                                      <  1:16079976/67‑1 (MQ=255)

                          gTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAg                                      <  1:12839773/67‑1 (MQ=255)

                          gTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAg                                      <  1:6578103/67‑1 (MQ=255)

                          gTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAg                                      <  1:9740594/67‑1 (MQ=255)

                          gTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAg                                      <  1:25002309/67‑1 (MQ=255)

                          gTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAg                                      <  1:189060/67‑1 (MQ=255)

                          gTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAg                                      <  1:8440037/67‑1 (MQ=255)

                          gTCGCGAGTTCGAGGCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAg                                      <  1:23803067/67‑1 (MQ=255)

                          gTCGCGAGTGCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAg                                      <  1:29070485/67‑1 (MQ=255)

                          gTCGCGAGTACGAGTCGCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAg                                      <  1:2458094/67‑1 (MQ=255)

                             gcgAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTT‑AAAGTGCTGTAAGGCACAGACCACCCAAGCggg                                  <  1:599685/68‑1 (MQ=255)

                               gatttCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGa                                 <  1:19431503/64‑1 (MQ=255)

                               gAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGa                                 <  1:24111642/67‑1 (MQ=255)

                               gAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGa                                 <  1:4022785/67‑1 (MQ=255)

                               gAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGa                                 <  1:8231315/67‑1 (MQ=255)

                               gAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGa                                 <  1:25154344/67‑1 (MQ=255)

                                 gTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATa                              >  1:5420394/1‑68 (MQ=255)

                                   tCGCGTCTGGTTTCCCGCTCCAAAATTTT‑AAAGTGCTGTAACGCACAGACCACCCAAGCGGGAATAGc                            <  1:8702874/68‑1 (MQ=14)

                                     gAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTc                          <  1:19662395/68‑1 (MQ=255)

                                       gTCTCGTTTCCCGCTCCAAAATTTT‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCa                         <  1:18239016/67‑1 (MQ=255)

                                         ctcGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAg                        <  1:26746962/66‑1 (MQ=255)

                                         ctcGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAg                        <  1:3394667/66‑1 (MQ=255)

                                         ctcGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAg                        <  1:13922805/66‑1 (MQ=255)

                                         ctcGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGt                       <  1:4581237/67‑1 (MQ=255)

                                         cgcGTTTCCCGCTCCAAAATATC‑AAAGTGCTGTAAGGAACAGACCACCCAAGCGGGAATAGCTCAGt                       <  1:23896204/65‑1 (MQ=17)

                                         cgcGTATCCCCCTCCAAAATATT‑AAAGTGCTCAAAGTCACAGACCACCCACGCGGGAATAGCTCAg                        <  1:21653723/64‑1 (MQ=9)

                                          tcGTTTCCCGGGCCAAAATTTT‑ATATTGCTGTAAGGCACATAACACCTAAGCGGGAATAGCTCAGTTg                     <  1:10063476/68‑1 (MQ=17)

                                               tCCCGGTCCAAAATTTT‑AAAGTGCTGAAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTaga                <  1:15043971/68‑1 (MQ=14)

                                                cccGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTaga                >  1:15566495/1‑67 (MQ=255)

                                                cccGCTCCAAAATTTG‑AAAGTGCTGAAAGGCAAAGACCACCCAAGAGGGAAAAGCTCAGTTGGTata                >  1:16372540/1‑65 (MQ=14)

                                                      ccAAAATTTT‑AAAGTGCTGAACGCCACACACCACCCAAGCGGGAATAGCTCAGTTCGTACAGCACGAc         <  1:20489592/68‑1 (MQ=14)

                                                      ccAAAATTTA‑AAATTGCTGAACGTCACATAACACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGAc         <  1:13446882/68‑1 (MQ=14)

                                                      caaaaaTTTG‑AAAGTGCTGTACGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGAc         <  1:22935246/66‑1 (MQ=255)

                                                      caaaaaTTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGAc         <  1:24760191/66‑1 (MQ=255)

                                                      caaaaaTTTC‑AAAGTTCTGTAATGCTCAGACCCCCCACGATGGAATAGATCAGTTGGTAGAGCACGAc         <  1:25597248/66‑1 (MQ=9)

                                                      acAAAATTTC‑AAAGTCCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGAc         <  1:29195809/67‑1 (MQ=17)

                                                       acaaaTTTGTAAAGTGCTCTATGTCAAATACCTCCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGAc         <  1:8357163/66‑1 (MQ=2)

                                                          attttt‑AAAGTGCTGTAATGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTg     <  1:12447201/66‑1 (MQ=11)

                                                           aTTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGc    <  1:25383502/68‑1 (MQ=255)

TAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCA  >  CP009273/4382308‑4382436

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