BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	692,035	A→T	5.6%	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	5.6%	2803.3 / 42.5	869	noncoding (38/75 nt)	glnV	tRNA‑Gln
Reads supporting (aligned to +/- strand): ref base A (574/236); new base T (15/44); total (589/280)
Fisher's exact test for biased strand distribution p-value = 5.44e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 905 total aligned reads displayed.

GCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATAT  >  CP009273/691968‑692102

gCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGa                                                                     >  1:36834799/1‑68 (MQ=255)

       cAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCACGGAATGCCGGTATCAgca                                                              >  1:7644771/1‑66 (MQ=21)

       cAAGGAACGCGTCTTTATTAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGaa                                                              <  1:4933170/68‑1 (MQ=14)

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

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

         aGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAAt                                                             >  1:34132705/1‑67 (MQ=255)

         aGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATc                                                            <  1:41175568/68‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

          ggTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCACGGAATGCCGGAATCaaaaac                                                            >  1:42890795/1‑65 (MQ=21)

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

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

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

              aGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGt                                                        >  1:45952521/1‑67 (MQ=255)

              aGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGt                                                        >  1:12413535/1‑67 (MQ=255)

                gagcttttttaatgTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTg                                                       >  1:26098513‑M2/15‑66 (MQ=255)

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

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

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

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

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

                ccgcTTTAATAATTTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGcc                                                     <  1:23071382/65‑1 (MQ=11)

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

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

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

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

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

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

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

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

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

                cGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCACGGAATGCCGGAATCAAAAACCTGtc                                                       >  1:37205032/1‑65 (MQ=17)

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

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

                 cgcTTGAATAAATTGGCTGGGGGACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGcc                                                     <  1:35760521/65‑1 (MQ=255)

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

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

                  tCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAAACCGGTGCCtt                                                   >  1:5507291/1‑68 (MQ=255)

                    ttttttaatgTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCtt                                                   <  1:58088261‑M2/56‑1 (MQ=255)

                     tttttaatgTGGCGGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTAcc                                                <  1:1692653‑M2/59‑1 (MQ=255)

                     tttttaaagTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTAcc                                                <  1:7215581‑M2/59‑1 (MQ=255)

                      ttttaatgTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAAt                                                                   <  1:3050452‑M2/40‑1 (MQ=255)

                          aatgTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTg                                           >  1:28055160‑M2/5‑68 (MQ=255)

                           aaTTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTg                                           >  1:8209154/1‑67 (MQ=35)

                            tttGGCTGGGGTACGAGGATTCGAACCTGGGAATGCCGGTATCATAATCCGGTGCCTTACCGCTTGGc                                         <  1:24219922/67‑1 (MQ=11)

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

                            tttGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGc                                         <  1:17977237/67‑1 (MQ=14)

                            tgTGGCTGGGGTCCGAGGATTCGAACCTCGGAATGCCGGAATCATAAACCGGTGCCTTACCGCTTGGc                                         <  1:56196886‑M2/66‑1 (MQ=255)

                            tgTGGCTGGGGTACGAGGATTCGTTCCTCGGAATGGCGGAATCAGAATCCGGTGCCTTACGGCTTGGc                                         <  1:49151455‑M2/66‑1 (MQ=255)

                            tgTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGATTAAGAATACGGTGCCTTACCGCTagga                                         >  1:52149688‑M2/3‑67 (MQ=255)

                            tgTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGCATCCGGTGCCTTACCGCttgtc                                         >  1:2489852‑M2/3‑66 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                             ttGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCaaaaa                                                             >  1:4646379/1‑46 (MQ=255)

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

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

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

                             gTGGCTGGGGTACGAGGAGTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCg                                        <  1:4768213‑M2/67‑1 (MQ=255)

                                          gAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCcaacagg                           >  1:48109780/1‑66 (MQ=14)

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

                                           aGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCcaacaggt                          >  1:51168825/1‑65 (MQ=14)

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

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

                                           aGGATTCGAACCTCGGAATGCCGGAATCAGAAAACTGTGCATTACCGCTGGGCGATACCCCAAGAAAt                          >  1:6300020/1‑68 (MQ=11)

                                           aGGATTCGAAATTCGGAATGCCGGAATCAGAATCCGCTGCCTTACCGCTTGGCGATACCCCAACaaaa                          >  1:29421112/1‑67 (MQ=11)

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

                                               ttCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTgga                      >  1:23691861/1‑67 (MQ=255)

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

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

                                                 cGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGttt                    >  1:4853655/1‑68 (MQ=255)

                                                        ctgAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAAttt             <  1:43368700/66‑1 (MQ=255)

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

                                                             tGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGGGTTGAATTTGCCGa        >  1:13600595/1‑68 (MQ=255)

                                                              gCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGccta        >  1:15271966/1‑65 (MQ=255)

                                                                cGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAAc      <  1:10046673/67‑1 (MQ=255)

                                                                cGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAAc      <  1:17223044/67‑1 (MQ=255)

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

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

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

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

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

GCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATAT  >  CP009273/691968‑692102

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