BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	692,043	T→A	6.2%	noncoding (30/75 nt)	glnV ←	tRNA‑Gln

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	692,043	0	T	A	6.2%	1526.0 / 30.9	488	noncoding (30/75 nt)	glnV	tRNA‑Gln
Reads supporting (aligned to +/- strand): ref base T (304/150); new base A (8/26); total (312/176)
Fisher's exact test for biased strand distribution p-value = 8.90e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 500 total aligned reads displayed.

GGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACA  >  CP009273/691978‑692110

ggTAAGCGTCTTGGATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATcc                                                                   >  1:20327249/1‑68 (MQ=255)

ggTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATcc                                                                   >  1:22588760/1‑68 (MQ=255)

ggTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATcc                                                                   >  1:24460822/1‑68 (MQ=255)

ggTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATcc                                                                   >  1:16722961/1‑68 (MQ=255)

ggTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATcc                                                                   >  1:16463535/1‑68 (MQ=255)

ggTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATcc                                                                   >  1:16186337/1‑68 (MQ=255)

ggTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATcc                                                                   >  1:15624376/1‑68 (MQ=255)

ggTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATcc                                                                   >  1:2696719/1‑68 (MQ=255)

ggTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATcc                                                                   >  1:28988294/1‑68 (MQ=255)

ggTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATcc                                                                   >  1:13329827/1‑68 (MQ=255)

ggTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATcc                                                                   >  1:12450481/1‑68 (MQ=255)

ggTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATcc                                                                   >  1:5683888/1‑68 (MQ=255)

  tAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCgg                                                                 >  1:23323236/1‑68 (MQ=255)

  tAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCgg                                                                 >  1:18395314/1‑68 (MQ=255)

  tAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCgg                                                                 >  1:24766568/1‑68 (MQ=255)

      gggcTTTAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGcc                                                             <  1:24866774/65‑1 (MQ=255)

      gagcttttttaatgTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGcc                                                             <  1:3204476‑M2/54‑1 (MQ=255)

      gagcttttttaatgTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGc                                                              <  1:20856550‑M2/53‑1 (MQ=255)

      gagcttttttaatgTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGc                                                              <  1:14013843‑M2/53‑1 (MQ=255)

      gagcttttttaatgTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGc                                                              <  1:22455602‑M2/53‑1 (MQ=255)

       gTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCt                                                            >  1:20161830/1‑68 (MQ=255)

       gTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCt                                                            >  1:29847686/1‑68 (MQ=255)

       gTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAAAAACCGGTGGCt                                                            >  1:15311929/1‑68 (MQ=17)

        tCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCtt                                                           >  1:22225917/1‑68 (MQ=255)

        gcttttttaatgTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCtt                                                           <  1:16963577‑M2/56‑1 (MQ=255)

        gcttttttaatgTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCt                                                            <  1:7937042‑M2/55‑1 (MQ=255)

           tttttaatgTGGCTTGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTAcc                                                        <  1:21232334‑M2/59‑1 (MQ=255)

             tttaatgTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGc                                                      <  1:17049008‑M2/61‑1 (MQ=255)

             tttaatgTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAAAAACCGGTGCCTTACCGc                                                      <  1:1802251‑M2/61‑1 (MQ=255)

             aaTAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTAcc                                                        >  1:3744606/1‑66 (MQ=255)

             aaTAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGc                                                      <  1:27381494/68‑1 (MQ=255)

              ttaatgTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCt                                                     <  1:8352663‑M2/62‑1 (MQ=255)

              ttAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCt                                                     <  1:4851573/67‑1 (MQ=39)

              aTAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCt                                                     >  1:19699467/1‑68 (MQ=40)

                aatgTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTg                                                   <  1:15539551‑M2/64‑1 (MQ=255)

                 atgTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTgg                                                  <  1:5282801‑M2/65‑1 (MQ=255)

                  tttGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGTATCATAAACCGGTGCCTTACCGCTTGGc                                                 <  1:27687580/67‑1 (MQ=11)

                  tttGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGTATCACAAACCGGTGCCTTACCGCTTGGc                                                 <  1:24235433/67‑1 (MQ=11)

                  tttGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGc                                                 <  1:24045608/67‑1 (MQ=11)

                  tgTGGCTGGGGTACGTGGATTCGAACCTCGGAATGCCGGAATCATAAACCGGTGCCTTACCGCTTGGc                                                 <  1:23707151‑M2/66‑1 (MQ=255)

                  tgTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGc                                                 <  1:10419184‑M2/66‑1 (MQ=255)

                  tgTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCATAAACCGGTGCCTTACCGCTTGGc                                                 <  1:9789118‑M2/66‑1 (MQ=255)

                  tgTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCgctttc                                                  >  1:1145015‑M2/3‑65 (MQ=255)

                  tgTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGc                                                 >  1:6891171‑M2/3‑68 (MQ=255)

                  tgTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGc                                                 <  1:12122198‑M2/66‑1 (MQ=255)

                  tgTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGc                                                 >  1:3336275‑M2/3‑68 (MQ=255)

                  tgTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGc                                                 <  1:29425573‑M2/66‑1 (MQ=255)

                  tgTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGc                                                 >  1:13552414‑M2/3‑68 (MQ=255)

                  tgTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGc                                                 >  1:6299463‑M2/3‑68 (MQ=255)

                  tgTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGc                                                 >  1:6517536‑M2/3‑68 (MQ=255)

                  tgTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGc                                                 >  1:13727310‑M2/3‑68 (MQ=255)

                  tgTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGc                                                 <  1:25464564‑M2/66‑1 (MQ=255)

                  tgTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGc                                                 >  1:21926625‑M2/3‑68 (MQ=255)

                   ttGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGATTCAGAATCCGGTGCCTTACCGCTTGGCg                                                >  1:4541832/1‑68 (MQ=14)

                   ttGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCg                                                >  1:19386379/1‑68 (MQ=31)

                   ttGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCg                                                >  1:25999132/1‑68 (MQ=31)

                   ttGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCg                                                >  1:26496752/1‑68 (MQ=31)

                   gTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCg                                                <  1:29349735‑M2/67‑1 (MQ=255)

                              aCGAGGATTCGAACCTCGGAATGCCGGTATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCcaaca                                     <  1:14527994/68‑1 (MQ=14)

                                 aGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCcaacagtg                                  >  1:29242448/1‑65 (MQ=14)

                                 aGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCcaacaggg                                  >  1:8635066/1‑65 (MQ=1)

                                 aGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCGACAAAt                                  >  1:6662726/1‑68 (MQ=255)

                                 aGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAAGaaag                                  >  1:9514621/1‑67 (MQ=9)

                                 aGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACatat                                  >  1:2882399/1‑68 (MQ=14)

                                 aGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACaaaa                                  >  1:17813781/1‑67 (MQ=25)

                                 aGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAGAt                                  >  1:10936620/1‑68 (MQ=14)

                                 aGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAAt                                  >  1:145967/1‑68 (MQ=37)

                                 aGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAAt                                  >  1:1482964/1‑68 (MQ=37)

                                 aGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAAt                                  >  1:8135573/1‑68 (MQ=37)

                                 aGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAAt                                  >  1:26331125/1‑68 (MQ=37)

                                  ggATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAAtt                                 >  1:11764155/1‑68 (MQ=39)

                                   gATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTg                                >  1:15989783/1‑68 (MQ=40)

                                     ttCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCATCAAATTTGt                              >  1:18674150/1‑68 (MQ=255)

                                     ttCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTgg                               <  1:25762796/67‑1 (MQ=255)

                                     ttCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGt                              <  1:5857078/68‑1 (MQ=255)

                                     ttCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGt                              <  1:4983774/68‑1 (MQ=255)

                                          aCCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACATTTTGGTTTTGa                         >  1:5490623/1‑68 (MQ=255)

                                          aCCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACATATTGGTTGTGa                         >  1:24655087/1‑68 (MQ=255)

                                          aCCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGa                         >  1:17273444/1‑68 (MQ=255)

                                          aCCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGa                         >  1:10697884/1‑68 (MQ=255)

                                          aCCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACACCAACAAAGTGGtattta                         >  1:2503241/1‑66 (MQ=14)

                                           ccTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGaa                        >  1:26112295/1‑68 (MQ=255)

                                              cGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAAttt                     <  1:10040903/68‑1 (MQ=255)

                                              cGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAAttt                     <  1:12972678/68‑1 (MQ=255)

                                                 aaTGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGc                   <  1:4168052/67‑1 (MQ=255)

                                                 aaTGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGc                   <  1:17524381/67‑1 (MQ=255)

                                                  aTGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCg                 >  1:19176622/1‑68 (MQ=255)

                                                   tGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGa                <  1:7748355/68‑1 (MQ=255)

                                                     ccGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAAc              >  1:28145325/1‑68 (MQ=255)

                                                       ggAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACat            <  1:8977156/68‑1 (MQ=255)

                                                         aaTCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACata           <  1:28324225/67‑1 (MQ=255)

                                                         aaTCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACata           <  1:23080178/67‑1 (MQ=255)

                                                          aTCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACatat          <  1:12706559/67‑1 (MQ=255)

                                                          aTCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACatat          <  1:10292553/67‑1 (MQ=255)

                                                           tCAGAATCCGGTGCCTT‑CCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCg       >  1:27921102/1‑68 (MQ=255)

                                                            cAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCg       <  1:9291769/68‑1 (MQ=255)

                                                                aTCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATAc   >  1:20493374/1‑68 (MQ=255)

                                                                aTCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATAc   >  1:21100483/1‑68 (MQ=255)

                                                                aTCCGGTGCCTTACCGCTTGGCGATACACCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATAc   >  1:7334549/1‑68 (MQ=255)

                                                                 tCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACa  >  1:26806630/1‑68 (MQ=255)

GGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACA  >  CP009273/691978‑692110

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