BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	692,249	T→A	8.7%	noncoding (38/75 nt)	glnW ←	tRNA‑Gln

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	692,249	0	T	A	8.7%	620.5 / 13.0	222	noncoding (38/75 nt)	glnW	tRNA‑Gln
Reads supporting (aligned to +/- strand): ref base T (114/87); new base A (2/19); total (116/106)
Fisher's exact test for biased strand distribution p-value = 3.45e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 225 total aligned reads displayed.

TGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCGAAG  >  CP009273/692183‑692316

tGAGCTACGTAGCCAGCTTGTTTATTCGAAGCCTGGTGAACCGGGAAGCGAACCAGGGATTGACTGTa                                                                    >  1:24272723/1‑68 (MQ=255)

tGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTa                                                                    >  1:725858/1‑68 (MQ=255)

tGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTa                                                                    >  1:33540112/1‑68 (MQ=255)

tGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTa                                                                    >  1:13819448/1‑68 (MQ=255)

 taGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTAt                                                                   <  1:30043599/67‑1 (MQ=255)

 gAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTAt                                                                   <  1:17492423/68‑1 (MQ=255)

 gAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTAt                                                                   <  1:26042286/68‑1 (MQ=255)

 gAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTAt                                                                   <  1:19698567/68‑1 (MQ=255)

 gAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTAt                                                                   <  1:6855641/68‑1 (MQ=255)

 gAGCTACGAAGCCAGATCTTTTCTACGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTAt                                                                   <  1:26365369/68‑1 (MQ=17)

   gCTACGTAGCCAGATTGTTTCTTCGATGGCTTGGTTACCTGGAATCGAACCAGGGAATGGAGGTAtaa                                                                 >  1:30761217/1‑66 (MQ=21)

   gCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCa                                                                 >  1:22801714/1‑68 (MQ=255)

   gCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCa                                                                 >  1:19769634/1‑68 (MQ=255)

   gCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCa                                                                 >  1:13518957/1‑68 (MQ=255)

    cTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCTATCAAGGAAAAGCCGCTAACaa                                                                >  1:28121686/1‑68 (MQ=21)

    cTACGTAGCCAGATTGTTTATTCGATGGCTGGGGTACCTGTATTCGAGCCATTGAATGCCTTTATCaa                                                                >  1:20770937/1‑68 (MQ=21)

     tACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCaaa                                                               >  1:32030699/1‑68 (MQ=255)

     tACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCa                                                                 <  1:6247759/66‑1 (MQ=255)

     tACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCGGGATTCGACCCAGGGAATGCCGGTATCaaa                                                               >  1:32533360/1‑68 (MQ=255)

     tACGTAGCCAGATTGTTTCTACGATTGCTGGGGTACCTGGATTCGAACCAGTGAATGCAGGTATCaaa                                                               >  1:2679823/1‑68 (MQ=18)

        gTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCaaaaa                                                             >  1:1569293/1‑67 (MQ=255)

        gTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCaaaaa                                                             >  1:30618900/1‑67 (MQ=255)

        gTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCaaaaa                                                             >  1:27279181/1‑67 (MQ=255)

        gTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCaaaaa                                                             >  1:34308116/1‑67 (MQ=255)

        gTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCaaaaa                                                             >  1:18667509/1‑67 (MQ=255)

        gTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCaaaaa                                                             >  1:17803150/1‑67 (MQ=255)

        gTAGCCAGATTGTTTATTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCaaaaa                                                             >  1:7736856/1‑67 (MQ=255)

          aGCCAGATATCTTCTTAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAAcc                                                           <  1:3776132/67‑1 (MQ=11)

           gCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCgg                                                         >  1:4909937/1‑68 (MQ=255)

            ccAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGt                                                        >  1:25688993/1‑68 (MQ=255)

             cAGATTGTTTCTTCGATGGCTGGGGTACCTGTATTCGAACCAGGGAATGCCGGTATCAAAAACCGGt                                                        >  1:5693733/1‑67 (MQ=255)

             cAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCaaaaa                                                             >  1:15178531/1‑62 (MQ=255)

             cAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGt                                                        <  1:29907780/67‑1 (MQ=255)

             cAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGt                                                        >  1:1280596/1‑67 (MQ=255)

             cAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGt                                                        >  1:30314092/1‑67 (MQ=255)

             cAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGt                                                        >  1:21390508/1‑67 (MQ=255)

             cAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGt                                                        >  1:17831982/1‑67 (MQ=255)

             cAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGt                                                        >  1:14418779/1‑67 (MQ=255)

             cAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGt                                                        >  1:10372631/1‑67 (MQ=255)

             cAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGACGGTAGCAAAAACCGGt                                                        >  1:31137614/1‑67 (MQ=255)

               ttttGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGcc                                                     <  1:20903675/66‑1 (MQ=255)

               gtttGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGc                                                      <  1:2264275/65‑1 (MQ=255)

                aTTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGt                                                        <  1:33847606/64‑1 (MQ=255)

                aTTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCt                                                    >  1:15430486/1‑68 (MQ=255)

                aTTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCt                                                    >  1:21613886/1‑68 (MQ=255)

                 ttGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGcc                                                     <  1:24116858/66‑1 (MQ=255)

                 ttGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGcc                                                     <  1:13978931/66‑1 (MQ=255)

                 ttGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGcc                                                     <  1:15545992/66‑1 (MQ=255)

                 ttGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCtt                                                   <  1:36369237/68‑1 (MQ=255)

                  tgattcttcGATGGCTGGGGTACCTGCATTCGAACCAGGGAATGCCGGAATCAAAAACCGCTGCCTTa                                                  >  1:36071198/4‑68 (MQ=14)

                  tGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCtt                                                   >  1:35632802/1‑67 (MQ=255)

                   gtgtcttcGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTAc                                                 <  1:19321674/65‑1 (MQ=255)

                    tttCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTAcc                                                <  1:10141091/68‑1 (MQ=255)

                     ttcttcGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTAcc                                                <  1:14653327/67‑1 (MQ=40)

                      tcttcGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTAcc                                                >  1:36547731/1‑66 (MQ=39)

                      tcttcGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTAcc                                                >  1:32803246/1‑66 (MQ=39)

                      tcttcGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGc                                              <  1:12356155/68‑1 (MQ=39)

                        ttcGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCtt                                            >  1:22942032/1‑68 (MQ=35)

                         tcGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTg                                           <  1:63415/68‑1 (MQ=33)

                          cGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTgg                                          <  1:14534227/68‑1 (MQ=31)

                          cGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTg                                           >  1:3983225/1‑67 (MQ=31)

                          cGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTg                                           >  1:28464756/1‑67 (MQ=31)

                                         gaggATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAAc                           <  1:35454037/66‑1 (MQ=255)

                                         gaggATTCGAACCACGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAAc                           <  1:12184391/66‑1 (MQ=255)

                                         cTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATac                            >  1:2030860/1‑66 (MQ=25)

                                          tGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATaggc                          >  1:23407489/1‑65 (MQ=25)

                                          tGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATATcc                          >  1:10118314/1‑68 (MQ=25)

                                             aTTCGAACCTGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCgg                        <  1:29250339/67‑1 (MQ=255)

                                              ttCGAACCTGGGAATGCCGGTATCATAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGc                      <  1:13336092/68‑1 (MQ=255)

                                              ttCGAACCTGGGAATGCCGGTATCAAAATCCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGc                      <  1:14849254/68‑1 (MQ=255)

                                              ttCGAACCTGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGc                      <  1:33273627/68‑1 (MQ=255)

                                              ttCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCGTTACCGCTTGGCGATACCCCAATAACCGGGc                      <  1:6947722/68‑1 (MQ=255)

                                              ttCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGCTACCCCAATAACCGGGc                      >  1:20528789/1‑68 (MQ=255)

                                              ttCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGc                      <  1:34952566/68‑1 (MQ=255)

                                              ttCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGc                      <  1:25116140/68‑1 (MQ=255)

                                              ttCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGc                      <  1:9400130/68‑1 (MQ=255)

                                              ttCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAAGAGCCGGGc                      >  1:4329137/1‑68 (MQ=255)

                                               tCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCg                     >  1:8888149/1‑68 (MQ=255)

                                                 gAACCTGGGAATGCCGGAATCACAATCCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGt                   <  1:37016/68‑1 (MQ=14)

                                                 gAACCTCGGAATGCCGGTATCATAATCCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGt                   <  1:22157870/68‑1 (MQ=14)

                                                 gAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGt                   <  1:6107561/68‑1 (MQ=255)

                                                 gAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGt                   <  1:21166275/68‑1 (MQ=255)

                                                   aCCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTgg                 >  1:4676652/1‑67 (MQ=255)

                                                     ctcggAATGCCGGTATCAAAATCCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAAc               <  1:32517123/65‑1 (MQ=255)

                                                      tcggAATGCCGGAATCATAATCCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAAcc              <  1:30979275/66‑1 (MQ=14)

                                                       cggAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCg             <  1:638325/67‑1 (MQ=255)

                                                       cggAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCg             <  1:8631886/67‑1 (MQ=14)

                                                        ggAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAAAAAACGGGCGGTGaaac              >  1:31360174/1‑64 (MQ=17)

                                                         gAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAAcc              <  1:545539/65‑1 (MQ=255)

                                                         gAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAAcc              <  1:8427200/65‑1 (MQ=255)

                                                            tGCCGGTATCACAAACCTGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTAc        <  1:31977997/68‑1 (MQ=255)

                                                            tGCCGGAATCATAATCCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTAc        <  1:35280169/68‑1 (MQ=18)

                                                             gCCGGTATCACAATCCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACt       <  1:36514434/68‑1 (MQ=255)

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

                                                             gCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAAGAACCGGGCGGTGACCCGATTagg       >  1:26594749/1‑66 (MQ=17)

                                                             gCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACt       <  1:195639/68‑1 (MQ=18)

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

                                                              ccGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTTAACCGCTTACTc      >  1:11124297/1‑68 (MQ=255)

                                                               cGGAATCATAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTc      <  1:34120688/67‑1 (MQ=255)

                                                                  tATCATAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCgaa   <  1:29602758/67‑1 (MQ=255)

TGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCGAAG  >  CP009273/692183‑692316

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