BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	1,020,968	A→T	5.2%	A329A (GCT→GCA)	uidB ←	glucuronide transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	1,020,968	0	A	T	5.2%	331.8 / 7.6	115	A329A (GCT→GCA)	uidB	glucuronide transporter
Reads supporting (aligned to +/- strand): ref base A (48/61); new base T (4/2); total (52/63)
Fisher's exact test for biased strand distribution p-value = 4.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.08e-01

Only 100 of 120 total aligned reads displayed.

TCACCGTATTCTACGGTATCAGCTTCCAGCGCCCACATCACGGTCATGGTAACGCCCTGACCAATTGAAGCGATGGCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCG  >  Exported/1020900‑1021037

tCACCGTATTCTACGGTATCAGCTTCCAGCGCCCACATCACGGTCATGGTAACGCCCTGACCAATTGAAGc                                                                     >  1:3072707/1‑71 (MQ=255)

 cACCGTATTCTACGGTATCAGCTTCCAGCGCCCACATCACGGTCATGGTAACGCCCTGACCAATTGAAGc                                                                     <  1:5291437/70‑1 (MQ=255)

 cACCGTATTCTACGGTATCAGCTTCCAGCGCCCACATCACGGTCATGGTAACGCCCTGACCAATTGAAGc                                                                     <  1:4264461/70‑1 (MQ=255)

 cACCGTATTCTACGGTATCAGCTTCCAGCGCCCACATCACGGTCATGGTAACGCCCTGACCAATTGAAGc                                                                     <  1:3985325/70‑1 (MQ=255)

 cACCGTATTCTACGGTATCAGCTTCCAGCGCCCACATCACGGTCATGGTAACGCCCTGACCAATTGAAGc                                                                     <  1:3980758/70‑1 (MQ=255)

 cACCGTATTCTACGGTATCAGCTTCCAGCGCCCACATCACGGTCATGGTAACGCCCTGACCAATTGAAGc                                                                     <  1:2300572/70‑1 (MQ=255)

   ccGTATTCTACGGTATCAGCTTCCAGCGCCCACATCACGGTCATGGTAACGCCCTGACCAATTGAAGCGAt                                                                  >  1:5129786/1‑71 (MQ=255)

   ccGTATTCTACGGTATCAGCTTCCAGCGCCCACATCACGGTCATGGTAACGCCCTGACCAATTGAAGCGAt                                                                  >  1:3917002/1‑71 (MQ=255)

    cGTATTCTACGGTATCAGCTTCCAGCGCCCACATCACGGTCATGGTAACGCCCTGACCAATTGAAg                                                                      >  1:1696935/1‑66 (MQ=255)

    cGTATTCTACGGTATCAGCTTCCAGCGCCCACATCACGGTCATGGTAACGCCCTGACCAATTGAAGc                                                                     <  1:1399457/67‑1 (MQ=255)

     gTATTCTACGGTATCAGCTTCCAGCGCCCACATCACGGTCATGGTAACGCCCTGACCAATTGAAGCGATgg                                                                >  1:5391364/1‑71 (MQ=255)

     gTATTCTACGGTATCAGCTTCCAGCGCCCACATCACGGTCATGGTAACGCCCTGACCAATTGAAGCGATgg                                                                >  1:4015475/1‑71 (MQ=255)

     gTATTCTACGGTATCAGCTTCCAGCGCCCACATCACGGTCATGGTAACGCCCTGACCAATTGAAGCGATgg                                                                >  1:2592269/1‑71 (MQ=255)

      tATTCTACGGTATCAGCTTCCAGCGCCCACATCACGGTCATGGTAACGCCCTGACCAATTGAAGCGATGGc                                                               >  1:3811872/1‑71 (MQ=255)

        ttCTACGGTATCAGCTTCCAGCGCCCACATCACGGTCATGGTAACGCCCTGACCAATTGAAGCGATGGcc                                                              >  1:993290/1‑70 (MQ=255)

        ttCTACGGTATCAGCTTCCAGCGCCCACATCACGGTCATGGTAACGCCCTGACCAATTGAAGCGATGGcc                                                              <  1:2702012/70‑1 (MQ=255)

        ttCTACGGTATCAGCTTCCAGCGCCCACATCACGGTCATGGTAACGCCCTGACCAATTGAAGCGATGGcc                                                              <  1:2033564/70‑1 (MQ=255)

        ttCTACGGTATCAGCTTCCAGCGCCCACATCACGGTCATGGTAACGCCCTGACCAATTGAAGCGATGGCCa                                                             <  1:3551493/71‑1 (MQ=255)

         tCTACGGTATCAGCTTCCAGCGCCCACATCACGGTCATGGTAACGCCCTGACCAATTGAAGCGATGGcc                                                              <  1:3844857/69‑1 (MQ=255)

          cTACGGTATCAGCTTCCAGCGCCCACATCACGGTCATGGTAACGCCCTGACCAATTGAAGCGATGGCCaa                                                            >  1:1431679/1‑70 (MQ=255)

          cTACGGTATCAGCTTCCAGCGCCCACATCACGGTCATGGTAACGCCCTGACCAATTGAAGCGATGGCCAAc                                                           >  1:4116470/1‑71 (MQ=255)

          cTACGGTATCAGCTTCCAGCGCCCACATCACGGTCATGGTAACGCCCTGACCAATTGAAGCGATGGCCAAc                                                           >  1:1131919/1‑71 (MQ=255)

           tACGGTATCAGCTTCCAGCGCCCACATCACGGTCATGGTAACGCCCTGACCAATTGAAGCGa                                                                   >  1:5691048/1‑62 (MQ=255)

           tACGGTATCAGCTTCCAGCGCCCACATCACGGTCATGGTAACGCCCTGACCAATTGAAGCGATGGCCAACg                                                          >  1:4539851/1‑71 (MQ=255)

             cGGTATCAGCTTCCAGCGCCCACATCACGGTCATGGTAACGCCCTGACCAATTGaa                                                                       <  1:4776965/56‑1 (MQ=255)

             cGGTATCAGCTTCCAGCGCCCACATCACGGTCATGGTAACGCCCTGACCAATTGaa                                                                       <  1:2288859/56‑1 (MQ=255)

             cGGTATCAGCTTCCAGCGCCCACATCACGGTCATGGTAACGCCCTGACCAATTGAAGCGATGGCCAACGc                                                         <  1:529998/70‑1 (MQ=255)

                 aTCAGCTTCCAGCGCCCACATCACGGTCATGGTAACGCCCTGACCAATTGAAGCGATGGCCAACGcaacc                                                     <  1:3700144/70‑1 (MQ=255)

                 aTCAGCTTCCAGCGCCCACATCACGGTCATGGTAACGCCCTGACCAATTGAAGCGATGGCCAACGcaacc                                                     >  1:4204222/1‑70 (MQ=255)

                  tCAGCTTCCAGCGCCCACATCACGGTCATGGTAACGCCCTGACCAATTGAAGCGATGGCCAACGcaaccaa                                                   >  1:4008438/1‑71 (MQ=255)

                  tCAGCTTCCAGCGCCCACATCACGGTCATGGTAACGCCCTGACCAATTGAAGCGATGGCCAACGcaacca                                                    <  1:904589/70‑1 (MQ=255)

                  tCAGCTTCCAGCGCCCACATCACGGTCATGGTAACGCCCTGACCAATTGAAGCGATGGCCAACGcaacca                                                    <  1:768352/70‑1 (MQ=255)

                   cAGCTTCCAGCGCCCACATCACGGTCATGGTAACGCCCTGACCAATTGAAGCGATGGCCAACGcaacca                                                    <  1:3104102/69‑1 (MQ=255)

                    aGCTTCCAGCGCCCACATCACGGTCATGGTAACGCCCTGACCAATTGaa                                                                       >  1:2462355/1‑49 (MQ=255)

                    aGCTTCCAGCGCCCACATCACGGTCATGGTAACGCCCTGACCAATTGAAGCGATGGCCAACGcaa                                                       >  1:3429065/1‑65 (MQ=255)

                    aGCTTCCAGCGCCCACATCACGGTCATGGTAACGCCCTGACCAATTGAAGCGATGGCCAACGCAACCAACg                                                 <  1:5576290/71‑1 (MQ=255)

                      cTTCCAGCGCCCACATCACGGTCATGGTAACGCCCTGACCAATTGAAGCGATGGCCAACGcaaccaa                                                   >  1:4215578/1‑67 (MQ=255)

                      cTTCCAGCGCCCACATCACGGTCATGGTAACGCCCTGACCAATTGAAGCGATGGCCAACGCAACCAACGcc                                               <  1:4788227/71‑1 (MQ=255)

                         ccAGCGCCCACATCACGGTCATGGTAACGCCCTGACCAATTGAAGCGa                                                                   >  1:744955/1‑48 (MQ=255)

                           aGCGCCCACATCACGGTCATGGTAACGCCCTGACCAATTGAAGCGATGGCCAACGCAACCAACGc                                                <  1:3608729/65‑1 (MQ=255)

                           aGCGCCCACATCACGGTCATGGTAACGCCCTGACCAATTGAAGCGATGGCCAACGCAACCAACGCCACCgg                                          <  1:3988408/71‑1 (MQ=255)

                           aGCGCCCACATCACGGTCATGGTAACGCCCTGACCAATTGAAGCGATGGCCAACGCAACCAACGCCACCgg                                          <  1:533918/71‑1 (MQ=255)

                           aGCGCCCACATCACGGTCATGGTAACGCCCTGACCAATTGAAGCGATGGCCAACGCAACCAACGCCACCgg                                          <  1:181057/71‑1 (MQ=255)

                           aGCGCCCACATCACGGTCATGGTAACGCCCTGACCAATTGAAGCGATGGCCAACGCAACCAACGCCACCgg                                          <  1:3473069/71‑1 (MQ=255)

                           aGCGCCCACATCACGGTCATGGTAACGCCCTGACCAATTGAAGCGATGGCCAACGCAACCAACGCCACCgg                                          <  1:4569368/71‑1 (MQ=255)

                           aGCGCCCACATCACGGTCATGGTAACGCCCTGACCAATTGAAGCGATGGCCAACGCAACCAACGCCACCgg                                          <  1:174188/71‑1 (MQ=255)

                            gcgcCCACATCACGGTCATGGTAACGCCCTGACCAATTGAAGCGATGGCCAACGCAACCAACGCCACCgg                                          >  1:2052206/1‑70 (MQ=255)

                             cgcCCACATCACGGTCATGGTAACGCCCTGACCAATTGAAGCGATGGCCAACGCAACCAACGCCACCGGCa                                        <  1:3637390/71‑1 (MQ=255)

                             cgcCCACATCACGGTCATGGTAACGCCCTGACCAATTGAAGCGATGGCCAACGCAACCAACGCCACCGGCa                                        <  1:5691560/71‑1 (MQ=255)

                             cgcCCACATCACGGTCATGGTAACGCCCTGACCAATTGAAGCGATGGCCAACGCAACCAACGCCACCGGCa                                        <  1:421009/71‑1 (MQ=255)

                             cgcCCACATCACGGTCATGGTAACGCCCTGACCAATTGAAGCGATGGCCAACGCAACCAACGCCACCGGCa                                        <  1:427759/71‑1 (MQ=255)

                              gcCCACATCACGGTCATGGTAGCGCCCTGACCAATTGAAGCGATGGCCAACGCAACCAACGCCACCGGCa                                        <  1:770949/70‑1 (MQ=255)

                                ccACATCACGGTCATGGTAACGCCCTGACCAATTGAAGCGATGGCCAACGCAACCAACGCCACCGGCAGTg                                     <  1:5261213/71‑1 (MQ=255)

                                   caTCACGGTCATGGTAACGCCCTGGCCAATTGATGCGATGGCCAACGCAACCAACGCCACCGGCAGTGa                                    >  1:2907770/1‑69 (MQ=255)

                                   caTCACGGTCATGGTAACGCCCTGGCCAATTGATGCGATGGCCAACGCAACCAACGCCACCGGCAGTGAc                                   >  1:1241314/1‑70 (MQ=255)

                                   caTCACGGTCATGGTAACGCCCTGGCCAATTGATGCGATGGCCAACGCAACCAACGCCACCGGCAGTGAc                                   >  1:87639/1‑70 (MQ=255)

                                   caTCACGGTCATGGTAACGCCCTGGCCAATTGATGCGATGGCCAACGCAACCAACGCCACCGGCAGTGAc                                   >  1:2775524/1‑70 (MQ=255)

                                    aTCACGGTCATGGTAACGCCCTGACCAATTGAAGCGATGGCCAACGcaacc                                                     >  1:2252618/1‑51 (MQ=255)

                                    aTCACGGTCATGGTAACGCCCTGACCAATTGAAGCGATGGCCAACGcaacc                                                     >  1:3631307/1‑51 (MQ=255)

                                    aTCACGGTCATGGTAACGCCCTGACCAATTGAAGCGATGGCCAACGcaacc                                                     >  1:3523967/1‑51 (MQ=255)

                                     taaCGGTCATGGTAACGCCCTGACCAATTGAAGCGATGGCCAACGCAACCAACGCCACCGGCAGTGAc                                   >  1:5367818/3‑68 (MQ=255)

                                     tCACGGTCATGGTAACGCCCTGACCAATTGAAGCGATGGCCAACGc                                                         <  1:1347321/46‑1 (MQ=255)

                                     tCACGGTCATGGTAACGCCCTGACCAATTGAAGCGATGGCCAACGc                                                         <  1:3562798/46‑1 (MQ=255)

                                     tCACGGTCATGGTAACGCCCTGACCAATTGAAGCGATGGCCAACGCAACCAACGCCACCGGCAGTGAc                                   >  1:1740024/1‑68 (MQ=255)

                                     tCACGGTCATGGTAACGCCCTGACCAATTGAAGCGATGGCCAACGCAACCAACGCCACCGGCAGTGAc                                   >  1:3053292/1‑68 (MQ=255)

                                     tCACGGTCATGGTAACGCCCTGACCAATTGAAGCGATGGCCAACGCAACCAACGCCACCGGCAGTGACCAg                                <  1:13581/71‑1 (MQ=255)

                                     tCACGGTCATGGTAACGCCCTGACCAATTGAAGCGATGGCCAACGCAACCAACGCCACCGGCAGTGACCAg                                <  1:2895597/71‑1 (MQ=255)

                                       aCGGTCATGGTAACGCCCTGACCAATTGAAGCGATGGCCaa                                                            >  1:4233316/1‑41 (MQ=255)

                                       aCGGTCATGGTAACGCCCTGACCAATTGAAGCGATGGCCAACGCAACCAACGCCACCGGCa                                        >  1:204941/1‑61 (MQ=255)

                                       aCGGTCATGGTAACGCCCTGACCAATTGAAGCGATGGCCAACGCAACCAACGCCACCGGCa                                        >  1:309281/1‑61 (MQ=255)

                                       aCGGTCATGGTAACGCCCTGACCAATTGAAGCGATGGCCAACGCAACCAACGCCACCGGCAGTGACCAGa                               <  1:894326/70‑1 (MQ=255)

                                       aCGGTCATGGTAACGCCCTGACCAATTGAAGCGATGGCCAACGCAACCAACGCCACCGGCAGTGACCAGAc                              >  1:4920963/1‑71 (MQ=255)

                                         ggTCATGGTAACGCCCTGACCAATTGAAGCGATGGCCAACGCAACCAACGCCACCgg                                          <  1:1073050/57‑1 (MQ=255)

                                          gTCATGGTAACGCCCTGTCCAATTGAAGCGATGGCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGa                           <  1:1540576/71‑1 (MQ=255)

                                           tCATGGTAACGCCCTGACCAATTGAAGCGATGGCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGa                           <  1:2585166/70‑1 (MQ=255)

                                           tCATGGTAACGCCCTGACCAATTGAAGCGATGGCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGa                           >  1:1155596/1‑70 (MQ=255)

                                                gTAACGCCCTGACCAATTGAAGCGATGGCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAAcc                       >  1:616881/1‑69 (MQ=255)

                                                   aCGCCCTGGCCAATTGATGCGATGGCCAACGCAACCAACGCCACCGGCAGt                                      <  1:3791714/51‑1 (MQ=255)

                                                   aCGCCCTGGCCAATTGATGCGATGGCCAACGCAACCAACGCCACCGGCAGt                                      <  1:418270/51‑1 (MQ=255)

                                                     gCCCTGACCAATTGAAGCGATGGCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCagaag                 >  1:1731338/1‑70 (MQ=255)

                                                     gCCCTGACCAATTGAAGCGATGGCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCagaa                  >  1:2003809/1‑69 (MQ=255)

                                                       ccTGACCAATTGAAGCGATGGCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAgaagaa               <  1:674143/70‑1 (MQ=255)

                                                         tGACCAATTGAAGCGATGGCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCagaaga                <  1:4216432/67‑1 (MQ=255)

                                                         tGACCAATTGAAGCGATGGCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAAcag            >  1:981342/1‑71 (MQ=255)

                                                         tGACCAATTGAAGCGATGGCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAAcag            >  1:967516/1‑71 (MQ=255)

                                                          gACCAATTGAAGCGATGGACAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAAcag            >  1:3706500/1‑70 (MQ=255)

                                                           aCCAATTGAAGCGATGGCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAAcagca          <  1:558847/71‑1 (MQ=255)

                                                              aaTTGAAGCGATGGCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGAt       <  1:4105449/71‑1 (MQ=255)

                                                              aaTTGAAGCGATGGCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGAt       <  1:5571749/71‑1 (MQ=255)

                                                               aTTGAAGCGATGGCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGAt       >  1:507000/1‑70 (MQ=255)

                                                                ttGAAGCGATGGCCAACGCAACCAACGCCACCGGCAGt                                      <  1:4613654/38‑1 (MQ=255)

                                                                ttGAAGCGATGGCCAACGCAACCAACGCCACCGGCAGt                                      <  1:2334524/38‑1 (MQ=255)

                                                                ttGAAGCGATGGCCAACGCAACCAACGCCACCGGCAGt                                      <  1:4646445/38‑1 (MQ=255)

                                                                ttGAAGCGATGGCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCaga                   <  1:750038/57‑1 (MQ=255)

                                                                ttGAAGCGATGGCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCaga                   <  1:2064551/57‑1 (MQ=255)

                                                                 tGAAGCGATGGCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCaga                   <  1:4627650/56‑1 (MQ=255)

                                                                  gAAGCGATGGCCAACGCAACCAACGCCACCGGCAGTGACCAGa                               >  1:3584359/1‑43 (MQ=255)

                                                                  gAAGCGATGGCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGAt       >  1:5570394/1‑67 (MQ=255)

                                                                   aaGCGATGGCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCg  >  1:1173408/1‑71 (MQ=255)

                                                                    aGCGATGGCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCg  >  1:4786607/1‑70 (MQ=255)

TCACCGTATTCTACGGTATCAGCTTCCAGCGCCCACATCACGGTCATGGTAACGCCCTGACCAATTGAAGCGATGGCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCG  >  Exported/1020900‑1021037

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