BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	692,040	G→A	5.8%	noncoding (33/75 nt)	glnV ←	tRNA‑Gln

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	692,040	0	G	A	5.8%	1444.2 / 39.5	487	noncoding (33/75 nt)	glnV	tRNA‑Gln
Reads supporting (aligned to +/- strand): ref base G (302/142); new base A (5/23); total (314/173)
Fisher's exact test for biased strand distribution p-value = 1.76e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 494 total aligned reads displayed.

TACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCG  >  CP009273/691973‑692105

tACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAg                                                                   >  1:15208492/1‑68 (MQ=255)

tACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAg                                                                   >  1:4187793/1‑68 (MQ=255)

  cAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCACaa                                                                 >  1:15736469/1‑68 (MQ=255)

   aaGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGaa                                                                 >  1:12437965/1‑67 (MQ=255)

     ggTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCaaaaa                                                                >  1:10939948/1‑65 (MQ=255)

     ggTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATcc                                                              >  1:24340353/1‑68 (MQ=255)

     ggTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATcc                                                              >  1:15624376/1‑68 (MQ=255)

     ggTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATcc                                                              >  1:166074/1‑68 (MQ=255)

     ggTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATcc                                                              >  1:18420664/1‑68 (MQ=255)

     ggTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATcc                                                              >  1:19264212/1‑68 (MQ=255)

     ggTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATcc                                                              >  1:13329845/1‑68 (MQ=255)

     ggTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATcc                                                              >  1:24688601/1‑68 (MQ=255)

     ggTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATcc                                                              >  1:2761275/1‑68 (MQ=255)

     ggTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATcc                                                              >  1:5468901/1‑68 (MQ=255)

     ggTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATcc                                                              >  1:8305772/1‑68 (MQ=255)

     ggTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAAAATcc                                                              >  1:1806193/1‑68 (MQ=255)

       tAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCgg                                                            >  1:578435/1‑68 (MQ=255)

       tAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAAACCgg                                                            >  1:227319/1‑68 (MQ=255)

           gcgcttttttaatgTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGcc                                                        <  1:10239554‑M2/54‑1 (MQ=255)

           gagcttttttaatgTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGcc                                                        <  1:1626973‑M2/54‑1 (MQ=255)

           gagcttttttaatgTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGcc                                                        <  1:21948532‑M2/54‑1 (MQ=255)

           gagcttttttaatgTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGc                                                         <  1:17130998‑M2/53‑1 (MQ=255)

           cGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTg                                                          >  1:15530386/1‑66 (MQ=255)

           cGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTg                                                          >  1:16138761/1‑66 (MQ=255)

           cGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTg                                                          >  1:26044957/1‑66 (MQ=255)

           cGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTg                                                          >  1:26171557/1‑66 (MQ=255)

           cGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTg                                                          >  1:21108128/1‑66 (MQ=255)

           cGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTg                                                          >  1:13743563/1‑66 (MQ=255)

            tgcttttttaattTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGcc                                                        <  1:10702930‑M2/54‑1 (MQ=255)

            gTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGATTGCCGGTATCAGAATCCGCTGCCt                                                       >  1:26435066/1‑68 (MQ=17)

            gTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGcc                                                        <  1:26949527/67‑1 (MQ=11)

             gcttttttaatgTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCt                                                       <  1:6421799‑M2/55‑1 (MQ=255)

               ttttttaattTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAACCAGAATCCGGTGCCTTAc                                                    <  1:25855883‑M2/58‑1 (MQ=255)

                tGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTAcc                                                   >  1:24805929/1‑68 (MQ=255)

                 ttttaatgTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCg                                                  <  1:22477063‑M2/60‑1 (MQ=255)

                  tttaattTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGTATCAAAAACCGGTGCCTTACCGc                                                 <  1:3875681‑M2/61‑1 (MQ=255)

                  tttaatgTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAAAATCCGGTGCCTTACCGc                                                 <  1:27883108‑M2/61‑1 (MQ=255)

                   ttaatgTGGCTGGGGTACGAGGATTCGAACCACGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCt                                                <  1:14076104‑M2/62‑1 (MQ=255)

                   ttAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCt                                                <  1:4851573/67‑1 (MQ=39)

                   aTAATTTGGCTGGGGTACGAGGATTCGAACCTGGGAATGCCGGTATCATAAACCGGTGCCTTACCGCt                                                <  1:6295183/68‑1 (MQ=11)

                   aTAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGc                                                 >  1:233922/1‑67 (MQ=40)

                   aTAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCt                                                >  1:20250977/1‑68 (MQ=40)

                   aTAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCt                                                >  1:20907029/1‑68 (MQ=40)

                   aTAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCt                                                >  1:25550077/1‑68 (MQ=40)

                      atgTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTgg                                             <  1:5282801‑M2/65‑1 (MQ=255)

                       tgTGGGTTGGGTACGAGGTTTGGATCCTCGGAATTCCGGAATCAGAATCCGGTGCCTTACCGCTtgga                                            <  1:28275185‑M2/66‑2 (MQ=255)

                       tgTGGCTGGGGTACGAGGATTCGACCCTCGGAATGCCGGATTCAGAATCCGGTGCCTTACCGCTtggg                                            >  1:7182941‑M2/3‑67 (MQ=255)

                       tgTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGc                                            <  1:10419184‑M2/66‑1 (MQ=255)

                       tgTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCATAAACCGGTGCCTTACCGCTTGGc                                            <  1:9789118‑M2/66‑1 (MQ=255)

                       tgTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCgctttc                                             >  1:1145015‑M2/3‑65 (MQ=255)

                       tgTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGc                                            >  1:7841643‑M2/3‑68 (MQ=255)

                       tgTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGc                                            >  1:2991540‑M2/3‑68 (MQ=255)

                       tgTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGc                                            >  1:27169077‑M2/3‑68 (MQ=255)

                       tgTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGc                                            >  1:19988864‑M2/3‑68 (MQ=255)

                       tgTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGc                                            <  1:26606642‑M2/66‑1 (MQ=255)

                       tgTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGc                                            >  1:2043159‑M2/3‑68 (MQ=255)

                       tgTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGc                                            >  1:17873100‑M2/3‑68 (MQ=255)

                       tgTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGc                                            >  1:5644404‑M2/3‑68 (MQ=255)

                       tgTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGc                                            >  1:24951422‑M2/3‑68 (MQ=255)

                       tgTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGc                                            <  1:12122198‑M2/66‑1 (MQ=255)

                       tgTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGc                                            >  1:7529026‑M2/3‑68 (MQ=255)

                       tgTGGCTGGAGTACGAGGATTCGAACCTCGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGc                                            <  1:16826731‑M2/66‑1 (MQ=255)

                       aTTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGc                                            <  1:17575096/68‑1 (MQ=33)

                       aTTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGc                                            >  1:13572115/1‑68 (MQ=33)

                        ttGGCTGGGGTACGAGGATTCGACCCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCg                                           >  1:2511429/1‑68 (MQ=14)

                        ttGGCTGGGGTACGAGGATTCGAACCTCTGAATGCCGGAATCATAATCCGGTGCCTTACCGCTTGGCg                                           >  1:19629461/1‑68 (MQ=11)

                        ttGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGATTCAGAATCCGGTGCCTTACCGCTTGGCg                                           >  1:4541832/1‑68 (MQ=14)

                        ttGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCg                                           >  1:14799590/1‑68 (MQ=31)

                        ttGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCg                                           >  1:23903001/1‑68 (MQ=31)

                        gTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCg                                           <  1:29349735‑M2/67‑1 (MQ=255)

                                   aCGAGGATTCGAACCTCGGAATGCCGGTATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCcaaca                                <  1:14527994/68‑1 (MQ=14)

                                      aGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCcaacatag                             >  1:29028345/1‑65 (MQ=14)

                                      aGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCcaacagtg                             >  1:29242448/1‑65 (MQ=14)

                                      aGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCcaacagtg                             >  1:18735614/1‑65 (MQ=14)

                                      aGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCcaacaggg                             >  1:6604013/1‑65 (MQ=1)

                                      aGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAAGaaag                             >  1:9514621/1‑67 (MQ=9)

                                      aGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAAt                             >  1:145967/1‑68 (MQ=37)

                                      aGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAAt                             >  1:28363311/1‑68 (MQ=37)

                                      aGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGAGATACCCCAACAAAt                             >  1:2808132/1‑68 (MQ=255)

                                       ggATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAAtt                            >  1:6804246/1‑68 (MQ=39)

                                       ggATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAAtt                            >  1:11764155/1‑68 (MQ=39)

                                        gATTCGAACCTCGGAATGCCGGAATCATAATCCGGTGCCTTACCGCTTGGCGATACCCCAGCAAGTTg                           >  1:3235147/1‑68 (MQ=11)

                                        gATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTg                           >  1:8068558/1‑68 (MQ=40)

                                          ttCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGt                         <  1:4983774/68‑1 (MQ=255)

                                               aCCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGa                    >  1:8613704/1‑68 (MQ=255)

                                                   cGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAAttt                <  1:12972678/68‑1 (MQ=255)

                                                   cGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAAttt                <  1:10040903/68‑1 (MQ=255)

                                                       aTGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCg            >  1:19176622/1‑68 (MQ=255)

                                                        tGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGa           <  1:22303064/68‑1 (MQ=255)

                                                         gCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGaa          <  1:1740563/68‑1 (MQ=255)

                                                         gCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGaa          <  1:28887457/68‑1 (MQ=255)

                                                         gCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGa           >  1:23180645/1‑67 (MQ=255)

                                                           cGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACa        <  1:244938/68‑1 (MQ=255)

                                                            ggAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACat       <  1:8977156/68‑1 (MQ=255)

                                                            ggAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACat       <  1:29468524/68‑1 (MQ=255)

                                                            ggAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACat       <  1:24151447/68‑1 (MQ=255)

                                                              aaTCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATGTGCCGAACatat     >  1:21406524/1‑68 (MQ=255)

                                                               aTCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACatat     <  1:12706559/67‑1 (MQ=255)

                                                                tCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTc   >  1:3393267/1‑68 (MQ=255)

                                                                 cAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCg  <  1:9291769/68‑1 (MQ=255)

TACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCG  >  CP009273/691973‑692105

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 14 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`