BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	2,811,524	G→T	5.3%	intergenic (‑30/+33)	argZ ← / ← argY	tRNA‑Arg/tRNA‑Arg

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	2,811,524	0	G	T	5.3%	2268.0 / 13.4	697	intergenic (‑30/+33)	argZ/argY	tRNA‑Arg/tRNA‑Arg
Reads supporting (aligned to +/- strand): ref base G (275/373); new base T (42/4); total (319/378)
Fisher's exact test for biased strand distribution p-value = 2.81e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 710 total aligned reads displayed.

TCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCT  >  CP009273/2811457‑2811591

    aGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTGATa                                                                  <  1:21142000/67‑1 (MQ=255)

     gCCGAGTACTCTATCCATCTGAGATAAGGAAGAATCCGCAAAT‑TAATTTTAAAGCTGATATTGATaa                                                                 >  1:10881650/1‑66 (MQ=12)

     gCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTGa                                                                    >  1:7607771/1‑64 (MQ=255)

     gCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTGATAca                                                                >  1:11162485/1‑67 (MQ=255)

       cGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTGATAcc                                                                <  1:39539238/66‑1 (MQ=255)

       cGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTGATAcc                                                                <  1:3457877/66‑1 (MQ=255)

       cGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTGATAcc                                                                >  1:10141686/1‑66 (MQ=255)

       cGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTGATAcc                                                                <  1:27429989/66‑1 (MQ=255)

       cGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTGATACCg                                                               <  1:28020196/67‑1 (MQ=255)

        gAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTGATACCGct                                                             <  1:40219711/68‑1 (MQ=255)

        gAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTGATACCGct                                                             <  1:29680175/68‑1 (MQ=255)

          gTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTGATACCGctac                                                           <  1:22943795/68‑1 (MQ=255)

          gTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTGATACCGcta                                                            <  1:42960510/67‑1 (MQ=255)

          gTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTGATACCGcta                                                            <  1:36697886/67‑1 (MQ=255)

          gTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTGATACCGct                                                             >  1:29389835/1‑66 (MQ=255)

            acgctATCCAGCGGAGCTACGGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTGATACCGctac                                                           <  1:26285109/63‑1 (MQ=255)

            aatctATCCAGATGTGCTTCTGTTGGTTTGGGATAT‑TTTTTTTAGTGATTATATTGATACCGctac                                                           <  1:16442079/64‑1 (MQ=255)

            aCTCTATCCAGCTGAGCTACGGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTGATACCGctac                                                           <  1:36213809/66‑1 (MQ=255)

            aCTCTATCCAGCTGAGCTACGGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTGATACCGctac                                                           <  1:12710742/66‑1 (MQ=255)

            aCTCTATCCAGCTGAGCTACGGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTGATACCGctac                                                           <  1:26909888/66‑1 (MQ=255)

            aCTCTATCCAGCTGAGCTACGGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTGATACCGctac                                                           <  1:34915660/66‑1 (MQ=255)

            aCTCTATCCAGCTGAGCTACGGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTGATACCGctac                                                           <  1:3421809/66‑1 (MQ=255)

            aCTCTATCCAGCTGAGCTACGGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTGATACCGctac                                                           <  1:6338893/66‑1 (MQ=255)

            aCTCTATCCAGCTGAGCTACGGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTGATACCGctac                                                           <  1:7984569/66‑1 (MQ=255)

            aCTCTATCCAGCTGAGCTACGGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTGATACCGctac                                                           <  1:14025511/66‑1 (MQ=255)

             ctctATCCAGCTGAGCTACGGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTGATACCGCTACTaa                                                        >  1:6024405/1‑68 (MQ=255)

              tctATCCAGCTGAGCTACGGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTTATACCGCTAATaaa                                                       >  1:24433431/1‑68 (MQ=255)

              tctATCCAGCTGAGCTACGGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTGATACCGCTACTaaa                                                       <  1:15417608/68‑1 (MQ=255)

              tctATCCAGCTGAGCTACGGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTGATACCGCTACTaaa                                                       >  1:40912605/1‑68 (MQ=255)

               ctATCCAGCTGAGCTACGGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTGATACCGCTACTAAAg                                                      <  1:31671537/68‑1 (MQ=255)

                tATCCAGCTGAGCTACGGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTGATACCGCTACTAAAGc                                                     <  1:32422348/68‑1 (MQ=255)

                tATCCAGCTGAGCTACGGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTGATACCGCTACTAAAGc                                                     <  1:37028350/68‑1 (MQ=255)

                tATCCAGCTGAGCTACGGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTGATACCGCTACTAAAGc                                                     <  1:26575506/68‑1 (MQ=255)

                 aTCCAGCTGAGCTACGGATGCATCGGGATAT‑TTATATTAGTGATGATATTGATAACCCGAATAAAGcc                                                    >  1:23498806/1‑68 (MQ=17)

                 aTCCAGCTGAGCTACGGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTGATACCGCTACTAAAGcc                                                    >  1:19718870/1‑68 (MQ=255)

                 aTCCAGCTGAGCTACGGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTGATACCGCTACTAAAGc                                                     <  1:24060812/67‑1 (MQ=255)

                  tCCAGCTGAGCTACGGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTGATACCGCTACTAAAGCCa                                                   <  1:14993006/68‑1 (MQ=255)

                   ccAGCTGAGCTACGGATGCATCTGGATTT‑TAATGTTCCTGGTGATAATGATACCGCTAGTAAGGCCag                                                  >  1:8361260/1‑67 (MQ=12)

                   ccAGCTGAGCTACGGATGCATCGGGAAATGTTTTTTTAGTGATGATATTGATACGGATAATAAAGcat                                                   >  1:31336602/1‑66 (MQ=2)

                   ccAGCTGAGCTACGGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTTATACCGCTACTAAAGCCat                                                  >  1:16636250/1‑68 (MQ=255)

                   ccAGCTGAGCTACGGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTGATACCGCTACTAAAGCCat                                                  >  1:38074665/1‑68 (MQ=255)

                   ccAGCTGAGCTACGGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTGATACCGCTACTAAAGCCat                                                  >  1:22289226/1‑68 (MQ=255)

                   ccAGCTGAGCTACGGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTGATACCGCTACTAAAGCCat                                                  >  1:6895076/1‑68 (MQ=255)

                   ccAGCTGAGCTACGGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTGATACCGCTAATAAAGCCat                                                  >  1:20867905/1‑68 (MQ=255)

                    cAGCTGAGCTACGGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTGATACCGCTAATAAAGCCata                                                 >  1:42156242/1‑68 (MQ=255)

                    cAGCTGAGCTACAGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTGATACCGCTACTAAAGCCata                                                 <  1:33870670/68‑1 (MQ=255)

                      gCTGAGCTACGGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTTATACCGCTACTAAAGCCATATc                                               >  1:3735397/1‑68 (MQ=255)

                      gCTGAGCTACGGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTTATACCGCTAATAAAGCCATATc                                               >  1:5614127/1‑68 (MQ=255)

                      gCTGAGCTACGGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTTATACCGCTAATAAAGCCATATc                                               >  1:21509462/1‑68 (MQ=255)

                      gCTGAGCTACGGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTTATAACGCTACTAAAGCCATATc                                               >  1:17720400/1‑68 (MQ=255)

                      gCTGAGCTACGGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTGATAACGCTACTAAAGCCATATc                                               >  1:8683995/1‑68 (MQ=255)

                      gCTGAGCTACGGATGCATCGGGAAAT‑TTATTTTACTGATGATATTGATAACCCTACTAAAGACAAATc                                               >  1:28885377/1‑68 (MQ=17)

                      gCTGAGCTACGGATGCATCGGGAAAT‑TTATATTACTGCTGATATAGATACGGCTACTTAAGCCAAATc                                               >  1:31991233/1‑68 (MQ=17)

                        ttaGCTACGTATGCATCGGTAAAT‑TTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCaa                                             <  1:9272853/66‑1 (MQ=17)

                        tGAGCTACGGATGCATGGGGAAAT‑TTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCaa                                             <  1:19073612/68‑1 (MQ=255)

                        tGAGCTACGGATGCATCGGTAAAT‑TTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCaa                                             <  1:22711749/68‑1 (MQ=255)

                        tGAGCTACGGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCaa                                             >  1:28328619/1‑68 (MQ=255)

                        tGAGCTACGGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCaa                                             >  1:25664842/1‑68 (MQ=255)

                        tGAGCTACGGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCaa                                             <  1:35732584/68‑1 (MQ=255)

                        tGAGCTACGGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCaa                                             <  1:41963828/68‑1 (MQ=255)

                        tGAGCTACGGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCaa                                             <  1:33518324/68‑1 (MQ=255)

                        tGAGCTACGGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCaa                                             >  1:29091736/1‑68 (MQ=255)

                        tGAGCTACGGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCaa                                             <  1:38583161/68‑1 (MQ=255)

                        tGAGCTACGGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTGATAACGCTACTAAAGCCCTATaaa                                             >  1:23154242/1‑68 (MQ=17)

                        tGAGCTACGGATGCATCGGGAAAT‑TTATTTTACTGAGGATATTGATACCCCTACTATACCCaaataa                                              >  1:20349789/1‑65 (MQ=17)

                        tGAGCTACGGATGCATAGGGAATT‑TTATATTACTGCAGATAATGATAACCCTAGTTAAACCaaataat                                             >  1:11359704/1‑65 (MQ=9)

                        tGAGCTACGGATGCATAGGGAAAT‑TTATATTACTTCTGATATTGATAACCCTAATAAAACCaaataa                                              >  1:7299951/1‑65 (MQ=12)

                         gAGCTACGGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAg                                            <  1:1796400/68‑1 (MQ=255)

                         gAGCTACGGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAg                                            <  1:24752988/68‑1 (MQ=255)

                           gCTACGGATGCATCGGTAAAT‑TTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGt                                           <  1:44432/67‑1 (MQ=255)

                           gCTACGGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGt                                           <  1:17459674/67‑1 (MQ=255)

                           gCTACGGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGt                                           <  1:43266521/67‑1 (MQ=255)

                             tACGGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTaa                                         <  1:15789157/67‑1 (MQ=255)

                              aCGGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAaga                                       <  1:28528560/68‑1 (MQ=255)

                              aCGGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAaga                                       <  1:4169551/68‑1 (MQ=255)

                              aCGGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAaga                                       <  1:4120404/68‑1 (MQ=255)

                                ggATGCATCGGGAAAT‑TTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAgaga                                     >  1:30645550/1‑68 (MQ=255)

                                ggATGCATCGGGAAAT‑TTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAgaga                                     >  1:30867249/1‑68 (MQ=255)

                                 gATGCATCGGGAAAT‑TTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGAt                                    >  1:17034710/1‑68 (MQ=255)

                                  aTGCATCGGGAAAT‑TTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATg                                   >  1:19430491/1‑68 (MQ=255)

                                   tGCTTCGGGAAAT‑TTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATgg                                  <  1:25226421/68‑1 (MQ=255)

                                   tGCATCGGGAAAT‑TTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATg                                   <  1:32894634/67‑1 (MQ=255)

                                     cATCGGGAAAT‑TTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGt                                 >  1:37610947/1‑67 (MQ=255)

                                      aTCGGGAAAT‑TTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGc                               >  1:18657941/1‑68 (MQ=255)

                                      aTCGGGAAAT‑TTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGc                               <  1:13278899/68‑1 (MQ=255)

                                         gggAAAT‑TTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATc                            <  1:11851427/68‑1 (MQ=255)

                                           gAAAT‑TTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCg                          >  1:27614709/1‑68 (MQ=255)

                                            aaaT‑TTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCg                          <  1:22071628/67‑1 (MQ=255)

                                                  tatttTACTGCTGATATTGATACCGCTACTACAGCCATATCACGTAAGAGATGGTGCATCCGggagg                     <  1:14586490/67‑1 (MQ=255)

                                                  tatttTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGggagga                    <  1:35384273/68‑1 (MQ=40)

                                                    ttttACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGAt                   >  1:13522533/1‑67 (MQ=40)

                                                         ctgctgATTATGATACAGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTCGAAc             <  1:4960317/68‑1 (MQ=11)

                                                         ctgctgATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAAc             <  1:39097533/68‑1 (MQ=255)

                                                           gctgATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAAcc            >  1:30290792/1‑67 (MQ=255)

                                                             tgATATTGATACCGCTACTAAAGCCATATCAAGTAAgaga                                     >  1:9080652/1‑40 (MQ=40)

                                                             tgATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTc          <  1:23765734/67‑1 (MQ=255)

                                                               atatTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCAt                             <  1:30030038/46‑1 (MQ=38)

                                                                  ttGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGAcc     >  1:42350004/1‑67 (MQ=255)

                                                                  ttGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGAcc     >  1:39824905/1‑67 (MQ=255)

                                                                   tGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGccc     >  1:12348869/1‑66 (MQ=255)

TCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAT‑TTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCT  >  CP009273/2811457‑2811591

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