BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	3,573,836	A→G	80.0%	V200A (GTG→GCG)	yihN ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	3,573,836	0	A	G	80.0%	98.8 / 19.2	45	V200A (GTG→GCG)	yihN	predicted transporter
Reads supporting (aligned to +/- strand): ref base A (4/5); new base G (17/19); total (21/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

GTTTTGTCTTTCAGCACCGAGGTCATACTGGTGTTTTTCGATGCTCCGTCTTCTTTTTTCACCTCAACGGACGGTGCCTCTTCTTTGTCATTCAGGACAAAGAAAATAATAATGCCTAC  >  W3110S.gb/3573776‑3573894

gTTTTGTCTTTCAGCACCGAGGTCATACTGGTGTTTTTCGATGCTCCGTCTTCTTTTTTCGCCTCAacg                                                    <  1:1488006/69‑1 (MQ=255)

   ttGTCTTTCAGCACCGAGGTCATACTGGTGTTTTTCGATGCTCCGTCTTCTTTTTTCACCTCAacggac                                                 <  1:3653437/69‑1 (MQ=255)

   ttGTCTTGCAGCACCGAGGTCATACTGGTGTTTTTCGATGCTCCGTCTTCTTTTTTCGCCTCAacg                                                    <  1:595311/66‑1 (MQ=255)

                  gAGGTCATACTGGTGTTTTTCGATGCTCCGTCTTCTTTTTTCGCCTCAACGGACGGTGc                                            <  1:1213163/59‑1 (MQ=255)

                  gAGGTCATACTGGTGTTTTTCGATGCTCCGTCTTCTTTTTTCGCCTCAACGGACGGTGc                                            <  1:2117967/59‑1 (MQ=255)

                   aGGTCATACTGGTGTTTTTCGATGCTCCGTCTTCTTTTTTCGCCTCAACGGACGGTGCCTCTTCttt                                   <  1:495919/67‑1 (MQ=255)

                        aTACTGGTGTTTTTCGATGCTCCGTCTTCTTTTTTCACCTCAACGGACGGTGCCtctt                                       <  1:546081/58‑1 (MQ=255)

                          aCTGGTGTTTTTCGATGCTCCGTCTTCTTTTTTCACCTCAACGGACGGTGCCTCTTCTTTGTCATTCa                           >  1:1443002/1‑68 (MQ=255)

                           cTGGTGTTTTTCGATGCTCCGTCTTCTTTTTTCGCCTCAACGGACGGTGCCTCTTCTTTGt                                 >  1:2534032/1‑61 (MQ=255)

                                 tttttCGATGCTCCGTCTTCTTTTTTCGCCTCAacggac                                                 <  1:3398494/39‑1 (MQ=255)

                                   tttCGATGCTCCGTCTTCTTTTTTCGCCTCAACGGACGGTGCCtctt                                       <  1:1633585/47‑1 (MQ=255)

                                   tttCGATGCTCCGTCTTCTTTTTTCGCCTCAACGGACGGTGCCTCTTCTTTGt                                 <  1:2742425/53‑1 (MQ=255)

                                   tttCGATGCTCCGTCTTCTTTTTTCGCCTCAACGGACGGTGCCTCTTCTTTGt                                 <  1:2025179/53‑1 (MQ=255)

                                   tttCGATGCTCCGTCTTCTTTTTTCGCCTCAACGGACGGTGCCTCTTCTTTGt                                 <  1:1672389/53‑1 (MQ=255)

                                   tttCGATGCTCCGTCTTCTTTTTTCGCCTCAACGGACGGTGCCTCTTCTTTGTCATTCAGGACAAAGaa                 >  1:1675972/1‑69 (MQ=255)

                                   tttCGATGCTCCGTCTTCTTTTTTCGCCTCAACGGACGGTGCCTCTTCTTTGTCATTCAGGACAAAGaa                 >  1:2142997/1‑69 (MQ=255)

                                   tttCGATGCTCCGTCTTCTTTTTTCGCCTCAACGGACGGTGCCTCTTCTTTGTCATTCAGGACAAAGaa                 >  1:3502368/1‑69 (MQ=255)

                                   tttCGATGCTCCGTCTTCTTTTTTCGCCTCAACGGACGGTGCCTCTTCTTTGTCATTCAGGACAAAGa                  >  1:32001/1‑68 (MQ=255)

                                   tttCGATGCTCCGTCTTCTTTTTTCGCCTCAACGGACGGTGCCTCTTCTTTGTCATTCAGGACAAAGa                  >  1:714399/1‑68 (MQ=255)

                                   tttCGATGCTCCGTCTTCTTTTTTCACCTCAACGGACGGTGCCTCTTCTTTGTCATTCAGGACAAAg                   >  1:2253322/1‑67 (MQ=255)

                                   tttCGATGCTCCGTCTTCTTTTTTCACCTCAACGGACGGTGCCTCTTCTTTGTCATTCAGGACAAAGaa                 <  1:842663/69‑1 (MQ=255)

                                   tttCGATGCTCCGTCTTCTTTTTTCACCTCAACGGACGGTGCCTCTTCTTTGTCATTCAGGACAAAGaa                 >  1:3460045/1‑69 (MQ=255)

                                     tCGATGCTCCGTCTTCTTTTTTCGCCTCAACGGACGGTGCCTCTTCTTTGTCATTCAGGACAAAGaaa                >  1:254329/1‑68 (MQ=255)

                                     tCGATGCTCCGTCTTCTTTTTTCGCCTCAACGGACGGTGCCTCTTCTTTGTCATTCAGGACAAAGaaa                >  1:3039359/1‑68 (MQ=255)

                                         tGCTCCGTCTTCTTTTTTCGCCTCAACGGACGGTGCCTCTTCTTTGTCATTCAGGACAAAGAAaataat           <  1:2175581/69‑1 (MQ=255)

                                         tGCTCCGTCTTCTTTTTTCGCCTCAACGGACGGTGCCTCTTCTTTGTCATTCAGGACAAAGAAaataat           <  1:1467119/69‑1 (MQ=255)

                                         tGCTCCGTCTTCTTTTTTCGCCTCAACGGACGGTGCCTCTTCTTTGTCATTCAGGACAAAGAAaataat           <  1:1403694/69‑1 (MQ=255)

                                         tGCTCCGTCTTCTTTTTTCGCCTCAACGGACGGTGCCTCTTCTTTGTCATTCAGGACAAAGAAaataa            >  1:453664/1‑68 (MQ=255)

                                         tGCTCCGTCTTCTTTTTTCGCCTCAACGGACGGTGCCTCTTCTTTGTCATTCAGGACAAAGAAaataa            >  1:3376923/1‑68 (MQ=255)

                                         tGCTCCGTCTTCTTTTTTCACCTCAACGGACGGTGCCTCTTCTTTGTCATTCAGGACAAAGAAaataa            <  1:1235266/68‑1 (MQ=255)

                                                tcttctTTTTTCGCCTCAACGGACGGTGCCTCTTCTTTGTCATTCAGGACAAAGAAAATAATAATGcc     <  1:710006/68‑1 (MQ=255)

                                                tcttctTTTTTCGCCTCAACGGACGGTGCCTCTTCTTTGTCATTCAGGACAAAGAAAATAATAATGcc     <  1:2008007/68‑1 (MQ=255)

                                                  ttcttTTTTCGCCTCAACGGACGGTGCCTCTTCTTTGTCATTCAGGACAAAGAAAATAATAATGCCTAc  >  1:3486061/1‑69 (MQ=255)

                                                  ttcttTTTTCGCCTCAACGGACGGTGCCTCTTCTTTGTCATTCAGGACAAAGAAAATAATAATGCCTAc  >  1:604664/1‑69 (MQ=255)

                                                  ttcttTTTTCACCTCAACGGACGGTGCCTCTTCTTTGTCATTCAGGACAAAGAAAATAATAATGCCTAc  >  1:179984/1‑69 (MQ=255)

                                                  ttcttTTTTCACCT‑AACGGACGGTGCCTCTTCTTTGTCATTCa                           <  1:302088/43‑1 (MQ=39)

                                                    cttTTTTCGCCTCAACGGACGGTGCCTCTTCTTTGt                                 <  1:3454785/36‑1 (MQ=255)

                                                    cttTTTTCGCCTCAACGGACGGTGCCTCTTCTTTGt                                 >  1:3396561/1‑36 (MQ=255)

                                                    cttTTTTCGCCTCAACGGACGGTGCCTCTTCTTTGt                                 >  1:2667982/1‑36 (MQ=255)

                                                    cttTTTTCGCCTCAACGGACGGTGCCTCTTCTTTGTCATTCAGGACAAAGAAAATAATAATGCCt    >  1:1297951/1‑65 (MQ=255)

                                                      tttttCGCCTCAACGGACGGTGCCTCTTCTTTGTCATTCAGGACAAAGAAAATAATAATGcc     <  1:880484/62‑1 (MQ=255)

                                                       ttttCGCCTCAACGGACGGTGCCTCTTCTTTGTCATTCAGGACAAAGAAaataat           <  1:663721/55‑1 (MQ=255)

                                                       ttttCGCCTCAACGGACGGTGCCTCTTCTTTGTCATTCAGGACAAAGAAaataat           >  1:787212/1‑55 (MQ=255)

                                                       ttttCGCCTCAACGGACGGTGCCTCTTCTTTGTCATTCAGGACAAAGAAaataat           >  1:2987701/1‑55 (MQ=255)

                                                       ttttCGCCTCAACGGACGGTGCCTCTTCTTTGTCATTCAGGACAAAGAAaataat           <  1:1616511/55‑1 (MQ=255)

GTTTTGTCTTTCAGCACCGAGGTCATACTGGTGTTTTTCGATGCTCCGTCTTCTTTTTTCACCTCAACGGACGGTGCCTCTTCTTTGTCATTCAGGACAAAGAAAATAATAATGCCTAC  >  W3110S.gb/3573776‑3573894

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