BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	2,511,629	A→T	5.4%	intergenic (‑34/+6)	alaX ← / ← alaW	tRNA‑Ala/tRNA‑Ala

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	2,511,629	0	A	T	5.4%	2387.1 / 10.9	823	intergenic (‑34/+6)	alaX/alaW	tRNA‑Ala/tRNA‑Ala
Reads supporting (aligned to +/- strand): ref base A (491/270); new base T (20/41); total (511/312)
Fisher's exact test for biased strand distribution p-value = 1.70e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 843 total aligned reads displayed.

CATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGC  >  CP009273/2511562‑2511696

  tgcAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCaaa                                                                   <  1:4235901/68‑1 (MQ=255)

    cAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTCCTGGGTGCAAAAt                                                                 >  1:4941652/1‑68 (MQ=255)

    cAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTATCAGTACCAAGTTTGCTGGGGGCaaaaa                                                                 >  1:5639156/1‑67 (MQ=25)

     aaGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAAt                                                                 <  1:45532089/67‑1 (MQ=255)

      aGCGCTCTCCCAGCTGAGCTATAGCCCCGCGAGGCGTTTAGGTTCGAAGTTTGCGGGGGGCAAAtttt                                                               >  1:36564510/1‑68 (MQ=18)

      aGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTAAGTACCAAGTTTGCTGCGTGCTAAtttt                                                               >  1:7157207/1‑68 (MQ=25)

        cgcTCTCCCAGCTGAGCTATAGCCCCAGGATGCGTTTACGTACCAAGTTTGGTGGGTGAAAAATTtgg                                                             >  1:54958567/1‑68 (MQ=37)

         gcTCTCCCAGCTGAGCTATAGCCCCCGGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTtggt                                                            >  1:44423149/1‑68 (MQ=255)

         gcTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTtggt                                                            >  1:15330723/1‑68 (MQ=255)

         gcTCTCCCAGCTGAGCTATAGCCCCACCATGCGTTTACGTACCAAGTTTGCGGGGTGCAAAATTtggt                                                            >  1:21915271/1‑68 (MQ=38)

          ctctCCCAGCTGAGCTATAGCCCCCCGATGCGTTTACGTACCAAGTTTGCTGGGTGAAAAATTtggtg                                                           >  1:20172612/1‑68 (MQ=38)

            ctcCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGa                                                         <  1:22890197/68‑1 (MQ=255)

            ctcCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGa                                                         <  1:33496026/68‑1 (MQ=255)

             tcCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAg                                                        <  1:11432393/68‑1 (MQ=255)

             tcCCAGCGGAGCTATAGCCCCACTAACCGTTTACGTACCAAGGTTGCTGGGTGCAAAATTTGGTGGAg                                                        <  1:46038551/68‑1 (MQ=37)

              cccAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGc                                                       <  1:22382316/68‑1 (MQ=255)

                 aGCTGAGCTATAGCCCCGAGAGGGGTTTACGTACCAATTTAGGTGCGTCCAAAATTTGGTTGAGCaaa                                                    >  1:30538525/1‑68 (MQ=255)

                 aGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTaa                                                    >  1:13018056/1‑68 (MQ=255)

                 aGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGAAAAATTTGGTGGAGCTaa                                                    >  1:24613791/1‑68 (MQ=255)

                 aGCTGAGCTATAGCCCCACGAGGCGTTTAGGTACCAAGTTTGCGGGGTCAAAATTTTGATTGAGCTaa                                                    >  1:10266100/1‑68 (MQ=255)

                  gCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTaa                                                    >  1:57729315/1‑67 (MQ=255)

                   cTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCTAAATTTGGTGGAGCTAAGc                                                  >  1:25816298/1‑68 (MQ=255)

                   cTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGc                                                  >  1:6412688/1‑68 (MQ=255)

                    tGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTaa                                                    <  1:23511510/65‑1 (MQ=255)

                    tGAGCTATAGCCCCACCCTGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTaa                                                    <  1:16094564/65‑1 (MQ=255)

                       gCTATAGCCCCACGATGCGGTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCggg                                               <  1:5545607/67‑1 (MQ=255)

                       gCTATACCCCCACGATGCCTTCACGTACCAAGTTTGCTGGGTGGTAAATTTGGTGGAGCTAAGCggg                                               <  1:36021635/67‑1 (MQ=37)

                           taGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGa                                          <  1:48358740/68‑1 (MQ=255)

                            aGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGTTGGAGCTAAGCGGGATCGa                                          >  1:4042863/1‑67 (MQ=255)

                            aGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGa                                          >  1:21100299/1‑67 (MQ=255)

                            aGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGa                                          >  1:48985411/1‑67 (MQ=255)

                            aGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGa                                          >  1:40665885/1‑67 (MQ=255)

                              ccccACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAAcc                                       >  1:25062423/1‑68 (MQ=255)

                              ccccACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAAcc                                       >  1:32415173/1‑68 (MQ=255)

                              ccccACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAAcc                                       >  1:46964670/1‑68 (MQ=255)

                              ccccACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAAcc                                       >  1:3096445/1‑68 (MQ=255)

                               cccACGCTGCGATACCGCACCAAGTTGCCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCg                                      <  1:57194673/68‑1 (MQ=18)

                               cccACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAccc                                       >  1:41934428/1‑67 (MQ=255)

                               cccACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAAcc                                       >  1:33172727/1‑67 (MQ=255)

                               cccACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAAcc                                       >  1:37678853/1‑67 (MQ=255)

                               cccACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAAcc                                       >  1:3449496/1‑67 (MQ=255)

                               cccACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAAcc                                       >  1:2804318/1‑67 (MQ=255)

                               cccACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAAcc                                       >  1:55975802/1‑67 (MQ=255)

                               cccACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAAcc                                       >  1:8680911/1‑67 (MQ=255)

                               cccACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAAcc                                       >  1:19474387/1‑67 (MQ=255)

                               cccACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAAcc                                       >  1:45085916/1‑67 (MQ=255)

                               cccACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAAcc                                       >  1:13837603/1‑67 (MQ=255)

                               cccACGAACCGATTACGGACCAAGTATGCTGGGTGCAAAATTTGGTGGANCNAAGCGGGATCGAACCg                                      <  1:52659406/68‑1 (MQ=21)

                                 cACGCTGCGTTTACGCACCAAGGTGGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCt                                    <  1:17596219/68‑1 (MQ=21)

                                 cACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAAcc                                       >  1:8316707/1‑65 (MQ=255)

                                 cACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAAcc                                       >  1:19227875/1‑65 (MQ=255)

                                 cACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCAAAGCGGGATCGAAcc                                       >  1:12263583/1‑65 (MQ=255)

                                   cGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTg                                   >  1:39301350/1‑67 (MQ=255)

                                     caccGTTTACGTACCAAGTTTGCTGGGTGGAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGAcc                                <  1:29094931/65‑1 (MQ=255)

                                      tGCTTTCACGTACCAAGTTTGCTGGGTGCAAAATTTGGTTGAGCTAAGCGGGATCGAACCGCTGACct                               <  1:31774203/68‑1 (MQ=21)

                                      tGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGAcc                                >  1:18870098/1‑67 (MQ=255)

                                      tGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGAcc                                >  1:46198597/1‑67 (MQ=255)

                                      agCGTATACGTACCAAGGTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGAcc                                <  1:32736143/66‑1 (MQ=21)

                                      accGTTTACGCACCAAGCTTGCTGGGTGCTAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACct                               <  1:35377452/66‑1 (MQ=18)

                                        cctttACGCACCAAGCTTGCTGGGTGGAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACct                               <  1:28548157/64‑1 (MQ=18)

                                        cGTTTACGTACCAAGTTTGCTGGGTTCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACctct                             <  1:11791684/68‑1 (MQ=255)

                                         gtatACGTACCAAGGTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACctct                             <  1:6749833/64‑1 (MQ=255)

                                         gattACGCACGAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACct                               <  1:4341963/63‑1 (MQ=21)

                                         gTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACctct                             <  1:33937448/67‑1 (MQ=255)

                                         gTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCtt                            >  1:17909050/1‑68 (MQ=255)

                                         gTTAACGTACCAAGTTTGCTGGGTTCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACctct                             <  1:5882657/67‑1 (MQ=255)

                                          tttACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCtt                            >  1:50589541/1‑67 (MQ=255)

                                            tACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGca                         >  1:54756868/1‑68 (MQ=255)

                                             aCGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGcat                        >  1:20776482/1‑68 (MQ=255)

                                              cGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGcat                        >  1:9544784/1‑67 (MQ=255)

                                              cGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGcat                        >  1:42704459/1‑67 (MQ=255)

                                              cGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGcat                        >  1:11161998/1‑67 (MQ=255)

                                              cGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGcat                        >  1:31277989/1‑67 (MQ=255)

                                              cGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGcat                        >  1:37013781/1‑67 (MQ=255)

                                                tcccAAGTTTGCTGGGTGCAATATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGcatgcc                     <  1:41541046/66‑1 (MQ=255)

                                                tACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGcatgcc                     >  1:5345777/1‑68 (MQ=255)

                                                 aCCAAGTTTGCTGGGTGGTAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGcatgcc                     <  1:17017214/67‑1 (MQ=255)

                                                 aCCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGcatgcca                    >  1:2697989/1‑68 (MQ=255)

                                                 aCCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGcatgcca                    >  1:16550667/1‑68 (MQ=255)

                                                 aCCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGcatgcc                     <  1:38577628/67‑1 (MQ=255)

                                                 aCCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGcatgcc                     <  1:9278632/67‑1 (MQ=255)

                                                 aCCAAGTTTGCAGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGcatgcc                     <  1:51155090/67‑1 (MQ=255)

                                                  ccAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGcatgccat                   <  1:29854606/68‑1 (MQ=255)

                                                       tttGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCaa               >  1:14811013/1‑67 (MQ=255)

                                                       tttGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCaa               >  1:27483891/1‑67 (MQ=255)

                                                       tttGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCaa               >  1:26641082/1‑67 (MQ=255)

                                                       tttGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCaa               >  1:19883153/1‑67 (MQ=255)

                                                        ttGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAgc             >  1:12655356/1‑68 (MQ=255)

                                                            tGGGTTCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGctc         <  1:37288003/68‑1 (MQ=255)

                                                            tGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAgcgc           >  1:38238961/1‑66 (MQ=255)

                                                            tGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAgcgc           >  1:1818636/1‑66 (MQ=255)

                                                            tGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAgcgc           >  1:50080931/1‑66 (MQ=255)

                                                            tGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGctc         <  1:44006123/68‑1 (MQ=255)

                                                             gggTTCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGctct        <  1:17245211/68‑1 (MQ=25)

                                                             gggTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGctct        <  1:23962817/68‑1 (MQ=255)

                                                             gggTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGctct        <  1:5390845/68‑1 (MQ=255)

                                                             gggTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGctct        <  1:5195175/68‑1 (MQ=255)

                                                               gTGGAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCtctc       <  1:7681697/67‑1 (MQ=18)

                                                               gTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCtctc       <  1:25988850/67‑1 (MQ=39)

                                                                 gCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTccc     <  1:47569677/67‑1 (MQ=35)

CATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGC  >  CP009273/2511562‑2511696

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