BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,020,959	A→G	5.1%	G332G (GGT→GGC)	uidB ←	glucuronide transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,020,959	0	A	G	5.1%	337.6 / 5.7	118	G332G (GGT→GGC)	uidB	glucuronide transporter
Reads supporting (aligned to +/- strand): ref base A (49/62); new base G (4/2); total (53/65)
Fisher's exact test for biased strand distribution p-value = 4.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01

Only 100 of 119 total aligned reads displayed.

CAGGTATTCACCGTATTCTACGGTATCAGCTTCCAGCGCCCACATCACGGTCATGGTAACGCCCTGACCAATTGAAGCGATGGCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCA  >  minE/1020893‑1021029

cAGGTATTCACCGTATTCTACGGTATCAGCTTCCAGCGCCCACATCACGGTCATGGTAACGCCCTGACCa                                                                     <  1:2650757/70‑1 (MQ=255)

 aGGTATTCACCGTATTCTACGGTATCAGCTTCCAGCGCCCACATCACGGTCATGGTAACGCCCTGACCaa                                                                    <  1:3426921/70‑1 (MQ=255)

 aGGTATTCACCGTATTCTACGGTATCAGCTTCCAGCGCCCACATCACGGTCATGGTAACGCCCTGACCaa                                                                    <  1:4952824/70‑1 (MQ=255)

 aGGTATTCACCGTATTCTACGGTATCAGCTTCCAGCGCCCACATCACGGTCATGGTAACGCCCTGACCAAt                                                                   >  1:4742366/1‑71 (MQ=255)

  ggTATTCACCGTATTCTACGGTATCAGCTTCCAGCGCCCACATCACGGTCATGGTAACGCCCTGACCAAtt                                                                  >  1:2690435/1‑71 (MQ=255)

  ggTATTCACCGTATTCTACGGTATCAGCTTCCAGCGCCCACATCACGGTCATGGTAACGCCCTGACCAAtt                                                                  <  1:315134/71‑1 (MQ=255)

  ggTATTCACCGTATTCTACGGTATCAGCTTCCAGCGCCCACATCACGGTCATAGTAACGCCCTGACCAAt                                                                   >  1:614081/1‑70 (MQ=255)

    tATTCACCGTATTCTACGGTATCAGCTTCCAGCGCCCACATCACGGTCATGGTAACGCCCTGAcc                                                                      <  1:5232966/65‑1 (MQ=255)

     aTTCACCGTATTCTACGGTATCAGCTTCCAGCGCCCACATCACGGTCATGGTAACGCCCTGACCAAtt                                                                  >  1:1029253/1‑68 (MQ=255)

     aTTCACCGTATTCTACGGTATCAGCTTCCAGCGCCCACATCACGGTCATGGTAACGCCCTGACCAATTg                                                                 >  1:5470164/1‑69 (MQ=255)

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

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

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

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

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

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

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

        cACCGTATTCTACGGTATCAGCTTCCAGCGCCCACATCACGGTCATGGTAACGCCCTGACCAATTGAAGc                                                             <  1:5291437/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:4015475/1‑71 (MQ=255)

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

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

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

               ttCTACGGTATCAGCTTCCAGCGCCCACATCACGGTCATGGTAACGCCCTGACCAATTGAAGc                                                             <  1:2904484/63‑1 (MQ=255)

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

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

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

               ttCTACGGTATCAGCTTCCAGCGCCCACATCACGGTCATGGTAACGCCCTGACCAATTGAAGCGATGGcc                                                      <  1:2033564/70‑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)

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

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

                   aCGGTATCAGCTTCCAGCGCCCACATCACGGTCATGGTAACGCCCTGACCAATTGAAGCg                                                            <  1:1248781/60‑1 (MQ=255)

                   aCGGTATCAGCTTCCAGCGCCCACATCACGGTCATGGTAACGCCCTGACCAATTGAAGCg                                                            <  1:878389/60‑1 (MQ=255)

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

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

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

                       tATCAGCTTCCAGCGCCCACATCACGGTCATGGTAACGCCCTGACCAAtt                                                                  <  1:986165/50‑1 (MQ=255)

                       tATCAGCTTCCAGCGCCCACATCACGGTCATGGTAACGCCCTGACCAATTg                                                                 <  1:1186709/51‑1 (MQ=255)

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

                         tCAGCTTCCAGCGCCCACATCACGGTCATGGTAACGCCCTGACCAATTg                                                                 <  1:121414/49‑1 (MQ=255)

                         tCAGCTTCCAGCGCCCACATCACGGTCATGGTAACGCCCTGACCAATTg                                                                 <  1:791384/49‑1 (MQ=255)

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

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

                         tCAGCTTCCAGCGCCCACATCACGGTCATGGTAACGCCCTGACCAATTGAAGCGATGGCCAACGcaacca                                            <  1:904589/70‑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:1602004/65‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

                                     gcCCACATCACGGTCATGGTAACGCCCTGACCAATTGAAGCGATGGCCAACGCAACCAACGCCACCGGCAg                               >  1:4296010/1‑71 (MQ=255)

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

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

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

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

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

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

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

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

                                           aTCACGGTCATGGTAACGCCCTGACCAATTGAAGCGATGGCCAACGcaacc                                             >  1:3631307/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:3053292/1‑68 (MQ=255)

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

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

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

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

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

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

                                                ggTCATGGTAACGCCCTGACCAATTGAAGCGATGGCCAACGCAACCAACGCCACCgg                                  <  1:1073050/57‑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)

                                                             cccTGACCAATTGAAGCGATGGCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAgaagaa       <  1:1891955/71‑1 (MQ=255)

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

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

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

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

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

CAGGTATTCACCGTATTCTACGGTATCAGCTTCCAGCGCCCACATCACGGTCATGGTAACGCCCTGACCAATTGAAGCGATGGCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCA  >  minE/1020893‑1021029

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