BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	1,547,470	C→T	64.1%	S91S (TCG→TCA)	yddL ←	predicted lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	1,547,470	0	C	T	64.1%	54.8 / 65.9	64	S91S (TCG→TCA)	yddL	predicted lipoprotein
Reads supporting (aligned to +/- strand): ref base C (15/8); new base T (16/25); total (31/33)
Fisher's exact test for biased strand distribution p-value = 6.77e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.69e-01

ATATAAGGTCAGTGATCATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGCGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGACCAAAACCAATCATTTGAT  >  W3110S.gb/1547410‑1547537

atatAAGGTCAGTGATCATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGTGAGCCGtt                                                             <  1:2527024/69‑1 (MQ=255)

atatAAGGTCAGTGATCATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGTGAGCCGtt                                                             <  1:1384918/69‑1 (MQ=255)

  ataAGGTCAGTGATCATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGTGAGCCGtt                                                             <  1:221544/67‑1 (MQ=255)

  ataAGGTCAGTGATCATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGCGAGCCGTTgg                                                           <  1:292971/69‑1 (MQ=255)

  ataAGGTCAGTGATCATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGCGAGCCGTTg                                                            <  1:842057/68‑1 (MQ=255)

   taAGGTCAGTGATCATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCTCGAGCCGTTGGc                                                          >  1:517163/1‑69 (MQ=255)

   taAGGTCAGTGATCATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGTGAGCCGTTGGc                                                          >  1:262463/1‑69 (MQ=255)

   taAGGTCAGTGATCATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGCGAGCCGTTGGc                                                          >  1:2554624/1‑69 (MQ=255)

   taAGGTCAGTGATCATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGCGAGCCGTTGGc                                                          >  1:1188667/1‑69 (MQ=255)

   taAGGTCAGTGATCATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGCGAGCCGTTGGc                                                          >  1:2981316/1‑69 (MQ=255)

   taAGGTCAGTGATCATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGCGAGCCGTTGGc                                                          >  1:369245/1‑69 (MQ=255)

   taAGGTCAGTGATCATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGCGAGCCGTTGGc                                                          >  1:3536808/1‑69 (MQ=255)

       gTCAGTGATCATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGTGAGCCGTTGGCTTc                                                       <  1:1747457/68‑1 (MQ=255)

            tGATCATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGTGAGCCGTTGGCTTCATCGt                                                  <  1:3421788/68‑1 (MQ=255)

            tGATCATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGCGAGCCGTTGGCTTCATCGt                                                  <  1:1975639/68‑1 (MQ=255)

            tGATCATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGCGAGCCGTTGGCTTCATCGt                                                  <  1:3070695/68‑1 (MQ=255)

            tGATCATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGCGAGCCGTTGGCTTCATCGt                                                  <  1:2379401/68‑1 (MQ=255)

               tCATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGCGAGCCGTTGGCTTCATCGTtata                                              >  1:3043993/1‑69 (MQ=255)

                cATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGTGAGCCGTTGGCTTCATCGt                                                  <  1:2660049/64‑1 (MQ=255)

                cATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGTGAGCCGTTGGCTTCATCGt                                                  <  1:2085659/64‑1 (MQ=255)

                cATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGCGAGCCGTTGGCTTCATCGTTATAg                                             >  1:2523587/1‑69 (MQ=255)

                 atatTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGTGAGCCGTTGGCTTCATCGTTATAGc                                            >  1:2868033/1‑69 (MQ=255)

                 atatTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGTGAGCCGTTGGCTTCATCGTTATAGc                                            >  1:2923282/1‑69 (MQ=255)

                 atatTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGTGAGCCGTTGGCTTCATCGTTATAGc                                            <  1:2834344/69‑1 (MQ=255)

                 atatTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGTGAGCCGTTGGCTTCATCGTTATAGc                                            >  1:3438028/1‑69 (MQ=255)

                 atatTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGTGAGCCGTTGGCTTCATCGTTATAGc                                            >  1:3354645/1‑69 (MQ=255)

                 atatTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGTGAGCCGTTGGCTTCATCGTTATAGc                                            >  1:1107879/1‑69 (MQ=255)

                 atatTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGTGAGCCGTTGGCTTCATCGTTATAGc                                            >  1:907045/1‑69 (MQ=255)

                 atatTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGTGAGCCGTTGGCTTCATCGTTATAGc                                            <  1:3685089/69‑1 (MQ=255)

                 atatTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGTGAGCCGTTGGCTTCATCGTTATAGc                                            >  1:457454/1‑69 (MQ=255)

                 atatTGTGATGGGCACATAACTCAATCAGAGATTCTTGCCGCGTGAGCCGTTGGCTTCATCGTTATAGc                                            >  1:186966/1‑69 (MQ=255)

                   atTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGCGAGCCGTTGGCTTCATCGt                                                  <  1:1226230/61‑1 (MQ=255)

                     tgtgATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGCGAGCCGt                                                              >  1:1395024/1‑47 (MQ=255)

                      gtgATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGCGAGCCGTTGGCTTc                                                       >  1:3534041/1‑53 (MQ=255)

                      gtgATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGCGAGCCGTTGGCTTCa                                                      >  1:2847740/1‑54 (MQ=255)

                      gtgATGGGCACTTAACTCAAT‑AGAGATTCTTGCCGCGCGAGCCGTTGGCTTCa                                                      >  1:2901573/1‑53 (MQ=255)

                           gggCACTTAACTCAATCAGAGATTCTTGCCGCGTGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAAt                                   <  1:2331936/68‑1 (MQ=255)

                           gggCACTTAACTCAATCAGAGATTCTTGCCGCGTGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAAt                                   <  1:1052073/68‑1 (MQ=255)

                           gggCACTTAACTCAATCAGAGATTCTTGCCGCGTGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAAt                                   <  1:2818435/68‑1 (MQ=255)

                           gggCACTTAACTCAATCAGAGATTCTTGCCGCGTGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAAt                                   <  1:3439620/68‑1 (MQ=255)

                             gCACTTAACTCAATCAGAGATTCTTGCCGCGTGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCa                                 <  1:3522057/68‑1 (MQ=255)

                             gCAATTAACTCAATCAGAGATTCTTGCCGCGCGAGc                                                                 >  1:2559718/1‑36 (MQ=255)

                              cACTTAACTCAATCAGAGATTCTTGCCGCGTGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCat                                <  1:1649726/68‑1 (MQ=255)

                              cACTTAACTCAATCAGAGATTCTTGCCGCGTGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCat                                <  1:1755926/68‑1 (MQ=255)

                              cACTTAACTCAATCAGAGATTCTTGCCGCGCGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCata                               <  1:2567271/69‑1 (MQ=255)

                                  tAACTCAATCAGAGATTCTTGCCGCGTGAGCCGTTGGCTTTATCGt                                                  <  1:2009618/46‑1 (MQ=38)

                                  tAACTCAATCAGAGATTCTTGCCGCGTGAGCCGTTGGCTTCATc                                                    >  1:2777642/1‑44 (MQ=255)

                                  tAACTCAATCAGAGATTCTTGCCGCGTGAGCCGTTGGCTTCATCGt                                                  <  1:973325/46‑1 (MQ=255)

                                  tAACTCAATCAGAGATTCTTGCCGCGTGAGCCGTTGGCTTCATCGt                                                  <  1:2664807/46‑1 (MQ=255)

                                  tAACTCAATCAGAGATTCTTGCCGCGTGAGCCGTTGGCTTCATCGt                                                  >  1:732922/1‑46 (MQ=255)

                                  tAACTCAATCAGAGATTCTTGCCGCGTGAGCCGTTGGCTTCATCGt                                                  <  1:1657724/46‑1 (MQ=255)

                                  tAACTCAATCAGAGATTCTTGCCGCGTGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATAtt                             >  1:1195040/1‑67 (MQ=255)

                                  tAACTCAATCAGAGATTCTTGCCGCGTGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTcc                           <  1:2906618/69‑1 (MQ=255)

                                  tAACTCAATCAGAGATTCTTGCCGCGCGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCa                                 <  1:923301/63‑1 (MQ=255)

                                  tAACTCAATCAGAGATTCTTGCCGCGCGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATAtt                             >  1:3112736/1‑67 (MQ=255)

                                   aaCTCAATCAGAGATTCTTGCCGCGTGAGCCGTTGGCTTCATCGt                                                  <  1:1139/45‑1 (MQ=255)

                                     cTCAATCAGAGATTCTTGCCGCGTGAGCCGTTGGCTTCATCGt                                                  >  1:461966/1‑43 (MQ=255)

                                     cTCAATCAGAGATTCTTGCCGCGTGAGCCGTTGGCTTCATCGt                                                  >  1:3133713/1‑43 (MQ=255)

                                       cAATCAGAGATTCTTGCCGCGCGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTg                      >  1:2988494/1‑69 (MQ=255)

                                                aTTCTTGCCGCGTGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGACGAAAAcc             <  1:147707/69‑1 (MQ=255)

                                                 ttCTTGCCGCGTGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGACCAAAAcc             <  1:1364881/68‑1 (MQ=255)

                                                 ttCTTGCCGCGTGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGACCAAAAcc             <  1:695664/68‑1 (MQ=255)

                                                  tCTTGCCGCGTGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGACCAAAACCa            >  1:1031057/1‑68 (MQ=255)

                                                   cTTGCCGCGTGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGACCAAAACCAAt          >  1:1427637/1‑69 (MQ=255)

                                                           gcgAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGACCAAAACCAATCATTTGAt  <  1:1736757/69‑1 (MQ=255)

ATATAAGGTCAGTGATCATATTGTGATGGGCACTTAACTCAATCAGAGATTCTTGCCGCGCGAGCCGTTGGCTTCATCGTTATAGCCTTTAAAATCATATTCCCAGTGACCAAAACCAATCATTTGAT  >  W3110S.gb/1547410‑1547537

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