BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	2,511,621	C→G	6.5%	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.5%	1213.5 / 9.7	413	intergenic (‑26/+14)	alaX/alaW	tRNA‑Ala/tRNA‑Ala
Reads supporting (aligned to +/- strand): ref base C (260/122); new base G (23/8); total (283/130)
Fisher's exact test for biased strand distribution p-value = 5.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 418 total aligned reads displayed.

TTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGC  >  CP009273/2511554‑2511687

   catgccatgcAAGCGCTCTCCCAGCTGAGCTATAGCCCCGCGAGGGGTTTAGGTACCAAGTTTGCggg                                                                 >  1:15664320/1‑68 (MQ=17)

   catgccatgcAAGCGCTCTCCCAGCTGAGCTATAGCCCCGCGAGGCGTTTACGTACCAAGTTTGCTgg                                                                 >  1:20448592/1‑68 (MQ=255)

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

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

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

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

      gccttgcAAGCGCTCTCCCAGCTGAGCTATAGCCCCGGGATGGGTTTACGTTCCAAGTTTGCTGGGt                                                               >  1:1441776/1‑67 (MQ=18)

      gccatgcAAGCGCTCTCCCAGCTGAGCTATAGCCCCCCGAGGCGTTTAGGTACCAAGTTTGCTGGGt                                                               >  1:10155026/1‑67 (MQ=21)

       ccatgcAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGc                                                             >  1:4938080/1‑68 (MQ=255)

            cAAGCGCTCTCCCAGCTGAGCTATAGCCCCGCGATGCTTATAGGTACCAGGTTTCGGGGGTGCAAAAt                                                        >  1:15048353/1‑68 (MQ=17)

            cAAGCGCTCTCCCAGCTGAGCTATAGCCCCGCGATGCGTTTACGTTCCAAGTGTGCTGGGTGCAAAAt                                                        >  1:2372928/1‑68 (MQ=25)

            cAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAGGTTTGCTGGGTGCAAAAt                                                        >  1:5480619/1‑68 (MQ=255)

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

              aGCGCTCTCCCAGCTGAGCTATAGCCCCCGGAGGGGTTTAGGTACCAATTTTCGTGGGGCCAAAAttt                                                      >  1:25061396/1‑68 (MQ=12)

              aGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAAttt                                                      >  1:22697064/1‑68 (MQ=255)

              aGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAAttt                                                      >  1:18230011/1‑68 (MQ=255)

                cgcTCTCCCAGCTGAGCTATAGCCCCGCGATGAGTTTAAGTTCCAGGTTTGCTGGGTGCAAAATTtgg                                                    >  1:8475626/1‑68 (MQ=25)

                cgcTCTCCCAGCTGAGCTATAGCCCCCCGATGGGTTTAAGTACCAAGTTTGGTGGGTGCAAAATTtgg                                                    >  1:19007701/1‑68 (MQ=37)

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

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

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

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

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

                    ctcACAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGa                                                <  1:5405251/68‑1 (MQ=255)

                     tcCCATCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAg                                               <  1:22550858/68‑1 (MQ=255)

                      cccAGCTGAGCTATAGCCCCGCGATGCGTTTACGTACCAATTTTCCTGGGTGCAAAAATTGGTGGAga                                              >  1:3831293/1‑67 (MQ=37)

                       ccAGCTGAGCTATAGCCCCGCGAGGGGTTTACGTTCGAATTTACGTGGGGGCAAATTTTTGTGGAGCt                                             >  1:25623750/1‑68 (MQ=12)

                       ccAGCTGAGCTATAGCCCCCCGATGCTTTTAGGTACCATGTTAGCTGGGGGCAAAATTTGGTGGAGCt                                             >  1:17245033/1‑68 (MQ=37)

                       ccAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCt                                             >  1:9563894/1‑68 (MQ=255)

                         aGCTGAGCTATAGCCCCGGGATGGGTTTAGGTTCCAATTTTGGTGGGTGCAAAATTTGTTGGAGCTaa                                           >  1:4349522/1‑68 (MQ=255)

                         aGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCaaaa                                                         >  1:26418427/1‑54 (MQ=255)

                         aGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCTAAATTTGGTGGAGCTaa                                           >  1:26172508/1‑68 (MQ=255)

                         aGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAAATTGGTGGAGCTaa                                           >  1:3477610/1‑68 (MQ=255)

                         aGCTGAGCTATAGCCCCACGAAGCGTTTACGTACCAAGTTTGCTGGGTGAAAATTTTGGTGGCGCTaa                                           >  1:11638115/1‑68 (MQ=255)

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

                           cTGAGCTATAGCCCCTCGAGGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGATaata                                         >  1:28678148/1‑66 (MQ=255)

                           cTGAGCTATAGCCCCGCGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGc                                         >  1:10707891/1‑68 (MQ=255)

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

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

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

                           cTGAGCTATAGCCCCACGATGCGTTTAAGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGc                                         >  1:309667/1‑68 (MQ=255)

                           cTGAGCTATAGCCCCACGAAGCGTTTACGTACCAAGTTTGCTGGGGGCTAAATTTGGTGGAGCTAAGc                                         >  1:14264186/1‑68 (MQ=39)

                            tgcgCTATAGCCACACCATGCTTATACGTCCCAAGTTTGCTGGGTTCTAAATTTGGTGGAGCTAAGc                                         <  1:15514974/64‑1 (MQ=21)

                            tGAGCTATAGCCCCACTATGCTTTTACGTACCAAGTTTGGTGGGTGCAAAATTTGGTGGAGCTaa                                           <  1:179862/65‑1 (MQ=255)

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

                            tGAGCTAAAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTaa                                           <  1:11898511/65‑1 (MQ=255)

                            tGAGCAAAACCCCCACGCTGCGTTTACGAACCAAGCTTGCTGGGTTCTAAATTTGGTGGAGCTaa                                           <  1:10987245/65‑1 (MQ=255)

                                 tataGCCCCACGATGCGTTTACGTACCAAGTTTGCAGGGTGCAAAATTTGGTGGAGCTAAGCGCGATc                                   >  1:18733929/1‑68 (MQ=255)

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

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

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

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

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

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

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

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

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

                                       cccACTCTGCGTTTACGTACGTAGTTTGCTGGGTGCTAAATTTGGTGGAGCTAAGCGGGATCGAAcc                              <  1:26982227/67‑1 (MQ=21)

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

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

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

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

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

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

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

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

                                       cccACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCg                             <  1:21988235/68‑1 (MQ=255)

                                        ccaccaTCCGATAACGTCCCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGc                            <  1:18421070/68‑1 (MQ=21)

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

                                         cACGAATCGATTACGCCCCAAGCATGCTGGGTTCTAAATTTGGTGGAGCTAAGCGGGATCGAACCGCt                           <  1:14818714/68‑1 (MQ=9)

                                          aCGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTg                          >  1:6193395/1‑68 (MQ=255)

                                           ctctGCGTTTACGTACAATGTTTGGTGGGTTCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGa                         <  1:21388996/65‑1 (MQ=18)

                                           cGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGa                         <  1:20686502/68‑1 (MQ=255)

                                           cGATGCGAATACGAACCAAGCTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGa                         <  1:7332263/68‑1 (MQ=21)

                                           cGATCCGTTAACGTCCCTAGCTTGCTGGGTGGTAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGa                         <  1:19845684/68‑1 (MQ=17)

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

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

                                               gCGTTTAGGTTCCAAGCTTCGTGGGTTCAAAATTTGGTGTAGCTAAGCGGGATCGAACCGCTGACctc                     <  1:131811/68‑1 (MQ=17)

                                               gCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGAcc                       >  1:28999372/1‑66 (MQ=255)

                                                 gattACGTACCAAGTTTGCTGGGTTCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACctct                    <  1:28382207/65‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

                                                         cccAAGTATGCTGGGTGCTAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGcatgcc            <  1:28181604/66‑1 (MQ=17)

                                                         aCCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGCCCTCTTGcatgcca           >  1:27128323/1‑68 (MQ=255)

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

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

                                                          ccAAGTTTCCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGcatgccat          <  1:16056306/68‑1 (MQ=255)

                                                          ccAAGCTTGCAGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGcatgccat          <  1:12374365/68‑1 (MQ=255)

                                                           cAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGcatgccatg         <  1:21860384/68‑1 (MQ=255)

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

                                                                ttGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAg     <  1:17797614/67‑1 (MQ=255)

                                                                 tGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAgc    >  1:27888533/1‑67 (MQ=255)

TTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGC  >  CP009273/2511554‑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: `‑`