BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	3,244,017	G→A	31.5%	Q86Q (CAG→CAA)	exuT →	hexuronate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	3,244,017	0	G	A	31.5%	183.4 / 114.7	130	Q86Q (CAG→CAA)	exuT	hexuronate transporter
Reads supporting (aligned to +/- strand): ref base G (48/41); new base A (24/17); total (72/58)
Fisher's exact test for biased strand distribution p-value = 7.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.72e-01

Only 100 of 133 total aligned reads displayed.

CCTGACGCGTAACACTGTGGCGGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCGGTAGCAGGTTATGTACTG  >  W3110S.gb/3243948‑3244086

ccTGACGCGTAACACTGTGGCGGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAg                                                                       <  1:3401665/70‑1 (MQ=255)

   gACGCGTAACACTGTGGCGGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTAtt                                                                   >  1:878681/1‑71 (MQ=255)

   gACGCGTAACACTGTGGCGGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTAt                                                                    >  1:294162/1‑70 (MQ=255)

   gACGCGTAACACTGTGGCGGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAATAtt                                                                   >  1:874371/1‑71 (MQ=255)

     cgcgTAACACTGTGGCGGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTcc                                                                 <  1:316133/71‑1 (MQ=255)

     cgcgTAACACTGTGGCGGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTcc                                                                 <  1:4449759/71‑1 (MQ=255)

     cgcgTAACACTGTGGCGGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTcc                                                                 >  1:3817862/1‑71 (MQ=255)

          aaCACTGTGGCGGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAATATTCCTatat                                                            >  1:2215230/1‑71 (MQ=255)

           acacTGTGGCGGCAGTTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTatat                                                            >  1:600515/1‑70 (MQ=255)

           acacTGTGGCGGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTatat                                                            >  1:3513330/1‑70 (MQ=255)

           acacTGTGGCGGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTatat                                                            >  1:3355403/1‑70 (MQ=255)

           acacTGTGGCGGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTatat                                                            >  1:3189091/1‑70 (MQ=255)

           acacTGTGGCGGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTatat                                                            >  1:4861316/1‑70 (MQ=255)

           acacTGTGGCGGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAATATTCCTatat                                                            >  1:3654019/1‑70 (MQ=255)

           acacTGTGGCGGCAGCTGCGCAAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAATATTCCTatat                                                            >  1:3281076/1‑70 (MQ=255)

              cTGTGGCGGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATatcatc                                                        <  1:3818612/71‑1 (MQ=255)

              cTGTGGCGGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATatcatc                                                        <  1:399532/71‑1 (MQ=255)

              cTGTGGCGGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATatcatc                                                        >  1:214630/1‑71 (MQ=255)

              cTGTGGCGGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATatcat                                                         >  1:2803635/1‑70 (MQ=255)

                gtgGCGGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCg                                                       >  1:990315/1‑70 (MQ=255)

                 tgGCGGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCg                                                       >  1:3270320/1‑69 (MQ=255)

                   gcggcAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTatat                                                            >  1:4318252/1‑62 (MQ=255)

                   gcggcAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGc                                                   <  1:2919377/71‑1 (MQ=255)

                      gcAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTa                                                <  1:3504210/71‑1 (MQ=255)

                      gcAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTa                                                <  1:1894662/71‑1 (MQ=255)

                       cAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTa                                                >  1:2451358/1‑70 (MQ=255)

                         gCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAATATTCCTATATCATCGCAGCCTATTc                                             <  1:3950673/71‑1 (MQ=255)

                          cTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTc                                                                  <  1:59636/49‑1 (MQ=255)

                          cTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTc                                             <  1:4983232/70‑1 (MQ=255)

                           tGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTcc                                                                 >  1:4907164/1‑49 (MQ=255)

                           tGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTcc                                                                 >  1:2962713/1‑49 (MQ=255)

                           tGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAATATTCCTATATCATCGCAGCCTATTctg                                           >  1:1410172/1‑71 (MQ=255)

                           tGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAATATTCCTATATCATCGCAGCCTATTctg                                           >  1:4805516/1‑71 (MQ=255)

                                cAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAATATTCCTATatca                                                          <  1:2236544/51‑1 (MQ=255)

                                cAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAATATTCCTATatca                                                          <  1:1516466/51‑1 (MQ=255)

                                cAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAATATTCCTATatca                                                          <  1:3712476/51‑1 (MQ=255)

                                  aCTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATatca                                                          <  1:229855/49‑1 (MQ=255)

                                  aCTCTGATGGAAGAGTTAAACATCTCCACCCAACAATa                                                                     >  1:311044/1‑38 (MQ=255)

                                    tctGATGGAAGAGTTAAACATCTCCACCCAACAatat                                                                    >  1:4797763/1‑37 (MQ=255)

                                    tctGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATatcatc                                                        >  1:1739257/1‑49 (MQ=255)

                                    tctGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTa                                                <  1:484943/57‑1 (MQ=255)

                                    tctGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTa                                                <  1:2310698/57‑1 (MQ=255)

                                    tctGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTtata                                  >  1:4399962/1‑71 (MQ=255)

                                    tctGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCACCCTa                                                <  1:4199083/57‑1 (MQ=255)

                                    tctGATGGAAGAGTTAAACATCTCCACCCAACAATATTcc                                                                 >  1:4132206/1‑40 (MQ=255)

                                    tctGATGGAAGAGTTAAACATCTCCACCCAACAATATTcc                                                                 >  1:2114577/1‑40 (MQ=255)

                                    tctGATGGAAGAGTTAAACATCTCCACCCAACAATATTcc                                                                 >  1:1644861/1‑40 (MQ=255)

                                    tctGATGGAAGAGTTAAACATCTCCACCCAACAATATTcc                                                                 >  1:4428795/1‑40 (MQ=255)

                                    tctGATGGAAGAGTTAAACATCTCCACCCAACAATATTcc                                                                 >  1:4319672/1‑40 (MQ=255)

                                     ctGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATAc                                 >  1:451829/1‑71 (MQ=255)

                                      tGATGGAAGAGTTAAACATCTCCACCCAACAATATTCCTATATCATCGCAGCCTATTCTGCTGCTtata                                  >  1:1866683/1‑69 (MQ=255)

                                       gATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTtata                                  <  1:472464/68‑1 (MQ=255)

                                        aTGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATAc                                 <  1:3067745/68‑1 (MQ=255)

                                        aTGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATAc                                 <  1:2654000/68‑1 (MQ=255)

                                        aTGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACgg                               <  1:4507689/70‑1 (MQ=255)

                                        aTGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGt                              >  1:4819746/1‑71 (MQ=255)

                                           gAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTAtt                                              >  1:161920/1‑52 (MQ=255)

                                           gAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCa                            >  1:2688796/1‑70 (MQ=255)

                                             agagTTAAACATCTCCACCCAACAATATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGc                         <  1:4536935/71‑1 (MQ=255)

                                                gTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAAc                      <  1:333266/71‑1 (MQ=255)

                                                 ttAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCAt                           <  1:4950045/65‑1 (MQ=255)

                                                 ttAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCAt                           <  1:728366/65‑1 (MQ=255)

                                                 ttAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAAcc                     <  1:2681978/71‑1 (MQ=255)

                                                   aaaCATCTCCACCCAACAGTATTCCTATATCATCGCAGc                                                   >  1:3043915/1‑39 (MQ=255)

                                                   aaaCATCTCCACCCAACAGTATTCCTATATCATCGCAGCCt                                                 >  1:1384672/1‑41 (MQ=255)

                                                   aaaCATCTCCACCCAACAGTATTCCTATATCATCGCAGCCt                                                 >  1:1026859/1‑41 (MQ=255)

                                                   aaaCATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTctgc                                          >  1:613495/1‑48 (MQ=255)

                                                   aaaCATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTctgc                                          >  1:4609970/1‑48 (MQ=255)

                                                   aaaCATCTCCACCCAACAATATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCgg                   >  1:221280/1‑71 (MQ=255)

                                                   aaaCATCTCCACCCAACAATATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCgg                   >  1:705394/1‑71 (MQ=255)

                                                    aaCATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCGGt                  <  1:3850455/71‑1 (MQ=255)

                                                    aaCATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCGGt                  <  1:3404588/71‑1 (MQ=255)

                                                    aaCATCTCCACCCAACAATATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCgg                   <  1:4203098/70‑1 (MQ=255)

                                                     aCATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCtt                                     >  1:1189561/1‑51 (MQ=255)

                                                     aCATCTCCACCCAACAATATTCCTATATCATCGCAGCCTATTctgct                                         >  1:3379956/1‑47 (MQ=255)

                                                      cATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCGGTAg                <  1:4275264/71‑1 (MQ=255)

                                                       aTCTCCACCCAACAGTATTCCTATTTCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCGGt                  <  1:3575007/68‑1 (MQ=255)

                                                       aTCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCGGt                  <  1:1713449/68‑1 (MQ=255)

                                                       aTCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCGGt                  <  1:732674/68‑1 (MQ=255)

                                                       aTCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCGGTAGc               <  1:3262479/71‑1 (MQ=255)

                                                         ctcCACCTAACAATATTCCTATATCATCGCAGCCTAt                                               <  1:2298565/37‑1 (MQ=38)

                                                         ctcCACCCAACAGTATTCCTATATCATCGCAGCCTATTctgc                                          >  1:1227112/1‑42 (MQ=255)

                                                         ctcCACCCAACAGTATTCCTATATCATCGCAGCCTATTctgc                                          >  1:1187792/1‑42 (MQ=255)

                                                         ctcCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATAc                                 <  1:3479563/51‑1 (MQ=255)

                                                          tcCACCCAACAGTATTCCTATATCATCGCAGCCTAtt                                              >  1:1403348/1‑37 (MQ=255)

                                                          tcCACCCAACAGTATTCCTATATCATCGCAGCCTAt                                               >  1:2956618/1‑36 (MQ=255)

                                                           ccACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCGGTAGCAGGt           >  1:2991279/1‑71 (MQ=255)

                                                           ccACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCGGTAGCAGGt           >  1:2752530/1‑71 (MQ=255)

                                                           ccACCCAACAATATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCGGTAGCAGGt           >  1:1706628/1‑71 (MQ=255)

                                                            cACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCGGTAGCa              <  1:2534283/67‑1 (MQ=255)

                                                            cACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCGGTAGCAGGt           <  1:4590794/70‑1 (MQ=255)

                                                            cACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCGGTAGCAGGt           <  1:4956288/70‑1 (MQ=255)

                                                               ccAACAATATTCCTATATCATCGCAGCCTATTctg                                           <  1:2555840/35‑1 (MQ=255)

                                                               ccAACAATATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCGGTAGCAGGTTATg       <  1:2779435/71‑1 (MQ=255)

                                                               ccAACAATATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCGGTAGCAGGTTATg       <  1:3806329/71‑1 (MQ=255)

                                                               ccAACAATATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCGGTAGCAGGTTATg       <  1:3182388/71‑1 (MQ=255)

                                                                cAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTtata                                  >  1:2744723/1‑43 (MQ=255)

                                                                  aCAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCGGTAGCAGGTTATGTAc    >  1:1909682/1‑71 (MQ=255)

                                                                  aCAATATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCGGTAGCAGGTTATGTAc    >  1:2107267/1‑71 (MQ=255)

                                                                     gTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCGGTAGCAGGTTATGTACTg  >  1:3933510/1‑70 (MQ=255)

CCTGACGCGTAACACTGTGGCGGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCGGTAGCAGGTTATGTACTG  >  W3110S.gb/3243948‑3244086

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