BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	3,177,041	T→C	10.3%	L91L (TTG→CTG)	tolC →	transport channel

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	3,177,041	0	T	C	10.3%	254.3 / 17.2	97	L91L (TTG→CTG)	tolC	transport channel
Reads supporting (aligned to +/- strand): ref base T (37/50); new base C (8/2); total (45/52)
Fisher's exact test for biased strand distribution p-value = 4.11e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

CAGATTACACCTATAGCAACGGCTACCGCGACGCGAACGGCATCAACTCTAACGCGACCAGTGCGTCCTTGCAGTTAACTCAATCCATTTTTGATATGTCGAAATGGCGTGCGTTAACGCTGCAGGAAAAAGCAG  >  W3110S.gb/3176973‑3177107

cAGATTACACCTATAGCAACGGCTACCGCGACGCGAACGGCATCAACTCTAACGCGACCAGTGCGTCCTTg                                                                  <  1:1969570/71‑1 (MQ=255)

cAGATTACACCTATAGCAACGGCTACCGCGACGCGAACGGCATCAACTCTAACGCGACCAGTGCGTCCTTg                                                                  <  1:3399446/71‑1 (MQ=255)

  gATTACACCTATAGCAACGGCTACCGCGACGCGAACGGCATCAACTCTAACGCGACCAGTGCGTCCTTg                                                                  >  1:3408446/1‑69 (MQ=255)

  gATTACACCTATAGCAACGGCTACCGCGACGCGAACGGCATCAACTCTAACGCGACCAGTGCGTCCTTg                                                                  >  1:3187572/1‑69 (MQ=255)

  gATTACACCTATAGCAACGGCTACCGCGACGCAAACGGCATCAACTCTAACGCGACCAGTGCGTCCTTg                                                                  >  1:3782987/1‑69 (MQ=255)

     tACACCTATAGCAACGGCTACCGCGACGCGAACGGCATCAACTCTAACGCGACCAGTGCGTCCTTGCAGt                                                              >  1:1672548/1‑70 (MQ=255)

     tACACCTATAGCAACGGCTACCGCGACGCGAACGGCATCAACTCTAACGCGACCAGTGCGTCCTTGCAGt                                                              >  1:367446/1‑70 (MQ=255)

      acacCTATAGCAACGGCTACCGCGACGCGAACGGCATCAACTCTAACGCGACCAGTGCGTCCTTGCAGtt                                                             >  1:2421575/1‑70 (MQ=255)

        acCTATAGCAACGGCTACCGCGACGCGAACGGCATCAACTCTAACGCGACCAGTGCGTCCTTGCAGTTAAc                                                          <  1:2953141/71‑1 (MQ=255)

        acCTATAGCAACGGCTACCGCGACGCGAACGGCATCAACTCTAACGCGACCAGTGCGTCCTTGCAGTTAAc                                                          <  1:3737837/71‑1 (MQ=255)

          cTATAGCAACGGCTACCGCGACGCGAACGGCATCAACTCTAACGCGACCAGTGCGTCCTTGCa                                                                >  1:1183699/1‑63 (MQ=255)

          cTATAGCAACGGCTACCGCGACGCGAACGGCATCAACTCTAACGCGACCAGTGCGTCCTTGCa                                                                >  1:1893208/1‑63 (MQ=255)

          cTATAGCAACGGCTACCGCGACGCGAACGGCATCAACTCTAACGCGACCAGTGCGTCCTTGCa                                                                >  1:3441396/1‑63 (MQ=255)

          cTATAGCAACGGCTACCGCGACGCGAACGGCATCAACTCTAACGCGACCAGTGCGTCCTTGCa                                                                >  1:3234179/1‑63 (MQ=255)

            ataGCAACGGCTACCGCGACGCGAACGGCATCAACTCTAACGCGACCAGTGCGTCCTTGCAGTTAACTCaa                                                      >  1:2518374/1‑71 (MQ=255)

             taGCAACGGCTACCGCGACGCGAACGGCATCAACTCTAACGCGACCAGTGCGTCCTTGCAGTTAACTCa                                                       <  1:1248573/69‑1 (MQ=255)

             taGCAACGGCTACCGCGACGCGAACGGCATCAACTCTAACGCGACCAGTGCGTCCTTGCAGTTAACTCAAt                                                     >  1:1869824/1‑71 (MQ=255)

             taGCAACGGCTACCGCGACGCGAACGGCATCAACTCTAACGCGACCAGTGCGTCCCTGCAGTTAACTCAAt                                                     >  1:3528972/1‑71 (MQ=255)

              aGCAACGGCTACCGCGACGCGAACGGCATCAACTCTAACGCGACCAGTGCGTCCTTGCAGTTAACTCAATc                                                    >  1:1956966/1‑71 (MQ=255)

                  aCGGCTACCGCGACGCGAACGGCATCAACTCTAACGCGACCAGTGCGTCCCTGCAGTTAACTCAATCCAtt                                                >  1:3304518/1‑71 (MQ=255)

                  aCGGCTACCGCGACGCGAACGGCATCAACTCTAACGCGACCAGTGCGTCCCTGCAGTTAACTCAATCCAtt                                                >  1:435382/1‑71 (MQ=255)

                    ggCTACCGCGACGCGAAAGGCATCAACTCTAACGCGACCAGTGCGTCCTTGCAGtt                                                             <  1:1535938/56‑1 (MQ=255)

                      cTACCGCGACGCGAACGGCATCAACTCTAACGCGACCAGTGCGTCCTTGCAGTTAACTCAATCCATTTTTg                                            <  1:1498855/71‑1 (MQ=255)

                       tACCGCGACGCGAACGGCATCAACTCTAACGCGACCAGTGCGTCCTTGCAGTTAACTCAATCCATTTTTGa                                           >  1:53776/1‑71 (MQ=255)

                       tACCGCGACGCGAACGGCATCAACTCTAACGCGACCAGTGCGTCCTTGCAGTTAACTCAATCCATTTTTGa                                           >  1:1217033/1‑71 (MQ=255)

                         ccGCGACGCGAACGGCATCAACTCTAACGCGACCAGTGCGTCCTTGCAGTTAACTCAATCCATTTTTGata                                         >  1:610586/1‑71 (MQ=255)

                         ccGCGACGCGAACGGCATCAACTCTAACGCGACCAGTGCGTCCTTGCAGTTAACTCAATCCATTTTTGata                                         >  1:351989/1‑71 (MQ=255)

                         ccGCGACGCGAACGGCATCAACTCTAACGCGACCAGTGCGTCCTTGCAGTTAACTCAATCCATTTTTGat                                          <  1:3573055/70‑1 (MQ=255)

                         ccGCGACGCGAACGGCATCAACTCTAACGCGACCAGTGCGTCCTTGCAGTTAACTCAATCCATTTTTGa                                           <  1:3455739/69‑1 (MQ=255)

                         ccGCGACGCGAACGGCATCAACTCTAACGCGACCAGTGCGTCCTTGCAGTTAACTCAATCCATTTTTGa                                           <  1:2965774/69‑1 (MQ=255)

                         ccGCGACGCGAACGGCATCAACTCTAACGCGACCAGTGCGTCCTTGCAGTTAACTCAATCCATTTTTGa                                           <  1:3642940/69‑1 (MQ=255)

                         ccGCGACGCGAACGGCATCAACTCTAACGCGACCAGTGCGTCCTTGCAGTTAACTCAATCCATTTTTGa                                           <  1:98842/69‑1 (MQ=255)

                         ccGCGACGCGAACGGCATCAACTCTAACGCGACCAGTGCGTCCTTGCAGTTAACTCAATCCATTTTTGa                                           <  1:1341327/69‑1 (MQ=255)

                         ccGCGACGCGAACGGCATCAACTCTAACGCGACCAGTGCGTCCTTGCAGTTAACTCAATCCATTTTTGa                                           <  1:3836101/69‑1 (MQ=255)

                          cgcgacgcgaACGGCATCAACTCTAACGCGACCAGTGCGTCCTTGCAGTTAACTCAATCCATTTTTGa                                           <  1:3147855/68‑1 (MQ=255)

                          cgcgacgcgaACGGCATCAACTCTAACGCGACCAGTGCGTCCCTGCAGTTAACTCAATCCATTTTTGatat                                        <  1:3144368/71‑1 (MQ=255)

                          cgcgacgcgaACGGCATCAACTCTAACGCGACCAGTGCGTCCCTGCAGTTAACTCAATCCATTTTTGatat                                        <  1:3164864/71‑1 (MQ=255)

                            cgacgcgaACGGCATCAACTCTAACGCGACCAGTGCGTCCTTGCAGTTAACTCAATCCATTTTTg                                            <  1:2102700/65‑1 (MQ=255)

                              acgcgaACGGCATCAACTCTAACGCGACCAGTGCGTCCTTGCAGTTAACTCAATCCATTTTTGATATGTCg                                    <  1:1331272/71‑1 (MQ=255)

                              acgcgaACGGCATCAACTCTAACGCGACCAGTGCGTCCTTGCAGTTAACTCAATCCATTTTTGATATGTCg                                    <  1:3862486/71‑1 (MQ=255)

                               cgcgaACGGCATCAACTCTAACGCGACCAGTGCGTCCTTGCAGTTAACTCAATCCATTTTTGATATGTCGa                                   <  1:1543/71‑1 (MQ=255)

                                gcgaACGGCATCAACTCTAACGCGACCAGTGCGTCCTTGCAGTTAACTCAATCCATTTTTGATATGTCGaa                                  <  1:3100563/71‑1 (MQ=255)

                                 cgaACGGCATCAACTCTAACGCGACCAGTGCGTCCTTGCAGTTAACTCAATCCATTTTTGATATGTCGaaa                                 <  1:3017490/71‑1 (MQ=255)

                                   aaCGGCATCAACTCTAACGCGACCAGTGCGTCCTTGCAGTTAACTCAATCCATTTTTGat                                          <  1:2385213/60‑1 (MQ=255)

                                   aaCGGCATCAACTCTAACGCGACCAGTGCGTCCCTGCAGTTAACTCAATCCATTTTTGatat                                        >  1:3211329/1‑62 (MQ=255)

                                   aaCGGCATCAACTCTAACGCGACCAGTGCGTCCCTGCAGTTAACTCAATCCATTTTTGatat                                        >  1:3695896/1‑62 (MQ=255)

                                   aaCGGCATCAACTCTAACGCGACCAGTGCGTCCCTGCAGTTAACTCAATCCATTTTTGatat                                        >  1:992836/1‑62 (MQ=255)

                                    aCGGCATCAACTCTAAGGCGACCAGTGCGTCCTTGCAGTTAACTCAATCCATTTTTGATGTGTAGAAATgg                              <  1:2795983/71‑1 (MQ=255)

                                    aCGGCATCAACTCTAACGCGACCAGTGCGTCCCTGCAGTTAACTCAATCCATTTTTGATATGTCGAAATgg                              >  1:305134/1‑71 (MQ=255)

                                    aCGGCATCAACTCTAACGCGACCAGTGCGTCCCTGCAGTTAACTCAATCCATTTTTGATATGTCGAAATgg                              >  1:416456/1‑71 (MQ=255)

                                       gCATCAACTCTAACGCGACCAGTGCGTCCTTGCAg                                                               >  1:1223725/1‑35 (MQ=255)

                                       gCATCAACTCTAACGCGACCAGTGCGTCCTTGCAg                                                               >  1:1903244/1‑35 (MQ=255)

                                         aTCAACTCTAACGCGACCAGTGCGTCCTTGCAGTTAACTCAATCCATTTTTGa                                           <  1:739102/53‑1 (MQ=255)

                                          tCAACTCTAACGCGACCAGTGCGTCCTTGCAGTTAACTCa                                                       <  1:530432/40‑1 (MQ=255)

                                          tCAACTCTAACGCGACCAGTGCGTCCTTGCAGTTAACTCa                                                       <  1:216331/40‑1 (MQ=255)

                                          tCAACTCTAACGCGACCAGTGCGTCCTTGCAGTTAACTCa                                                       <  1:2332855/40‑1 (MQ=255)

                                          tCAACTCTAACGCGACCAGTGCGTCCTTGCAGTTAACTCAATCCATTTTTGat                                          <  1:3568471/53‑1 (MQ=255)

                                          tCAACTCTAACGCGACCAGTGCGTCCTTGCAGTTAACTCAATCCATTTTTGa                                           <  1:1151590/52‑1 (MQ=255)

                                            aaCTCTAACGCGACCAGTGCGTCCTTGCAGTTAACTCAATCCATTTTTGatat                                        >  1:2343725/1‑53 (MQ=255)

                                            aaCTCTAACGCGACCAGTGCGTCCTTGCAGTTAACTCAATCCATTTTTGATATGTCGaaa                                 >  1:349591/1‑60 (MQ=255)

                                            aaCTCTAACGCGACCAGTGCGTCCTTGCAGTTAACTCAATCCATTTTTGATATGTCGAAATGgcgtgcgt                       >  1:3061099/1‑70 (MQ=255)

                                            aaCTCTAACGCGACCAGTGCGTCCTTGCAGTTAACTCAATCCATTTTTGATATGTCGAAATGgcgtgcgt                       >  1:1553963/1‑70 (MQ=255)

                                            aaCTCTAACGCGACCAGTGCGTCCTTGCAGTTAACTCAATCCATTTTTGATATGTCGAAATGgcgtgcgt                       >  1:3180562/1‑70 (MQ=255)

                                            aaCTCTAACGCGACCAGTGCGTCCTTGCAGTTAACTCAATCCATTTTTGATATGTCGAAATGGCGTGCGtt                      <  1:1246481/71‑1 (MQ=255)

                                            aaCTCTAACGCGACCAGTGCGTCCTTGCAGTTAACTCAATCCATTTTTGATATGTCGAAATGGCGTGCGtt                      <  1:1507302/71‑1 (MQ=255)

                                            aaCTATAACGCGACCAGTGCGTCCTTGCAGTTAACTCAATCCATTTTTGATATGTCGAAATGGCGTGCGtt                      <  1:3033276/71‑1 (MQ=255)

                                              ctctAACGCGACCAGTGCGTCCTTGCAGTTAACTCAATCCATTTTTGATATGTCGaa                                  >  1:1990848/1‑57 (MQ=255)

                                              ctctAACGCGACCAGTGCGTCCTTGCAGTTAACTCAATCCATTTTTGATATGTCGa                                   >  1:2662868/1‑56 (MQ=255)

                                               tctAACGCGACCAGTGCGTCCTTGCAGTTAACTCAATCCATTTTTGATATGTCGAAATGGCGTGCGTTaa                    >  1:1235613/1‑70 (MQ=255)

                                               tctAACGCGACCAGTGCGTCCTTGCAGTTAACTCAATCCATTTTTGATATGTCGAAATGGCGTGCGTTAAc                   >  1:3794189/1‑71 (MQ=255)

                                                ctAACGCGACCAGTGCGTCCTTGCAGTTAACTCAATCCATTTTTGATATGTCGaaa                                 >  1:1282272/1‑56 (MQ=255)

                                                 tAACGCGACCAGTGCGTCCTTGCAGTTAACTCAATCCATTTTTGa                                           <  1:285852/45‑1 (MQ=255)

                                                 tAACGCGACCAGTGCGTCCTTGCAGTTAACTCAATCCATTTTTGa                                           <  1:1490432/45‑1 (MQ=255)

                                                 tAACGCGACCAGTGCGTCCTTGCAGTTAACTCAATCCATTTTTGATATGTCGAAATGGCGTGCGTTAAc                   <  1:512106/69‑1 (MQ=255)

                                                 tAACGCGACCAGTGCGTCCTTGCAGTTAACTCAATCCATTTTTGATATGTCGAAATGGCGTGCGTTAAc                   <  1:3439622/69‑1 (MQ=255)

                                                  aaCGCGACCAGTGCGTCCTTGCAGTTAACTCAATCCATTTTTGata                                         >  1:2828837/1‑46 (MQ=255)

                                                   aCGCGACCAGTGCGTCCTTGCAGTTAACTCAATcc                                                   >  1:2915352/1‑35 (MQ=255)

                                                   aCGCGACCAGTGCGTCCTTGCAGTTAACTCAATCCATTTTTGata                                         >  1:3753599/1‑45 (MQ=255)

                                                   aCGCGACCAGTGCGTCCTTGCAGTTAACTCAATCCATTTTTGata                                         >  1:329641/1‑45 (MQ=255)

                                                   aCGCGACCAGTGCGTCCTTGCAGTTAACTCAATCCATTTTTGata                                         >  1:2929054/1‑45 (MQ=255)

                                                   aCGCGACCAGTGCGTCCTTGCAGTTAACTCAATCCATTTTTGata                                         >  1:3386687/1‑45 (MQ=255)

                                                   aCGCGACCAGTGCGTCCTTGCAGTTAACTCAATCCATTTTTGATATGTCGAAATGgcgtgcgt                       >  1:1189867/1‑63 (MQ=255)

                                                     gcgACCAGTGCGTCCTTGCAGTTAACTCAATCCATTTTTGATATGTCg                                    <  1:3596516/48‑1 (MQ=255)

                                                     gcgACCAGTGCGTCCTTGCAGTTAACTCAATCCATTTTTGATATGTCGAAATGGCGTGCGTTAACGCTGc              <  1:2143040/70‑1 (MQ=255)

                                                         ccAGTGCGTCCTTGCAGTTAACTCAATCCATTTTTGATATGTCGAAATGGCGTGCGTTAACGCTGCAGGaa         <  1:3276566/71‑1 (MQ=255)

                                                         ccAGTGCGTCCTTGCAGTTAACTCAATCCATTTTTGATATGTCGAAATGGCGTGCGTTAACGCTGCAGGaa         <  1:573047/71‑1 (MQ=255)

                                                            gTGCGTCCTTGCAGTTAACTCAATCCATTTTTGATATGTCGAAATGgcgtgcg                        <  1:1337890/53‑1 (MQ=255)

                                                            gTGCGTCCTTGCAGTTAACTCAATCCATTTTTGATATGTCGAAATGgcgtgcg                        <  1:639271/53‑1 (MQ=255)

                                                             tGCGTCCTTGCAGTTAACTCAATCCATTTTTGATATGTCGa                                   <  1:2814642/41‑1 (MQ=255)

                                                             tGCGTCCTTGCAGTTAACTCAATCCATTTTTGATATGTCGAAATGgc                             <  1:2188825/47‑1 (MQ=255)

                                                             tGCGTCCTTGCAGTTAACTCAATCCATTTTTGATATGTCGAAATGgc                             <  1:3220172/47‑1 (MQ=255)

                                                                gTCCTTGCAGTTAACTCAATCCATTTTTGATATGTCGAAATGGCGTGCGTTAACGCTGCAGGAAAAagcag  <  1:1377044/71‑1 (MQ=255)

                                                                gTCCTTGCAGTTAACTCAATCCATTTTTGATATGTCGAAATGGCGTGCGTTAACGCTGCAGGAAAAagcag  <  1:1650309/71‑1 (MQ=255)

                                                                gTCCTTGCAGTTAACTCAATCCATTTTTGATATGTCGAAATGGCGTGCGTTAACGCTGCAGGAAAAagcag  <  1:84866/71‑1 (MQ=255)

                                                                gTCCTTGCAGTTAACTCAATCCATTTTTGATATGTCGAAATGGCGTGCGTTAACGCTGCAGGAAAAagcag  <  1:925682/71‑1 (MQ=255)

                                                                 tCCTTGCAGTTAACTCAATCCATTTTTGATATGTCGa                                   <  1:761635/37‑1 (MQ=255)

                                                                   cTTGCAGTTAACTCAATCCATTTTTGATATGTCGAAATGGCGTGCGTTAACGCTGCAGGAAAAagcag  >  1:57713/1‑68 (MQ=255)

CAGATTACACCTATAGCAACGGCTACCGCGACGCGAACGGCATCAACTCTAACGCGACCAGTGCGTCCTTGCAGTTAACTCAATCCATTTTTGATATGTCGAAATGGCGTGCGTTAACGCTGCAGGAAAAAGCAG  >  W3110S.gb/3176973‑3177107

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