BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	691,986	T→G	5.1%	intergenic (‑26/+12)	glnX ← / ← glnV	tRNA‑Gln/tRNA‑Gln

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	691,986	0	T	G	5.1%	1302.6 / 4.1	391	intergenic (‑26/+12)	glnX/glnV	tRNA‑Gln/tRNA‑Gln
Reads supporting (aligned to +/- strand): ref base T (289/73); new base G (17/9); total (309/82)
Fisher's exact test for biased strand distribution p-value = 8.67e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 412 total aligned reads displayed.

CGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTT  >  CP009273/691920‑692053

cGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGt                                                                     <  1:6269284/67‑1 (MQ=255)

cGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGt                                                                     <  1:33754493/67‑1 (MQ=255)

cGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGt                                                                     <  1:11399972/67‑1 (MQ=255)

          aTCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATaaa                                                           <  1:26450826/67‑1 (MQ=255)

          aTCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATaa                                                            <  1:34242008/66‑1 (MQ=255)

           tCCGGTGCCTTACCGCTTGGCGATACCCCAACTAGGTGCAATTAAATGTGAAGAGTATTGAATaaaaa                                                         >  1:1235043/1‑66 (MQ=12)

           tCCGGTGCCTTACCGCTTGGCGATACCCCAACAGGGTGCTGTTAAATGGTAAGAGGATTGAATAAAta                                                         >  1:31594550/1‑67 (MQ=12)

            ccGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTg                                                        <  1:12127817/68‑1 (MQ=255)

                  ccTTACCGCTTGGCGATACCCCATCAGGGTGCAGTTGCAAGGTGAGAGGCTAGAATATCTTGGCTggg                                                  >  1:26212478/1‑68 (MQ=12)

                   cTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTgggg                                                 <  1:16692803/68‑1 (MQ=255)

                   cTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTgggg                                                 <  1:1871142/68‑1 (MQ=255)

                   cTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTgggg                                                 <  1:7572576/68‑1 (MQ=255)

                   cTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTgggg                                                 <  1:10017503/68‑1 (MQ=255)

                   cTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTggg                                                  <  1:11882228/67‑1 (MQ=255)

                   cTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTggg                                                  <  1:32029855/67‑1 (MQ=255)

                    ttACCGCTTGGCGATACCCCAACTGGGTGCACTTACTAGGTAAGCGGCTTGAATAAATTGGCTGGGGt                                                <  1:29224483/68‑1 (MQ=255)

                    ttACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGt                                                <  1:24067581/68‑1 (MQ=255)

                      aCCGCTTGGCGATACCCCATCAGGGTGCACTTAGAAGGTAATCGGAGTGAGTAAATTGGCTTGGGTAc                                              >  1:22682157/1‑68 (MQ=19)

                      aCCGCTTGGCGATACCCCAACAGGGTGCAGTTTAAAGGTAAGAGGATTGAATATAGAGGCTGGAGTAc                                              >  1:26064542/1‑68 (MQ=12)

                        cGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGa                                            >  1:16081657/1‑68 (MQ=255)

                        cGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATAGGCTGGGGTACGa                                            >  1:8019109/1‑68 (MQ=255)

                           ttGGCTATACCCCAACAGGGTGAACTTACAAGGTAAGCTTCTTTAATAAATTGGCTGGGGTACGAGGa                                         <  1:3305654/68‑1 (MQ=255)

                           ttGGCGATACCCCAACTGGGTGCACTTACCAGGTACGCGGCTTGAATAAATTGGCTGGGGTACGAGGa                                         <  1:8653633/68‑1 (MQ=255)

                           ttGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGa                                         >  1:10670750/1‑68 (MQ=255)

                           ttGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGa                                         <  1:6452513/68‑1 (MQ=255)

                           ttGGCGATACCCCAACTGGGTGCACTTACAACGTAAGCTTCTTGAATAAATTGGCTGGGGTACGAGGa                                         <  1:21120945/68‑1 (MQ=255)

                              gCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTTCGAGGATTc                                      >  1:27467811/1‑68 (MQ=255)

                              gCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGGCTTGAATAAATTGGCTGGGGTACGAGGATTc                                      >  1:2710583/1‑68 (MQ=255)

                                gATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGa                                    >  1:11523232/1‑68 (MQ=255)

                                gATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGa                                    >  1:21623419/1‑68 (MQ=255)

                                 aTACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGa                                    >  1:20882309/1‑67 (MQ=255)

                                 aTACCCCAACTGGGTGCACTTAAAAGGTAAGCGTATTGAATAAATTGGCT‑GGGTAAGAGGATTCGaa                                   >  1:17043162/1‑67 (MQ=255)

                                        aaCTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGa                            >  1:3705040/1‑68 (MQ=255)

                                        aaCTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGa                            >  1:2961488/1‑68 (MQ=255)

                                        aaCTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGa                            >  1:6765918/1‑68 (MQ=255)

                                         aCTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGaa                           >  1:19784849/1‑68 (MQ=255)

                                         aCTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGaa                           >  1:19561526/1‑68 (MQ=255)

                                         aCTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGaa                           >  1:31170499/1‑68 (MQ=255)

                                         aCTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGaa                           >  1:116547/1‑68 (MQ=255)

                                         aCTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCACGGaa                           >  1:14049945/1‑68 (MQ=255)

                                          cTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAAt                          >  1:5797896/1‑68 (MQ=255)

                                          cTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAAt                          >  1:4364119/1‑68 (MQ=255)

                                           tGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAAt                          >  1:10512769/1‑67 (MQ=255)

                                           tGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAAt                          >  1:18093499/1‑67 (MQ=255)

                                           tGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAAt                          >  1:18477003/1‑67 (MQ=255)

                                           tGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAAt                          >  1:9599901/1‑67 (MQ=255)

                                           tGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAAt                          >  1:28995612/1‑67 (MQ=255)

                                           tGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAAt                          >  1:31544377/1‑67 (MQ=255)

                                           tGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAAt                          >  1:2035384/1‑67 (MQ=255)

                                            gggTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGaaag                         >  1:26963355/1‑65 (MQ=255)

                                                  aCTTCCAAGGCAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAAt                  <  1:23402087/68‑1 (MQ=255)

                                                  aCTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGaa                   >  1:33437602/1‑67 (MQ=255)

                                                  aCTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGaa                   >  1:18986892/1‑67 (MQ=255)

                                                  aCTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGaa                   >  1:20641410/1‑67 (MQ=255)

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

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

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

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

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

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

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

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

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

                                                      aCAAGGTAATCTTCTTGAATAACTTGGCTGGGGTACGAGGATTCGAACCAGGGAATGCCGGTAt                  <  1:15586489/64‑1 (MQ=17)

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

                                                        aaGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAAt            <  1:27810506/68‑1 (MQ=255)

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

                                                          ggTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCATAATcc          >  1:24312386/1‑68 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                             aaGCTTCTTTATTAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGt       <  1:18878641/68‑1 (MQ=14)

                                                             aaGCGTCTTTTATAAATTTGCTGGGGTACGATGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGt       <  1:4273567/68‑1 (MQ=14)

                                                             aaGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGt       >  1:1251963/1‑68 (MQ=255)

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

                                                                gttCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTGGGAATGCCGGTATCAAAAACCGGTGcc    <  1:3491427/66‑1 (MQ=14)

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

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

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

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

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

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

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

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

                                                                cGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCACGGAATGCCGGTATCAAAAACCTGtc      >  1:28503457/1‑65 (MQ=17)

                                                                cGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCACGGAATGCCGGAATCAGAATCCGGTg      >  1:24783141/1‑66 (MQ=255)

                                                                cGTCTTGAATAAATTGGCTGGGGAACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTg      >  1:13765562/1‑66 (MQ=255)

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

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

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

CGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTT  >  CP009273/691920‑692053

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