BRESEQ :: Evidence

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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	692,358	0	T	A	5.9%	1060.0 / 13.6	348	noncoding (38/75 nt)	glnU	tRNA‑Gln
Reads supporting (aligned to +/- strand): ref base T (184/141); new base A (2/21); total (186/162)
Fisher's exact test for biased strand distribution p-value = 4.72e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 360 total aligned reads displayed.

CCGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCATCCGTACAACGCTTTCTGGTGAATGGTGCG  >  CP009273/692291‑692425

ccGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGt                                                                     >  1:12266975/1‑68 (MQ=255)

ccGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGt                                                                     >  1:9113280/1‑68 (MQ=255)

ccGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGt                                                                     >  1:13747235/1‑68 (MQ=255)

 cGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTa                                                                    >  1:22762360/1‑68 (MQ=255)

 cGGGCGGTGAACCGCTTACTCGAAGAAGATGGATGGGGTACCTGGATTCGAACCAGGGAATGCCGGTa                                                                    >  1:25869671/1‑68 (MQ=255)

   ggCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATG‑CGGTATCa                                                                 >  1:23285288/1‑68 (MQ=255)

          aaCCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAAcc                                                           >  1:30209984/1‑68 (MQ=255)

          aaCCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAAcc                                                           >  1:54511/1‑68 (MQ=255)

          aaCCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAAcc                                                           >  1:6227171/1‑68 (MQ=255)

              gCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGt                                                        >  1:11610635/1‑67 (MQ=255)

              gCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGt                                                        >  1:19875187/1‑67 (MQ=255)

              gCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGt                                                        >  1:32702505/1‑67 (MQ=255)

              gCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGt                                                        >  1:15397188/1‑67 (MQ=255)

              gCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGt                                                        >  1:6501537/1‑67 (MQ=255)

              gCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGt                                                        >  1:14653202/1‑67 (MQ=255)

                  aCTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCtt                                                   >  1:26446175/1‑68 (MQ=255)

                  aCTAGAATAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCtt                                                   <  1:34154254/68‑1 (MQ=255)

                   cTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCtt                                                   >  1:24528233/1‑67 (MQ=255)

                   cTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCtt                                                   >  1:23843665/1‑67 (MQ=255)

                   cTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCtt                                                   >  1:1475792/1‑67 (MQ=255)

                    tCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTa                                                  >  1:15525575/1‑67 (MQ=255)

                     cGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTAcc                                                >  1:18695007/1‑68 (MQ=255)

                          aagcTGGCTGGGGTACCTGGATTCGAACCACGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTg                                           <  1:21335450/68‑1 (MQ=255)

                          aagaTGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTg                                           >  1:37037754/1‑68 (MQ=33)

                                      gTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAt                               >  1:17228846/1‑68 (MQ=31)

                                       tACCTGGATTGGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAt                               >  1:9296498/1‑67 (MQ=14)

                                       tACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAtt                              >  1:18398066/1‑67 (MQ=14)

                                        aCCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCATc                              <  1:8353609/67‑1 (MQ=33)

                                          cTGGATTCGAACCTCGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCATCCGt                           <  1:200340/68‑1 (MQ=255)

                                          cTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAtgag                            >  1:3037978/1‑64 (MQ=14)

                                          cTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCATCCGt                           >  1:11000357/1‑68 (MQ=39)

                                          cTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCATCCGt                           <  1:32195716/68‑1 (MQ=39)

                                          cTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCATCCGt                           >  1:36750739/1‑68 (MQ=39)

                                           tGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCATCCGt                           >  1:13569301/1‑67 (MQ=39)

                                           tGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCATCCGt                           >  1:4980673/1‑67 (MQ=39)

                                           aggATACGAACCTCGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCATCCGt                           >  1:33672382/2‑67 (MQ=14)

                                            ggATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCATACgga                          >  1:20766237/1‑65 (MQ=9)

                                               ttCGAACCTCGGAATGCCGGAATCATAAACCGGTGCCTTACCGCTTGGCGATACCCCATCCGTACaa                       <  1:328606/67‑1 (MQ=11)

                                               ttCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCATCCGTACAAc                      >  1:18888784/1‑68 (MQ=255)

                                               ttCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCATCCGTACAAc                      >  1:36429303/1‑68 (MQ=255)

                                               ttCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCATCCGTACAAc                      >  1:38734440/1‑68 (MQ=255)

                                               ttCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCATCCGTACAAc                      >  1:2734187/1‑68 (MQ=255)

                                                 cGAACCTGGGAATGCCGGTATCATAAACCGGTGCCTTACCGCTTGGCGATACCCCATCCGTACAACGc                    <  1:8112565/68‑1 (MQ=255)

                                                 cGAACCAGGGAATGCCGGTATCATAAACCGGTGCCTTACCGCTTGGCGATACCCCATCCGTACAACGc                    <  1:5699135/68‑1 (MQ=255)

                                                 cGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCATCCGTACAACGc                    <  1:20357003/68‑1 (MQ=255)

                                                 cGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCATCCGTACAACGc                    <  1:31917880/68‑1 (MQ=255)

                                                 cGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCATCCGTACAACGc                    <  1:21237653/68‑1 (MQ=255)

                                                   aaCCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCATCCGTACAACGCtt                  >  1:37929061/1‑68 (MQ=255)

                                                   aaCCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCATCCGTACAACGCtt                  >  1:36227240/1‑68 (MQ=255)

                                                   aaCCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCATCAGAACAACGGtt                  >  1:6712024/1‑68 (MQ=14)

                                                    aCCTGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCATCCGTACAACGCttt                 <  1:2884828/68‑1 (MQ=255)

                                                    aCCTGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCATCCGTACAACGCtt                  <  1:10177594/67‑1 (MQ=255)

                                                    aCCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCATCCGTACAACGCtt                  >  1:37392093/1‑67 (MQ=255)

                                                    aCCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCATCCGTACAACGCtt                  >  1:31508727/1‑67 (MQ=255)

                                                    aCCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCATCCGTACAACGCtt                  >  1:35353775/1‑67 (MQ=255)

                                                    aCCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCATCCGTACAACGCtt                  >  1:26914368/1‑67 (MQ=255)

                                                    aCCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCATCCGGACAACGCttt                 >  1:24178916/1‑68 (MQ=255)

                                                    aCCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAGCCGTACAACGCtt                  >  1:18211052/1‑67 (MQ=255)

                                                     ccAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCATCCGTACAACGCtt                  >  1:29959510/1‑66 (MQ=255)

                                                     ccAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCATCCGTACAACGCtt                  >  1:25173965/1‑66 (MQ=255)

                                                      ctcggAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCATCCGTACAACGCTTTCt               <  1:29010925/65‑1 (MQ=14)

                                                      ctcggAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCATCCGTACAACGCTTTCt               <  1:141247/65‑1 (MQ=14)

                                                      cacggAATGCCGGTATCACAAACCGGTGCCTTACCGCTTGGCGATACCCCATCCGTACAACGCTTTCt               <  1:29686620/65‑1 (MQ=255)

                                                      cacggAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCATCCGTACAACGCTTTCt               <  1:1616577/65‑1 (MQ=255)

                                                      cAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCATCCGTACAACGCTTTCt               <  1:22404482/68‑1 (MQ=255)

                                                      cAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCATCCGTACAACGCTTTCt               <  1:9633081/68‑1 (MQ=255)

                                                      cAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCATCCGTACAACGCTTTCt               <  1:16706674/68‑1 (MQ=255)

                                                      cAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCATCCGTACAACGCTTTCt               <  1:16958314/68‑1 (MQ=255)

                                                      cAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCATCCGTACAACGCTTTCt               <  1:4718530/68‑1 (MQ=255)

                                                      cAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCATCCGTACAACGCTTTCt               <  1:26158641/68‑1 (MQ=255)

                                                        gggAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCATCCGTACAACGCTTTCt               <  1:12079583/66‑1 (MQ=255)

                                                        cggAATGCCGGTATCATAATCCGGTGCCTTACCGCTTGGCGATACCCCATCCGTACAACGCTTTCTgg             <  1:27989438/67‑1 (MQ=18)

                                                         ggAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATCCCCCATCCGTACAACGCTTTCTGGt            >  1:12903954/1‑68 (MQ=255)

                                                         ggAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCATCCgga                          >  1:10288947/1‑52 (MQ=255)

                                                         ggAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCATCCGTACAACGCTTTCTGGt            >  1:18040630/1‑68 (MQ=255)

                                                         ggAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCATCCGTACAACGCTTTCTGGt            >  1:13139858/1‑68 (MQ=255)

                                                         ggAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCATCCGTACAACGCTTTCTGGt            >  1:6950876/1‑68 (MQ=255)

                                                         ggAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCATCCGTACAACGCTTTCTGGt            >  1:14465722/1‑68 (MQ=255)

                                                         ggAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCATCCGTACAACGCTTTCTGGt            >  1:37703740/1‑68 (MQ=255)

                                                          gAATGCCGGTATCAAAAACGGGTGCCTTACCGCTTGGCGATACCCCATCCGTACAACGCTTTCTgg             <  1:6744395/66‑1 (MQ=255)

                                                          gAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCATCCGTACAACGCTTTCTGGTg           >  1:35180071/1‑68 (MQ=255)

                                                            aTGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCATCCGTACAACGCTTTCTGGTGa          <  1:29377138/67‑1 (MQ=255)

                                                            aTGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCATCCGTACAACGCTTTCTGGTGa          <  1:28305905/67‑1 (MQ=255)

                                                              gCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCATCCGTACAACGCTTTCTGGTGAATg       >  1:19501163/1‑68 (MQ=255)

                                                               ccGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCATCCGTACAACGCTTTCTGGTGAATgg      >  1:4402232/1‑68 (MQ=255)

                                                               ccGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCATCCGTACAACGCTTTCTGGTGAATgg      >  1:35808887/1‑68 (MQ=255)

                                                               ccGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCATCCGTACAACGCTTTCTGGTGAATgg      >  1:6018771/1‑68 (MQ=255)

                                                               ccGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCATCCGTACAACGCTTTCTGGTGAATgg      >  1:7614445/1‑68 (MQ=255)

                                                                cGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCATCCGTACAACGCTTTCTGGTGAATgg      <  1:31053626/67‑1 (MQ=255)

                                                                cGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCATCCGTACAACGCTTTCTGGTGAATGGt     <  1:26068395/68‑1 (MQ=255)

                                                                cGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCATCCGTACAACGCTTTCTGGTGAATGGt     <  1:35035236/68‑1 (MQ=255)

                                                                cGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCATCCGTACAACGCTTTCTGGTGAATGGt     <  1:7154813/68‑1 (MQ=255)

                                                                cGGTATCAAAAACCGGTGCCTTAACGCTTGGCGATACCCCATCCGTACAACGCTTTCTGGTGAATGGt     <  1:21640419/68‑1 (MQ=255)

                                                                 ggTATCATAAACCGGTGCCTTACCGCTTGGCGATACCCCATCCGTACAACGCTTTCTGGTGAATGGTg    <  1:34476156/68‑1 (MQ=255)

                                                                 ggTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCATCCGTACAACGCTTTCTGGTGAATGGt     >  1:3124702/1‑67 (MQ=255)

                                                                 ggTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCATCCGTACAACGCTTTCTGGTGAATGGt     >  1:13945815/1‑67 (MQ=255)

                                                                 ggTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCATCCGTACAACGCTTTCTGGTGAATGGTg    <  1:17440036/68‑1 (MQ=255)

                                                                 ggTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCATCCGTACAACGCTTTCTGGTGAATGGTg    <  1:12497366/68‑1 (MQ=255)

                                                                  gTATCAAAAACCGTTGCCTTTCCGCTTGGCGATACCCCATGCGTACATCGCTTTCTGGTGAATGGTg    <  1:11146758/67‑1 (MQ=21)

                                                                  gTAGCAAAAACCGTTGCCTTACAGCTTGGCGATACCCCATCCGTACAACGCTTTCTGGTGAATGGTg    <  1:39127654/67‑1 (MQ=25)

CCGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCATCCGTACAACGCTTTCTGGTGAATGGTGCG  >  CP009273/692291‑692425

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