BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,596,014	G→T	8.8%	intergenic (‑32/+3)	leuV ← / ← leuP	tRNA‑Leu/tRNA‑Leu

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	4,596,014	0	G	T	8.8%	1333.2 / 8.0	517	intergenic (‑32/+3)	leuV/leuP	tRNA‑Leu/tRNA‑Leu
Reads supporting (aligned to +/- strand): ref base G (256/176); new base T (2/59); total (259/258)
Fisher's exact test for biased strand distribution p-value = 1.03e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 536 total aligned reads displayed.

TGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCG  >  CP009273/4595947‑4596080

tGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCg                                                                    <  1:36958403/68‑1 (MQ=255)

tGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCg                                                                    <  1:42298429/68‑1 (MQ=255)

tGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCg                                                                    <  1:50068470/68‑1 (MQ=255)

tGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCg                                                                    <  1:24523996/68‑1 (MQ=255)

 gAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGt                                                                   <  1:7450378/68‑1 (MQ=255)

   aGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGttt                                                                 >  1:12714354/1‑68 (MQ=255)

   aGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGtt                                                                  <  1:25143021/67‑1 (MQ=255)

     cTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTgg                                                               >  1:33702392/1‑68 (MQ=255)

      tAGCGCGTCTACCAATTCCGCCCCCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTgg                                                               <  1:53519344/67‑1 (MQ=255)

      tAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGt                                                              <  1:9074238/68‑1 (MQ=255)

      tAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGt                                                              <  1:8150364/68‑1 (MQ=255)

      tAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGt                                                              <  1:23515831/68‑1 (MQ=255)

      aaGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGt                                                              <  1:6015817/67‑1 (MQ=255)

          gcgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGa                                                          <  1:39651395/68‑1 (MQ=255)

            gTCTAGCAATTCCGCCACCTTCGCACAGTCATGTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAg                                                         <  1:8388034/67‑1 (MQ=255)

              cTACCAATTCCGCCACCTTCGCGCAGTCATCTTACTTTTTTTGATATCGCCTCCTTTGGTGCGAgggg                                                      <  1:54997470/68‑1 (MQ=38)

              cTACCAATTCCGCCACCTTCGCACAGGCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAggg                                                       <  1:5605289/67‑1 (MQ=255)

              cTACCAATTCCGCCACCTTAGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAggg                                                       <  1:11091958/67‑1 (MQ=255)

                  cAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAgggggggg                                                  >  1:49464505/1‑68 (MQ=255)

                   aaTTCCGCCACCAAAACACATGCATCTTACATTTTTTGATATCCCCACCTTTGGTGCGAGGGGGGGGa                                                 <  1:40859039/68‑1 (MQ=16)

                    aTTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGACTCGTTTGGTGCGAgggggggg                                                  >  1:2932360/1‑66 (MQ=255)

                    aTTCCGCCACCTTCGCACACACTGATTACTTTAAATGATATCGGCTCTGGTGGTGCGAgggggggg                                                  >  1:47231370/1‑66 (MQ=255)

                    aTTCCGCCACCTTCGCAAAGTCAGCTTACTTTTTTTGATATCGGCGCGGTTGGTGCGAgggggggg                                                  >  1:37414482/1‑66 (MQ=255)

                     ttCCGCCACCTTCGCAAAGTCAGCTTACTTTTTATGATATCGGCTCGTTTGGTGCGAGGGGTGGGCCt                                               >  1:20094806/1‑68 (MQ=255)

                     ttCCGCCACCTACGCACAGTCAACTAACTATTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACt                                               <  1:51470029/68‑1 (MQ=38)

                      tCCGCCACCTTCGCACAGTCATCTTACATTTTTTGATATCGCCACCTTTGGTGCGAGGGGGGGGACtt                                              <  1:13839330/68‑1 (MQ=38)

                        cGCCACCTTCCCACATGCATCTAACTTTTTTTGATATCCCCTCCTTTGGTGCGAGGGGGGGGACTTGa                                            <  1:27117768/68‑1 (MQ=25)

                           cACCTTCGCACAGTCAGATTACTTTTTTTGATATCGACTCGTTTGGTGCGAGGGGGGGGACATGAAc                                          >  1:32382302/1‑67 (MQ=37)

                            aCCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAAccc                                        <  1:430308/68‑1 (MQ=255)

                                   cacaGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGt                                  >  1:31579029/1‑67 (MQ=255)

                                   cacaGGCATCAAACATTATTTGATATCCCCACTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTc                                 <  1:31167336/68‑1 (MQ=14)

                                   aacaGTCATCTTACTTTTATTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTc                                 <  1:43652444/67‑1 (MQ=255)

                                   aacaGTCAAGTTACTTTTTTTGATATAGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTc                                 <  1:11990272/67‑1 (MQ=21)

                                    acaGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAAcc                                         >  1:562836/1‑59 (MQ=255)

                                    acaGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTcc                                >  1:44303663/1‑68 (MQ=255)

                                    acaGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTcc                                >  1:18404724/1‑68 (MQ=255)

                                    acaGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTcc                                >  1:3644661/1‑68 (MQ=255)

                                    acaGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTcc                                >  1:44513133/1‑68 (MQ=255)

                                    acaGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTcc                                >  1:13466958/1‑68 (MQ=255)

                                    acaGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTcc                                >  1:54262037/1‑68 (MQ=255)

                                    acaGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTcc                                >  1:1673107/1‑68 (MQ=255)

                                    acaGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTcc                                >  1:21377052/1‑68 (MQ=255)

                                    acaGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTcc                                >  1:57089254/1‑68 (MQ=255)

                                    acaGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGGGCGAGGGGGGGGACTTGAACCCCCACGTcc                                <  1:53110987/68‑1 (MQ=255)

                                    acaGTCATCTTACATTATTTGAAATCCCCTATTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTcc                                <  1:14752280/68‑1 (MQ=17)

                                     cattcatcaTACTTTTTTTGATATCCCCTAGATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCg                               <  1:27821120/65‑1 (MQ=18)

                                     cattCATCTAACATTTATTGATATCCCCTCTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCg                               <  1:29899207/65‑1 (MQ=17)

                                     cattCAACTAACATTTTTTGAAATCGCCTAGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCg                               <  1:19161186/65‑1 (MQ=17)

                                      aGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGt                              >  1:41999830/1‑68 (MQ=255)

                                      aGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGt                              >  1:9477943/1‑68 (MQ=255)

                                       gTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTa                             >  1:10638471/1‑68 (MQ=255)

                                       gTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTa                             >  1:7141045/1‑68 (MQ=255)

                                       gTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTa                             >  1:18844841/1‑68 (MQ=255)

                                        tCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTaa                            >  1:41671326/1‑68 (MQ=255)

                                         gaacTTACTTTTTTTGATATCGCCTCCTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCg                               <  1:28156521/61‑1 (MQ=255)

                                         cATCTTACTTTTTTTGATATCTCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCg                               <  1:54048236/64‑1 (MQ=255)

                                         cATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTa                             <  1:13060677/66‑1 (MQ=255)

                                         cATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAg                           >  1:56914211/1‑68 (MQ=255)

                                         cATCTTACTTTTTTTGATATCCCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCg                               <  1:52724331/64‑1 (MQ=255)

                                         cATCTAACTATATTTGATATCCCCACGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTaa                            <  1:52042654/67‑1 (MQ=17)

                                          aTCTTACATTTTTTGATATACCCTTTATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGa                          <  1:2372147/68‑1 (MQ=14)

                                          aTCTAACTTTTTTTGATATACCCACTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGa                          <  1:24703399/68‑1 (MQ=17)

                                            cTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAac                        >  1:22113147/1‑68 (MQ=255)

                                             ttACTTTTTTTGATATCGCCACTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAaca                       <  1:35814290/68‑1 (MQ=255)

                                             atACTTTTTTTGATATCGCCTCTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAaca                       <  1:28760250/67‑1 (MQ=255)

                                             atACTTTTTTTGAAATCCCCTCCTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAaca                       <  1:5107977/67‑1 (MQ=18)

                                              tACTTTTTTTGATATCGCCTCGTTTGGTGCGTGGGGGGGGACTTGAACCCCCACGTCCGTAAGAac                        >  1:39263022/1‑66 (MQ=255)

                                              tACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCCCGTCCGTAAGaaa                        >  1:50688879/1‑65 (MQ=255)

                                              tACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGac                        >  1:45439691/1‑66 (MQ=255)

                                              tACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAac                        >  1:48149054/1‑66 (MQ=255)

                                              tACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAac                        >  1:6734214/1‑66 (MQ=255)

                                              tACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAac                        >  1:6325628/1‑66 (MQ=255)

                                              tACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAac                        >  1:57709758/1‑66 (MQ=255)

                                              tACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAac                        >  1:33068311/1‑66 (MQ=255)

                                              tACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAac                        >  1:25928846/1‑66 (MQ=255)

                                              tACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAac                        >  1:45818326/1‑66 (MQ=255)

                                              tACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAac                        >  1:26552273/1‑66 (MQ=255)

                                              tACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAac                        >  1:47660577/1‑66 (MQ=255)

                                              tACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAac                        >  1:49081609/1‑66 (MQ=255)

                                              tACTTTTTTTGAAATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAacac                      <  1:26176114/68‑1 (MQ=255)

                                                cTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTa                    >  1:17365043/1‑68 (MQ=255)

                                                 tttttttGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTa                    >  1:16214599/1‑67 (MQ=255)

                                                 tttttttGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTa                    >  1:30415314/1‑67 (MQ=255)

                                                  ttttttGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAac                  >  1:56579912/1‑68 (MQ=255)

                                                   tttttGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAaca                 <  1:15317220/68‑1 (MQ=255)

                                                   tttttGATATCGCCTAGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAaca                 <  1:3181480/68‑1 (MQ=255)

                                                   tttttGATATCCCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAaca                 <  1:20580872/68‑1 (MQ=255)

                                                   tttttGAAATCGCGTCTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAaca                 <  1:21814095/68‑1 (MQ=17)

                                                   tttttGAAATACCCACGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAaca                 <  1:23005012/68‑1 (MQ=14)

                                                   tatttGATATCCCCTCTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAaca                 <  1:46616237/66‑1 (MQ=17)

                                                     tttGATCTCGACTATTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACAcc               <  1:54695051/68‑1 (MQ=14)

                                                     tttGATATCTCCTTGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACAcc               <  1:15815335/68‑1 (MQ=255)

                                                     tttGATATCGCCTCTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACAcc               <  1:38513366/68‑1 (MQ=255)

                                                     tttGAATTCGCCTGTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACAcc               <  1:46094678/68‑1 (MQ=14)

                                                      ttGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAaca                 >  1:7922141/1‑65 (MQ=255)

                                                       tGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTg             >  1:10000077/1‑68 (MQ=255)

                                                          aatCCCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGaa           <  1:35052709/66‑1 (MQ=255)

                                                               ccTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCCCGTCCGTAAGAACACTAACACCTGAAGCTAgc     >  1:34407464/1‑68 (MQ=255)

                                                               ccTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAgc     >  1:48403558/1‑68 (MQ=255)

                                                               ccTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAgc     >  1:38813335/1‑68 (MQ=255)

TGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCG  >  CP009273/4595947‑4596080

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