BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	2,511,621	C→G	6.0%	intergenic (‑26/+14)	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,621	0	C	G	6.0%	1898.7 / 15.5	627	intergenic (‑26/+14)	alaX/alaW	tRNA‑Ala/tRNA‑Ala
Reads supporting (aligned to +/- strand): ref base C (383/200); new base G (34/10); total (417/210)
Fisher's exact test for biased strand distribution p-value = 1.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 635 total aligned reads displayed.

GCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGC  >  CP009273/2511556‑2511687

 catgccatgcAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCt                                                                   >  1:31241661/1‑66 (MQ=255)

 catgccatgcAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTgg                                                                 >  1:38843918/1‑68 (MQ=255)

   tgccatgcAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGt                                                               <  1:34500612/68‑1 (MQ=255)

   tgccatgcAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGt                                                               <  1:8481119/68‑1 (MQ=255)

   tgccatgcAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGt                                                               <  1:25284052/68‑1 (MQ=255)

    gccatgcAAGCGCTCTCCCAGCTGAGCTATAGCCCCGGGATGCGTTTACGTACGATGTTTGGTGGGt                                                               >  1:17614854/1‑67 (MQ=18)

    gccatgcAAGCGCTCTCCCAGCTGAGCTATAGCCCCGCGATGCGTTTACGTACCAAGTTTGGTGGGt                                                               >  1:19465095/1‑67 (MQ=255)

    gccatgcAAGCGCTCTCCCAGCTGAGCTATAGCCCCGCGATGCGTTTACGTACCAAGTTTGCTGGGt                                                               >  1:21937476/1‑67 (MQ=255)

    gccatgcAAGCGCTCTCCCAGCTGAGCTATAGCCCCCGGAGGCGTTTACGTACCAAGTTTggtgggtg                                                              >  1:9037280/1‑68 (MQ=21)

    gccatgcAAGCGCTCTCCCAGCTGAGCTATAGCCCCCCGAGGCGTTTACGTACCAAGTTTCCTGGGt                                                               >  1:24284304/1‑67 (MQ=21)

    gccatgcAAGCGCTCTCCCAGCTGAGCTATAGCCCCAGGAGGGGTTTAGGTACGAAGTTTCGTGGGTg                                                              >  1:12356690/1‑68 (MQ=17)

    gccatgcAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGGGTTTACGTACGAAGTTTGCTGGGt                                                               >  1:37200282/1‑67 (MQ=255)

    gccatgcAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAATTTTGCTGGGt                                                               >  1:20325417/1‑67 (MQ=255)

    gccatgcAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGt                                                               >  1:36915963/1‑67 (MQ=255)

    gccatgcAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTg                                                              <  1:43569391/68‑1 (MQ=255)

    gccatgcAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCGGGGt                                                               >  1:11633793/1‑67 (MQ=255)

    gccatgcAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGAGGCGTTTAAGTACCATTTTTGGTgggg                                                               >  1:16939986/1‑66 (MQ=18)

    gccatgcAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGAGGCGATTACGGACCAAGTTTGCTGGGt                                                               >  1:32873379/1‑67 (MQ=21)

      catgcAAGCGCTCTCCCAGCTGAGCTATAGCCCCGGGATGAGTTTACTTTCCATGTTAGGTGGGTGCa                                                            >  1:22497936/1‑68 (MQ=14)

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

            aGCGCTCTCCCAGCTGAGCTATAGCCCCCCGATGCGTTTACGTTCCAGGTTTCCTGGGTGCAAAAttt                                                      >  1:7396828/1‑68 (MQ=25)

            aGCGCTCTCCCAGCTGAGCTATAGCCCCCCGATGATATTACGTACCATGTTTCGTGGGTGCTAAAttt                                                      >  1:33151405/1‑68 (MQ=18)

            aGCGCTCTCCCAGCTGAGCTATAGCCCCCCGAGGCGTTTAAGTACCAAGTTTGCTGGGTGCAAAAttt                                                      >  1:4760553/1‑68 (MQ=37)

            aGCGCTCTCCCAGCTGAGCTATAGCCCCCCGAAGCGTTTACGTACCAAGTTACCTGGGTGCAAAAttt                                                      >  1:16421484/1‑68 (MQ=25)

              cgcTCTCCCAGCTGAGCTATAGCCCCACGATGGGTTTACGTACCAAGTTTGCTGGGGGCAAAATTtgg                                                    >  1:10536512/1‑68 (MQ=255)

              cgcTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTCCTGCGTGCAAAATTtgg                                                    >  1:5249853/1‑68 (MQ=255)

               gcTCTCCCAGCTGAGCTATAGCCCCACGATTCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTtggt                                                   <  1:40615936/68‑1 (MQ=255)

               gcTCTCCCAGCTGAGCTATAGCCCCACGATGCATTTAATTACGAAGTTTGCGGCGTGCAAAATTtggt                                                   >  1:37563525/1‑68 (MQ=21)

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

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

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

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

                       atcTGAGCTATAGCCCCACGATTCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTaa                                           <  1:32288970/66‑1 (MQ=255)

                       aGCTGAGCTATAGCCCCCCCATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCCaa                                           >  1:23194874/1‑68 (MQ=255)

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

                       aGCTGAGCTATAGCCCCACGAGGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTa                                            >  1:26104681/1‑67 (MQ=255)

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

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

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

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

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

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

                             gCTATATCCCCACTATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCggg                                      <  1:20039097/67‑1 (MQ=255)

                             gCTATACCCCCACGATCCGTTTACGTACCAAGTATGCTGGGTGCAAAATTTGGTGGAGCTAAGCggg                                      <  1:30267382/67‑1 (MQ=38)

                               tataGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATc                                   >  1:27341092/1‑68 (MQ=255)

                               tataGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATc                                   >  1:13479868/1‑68 (MQ=255)

                               tataGCCCCACCATTCTTTTACGAACGAAGTTTGCTGGGTTCAAAATTTGGTGGAGCTAAGCggg                                      <  1:20657427/65‑1 (MQ=25)

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

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

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

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

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

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

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

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

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

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

                                    ccccACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAAc                               >  1:15514427/1‑67 (MQ=255)

                                     cccACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGaaac                              >  1:23890415/1‑65 (MQ=255)

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

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

                                     cccACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAAcc                              <  1:14425318/67‑1 (MQ=255)

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

                                      ccACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGc                            <  1:38941257/68‑1 (MQ=255)

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

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

                                         cGATTCGTTTACGTACCAAGGTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGa                         <  1:4160211/68‑1 (MQ=255)

                                          gATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGa                         >  1:22700645/1‑67 (MQ=255)

                                          gATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGa                         >  1:34302316/1‑67 (MQ=255)

                                          gATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGa                         >  1:23041694/1‑67 (MQ=255)

                                           aagCGTTTACGCACCAAGGTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGAcc                       <  1:19864419/66‑1 (MQ=21)

                                           aTGCCTTTACGTACCAAGTTTGGTGGGTGCTAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGAcc                       <  1:31628557/68‑1 (MQ=21)

                                            tccGTTTACGGACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGAcc                       <  1:9416271/65‑1 (MQ=255)

                                            tGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACct                      <  1:37809998/68‑1 (MQ=255)

                                            tGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGATAAGCGGGATCGAACCGCTGACct                      >  1:42227872/1‑68 (MQ=255)

                                            tGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGCTGGAGCTAAGCGGGCTCGAACCGCTGAcc                       >  1:30860206/1‑67 (MQ=255)

                                             gCGTTTACGTACCAAGTTTGCTGGGTTCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACctc                     <  1:29115779/68‑1 (MQ=255)

                                              cctttACGTACCTAGCATGCTGGGTGGAAAATTTGGTGGAGCTAAGCGGGATTGAACCGCTGACctct                    <  1:18016582/66‑1 (MQ=17)

                                               ctttACCCACCAAGCTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCtt                   <  1:42748383/67‑1 (MQ=18)

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

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

                                                    cGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGcatg              <  1:32575580/68‑1 (MQ=255)

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

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

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

                                                       aCCAAGTTTGCTGGGTTCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGcatgcc            <  1:27072326/67‑1 (MQ=255)

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

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

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

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

                                                       aCCAAGTTTCCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGcatgcc            <  1:29738507/67‑1 (MQ=255)

                                                       aCCAAGTATGCTGGGTGCTAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGcatgcc            <  1:34826888/67‑1 (MQ=255)

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

                                                         cAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGGTCGAACCGCTGACCTCTTGcatgccat          >  1:6072846/1‑67 (MQ=255)

                                                           aGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCa       >  1:42883964/1‑68 (MQ=255)

                                                           aGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCa       >  1:34119196/1‑68 (MQ=255)

                                                            gTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCaa      >  1:28512115/1‑68 (MQ=255)

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

                                                              ttGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAgc    <  1:28010619/68‑1 (MQ=255)

                                                                 cTGGGTTCAAAATTTGGTGGAGCTAAGCGGGAGCGAACCGCTGACCTCTTGCATGCCATGCAAgc    <  1:15255487/65‑1 (MQ=255)

GCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGC  >  CP009273/2511556‑2511687

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