BRESEQ :: Evidence

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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	692,035	0	A	T	7.1%	1420.3 / 41.6	461	noncoding (38/75 nt)	glnV	tRNA‑Gln
Reads supporting (aligned to +/- strand): ref base A (283/141); new base T (6/31); total (289/172)
Fisher's exact test for biased strand distribution p-value = 1.95e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 479 total aligned reads displayed.

GCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATA  >  CP009273/691968‑692101

 cACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGaa                                                                   >  1:31897897/1‑68 (MQ=255)

  aCTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAAt                                                                  <  1:17530599/68‑1 (MQ=255)

   cTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATc                                                                 >  1:15785741/1‑68 (MQ=255)

     tACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAg                                                               >  1:15919173/1‑68 (MQ=255)

        aaGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGaa                                                             >  1:26538802/1‑67 (MQ=255)

        aaGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGaa                                                             >  1:16393003/1‑67 (MQ=255)

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

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

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

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

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

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

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

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

          ggTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAAAcc                                                          >  1:31696514/1‑68 (MQ=255)

          ggTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCACAATcc                                                          >  1:24973911/1‑68 (MQ=255)

            tAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGTATCAAAATCCgg                                                        >  1:13569460/1‑68 (MQ=255)

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

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

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

             aaGCTTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGt                                                       <  1:22766122/68‑1 (MQ=14)

             aaGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGt                                                       <  1:25736194/68‑1 (MQ=255)

              aGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCACGGAATGCCGGTATCAAAAACCTGtc                                                      >  1:22558721/1‑67 (MQ=17)

              aGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCACGGAATGCCGGAATCAAAAACCTGtc                                                      >  1:28419190/1‑67 (MQ=18)

                cGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCaaaaa                                                            >  1:30721618/1‑59 (MQ=255)

                cGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTg                                                      >  1:4233223/1‑66 (MQ=255)

                cGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTg                                                      >  1:10827051/1‑66 (MQ=255)

                cGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTg                                                      >  1:20905599/1‑66 (MQ=255)

                cGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTg                                                      >  1:1112218/1‑66 (MQ=255)

                cGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTg                                                      >  1:11560883/1‑66 (MQ=255)

                cGTCTTGAATAAATTGGCTGGGGTACAAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTg                                                      >  1:27443896/1‑66 (MQ=255)

                 ttattGAATAAATTGGCTGGGGTACGAGGATTCGAACCTGGGAATGCCGGTATCAAAAACCGGTGcc                                                    <  1:27569790/64‑1 (MQ=14)

                 gTCTTGAATAACTTGGCTGGGGTACGAGGATTCGAACCTGGGAATGCCGGTATCAAAAACCGGTGc                                                     <  1:25172247/66‑1 (MQ=14)

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

                 gTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGc                                                     <  1:12532765/66‑1 (MQ=11)

                 gTCTTGAATAAATTGGCTGGGGTAAGAGGATTCGAACCTCGGAATGCCGGTATCaaaaa                                                            >  1:21713802/1‑58 (MQ=18)

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

                    ttGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTAc                                                <  1:1915041/68‑1 (MQ=255)

                     tttttaatgTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTAcc                                               <  1:24328662‑M2/59‑1 (MQ=255)

                     ttaaTAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTAcc                                               <  1:10156623/66‑1 (MQ=255)

                      gAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTAcc                                               >  1:23674838/1‑67 (MQ=255)

                       tttaatgTGGGTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGc                                             <  1:23319382‑M2/61‑1 (MQ=255)

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

                       tttaaagTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGc                                             <  1:28772485‑M2/61‑1 (MQ=255)

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

                        ttgatgTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCt                                            <  1:17374331‑M2/62‑1 (MQ=255)

                        ttaattTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCt                                            <  1:28925582‑M2/62‑1 (MQ=255)

                        ttaatgTGGCTGGTGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGc                                             <  1:21864458‑M2/61‑1 (MQ=255)

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

                        ttaatgTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCt                                            >  1:8001349‑M2/7‑68 (MQ=255)

                        ttaatgTGGCTGGGGTAAGAGGATTCGAACCTCGGAATGCAGGAATCAGAATAAGGTGCATTACCGCt                                            >  1:22151039‑M2/7‑68 (MQ=255)

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

                        aTAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAAACCGGTGCCTTAcc                                               >  1:1679623/1‑65 (MQ=255)

                        aTAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGAAGGAATCAGAATCCGGTGCCTTACCGCt                                            >  1:2683736/1‑68 (MQ=255)

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

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

                           atttGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCACAAACCGGTGCCTTACCGCTTgg                                         <  1:27835500/66‑1 (MQ=11)

                           atgTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTgg                                         <  1:25547885‑M2/65‑1 (MQ=255)

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

                           atgTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCCGTGCCTTACCGCTTgg                                         <  1:16021277‑M2/65‑1 (MQ=255)

                            tttGGCTGGGGTACGAGGATTCGAACCACGGAATGCCGGTATCATAAACCGGTGCCTTACCGCTTGGc                                        <  1:31097067/67‑1 (MQ=11)

                            tgTGGCTGGGGTACGTGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGc                                        >  1:843956‑M2/3‑68 (MQ=255)

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

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

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

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

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

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

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

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

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

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

                            tgTGGCTGGGGTACGAGGATTCGAACCACGGAATGCCGGAATCACAAACCGGTGCCTTACCGCTTGGc                                        <  1:21538601‑M2/66‑1 (MQ=255)

                            tgTGGCTGGGGTACGAGGATTCGAAACTCGGAATGCCGGAATCAGAATCCGGTGCCTTAACGCTTCGc                                        >  1:13273172‑M2/3‑68 (MQ=255)

                            agTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGc                                        <  1:14078194‑M2/66‑1 (MQ=255)

                            aTTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGc                                        >  1:7486968/1‑68 (MQ=33)

                            aTTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGc                                        >  1:1531929/1‑68 (MQ=33)

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

                                          gAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCcaaccgg                          >  1:28071058/1‑65 (MQ=1)

                                          gAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACaaa                          <  1:2717564/68‑1 (MQ=35)

                                          gAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACaaa                          <  1:535620/68‑1 (MQ=35)

                                           aGGATTCGAACCTCGGAATGCCGGAATCATAATCCGGTGCATAACCGCTTGGCGATAACCccacaggg                         >  1:29562472/1‑65 (MQ=1)

                                           aGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCcaacaagt                         >  1:12855059/1‑66 (MQ=18)

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

                                           aGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACAGCTTGGCGATACCCCAACaaaa                         >  1:5534428/1‑67 (MQ=255)

                                           aGGATTCGAACCTCGGAATCCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACaaaa                         >  1:1793687/1‑67 (MQ=255)

                                             gATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAGTTg                       >  1:30906672/1‑68 (MQ=25)

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

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

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

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

                                                   aaCCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGGTTTg                 >  1:1464057/1‑68 (MQ=255)

                                                    aCCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGGTTTGa                >  1:26043939/1‑68 (MQ=255)

                                                             tGCCGGAATCAGAATCCGGGGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGa       <  1:30178111/68‑1 (MQ=255)

                                                              gCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGaa      <  1:26762259/68‑1 (MQ=255)

                                                              gCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGa       >  1:10705674/1‑67 (MQ=255)

                                                                cGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACa    <  1:7317033/68‑1 (MQ=255)

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

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

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

GCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATA  >  CP009273/691968‑692101

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