BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	1,861,208	C→T	12.7%	G246S (GGT→AGT)	ydjK ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	1,861,208	0	C	T	12.7%	151.6 / 11.9	55	G246S (GGT→AGT)	ydjK	predicted transporter
Reads supporting (aligned to +/- strand): ref base C (17/31); new base T (3/4); total (20/35)
Fisher's exact test for biased strand distribution p-value = 6.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.88e-01

AATATCACGCGTTTCAGTAATACTCCTGTCAGTAAGGCTGAATACGGCACCGCTTGTGGCGCTTTACCGTCATCAGCAATAACCACAGGCGGCAAAGGTTTACCCGTCTGGCGTATGACGCCTTCTTCTA  >  W3110S.gb/1861141‑1861270

aaTATCACGCGTTTCAGTAATACTCCTGTCAGTAAGGCTGAATACGGCACCGCTTGTGGCGCTTTACTGt                                                              >  1:871562/1‑70 (MQ=255)

 atatCACGCGTTTCAGTAATACTCCTGTCAGTAAGGCTGAATACGGCACCGCTTGTGGCGCTTTACCGtca                                                            >  1:2408286/1‑71 (MQ=255)

 atatCACGCGTTTCAGTAATACTCCTGTCAGTAAGGCTGAATACGGCACCGCTTGTGGCGCTTTACCGtca                                                            >  1:3065125/1‑71 (MQ=255)

     cACGCGTTTCAGTAATACTCCTGTCAGTAAGGCTGAATACGGCACCGCTTGTGGCGCTTTACCGTCATCAg                                                        >  1:1685246/1‑71 (MQ=255)

        gcgTTTCAGTAATACTCCTGTCAGTAAGGCTGAATACGGCACCGCTTGTGGCGCTTTACCGTCATCAGCa                                                      <  1:3311741/70‑1 (MQ=255)

         cgTTTCAGTAATACTCCTGTCAGTAAGGCTGAATACGGCACCGCTTGTGGCGCTTTACTGTCATCAGCAAt                                                    <  1:328959/71‑1 (MQ=255)

         cgTTTCAGTAATACTCCTGTCAGTAAGGCTGAATACGGCACCGCTTGTGGCGCTTTACTGTCATCAGCAAt                                                    <  1:870776/71‑1 (MQ=255)

         cgTTTCAGTAATACTCCTGTCAGTAAGGCTGAATACGGCACCGCTTGTGGCGCTTTACCGTCATCAGCAAt                                                    <  1:3860867/71‑1 (MQ=255)

         cgTTTCAGTAATACTCCTGTCAGTAAGGCTGAATACGGCACCGCTTGTGGCGCTTTACCGTCATCAGCAAt                                                    <  1:2295949/71‑1 (MQ=255)

         cgTTTCAGTAATACTCCTGTCAGTAAGGCTGAATACGGCACCGCTTGTGGCGCTTTACCGTCATCAGCAAt                                                    <  1:1863645/71‑1 (MQ=255)

          gTTTCAGTAATACTCCTGTCAGTAAGGCTGAATACGGCACCGCTTGTGGCGCTTTACCGTCATCAGCAAt                                                    <  1:1913701/70‑1 (MQ=255)

          gTTTCAGTAATACTCCTGTCAGTAAGGCTGAATACGGCACCGCTTGTGGCGCTTTACCGTCATCAGCAATa                                                   <  1:266748/71‑1 (MQ=255)

             tCAGTAATACTCCTGTCAGTAAGGCTGAATACGGCACCGCTTGTGGCGCTTTACCGt                                                              >  1:3673192/1‑57 (MQ=255)

             tCAGTAATACTCCTGTCAGTAAGGCTGAATACGGCACCGCTTGTGGCGCTTTACCGt                                                              >  1:2437716/1‑57 (MQ=255)

             tCAGTAATACTCCTGTCAGTAAGGCTGAATACGGCACCGCTTGTGGCGCTTTACCGTCATCAGCAATAAc                                                 >  1:89915/1‑70 (MQ=255)

             tCAGTAATACTCCTGTCAGTAAGGCTGAATACGGCACCGCTTGTGGCGCTTTACCGTCATCAGCAATAAc                                                 >  1:2308268/1‑70 (MQ=255)

             tCAGTAATACTCCTGTCAGTAAGGCTGAATACGGCACCGCTTGTGGCGCTTTACCGTCATCAGCAATAAc                                                 >  1:3609692/1‑70 (MQ=255)

             tCAGTAATACTCCTGTCAGTAAGGCTGAATACGGCACCGCTTGTGGCGCTTTACC‑TCATCAGCAATAACca                                               >  1:2221686/1‑71 (MQ=255)

                   aTACTCCTGTAAGTAAGGCTGAATAAGGCACCGCTTGTGGCGCTTTACCGTGATCAGCAATAACCACAggc                                          >  1:4026397/1‑71 (MQ=255)

                      cTCCTGTCAGTAAGGCTGAATACGGCACCGCTTGTGGCGCTTTACCGTCATCAGCAATAACCACAggcgg                                        >  1:3936319/1‑70 (MQ=255)

                      cTCCTGTCAGTAAGGCTGAATACGGCACCGCTTGTGGCGCTTTACCGTCATCAGCAATAACCACAggcgg                                        >  1:688452/1‑70 (MQ=255)

                      cTCCTGTCAGTAAGGCTGAATACGGCACCGCTTGTGGCGCTTTACCGTCATCAGCAATAACCACAggcgg                                        >  1:1107772/1‑70 (MQ=255)

                         cTGTCAGTAAGGCTGAATACGGCACCGCTTGTGGCGCTTTACCGTCATCAGCAATAACCACAGGCGGCa                                      >  1:1923088/1‑69 (MQ=255)

                              aGTAAGGCTGAATACGGCACCGCTTGTGGCGCTTTACCGTCATCAGCAATAACCACAGGCGGCAAAGGtt                                <  1:578791/70‑1 (MQ=255)

                              aGTAAGGCTGAATACGGCACCGCTTGTGGCGCTTTACCGTCATCAGCAATAACCACAGGCGGCAAAGGtt                                <  1:2659322/70‑1 (MQ=255)

                              aGTAAGGCTGAATACGGCACCGCTTGTGGCGCTTTACCGTCATCAGCAATAACCACAGGCGGCAAAGGtt                                <  1:840473/70‑1 (MQ=255)

                              aGTAAGGCTGAATACGGCACCGCTTGTGGCGCTTTACCGTCATCAGCAATAACCACAGGCGGCAAAGGtt                                <  1:1127580/70‑1 (MQ=255)

                              aGTAAGGCTGAATACGGCACCGCTTGTGGCGCTTTACCGTCATCAGCAATAACCACAGGCGGCAAAGGtt                                <  1:1084095/70‑1 (MQ=255)

                               gTAAGGCTGAATACGGCACCGCTTGTGGCGCTTTACCGTCATCAGCAATAACCACAGGCGGCAAAGGTTTa                              <  1:276979/71‑1 (MQ=255)

                               gTAAGGCTGAATACGGCACCGCTTGTGGCGCTTTACCGTCATCAGCAATAACCACAGGCGGCAAAGGTTTa                              <  1:1410205/71‑1 (MQ=255)

                                  aGGCTGAATACGGCACCGCTTGTGGCGCTTTACCGTCATCAGCAATAACCACAGGCGGCAAAGGTTTAccc                           >  1:3956208/1‑71 (MQ=255)

                                    gCTGAATACGGCACCGCTTGTGGCGCTTTACCGTCATCAGCAATAACca                                               <  1:1183760/49‑1 (MQ=255)

                                        aaTACGGCACCGCTTGTGGCGCTTTACCGTCATCAGCAATAACCACAGGCGGCAAAGGTTTACCCGTCt                       <  1:91086/69‑1 (MQ=255)

                                             ggCACCGCTTGTGGCGCTTTACTGTCATCAGCAATAACCACAGGCGGCAAAGGTTTACCCGTCTGGCGTAt                <  1:3235290/71‑1 (MQ=255)

                                             ggCACCGCTTGTGGCGCTTTACTGTCATCAGCAATAACCACAGGCGGCAAAGGTTTACCCGTCTGGCGTAt                <  1:830669/71‑1 (MQ=255)

                                             ggCACCGCTTGTGGCGCTTTACCGTCATCAGCAATAACCACAGGCGGCAAAGGTTTACCCGTCTGGCGTAt                <  1:3942387/71‑1 (MQ=255)

                                             ggCACCGCTTGTGGCGCTTTACCGTCATCAGCAATAACCACAGGCGGCAAAGGTTTACCCGTCTGGCGTAt                <  1:1141567/71‑1 (MQ=255)

                                             ggCACCGCTTGTGGCGCTTTACCGTCATCAGCAATAACCACAGGCGGCAAAGGTTTACCCGTCTGGCGTAt                <  1:701500/71‑1 (MQ=255)

                                             ggCACCGCTTGTGGCGCTTTACCGTCATCAGCAATAACCACAGGCGGCAAAGGTTTACCCATCTGGCGTAt                <  1:1301412/71‑1 (MQ=255)

                                              gCACCGCTTGTGGCGCTTTACCGTCATCAGCAATAACCACAGGCGGCAAAGGTTTACCCGTCTGGCGTAt                <  1:2814547/70‑1 (MQ=255)

                                              gCACCGCTTGTGGCGCTTTACCGTCATCAGCAATAACCACAGGCGGCAAAGGTTTACCCGTCTGGCGTAt                >  1:874420/1‑70 (MQ=255)

                                                 ccGCTTGTGGCGCTTTACCGTCATCAGCAATAACca                                               <  1:1307718/36‑1 (MQ=255)

                                                 ccGCTTGTGGCGCTTTACCGTCATCAGCAATAACca                                               <  1:2897140/36‑1 (MQ=255)

                                                 ccGCTTGTGGCGCTTTACCGTCATCAGCAATAACca                                               <  1:834539/36‑1 (MQ=255)

                                                 ccGCTTGTGGCGCTTTACCGTCATCAGCAATAACca                                               <  1:493096/36‑1 (MQ=255)

                                                 ccGCTTGTGGCGCTTTACCGTCATCAGCAATAACCACAGGCGGCAAAGGttt                               <  1:221477/52‑1 (MQ=255)

                                                 ccGCTTGTGGCGCTTTACCGTCATCAGCAATAACCACAGGCGGCAAAGGTTTACCCGTCTGGCGTATGACg            <  1:1934962/71‑1 (MQ=255)

                                                   gCTTGTGGCGCTTTACCGTCATCAGCAATAACCACAGGCGGCAAAGGTTTACCCGTCTGGCGTATGACGcc          >  1:3644469/1‑71 (MQ=255)

                                                       gtgGCGCTTTACCGTCATCAGCAATAACCACAGGCGGCAAAGGTTTACCCGTCTGGCGTATGACGCcttc       <  1:435144/70‑1 (MQ=255)

                                                          gcgcTTTACCGTCATCAGCAATAACCACAGGCGGCAAAGGTTTACCCGTCTGGCGTATGACGCCttcttc    <  1:1070385/70‑1 (MQ=255)

                                                          gcgcTTTACCGTCATCAGCAATAACCACAGGCGGCAAAGGTTTACCCGTCTGGCGTATGACGCCttcttc    <  1:264278/70‑1 (MQ=255)

                                                           cgcTTTACCGTCATCAGCAATAACCACAGGCGGCAAAGGTTTACCCGTCTGGCGTATGACGCCTTCTTCta  <  1:1792162/71‑1 (MQ=255)

                                                           cgcTTTACCGTCATCAGCAATAACCACAGGCGGCAAAGGTTTACCCGTCTGGCGTATGACGCCTTCTTCta  <  1:600240/71‑1 (MQ=255)

                                                              tttACTGTCATCAGCAATAACCACAGGCGGCAAAGGTTTAc                             >  1:1145125/1‑41 (MQ=255)

                                                              tttACTGTCATCAGCAATAACCACAGGCGGCAAAGGTTTAc                             >  1:2185266/1‑41 (MQ=255)

AATATCACGCGTTTCAGTAATACTCCTGTCAGTAAGGCTGAATACGGCACCGCTTGTGGCGCTTTACCGTCATCAGCAATAACCACAGGCGGCAAAGGTTTACCCGTCTGGCGTATGACGCCTTCTTCTA  >  W3110S.gb/1861141‑1861270

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