BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	692,307	T→G	5.2%	intergenic (‑21/+14)	glnW ← / ← glnU	tRNA‑Gln/tRNA‑Gln

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	692,307	0	T	G	5.2%	1430.6 / 6.5	463	intergenic (‑21/+14)	glnW/glnU	tRNA‑Gln/tRNA‑Gln
Reads supporting (aligned to +/- strand): ref base T (261/168); new base G (17/15); total (280/183)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 471 total aligned reads displayed.

GAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGC  >  CP009273/692240‑692373

  aTGCCGGAATCACAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTa                                                                  <  1:48704763/68‑1 (MQ=255)

   tGCCGGAATCATAAACCGTTGCCTTACCTCTTGGCGATACCCCAATAACCTGGCGGTGAACCGCTTAc                                                                 <  1:33085739/68‑1 (MQ=21)

   tGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTAc                                                                 <  1:35240785/68‑1 (MQ=17)

   tGCCGGAATCACAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTAc                                                                 <  1:22584964/68‑1 (MQ=255)

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

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

    gCCGGAATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACt                                                                <  1:2905026/68‑1 (MQ=255)

    gCCGGAATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACt                                                                <  1:47492992/68‑1 (MQ=255)

     ccGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAAAAAGCGGTCGGTGACCCGATTACTc                                                               >  1:49561088/1‑68 (MQ=17)

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

      cGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCg                                                              <  1:35006510/68‑1 (MQ=255)

       ggaaTCAGAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCga                                                             <  1:26129367/65‑1 (MQ=255)

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

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

            cAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCGaagaag                                                        <  1:27801131/68‑1 (MQ=255)

               aatccGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCGAAGAAGATg                                                     <  1:39797872/65‑1 (MQ=255)

               aatccGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCGAAGAAGATg                                                     <  1:43060996/65‑1 (MQ=255)

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

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

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

                 tccGGTGCCTTACCGCTTGGCGATACCCCAATAAAAGGGCGGTGAAGCGGGTGCTCGAAGAAGATg                                                     >  1:3748046/2‑66 (MQ=18)

                 tccGGTGCCTTACCGCTTGGCGATACCCCAAGAACCTGGCGGTGAACGGCGTACTCGAAGAAGATg                                                     >  1:31394702/2‑66 (MQ=21)

                 aCCGGTGCCTTACCGCTTGGCGATACCCCAATATCCGGGCGGTGAACCGCTTACTCGAAGAAGATg                                                     >  1:12730080/1‑66 (MQ=255)

                 aCCGGTGCCTTACCGCTTGGCGATACCCCAATAGCCGGGCGGTGACCCGATTAATCGAAGACCATGGc                                                   >  1:40934059/1‑68 (MQ=21)

                 aCCGGTGCCTTACCGCTTGGCGATACCCCAATAGCCGGGCGGTGAACCGATTACGCGAAGAAGATCGc                                                   >  1:24470786/1‑68 (MQ=25)

                 aCCGGTGCCTTACCGCTTGGCGATACCCCAATAGCACGGCGGTGAACCGCTAACTCGAAGACGAtc                                                     >  1:24889706/1‑65 (MQ=21)

                 aCCGGTGCCTTACCGCTTGGCGATACCCCAATACCCGGGCGGTGAACCGCTTACTCGTAGCAGGTg                                                     >  1:19984597/1‑66 (MQ=25)

                 aCCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTAACTCGAAGAAGATg                                                     >  1:41634206/1‑66 (MQ=255)

                  ccGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCGAAGAAGATGGCt                                                  <  1:23709023/68‑1 (MQ=255)

                   cGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTg                                                 >  1:30866801/1‑68 (MQ=255)

                   cGGTGCCTTACCGCTTGGCGATACCCCAAGAGACGGGCGGTGACCGGGTTAGTCGCAGATGATggagg                                                 >  1:20427496/1‑65 (MQ=9)

                    ggTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTgg                                                <  1:48332748/68‑1 (MQ=255)

                      tGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTggg                                               <  1:18381736/67‑1 (MQ=255)

                      tGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTggg                                               <  1:5618098/67‑1 (MQ=255)

                         cTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCgaaga                                                          <  1:16976907/53‑1 (MQ=255)

                         cTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTAc                                           <  1:22167405/68‑1 (MQ=255)

                         cTTACCGCTTGGCGATACCCCAATAAACGGGAGCTGAACCGGTTACTCGCAGAAGATGGCTGGGGTa                                            >  1:30329328/1‑67 (MQ=25)

                         cTTACCGCTTGGCGATACCCCAAAAACCGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTa                                            <  1:44805926/67‑1 (MQ=255)

                            aCCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTAAGCGAAGAAGATGGCTGGGGTACCTg                                        >  1:28654660/1‑68 (MQ=255)

                              cgcgTGGCGATACCCCAACTACCTGTCGGTTAACCGCGTACTCGAATAATATGGCTGGGGTACCTGGa                                      <  1:16721549/68‑1 (MQ=255)

                              cGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGa                                      <  1:42050198/68‑1 (MQ=255)

                              cGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTGCGGTACCTGGa                                      >  1:14346021/1‑68 (MQ=255)

                              cGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCGAAGAAGATCGCTGGGGTACCTGGa                                      >  1:50952508/1‑68 (MQ=255)

                              cGCTTGGCGATACCCCAAGAGCCGGGCGGTGAACCGCGTACTCGAAGAAGATTGCGTGGTTACCTCGa                                      >  1:6549827/1‑68 (MQ=255)

                               gCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGAt                                     >  1:7427219/1‑68 (MQ=255)

                               gCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGAt                                     >  1:36095819/1‑68 (MQ=255)

                               gCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGAt                                     >  1:5940096/1‑68 (MQ=255)

                                 ttGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGCTTc                                   >  1:8097449/1‑68 (MQ=255)

                                   ggCGATACCGCAATAACCGGGCGGTGAACCGCTTACTGGTAGAAGATGGCTGGGGTACCTGGATTCGa                                 <  1:27571883/68‑1 (MQ=38)

                                   ggCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGa                                 <  1:25762947/68‑1 (MQ=255)

                                   ggCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGa                                 >  1:6961205/1‑68 (MQ=255)

                                    gCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGATACGa                                 >  1:32321170/1‑67 (MQ=255)

                                     cGATCCCCCAATAACCGGGCGGTGCACCGCGTTCTCTATGAAGATGGCTGGGGTACCTGGATTCGaa                                <  1:19645553/67‑1 (MQ=25)

                                     cGATACCCCAATACCCGGGCGTTGAACCGCTTACTCGAATAAGATGGCTGGGGTACCTGGATTCGaa                                <  1:17590479/67‑1 (MQ=37)

                                     cGATACCCCAATAACCGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGaa                                <  1:51313298/67‑1 (MQ=255)

                                     cGATACCCCAATAACCGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGaa                                <  1:34067619/67‑1 (MQ=255)

                                     cGATACCCCAATAACCGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGaa                                <  1:24584983/67‑1 (MQ=255)

                                     cGATACCAAAATAACCGCGCTTCGAACCCCGTACTCTACGAAGATGGCTGGGGTACCTGGATTCGaa                                <  1:33386485/67‑1 (MQ=12)

                                        tACCCAAATCACCCGGCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAg                            <  1:38580603/68‑1 (MQ=37)

                                         aCCCCAATAACCGGGCGGTGAACCGCTTTCTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAgg                           <  1:36585784/68‑1 (MQ=255)

                                         aCCCCAATAACCGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAgg                           >  1:4636062/1‑68 (MQ=255)

                                         aCCCCAATAACCGGGCGCTGAACCGCTTTCTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAgg                           <  1:6340745/68‑1 (MQ=255)

                                          ccccAATAACCGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAggg                          >  1:37175971/1‑68 (MQ=255)

                                          ccccAATAACCGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAggg                          >  1:45607843/1‑68 (MQ=255)

                                           cccAATAACCGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGa                         >  1:47833009/1‑68 (MQ=255)

                                           cccAATAACCGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGa                         >  1:13678697/1‑68 (MQ=255)

                                                 aaCCGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGCTCCAGGGAAGGCCg                   >  1:16176661/1‑68 (MQ=37)

                                                 aaCCGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCg                   >  1:45931817/1‑68 (MQ=255)

                                                 aaCCGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCg                   >  1:8253397/1‑68 (MQ=255)

                                                 aaCCGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCg                   >  1:8922808/1‑68 (MQ=255)

                                                 aaCCGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCg                   >  1:10480268/1‑68 (MQ=255)

                                                 aaCCGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCg                   >  1:40126557/1‑68 (MQ=255)

                                                 aaCCGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCg                   >  1:3783514/1‑68 (MQ=255)

                                                 aaCCGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCg                   >  1:35742831/1‑68 (MQ=255)

                                                 aaCCGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCg                   >  1:12232664/1‑68 (MQ=255)

                                                 aaCCGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCg                   >  1:29314765/1‑68 (MQ=255)

                                                 aaCCGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCg                   >  1:25452534/1‑68 (MQ=255)

                                                 aaCCGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCg                   >  1:17269036/1‑68 (MQ=255)

                                                 aaCCGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCg                   >  1:18903055/1‑68 (MQ=255)

                                                 aaCCGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCg                   >  1:42809348/1‑68 (MQ=255)

                                                 aaCCGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCg                   >  1:18598577/1‑68 (MQ=255)

                                                  aCCGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCgg                  >  1:9192841/1‑68 (MQ=255)

                                                  aCCGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCgg                  >  1:46873858/1‑68 (MQ=255)

                                                   ccGGGCTGCGAACCGCGTAATCTAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGt                 <  1:26524804/68‑1 (MQ=17)

                                                   ccGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGt                 >  1:50496725/1‑68 (MQ=255)

                                                   ccGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGt                 >  1:11205683/1‑68 (MQ=255)

                                                   ccGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGt                 >  1:11634678/1‑68 (MQ=255)

                                                    cGGGCGGTGCACCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTa                >  1:38165130/1‑68 (MQ=255)

                                                    cGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTa                >  1:15664234/1‑68 (MQ=255)

                                                    cGGGCGGTGAACCGCTTACTAGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTa                >  1:43488645/1‑68 (MQ=255)

                                                     gggCGGTGTCCCGCTTACTAGAAGAACATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTAt               <  1:30701687/68‑1 (MQ=18)

                                                      ggCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATc              >  1:4128585/1‑68 (MQ=255)

                                                       gCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCa             >  1:21328552/1‑68 (MQ=255)

                                                             aaCCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAAcc       >  1:14676979/1‑68 (MQ=255)

                                                             aaCCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAAcc       >  1:7695646/1‑68 (MQ=255)

                                                              aCCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCg      >  1:1500152/1‑68 (MQ=255)

                                                              aCCGCGTACTTGAATAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAAc        <  1:13911580/66‑1 (MQ=14)

                                                                 gCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGt    >  1:12965943/1‑67 (MQ=255)

                                                                   ttACTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTg   >  1:8516997/1‑66 (MQ=255)

                                                                   ttACTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTg   >  1:10722534/1‑66 (MQ=255)

GAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGC  >  CP009273/692240‑692373

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