BRESEQ :: Evidence

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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	692,257	0	A	T	5.9%	813.9 / 7.8	287	noncoding (30/75 nt)	glnW	tRNA‑Gln
Reads supporting (aligned to +/- strand): ref base A (149/113); new base T (2/20); total (152/135)
Fisher's exact test for biased strand distribution p-value = 9.61e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 364 total aligned reads displayed.

CGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCGAAGAAGATGGC  >  CP009273/692190‑692324

cGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCaaaaa                                                                     <  1:25815289/68‑1 (MQ=255)

 gTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGTGAATGCAGGTATCaaaaa                                                                     >  1:37351703/1‑67 (MQ=255)

 gTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCaaaaa                                                                     >  1:17099161/1‑67 (MQ=255)

 gTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCaaaaa                                                                     >  1:39362800/1‑67 (MQ=255)

 gTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCaaaaa                                                                     >  1:44454278/1‑67 (MQ=255)

 gTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCaaaaa                                                                     >  1:4575134/1‑67 (MQ=255)

 gTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCaaaaa                                                                     >  1:460970/1‑67 (MQ=255)

 gTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCaaaaa                                                                     >  1:26807775/1‑67 (MQ=255)

 gTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCaaaaa                                                                     >  1:14249553/1‑67 (MQ=255)

 gTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCaaaaa                                                                     >  1:10660727/1‑67 (MQ=255)

 gTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCAGGTATCaaaaa                                                                     >  1:36169953/1‑67 (MQ=255)

  tAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAAcc                                                                   >  1:57532255/1‑68 (MQ=255)

     ccAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGt                                                                >  1:55950587/1‑68 (MQ=255)

      cAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGt                                                                >  1:16002244/1‑67 (MQ=255)

      cAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGt                                                                >  1:37500649/1‑67 (MQ=255)

      cAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGt                                                                >  1:17818019/1‑67 (MQ=255)

      cAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGt                                                                >  1:38495235/1‑67 (MQ=255)

      cAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGt                                                                >  1:31280419/1‑67 (MQ=255)

      cAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGt                                                                >  1:47266349/1‑67 (MQ=255)

      cAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGt                                                                >  1:13322798/1‑67 (MQ=255)

      cAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATC‑AAAACCGGTg                                                               >  1:54960770/1‑67 (MQ=255)

      cAGATTGTTTCTTCGATGGCTGGGGTACATGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGt                                                                >  1:48933972/1‑67 (MQ=255)

        gATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTg                                                               >  1:882231/1‑66 (MQ=255)

         aTTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCt                                                            >  1:46875871/1‑68 (MQ=255)

          ttttttCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGcc                                                             <  1:5073584/63‑1 (MQ=255)

          ttGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGcc                                                             <  1:37183475/66‑1 (MQ=255)

          ttGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGc                                                              <  1:11635154/65‑1 (MQ=255)

          ttGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTAT‑AAAAACCGGTGcc                                                             <  1:27956283/65‑1 (MQ=255)

          ttGGTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGc                                                              <  1:37593979/65‑1 (MQ=255)

          tcgTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGcc                                                             <  1:32974278/64‑1 (MQ=255)

               tcttcGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGc                                                      <  1:30206894/68‑1 (MQ=39)

               tcttcGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGc                                                      <  1:1314836/68‑1 (MQ=39)

                   cGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTg                                                   >  1:51599642/1‑67 (MQ=31)

                              gTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCaa                                       >  1:18702255/1‑68 (MQ=31)

                                aCCAGGATTCGAACCAGGGAATGCCGGAATCATAAACCGGTGCCTTACCGCTTGGCGATACCCCAATa                                     <  1:28716217/68‑1 (MQ=11)

                                  cTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATag                                    >  1:25352018/1‑66 (MQ=25)

                                   tGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAaaaaa                                   >  1:19830756/1‑66 (MQ=9)

                                   tGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAAcc                                  >  1:21858290/1‑68 (MQ=40)

                                     gATTCGAACCAGGGAATGCCGGTATCAAAAACCGGAGCCTTACCGCTTGGCGATACCCCAATAACCgg                                >  1:8120626/1‑68 (MQ=255)

                                       ttCGAACCTGGGTATGAGTGTATCATAAGACTGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGc                              <  1:49466664/68‑1 (MQ=12)

                                       ttCGAACCTGGGAATGCCGGTATCTTAATCCGGTGCCTTACCGGTTGGCGATACCCCAATAACCGGGc                              <  1:9680712/68‑1 (MQ=14)

                                       ttCGAACCTGGGAATGCCGGTATCATAATCCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGc                              <  1:48273289/68‑1 (MQ=14)

                                       ttCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCATTAACCGGGc                              >  1:21293214/1‑68 (MQ=255)

                                       ttCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGc                              <  1:51260984/68‑1 (MQ=255)

                                       ttCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCCGGc                              >  1:5575932/1‑68 (MQ=255)

                                       ttCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAAAAGCCGGGc                              >  1:20093786/1‑68 (MQ=255)

                                       ttCGAACCACGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGc                              <  1:40859182/68‑1 (MQ=255)

                                        tCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCg                             >  1:23104425/1‑68 (MQ=255)

                                          gAACCTCGGAATGCCGGTATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGt                           <  1:4274858/68‑1 (MQ=11)

                                          gAACCTCGGAATGCCGGTATCACAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGt                           <  1:5440599/68‑1 (MQ=17)

                                          gAACCTCGGAATGCCGGAATCACAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGt                           <  1:8651189/68‑1 (MQ=14)

                                          gAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGt                           <  1:11484263/68‑1 (MQ=255)

                                           aaCCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTg                          >  1:38953444/1‑68 (MQ=255)

                                            aCCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAAGCGGGCGGTg                          >  1:48733545/1‑67 (MQ=255)

                                             ccAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCGCAAAAACCGGGCGGttga                        >  1:28296014/1‑65 (MQ=17)

                                             ccAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGa                         >  1:49984185/1‑67 (MQ=255)

                                             ccAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAAGAAACGGGCGGTgt                         >  1:48179662/1‑66 (MQ=14)

                                               acggAATGCCGGTATCAAAACCCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAAcc                      <  1:43035851/66‑1 (MQ=255)

                                                cggAATGCCGGTATCACAATCCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCg                     <  1:52664030/67‑1 (MQ=18)

                                                cggAATGCCGGAATCATAATCCTGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCg                     <  1:41463327/67‑1 (MQ=17)

                                                 ggAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATATCCGGGCGGTGAACCGc                    >  1:41958767/1‑68 (MQ=255)

                                                 ggAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAAcc                      >  1:46201581/1‑66 (MQ=255)

                                                  gAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCg                     <  1:1542920/66‑1 (MQ=255)

                                                  gAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCt                   >  1:31501708/1‑68 (MQ=255)

                                                  gAATGCCGGTATC‑AAAACCGGTGCCTTACCGCTTGGCGATACCCCAATACCAGGGCGGTGAACCGGtt                  >  1:12832235/1‑68 (MQ=17)

                                                    aTGCCGGAATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTa                 <  1:24234919/68‑1 (MQ=255)

                                                     tGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTAc                <  1:24586837/68‑1 (MQ=255)

                                                     tGCCGGAATCATAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTAc                <  1:5024453/68‑1 (MQ=255)

                                                      gCCGGTATCACAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACt               <  1:6328045/68‑1 (MQ=255)

                                                      gCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACt               <  1:43970141/68‑1 (MQ=255)

                                                      gCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACt               <  1:7125179/68‑1 (MQ=255)

                                                      gCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACt               <  1:56335568/68‑1 (MQ=18)

                                                       ccGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGACGCGCTTAATc              >  1:55622022/1‑68 (MQ=18)

                                                       ccGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTc              >  1:44570310/1‑68 (MQ=255)

                                                        cGGTATCACAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTc              <  1:33631315/67‑1 (MQ=255)

                                                        cGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTc              <  1:19624514/67‑1 (MQ=255)

                                                        cGGAATCATAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCg             <  1:45198704/68‑1 (MQ=255)

                                                         ggTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCga            <  1:5553986/68‑1 (MQ=255)

                                                         ggTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCga            <  1:2618246/68‑1 (MQ=255)

                                                         ggTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCga            <  1:10181164/68‑1 (MQ=255)

                                                           tATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCgaa           <  1:14840002/67‑1 (MQ=255)

                                                           aaTCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCgaa           <  1:23702366/66‑1 (MQ=255)

                                                            aTCATAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCgaaga         <  1:32363620/68‑1 (MQ=255)

                                                            aTCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCgaaga         <  1:18502381/68‑1 (MQ=255)

                                                            aTCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCgaaga         <  1:2708987/68‑1 (MQ=255)

                                                            aTCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCgaaga         <  1:35668967/68‑1 (MQ=255)

                                                             tCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGaa                        >  1:29144181/1‑52 (MQ=17)

                                                                 aaaCCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCGAAGAAGATg    <  1:14321862/68‑1 (MQ=255)

                                                                 aaaCCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCGAAGAAGATg    <  1:30769565/68‑1 (MQ=255)

                                                                  aaCCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCGAAGAAGATgg   <  1:1519141/68‑1 (MQ=255)

                                                                   aCCGGTGCCTTACCGCTTGGCGATACCCCATTAGCCGGGCGGTGAACCGCTTAATCGAAGAAGAtc    >  1:29429812/1‑65 (MQ=25)

                                                                   aCCGGTGCCTTACCGCTTGGCGATACCCCATGAACCGGGCGGTTTACCGGTTACTCGGAGTTGATg    >  1:35171722/1‑66 (MQ=12)

                                                                   aCCGGTGCCTTACCGCTTGGCGATACCCCATACACCGGGCGGTGAACCGCGTACTCGGAGAAGATGGc  >  1:22616230/1‑68 (MQ=21)

                                                                   aCCGGTGCCTTACCGCTTGGCGATACCCCAATACCCGGGCGGTGACCCGCTTACTGGAAGatgatg    >  1:1217750/1‑66 (MQ=25)

                                                                   aCCGGTGCCTTACCGCTTGGCGATACCCCAATACACGGGCGGTGAACCGCGTAATCGAAGAAGATg    >  1:46501562/1‑66 (MQ=25)

                                                                   aCCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGGGTACTCGGAGAAGAGGGc  >  1:20623722/1‑68 (MQ=25)

                                                                   aCCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCGTAGAAGATg    >  1:15878266/1‑66 (MQ=255)

                                                                   aCCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCGAAGAAGATg    >  1:36351397/1‑66 (MQ=255)

                                                                   aCCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGCTGACACGATTAATCGAAGAAAAtc    >  1:34220605/1‑65 (MQ=18)

                                                                   aCCGGTGCCTTACCGCTTGGCGATACCCAAAAAAACGGCCGGTGAAGCGATAACTCGAAGAAGAtcga  >  1:30058731/1‑65 (MQ=18)

CGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCGAAGAAGATGGC  >  CP009273/692190‑692324

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