BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	1,210,013	T→C	10.3%	L102P (CTG→CCG)	ymfS →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	1,210,013	0	T	C	10.3%	178.7 / 11.9	68	L102P (CTG→CCG)	ymfS	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (26/35); new base C (4/3); total (30/38)
Fisher's exact test for biased strand distribution p-value = 6.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

TGATTTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCG  >  W3110S.gb/1209948‑1210077

tGATTTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGaa                                                               <  1:169902/69‑1 (MQ=255)

tGATTTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGaa                                                               <  1:1388655/69‑1 (MQ=255)

tGATTTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGaa                                                               <  1:3301529/69‑1 (MQ=255)

tGATTTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGa                                                                >  1:2198293/1‑68 (MQ=255)

            tAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAAtt                                                    <  1:1897864/68‑1 (MQ=255)

            tAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAAtt                                                    <  1:2515022/68‑1 (MQ=255)

            tAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAAtt                                                    <  1:3138732/68‑1 (MQ=255)

            tAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAAtt                                                    <  1:3134390/68‑1 (MQ=255)

                aGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAAtt                                                    <  1:3101241/64‑1 (MQ=255)

                     gTGGCAGGGTAAGGCTGCATT‑GGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGAt                                         <  1:527818/69‑1 (MQ=255)

                        gCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTg                                            >  1:1634873/1‑64 (MQ=255)

                        gCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTa                                       <  1:2223429/69‑1 (MQ=255)

                           gggTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAatat                                                        >  1:918885/1‑49 (MQ=255)

                             gTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAAGAATATAATTTGTGGCTGGATTATCt                                    >  1:2090922/1‑67 (MQ=255)

                             gTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCt                                    >  1:2399818/1‑67 (MQ=255)

                              tAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACa                                                            <  1:208281/42‑1 (MQ=255)

                              tAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTg                                            >  1:1437464/1‑58 (MQ=255)

                              tAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTg                                            <  1:3359581/58‑1 (MQ=255)

                              tAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGa                                 >  1:3243140/1‑69 (MQ=255)

                               aaGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGa                                 >  1:2794653/1‑68 (MQ=255)

                                aGGCTGCATTGGGAAGACTTAAGGAAGATGAGCTGAAACAATATAATTTGTGGCTg                                            <  1:1523407/56‑1 (MQ=255)

                                aGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTa                                       >  1:980331/1‑61 (MQ=255)

                                  gCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTtgtg                                                >  1:1596926/1‑50 (MQ=255)

                                  gCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTtgtg                                                >  1:981153/1‑50 (MQ=255)

                                  gCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCa                             >  1:3498303/1‑69 (MQ=255)

                                  gCTGCATTGGGAAGACTTAAAGAAGATGAGCCGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCa                             >  1:1042266/1‑69 (MQ=255)

                                  gCTGCATTGGGAAGACTTAAAGAAGATGAGCCGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCa                             >  1:188469/1‑69 (MQ=255)

                                  gCTGCATTGGGAAGACTTAAAGAAGATGAGCCGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCa                             >  1:1472333/1‑69 (MQ=255)

                                        tttggAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTgg                         <  1:280773/64‑1 (MQ=255)

                                        ttGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTgg                         <  1:308764/67‑1 (MQ=255)

                                          gggAAGTCTTAAAGAAGATGAGCTGAAACAATATAATTtg                                                  <  1:3182078/40‑1 (MQ=255)

                                          gggAAGACTTAAAGAAGATGAGCTGAAACAATATAATTtgtg                                                >  1:2013851/1‑42 (MQ=255)

                                          gggAAGACTTAAAGAAGATGAGCTGAAACAATATAATTtgtg                                                >  1:1179739/1‑42 (MQ=255)

                                          gggAAGACTTAAAGAAGATGAGCTGAAACAATATAATTtgtg                                                >  1:3594285/1‑42 (MQ=255)

                                          gggAAGACTTAAAGAAGATGAGCTGAAACAATATAATTtgtg                                                >  1:3069237/1‑42 (MQ=255)

                                          gggAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTg                                            <  1:2657071/46‑1 (MQ=255)

                                          gggAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTg                                            <  1:1691359/46‑1 (MQ=255)

                                          gggAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTg                                            <  1:2371480/46‑1 (MQ=255)

                                          gggAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTg                                            <  1:3027933/46‑1 (MQ=255)

                                          gggAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCa                             <  1:3600831/61‑1 (MQ=255)

                                          gggAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACt                     >  1:3222143/1‑69 (MQ=255)

                                          gggAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACt                     >  1:3072145/1‑69 (MQ=255)

                                          gggAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACt                     >  1:2852546/1‑69 (MQ=255)

                                          gggAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACt                     >  1:2534863/1‑69 (MQ=255)

                                          gggAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACt                     >  1:2761062/1‑69 (MQ=255)

                                          gggAAGAATTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACt                     >  1:2498749/1‑69 (MQ=255)

                                           ggAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTg                                            <  1:1332269/45‑1 (MQ=255)

                                             aaGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTg                                            <  1:1621055/43‑1 (MQ=255)

                                              aGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCCCTGGAACTGGtt                 <  1:2032360/69‑1 (MQ=255)

                                              aGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGtt                 <  1:1503438/69‑1 (MQ=255)

                                              aGACTTAAAGAAGATGAGCCGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGtt                 <  1:2987929/69‑1 (MQ=255)

                                              aGACTTAAAGAAGATGAGCCGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGtt                 <  1:1587376/69‑1 (MQ=255)

                                              aGACTTAAAGAAGATGAGCCGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGtt                 <  1:778098/69‑1 (MQ=255)

                                                aCTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGa               <  1:254911/69‑1 (MQ=255)

                                                 cTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTgg                         >  1:3311032/1‑58 (MQ=255)

                                                 cTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGa               >  1:107664/1‑68 (MQ=255)

                                                 cTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGAt              <  1:1453663/69‑1 (MQ=255)

                                                   tAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTg                                            >  1:1313217/1‑37 (MQ=255)

                                                    aaaGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGa               <  1:1482148/65‑1 (MQ=255)

                                                    aaaGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGa               <  1:642095/65‑1 (MQ=255)

                                                       gaagaTGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAGGCACTGGAACTGGTTGATACATc         <  1:1602686/68‑1 (MQ=255)

                                                       gaagaTGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGCAAGCACTGGAACTGGTTGATACATc         <  1:3283685/68‑1 (MQ=255)

                                                           aTGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACt                     <  1:2562872/52‑1 (MQ=255)

                                                           aTGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACt                     <  1:1831533/52‑1 (MQ=255)

                                                             gAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTgg                         <  1:3216612/46‑1 (MQ=255)

                                                             gAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTgcg  <  1:1307393/69‑1 (MQ=255)

                                                             gAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTgc   >  1:1032998/1‑68 (MQ=255)

                                                             gAGCCGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTgc   >  1:3070114/1‑68 (MQ=255)

TGATTTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCG  >  W3110S.gb/1209948‑1210077

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