BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	742,687	T→A	I98F (ATC→TTC)	ybgH ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	742,687	0	T	A	98.3%	164.9 / ‑5.7	60	I98F (ATC→TTC)	ybgH	predicted transporter
Reads supporting (aligned to +/- strand): ref base T (1/0); new base A (52/7); total (53/7)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACAGGCCATAGCCGCAGACGATAATCGCCAGGGACAGATAGAGGAATGACGGATGGATCTCACTGGCACCCAGCACCACATGACCGATCGCCATCAACAACGCCCCCAGCATCACCGCCAT  >  W3110S.gb/742629‑742750

aaCAGGCCATAGCCGCAGACGATAATCGCCAGGGACAGATAGAGGAATGACGGATGGAACTCACTGGCAc                                                      <  1:3086309/70‑1 (MQ=255)

           gCCGCAGACGATAATCGCCAGGGACAGATAGAGGAATGACGGATGGAACTCACTGGCAcc                                                     >  1:2967772/1‑60 (MQ=255)

           gCCGCAGACGATAATCGCCAGGGACAGATAGAGGAATGACGGATGGAACTCACTGGCAcc                                                     >  1:703372/1‑60 (MQ=255)

           gCCGCAGACGATAATCGCCAGGGACAGATAGAGGAATGACGGATGGAACTCACTGGCAcc                                                     >  1:2444265/1‑60 (MQ=255)

           gCCGCAGACGATAATCGCCAGGGACAGATAGAGGAATGACGGATGGAACTCACTGGCAcc                                                     >  1:1278513/1‑60 (MQ=255)

           gCCGCAGACGATAATCGCCAGGGACAGATAGAGGAATGACGGATGGAACTCACTGGCAcc                                                     >  1:1343022/1‑60 (MQ=255)

           gCCGCAGACGATAATCGCCAGGGACAGATAGAGGAATGACGGATGGAACTCACTGGCAcc                                                     >  1:3225024/1‑60 (MQ=255)

           gCCGCAGACGATAATCGCCAGGGACAGATAGAGGAATGACGGATGGAACTCACTGGCAcc                                                     >  1:3633112/1‑60 (MQ=255)

           gCCGCAGACGATAATCGCCAGGGACAGATAGAGGAATGACGGATGGAACTCACTGGCAcc                                                     >  1:4042423/1‑60 (MQ=255)

           gCCGCAGACGATAATCGCCAGGGACAGATAGAGGAATGACGGATGGAACTCACTGGCAcc                                                     >  1:4069588/1‑60 (MQ=255)

           gCCGCAGACGATAATCGCCAGGGACAGATAGAGGAATGACGGATGGAACTCACTGGCAc                                                      >  1:293999/1‑59 (MQ=255)

             cGCAGACGATAATCGCCAGGGACAGATAGATGAATGACGGATGGAACTCACTGGCACCCAGCACCACATGa                                        >  1:2820753/1‑71 (MQ=255)

             cGCAGACGATAATCGCCAGGGACAGATAGAGGAATGACGGATGGAACTCACTGGCACCCAGCACCACATGa                                        >  1:4082277/1‑71 (MQ=255)

             cGCAGACGATAATCGCCAGGGACAGATAGAGGAATGACGGATGGAACTCACTGGCACCCAGCACCACATGa                                        >  1:2254127/1‑71 (MQ=255)

             cGCAGACGATAATCGCCAGGGACAGATAGAGGAATGACGGATGGAACTCACTGGCACCCAGCACCACATGa                                        >  1:2013336/1‑71 (MQ=255)

             cGCAGACGATAATCGCCAGGGACAGATAGAGGAATGACGGATGGAACTCACTGGCACCCAGCACCACATGa                                        >  1:2645578/1‑71 (MQ=255)

             cGCAGACGATAATCGCCAGGGACAGATAGAGGAATGACGGATGGAACTCACTGGCACCCAGCACCACATGa                                        >  1:3503828/1‑71 (MQ=255)

             cGCAGACGATAATCGCCAGGGACAGATAGAGGAATGACGGATGGAACTCACTGGCACCCAGCACCACATGa                                        >  1:1381034/1‑71 (MQ=255)

             cGCAGACGATAATCGCCAGGGACAGATAGAGGAATGACGGATGGAACTCACTGGCACCCAGCACCACATGa                                        >  1:3672515/1‑71 (MQ=255)

             cGCAGACGATAATCGCCAGGGACAGATAGAGGAATGACGGATGGAACTCACTGGCACCCAGCACCACATGa                                        >  1:41425/1‑71 (MQ=255)

             cGCAGACGATAATCGCCAGGGACAGATAGAGGAATGACGGATGGAACTCACTGGCACCCAGCACCACATGa                                        >  1:1457699/1‑71 (MQ=255)

                   cGATAATCGCCAGGGACAGATAGAGGAATGACGGATGGAACTCACTGGCACCCAGCACCACATGACCGATc                                  <  1:1668413/71‑1 (MQ=255)

                    gATAATCGCCAGGGACAGATAGAGGAATGACGGATGGAACTCACTGGCACCCAGCACCACATGACCGATc                                  <  1:458870/70‑1 (MQ=255)

                     aTAATCGCCAGGGACAGATAGAGGAATGACGGATGGAACTCACTGGCACCCAGCACCACATGACCGATCGc                                <  1:1173857/71‑1 (MQ=255)

                          cGCCAGGGACAGATAGAGGAATGACGGATGGATCTCACTGGCACCCAGCACCACATGAc                                       >  1:210631/1‑59 (MQ=255)

                          cGCCAGGGACAGATAGAGGAATGACGGATGGAACTCACTGGGACCCAGCACCACATGACCGATCGCCAt                             >  1:3519542/1‑69 (MQ=255)

                          cGCCAGGGACAGATAGAGGAATGACGGATGGAACTCACTGGCACCCAGCACCACATGACCGATCGc                                >  1:1381567/1‑66 (MQ=255)

                          cGCCAGGGACAGATAGAGGAATGACGGATGGAACTCACTGGCACCCAGCACCACATGACCGATCGCCAt                             >  1:1087231/1‑69 (MQ=255)

                          cGCCAGGGACAGATAGAGGAATGACGGATGGAACTCACTGGCACCCAGCACCACATGACCGATCGCCAt                             >  1:2638054/1‑69 (MQ=255)

                          cGCCAGGGACAGATAGAGGAATGACGGATGGAACTCACTGGCACCCAGCACCACATGACCGATCGCCAt                             >  1:421632/1‑69 (MQ=255)

                          cGCCAGGGACAGATAGAGGAATGACGGATGGAACTCACTGGCACCCAGCACCACATGACCGATCGCCAt                             >  1:515087/1‑69 (MQ=255)

                          cGCCAGGGACAGATAGAGGAATGACGGATGGAACTCACTGGCACCCAGCACCACATGACCGATCGCCAt                             >  1:970673/1‑69 (MQ=255)

                          cGCCAGGGACAGATAGAGGAATGACGGATGGAACTCACTGGCACCCAGCACCACATGACCGATCGCCAt                             >  1:1515287/1‑69 (MQ=255)

                          cGCCAGGGACAGATAGAGGAATGACGGATGGAACTCACTGGCACCCAGCACCACATGACCGATCGCCAt                             >  1:1366049/1‑69 (MQ=255)

                          cGCCAGGGACAGATAGAGGAATGACGGATGGAACTCACTGGCACCCAGCACCACATGACCGATCGCCAt                             >  1:1279130/1‑69 (MQ=255)

                          cGCCAGGGACAGATAGAGGAATGACGGATGGAACTCACTGGCACCCAGCACCACATGACCGATCGCCAt                             >  1:160080/1‑69 (MQ=255)

                          cGCCAGGGACAGATAGAGGAATGACGGATGGAACTCACTGGCACCCAGCACCACATGACCGATCGCCAt                             >  1:165944/1‑69 (MQ=255)

                          cGCCAGGGACAGATAGAGGAATGACGGATGGAACTCACTGGCACCCAGCACAACATGACCGATCGCCAt                             >  1:1206915/1‑69 (MQ=255)

                                         gagGAATGACGGATGGAACTCACTGGCACCCAGCACCACATGACCGATCGcc                               >  1:2972628/1‑52 (MQ=255)

                                         gagGAATGACGGATGGAACTCACTGGCACCCAGCACCACATGACCGATCGCCATCAACAACGCCCCCAGc             >  1:1506502/1‑70 (MQ=255)

                                         gagGAATGACGGATGGAACTCACTGGCACCCAGCACCACATGACCGATCGCCATCAACAACGCCCCCAGCa            >  1:2743051/1‑71 (MQ=255)

                                         gagGAATGACGGATGGAACTCACTGGCACCCAGCACCACATGACCGATCGCCATCAACAACGCCCCCAGCa            >  1:845775/1‑71 (MQ=255)

                                         gagGAATGACGGATGGAACTCACTGGCACCCAGCACCACATGACCGATCGCCATCAACAACGCCCCCAGCa            >  1:1192474/1‑71 (MQ=255)

                                         gagGAATGACGGATGGAACTCACTGGCACCCAGCACCACATGACCGATCGCCATCAACAACGCCCCCAGCa            >  1:4026188/1‑71 (MQ=255)

                                         gagGAATGACGGATGGAACTCACTGGCACCCAGCACCACATGACCGATCGCCATCAACAACGCCCCCAGCa            >  1:1988230/1‑71 (MQ=255)

                                           ggAATGACGGATGGAACTCACTGGCACCCAGCACCACATGACCGATCGCCATCAACAACGCCCCCAGCAt           >  1:1727664/1‑70 (MQ=255)

                                           ggAATGACGGATGGAACTCACTGGCACCCAGCACCACATGACCGATCGCCATCAACAACGCCCCCAGCAt           >  1:1455464/1‑70 (MQ=255)

                                           ggAATGACGGATGGAACTCACTGGCACCCAGCACCACATGACCGATCGCCATCAACAACGCCCCCAGCAt           >  1:3813260/1‑70 (MQ=255)

                                           ggAATGACGGATGGAACTCACTGGCACCCAGCACCACATGACCGATCGCCATCAACAACGCCCCCAGCAt           >  1:1647216/1‑70 (MQ=255)

                                           ggAATGACGGATGGAACTCACTGGCACCCAGCACCACATGACCGATCGCCATCAACAACGCCCCCAGCAt           >  1:3517179/1‑70 (MQ=255)

                                           ggAATGACGGATGGAACTCACTGGCACCCAGCACCACATGACCGATCGCCATCAACAACGCCCCCAGCAt           >  1:3335710/1‑70 (MQ=255)

                                           ggAATGACGGATGGAACTCACTGGCACCCAGCACCACATGACCGATCGCCATCAACAACGCCCCCAGCAt           >  1:551650/1‑70 (MQ=255)

                                           ggAATGACGGATGGAACTCACTGGCACCCAGCACCACATGACCGATCGCCATCAACAACGCCCCCAGCAt           >  1:1944085/1‑70 (MQ=255)

                                           ggAATGACGGATGGAACTCACTGGCACCCAGCACCACATGACCGATCGCCATCAACAACGCCCCCAGCAt           >  1:559523/1‑70 (MQ=255)

                                           ggAATGACGGATGGAACTCACTGGCACCCAGCACCACATGACCGATCGCCATCAACAACGCCCCCAGCAt           >  1:660076/1‑70 (MQ=255)

                                           ggAATGACGGATGGAACTCACTGGCACCCAGCACCACATGACCGATCGCCATCAACAACGCCCCCAGCAt           >  1:2335241/1‑70 (MQ=255)

                                           ggAATGACGGATGGAACTCACTGGCACCCAGCACCACATGACCGATCGCCATCAACAACGCCCCCAGCAt           >  1:2580344/1‑70 (MQ=255)

                                             aaTGACGTATGGAACTCACTGGCACCCAGCACCACATGACCGATCGCCATCAACAACGCCCCCAGCATCAc        <  1:1033708/71‑1 (MQ=255)

                                             aaTGACGGATGGAACTCACTGGCACCCAGCACCACATGACCGATCGCCATCAACAACGCCCCCAGCATCAc        <  1:3707789/71‑1 (MQ=255)

                                             aaTGACGGATGGAACTCACTGGCACCCAGCACCACATGACCGATCGCCATCAACAACGCCCCCAGCATCAc        <  1:596944/71‑1 (MQ=255)

                                                   ggatggaaCTCACTGGCACCCAGCACCACATGACCGATCg                                 >  1:3850599/1‑40 (MQ=255)

                                                   ggatggaaCTCACTGGCACCCAGCACCACATGACCGATCGCCATcaacaa                       >  1:364824/1‑50 (MQ=255)

                                                   ggatggaaCTCACTGGCACCCAGCACCACATGACCGATCGCCATcaacaa                       >  1:2042179/1‑50 (MQ=255)

                                                   ggatggaaCTCACTGGCACCCAGCACCACATGACCGATCGCCATcaacaa                       >  1:2248288/1‑50 (MQ=255)

                                                   ggatggaaCTCACTGGCACCCAGCACCACATGACCGATCGCCATcaacaa                       >  1:2936824/1‑50 (MQ=255)

                                                   ggatggaaCTCACTGGCACCCAGCACCACATGACCGATCGCCATcaacaa                       >  1:900064/1‑50 (MQ=255)

                                                   ggatggaaCTCACTGGCACCCAGCACCACATGACCGATCGCCATcaacaa                       >  1:2628018/1‑50 (MQ=255)

                                                   ggatggaaCTCACTGGCACCCAGCACCACATGACCGATCGCCATCAACAACGCCCCCAGCATCACCGCCAt  <  1:1496395/71‑1 (MQ=255)

                                                   ggatggaaCTCACTGGCACCCAGCACCACATGACCGATCGCCATCAACAACGCCCCCAGCATCACCGCCAt  <  1:3681887/71‑1 (MQ=255)

                                                   ggatggaaCTCACTGGCACCCAGCACCACATGACCGATCGCCATCAACAACGCCCCCAGCATCACCGCCAt  <  1:3285143/71‑1 (MQ=255)

                                                   ggatggaaCTCACTGGCACCCAGCACCACATGACCGATCGCCATCAACAACGCCCCCAGCATCACCGCCAt  <  1:2317497/71‑1 (MQ=255)

AACAGGCCATAGCCGCAGACGATAATCGCCAGGGACAGATAGAGGAATGACGGATGGATCTCACTGGCACCCAGCACCACATGACCGATCGCCATCAACAACGCCCCCAGCATCACCGCCAT  >  W3110S.gb/742629‑742750

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