BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	2,759,336	C→T	65.8%	T45I (ACA→ATA)	yfjJ →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,759,336	0	C	T	65.8%	75.1 / 71.8	73	T45I (ACA→ATA)	yfjJ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (5/20); new base T (17/31); total (22/51)
Fisher's exact test for biased strand distribution p-value = 1.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.97e-01

GACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCTATATCTCGCTTTACCAGTG  >  W3110S.gb/2759272‑2759404

gaCTCCCTAAACGAACATCTTAAGTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGa                                                                   <  1:2891324/68‑1 (MQ=255)

gaCTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGaa                                                                  >  1:3351996/1‑69 (MQ=255)

gaCTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGa                                                                   <  1:2407954/68‑1 (MQ=255)

gaCTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGa                                                                   <  1:2612473/68‑1 (MQ=255)

gaCTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGa                                                                   <  1:43886/68‑1 (MQ=255)

gaCTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGa                                                                   <  1:1999730/68‑1 (MQ=255)

gaCTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGa                                                                   <  1:437731/68‑1 (MQ=255)

   tCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGaaaa                                                                <  1:2145675/68‑1 (MQ=255)

   tCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGaaaa                                                                <  1:1987877/68‑1 (MQ=255)

   tCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGaaaa                                                                <  1:3487387/68‑1 (MQ=255)

   tCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGaaaa                                                                <  1:3013712/68‑1 (MQ=255)

   tCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGaaaa                                                                <  1:1032468/68‑1 (MQ=255)

   tCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGaaaa                                                                <  1:63010/68‑1 (MQ=255)

       tAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATg                                                            >  1:396062/1‑68 (MQ=255)

       tAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATg                                                            >  1:2169598/1‑68 (MQ=255)

         aaCGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGAt                                                         >  1:1547459/1‑69 (MQ=255)

            gAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAAc                                                      <  1:2435310/69‑1 (MQ=255)

             aaCATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGATAAcc                                                     <  1:3355923/69‑1 (MQ=255)

             aaCATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGATAAcc                                                     <  1:3052300/69‑1 (MQ=255)

              aCATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATgg                                                           >  1:2623656/1‑62 (MQ=255)

              aCATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGa                                                          <  1:359618/63‑1 (MQ=255)

              aCATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGa                                                          <  1:2077910/63‑1 (MQ=255)

              aCATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGa                                                          <  1:2082149/63‑1 (MQ=255)

              aCATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGa                                                          <  1:271372/63‑1 (MQ=255)

              aCATCTTAACTCAATCATTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGATAAc                                                      >  1:1572969/1‑67 (MQ=255)

                aTCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCAt                                                  <  1:3319169/69‑1 (MQ=255)

                aTCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCAt                                                  <  1:214163/69‑1 (MQ=255)

                 tCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGaa                                                                  <  1:1344790/52‑1 (MQ=255)

                  cTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTc                                                >  1:1207928/1‑69 (MQ=255)

                   ttAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGaaaa                                                                >  1:1964473/1‑52 (MQ=255)

                   ttAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGaaaa                                                                >  1:468045/1‑52 (MQ=255)

                   ttAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGAt                                                         >  1:865753/1‑59 (MQ=255)

                       cTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCt                                               <  1:1256097/65‑1 (MQ=255)

                       cTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCt                                               <  1:317359/65‑1 (MQ=255)

                        tCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCt                                               <  1:1044201/64‑1 (MQ=255)

                            tGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCAAccc                                       <  1:3593602/68‑1 (MQ=255)

                              aTTCTACGTGTTGATCTGCATGACCCAATTGATATAGaaa                                                                 <  1:1841389/40‑1 (MQ=255)

                              aTTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAgg                                    >  1:2042494/1‑69 (MQ=255)

                              aTTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAgg                                    <  1:3139313/69‑1 (MQ=255)

                                tCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGATAAcc                                                     <  1:2955351/50‑1 (MQ=255)

                                tCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGATAAcc                                                     <  1:3532062/50‑1 (MQ=255)

                                tCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGATAAcc                                                     <  1:3632666/50‑1 (MQ=255)

                                tCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCAACCCAggg                                   >  1:2124893/1‑68 (MQ=255)

                                 cTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGtt                                 <  1:1147013/69‑1 (MQ=255)

                                 cTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGtt                                 <  1:2559565/69‑1 (MQ=255)

                                    cGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGAc                              >  1:3545224/1‑69 (MQ=255)

                                       gtTGATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGa                               <  1:2274068/65‑1 (MQ=255)

                                          gATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAAccc                                       >  1:3627817/1‑54 (MQ=255)

                                          gATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAAccc                                       >  1:1228833/1‑54 (MQ=255)

                                          gATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAAccc                                       >  1:1809001/1‑54 (MQ=255)

                                          gATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAAc                                         >  1:1365387/1‑52 (MQ=255)

                                          gATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGt                        >  1:2864402/1‑69 (MQ=255)

                                            tCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCa                                           <  1:2770230/48‑1 (MQ=255)

                                            tCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCa                                           <  1:228255/48‑1 (MQ=255)

                                            tCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCa                                           <  1:184524/48‑1 (MQ=255)

                                             cTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGc                      <  1:2724296/68‑1 (MQ=255)

                                               gcttGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCta                    <  1:899025/65‑1 (MQ=255)

                                               gCATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCta                    <  1:2718160/68‑1 (MQ=255)

                                               gCATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCta                    <  1:2736730/68‑1 (MQ=255)

                                                cATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCtat                   >  1:1840183/1‑68 (MQ=255)

                                                cATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCta                    <  1:328902/67‑1 (MQ=255)

                                                cATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCta                    <  1:46571/67‑1 (MQ=255)

                                                cATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCta                    <  1:2128985/67‑1 (MQ=255)

                                                  tGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCta                    <  1:3320208/65‑1 (MQ=255)

                                                  tGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCta                    <  1:1400037/65‑1 (MQ=255)

                                                      ccAATTGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCTATATCTCGc            >  1:1171247/1‑69 (MQ=255)

                                                       cAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCa                                      <  1:589406/42‑1 (MQ=255)

                                                       cAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCa                                      <  1:2312145/42‑1 (MQ=255)

                                                       cAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCa                                      <  1:2316463/42‑1 (MQ=255)

                                                           tGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGtt                                 >  1:2978245/1‑43 (MQ=255)

                                                           tGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGtt                                 <  1:3179507/43‑1 (MQ=255)

                                                           tGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCTATATCTCGCttt         >  1:2739287/1‑67 (MQ=255)

                                                           tGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCTATATCTCGCTTTAc       <  1:1178083/69‑1 (MQ=255)

                                                                cAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCTATATCTCGCTTTACCAGTg  >  1:1659391/1‑69 (MQ=255)

GACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCTATATCTCGCTTTACCAGTG  >  W3110S.gb/2759272‑2759404

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