BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	1,055,413	A→G	D396D (GAT→GAC)	torS ←	hybrid sensory histidine kinase in two‑component regulatory system with TorR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	1,055,413	0	A	G	100.0%	147.7 / NA	40	D396D (GAT→GAC)	torS	hybrid sensory histidine kinase in two‑component regulatory system with TorR
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (31/9); total (31/9)

GGTGGCGATTCAGTGCATGAACATTGCTGCGAAACGCATCCATCAGCCGCCCGATGGTATCCAGCTCCCGTACGCCAGCGGTTTCCGGGAAAG  >  W3110S.gb/1055376‑1055468

ggTGGCGATTCAGTGCATGAACATTGCTGCGAAACGCGTCCATCAGCCGCCCGATGGTATc                                  >  1:2125666/1‑61 (MQ=255)

              gCATGAACATTGCTGCGAAACGCGTCCATCAGCCCCCCGATGGTATCCAGCTCCCGTACGcc                   <  1:2856481/62‑1 (MQ=255)

              gCATGAACATTGCTGCGAAACGCGTCCATCAGCCCCCCGATGGTATCCAGCTCCCGTACGcc                   <  1:2781746/62‑1 (MQ=255)

                      aTTGCTGCGAAACGCGTCCATCAGCCGCCCGATGGTATCCAGCTCCCGTACGCCAGCGGttt           >  1:3986273/1‑62 (MQ=255)

                      aTTGCTGCGAAACGCGTCCATCAGCCGCCCGATGGTATCCAGCTCCCGTACGCCAGCGGttt           >  1:2281073/1‑62 (MQ=255)

                       ttGCTGCGAAACGCGTCCATCAGCCGCCCGATGGTATCCa                                <  1:4521790/40‑1 (MQ=255)

                       ttGCTGCGAAACGCGTCCATCAGCCGCCCGATGGTATCCa                                <  1:3415034/40‑1 (MQ=255)

                       ttGCTGCGAAACGCGTCCATCAGCCGCCCGATGGTATCCa                                <  1:358474/40‑1 (MQ=255)

                       ttGCTGCGAAACGCGTCCATCAGCCGCCCGATGGTATCCa                                <  1:2687592/40‑1 (MQ=255)

                       ttGCTGCGAAACGCGTCCATCAGCCGCCCGATGGTATCCa                                <  1:2664353/40‑1 (MQ=255)

                       ttGCTGCGAAACGCGTCCATCAGCCGCCCGATGGTATCCa                                <  1:503524/40‑1 (MQ=255)

                        tgctgcGAAACGCGTCCATCAGCCGCCCGATGGTATCCa                                <  1:686623/39‑1 (MQ=255)

                        tgctgcGAAACGCGTCCATCAGCCGCCCGATGGTATCCAGCTCCCGTACGCCAGCGGTTTcc         >  1:4521996/1‑62 (MQ=255)

                        tgctgcGAAACGCGTCCATCAGCCGCCCGATGGTATCCAGCTCCCGTACGCCAGCGGTTTcc         >  1:2266661/1‑62 (MQ=255)

                        tgctgcGAAACGCGTCCATCAGCCGCCCGATGGTATCCAGCTCCCGTACGCCAGCGGTTTcc         >  1:4174700/1‑62 (MQ=255)

                        tgctgcGAAACGCGTCCATCAGCCGCCCGATGGTATCCAGCTCCCGTACGCCAGCGGTTTcc         >  1:953614/1‑62 (MQ=255)

                        tgctgcGAAACGCGTCCATCAGCCGCCCGATGGTATCCAGCTCCCGTACGCCAGCGGTTTcc         >  1:3881920/1‑62 (MQ=255)

                        tgctgcGAAACGCGTCCATCAGCCGCCCGATGGTATCCAGCTCCCGTACGCCAGCGGTTTcc         >  1:1035496/1‑62 (MQ=255)

                        tgctgcGAAACGCGTCCATCAGCCGCCCGATGGTATCCAGCTCCCGTACGCCAGCGGTTTca         >  1:1855823/1‑61 (MQ=255)

                              gAAACGCGTCCATCAGCCGCCCGATGGTATCCAGCTCCCg                         >  1:1092156/1‑40 (MQ=255)

                              gAAACGCGTCCATCAGCCGCCCGATGGTATCCAGCTCCCg                         >  1:4376279/1‑40 (MQ=255)

                              gAAACGCGTCCATCAGCCGCCCGATGGTATCCAGCTCCCg                         >  1:721091/1‑40 (MQ=255)

                              gAAACGCGTCCATCAGCCGCCCGATGGTATCCAGCTCCCg                         >  1:3546293/1‑40 (MQ=255)

                              gAAACGCGTCCATCAGCCGCCCGATGGTATCCAGCTCCCg                         >  1:189919/1‑40 (MQ=255)

                              gAAACGCGTCCATCAGCCGCCCGATGGTATCCAGCTCCCg                         >  1:3240752/1‑40 (MQ=255)

                              gAAACGCGTCCATCAGCCGCCCGATGGTATCCAGCTCCCg                         >  1:3170945/1‑40 (MQ=255)

                              gAAACGCGTCCATCAGCCGCCCGATGGTATCCAGCTCCCg                         >  1:311966/1‑40 (MQ=255)

                              gAAACGCGTCCATCAGCCGCCCGATGGTATCCAGCTCCCg                         >  1:3534329/1‑40 (MQ=255)

                              gAAACGCGTCCATCAGCCGCCCGATGGTATCCAGCTCCCg                         >  1:2275037/1‑40 (MQ=255)

                              gAAACGCGTCCATCAGCCGCCCGATGGTATCCAGCTCCCg                         >  1:2479480/1‑40 (MQ=255)

                              gAAACGCGTCCATCAGCCGCCCGATGGTATCCAGCTCCCg                         >  1:2341452/1‑40 (MQ=255)

                               aaaCGCGTCCATCAGCCGCCCGATGGTATCCAGCTCCCGTACGCCAGCGGTTTCCGGGAAAg  >  1:4361502/1‑62 (MQ=255)

                               aaaCGCGTCCATCAGCCGCCCGATGGTATCCAGCTCCCGTACGCCAGCGGTTTCCGGGAAAg  >  1:981570/1‑62 (MQ=255)

                               aaaCGCGTCCATCAGCCGCCCGATGGTATCCAGCTCCCGTACGCCAGCGGTTTCCGGGAAAg  >  1:2092157/1‑62 (MQ=255)

                               aaaCGCGTCCATCAGCCGCCCGATGGTATCCAGCTCCCGTACGCCAGCGGTTTCCGGGAAAg  >  1:2140613/1‑62 (MQ=255)

                               aaaCGCGTCCATCAGCCGCCCGATGGTATCCAGCTCCCGTACGCCAGCGGTTTCCGGGAAAg  >  1:2620929/1‑62 (MQ=255)

                               aaaCGCGTCCATCAGCCGCCCGATGGTATCCAGCTCCCGTACGCCAGCGGTTTCCGGGAAAg  >  1:2962832/1‑62 (MQ=255)

                               aaaCGCGTCCATCAGCCGCCCGATGGTATCCAGCTCCCGTACGCCAGCGGTTTCCGGGAAAg  >  1:3591035/1‑62 (MQ=255)

                               aaaCGCGTCCATCAGCCGCCCGATGGTATCCAGCTCCCGTACGCCAGCGGTTTCCGGGAAAg  >  1:3548335/1‑62 (MQ=255)

                               aaaCGCGTCCATCAGCCGCCCGATGGTATCCAGCTCCCGTACGCCAGCGGTTTCCGGGAAAg  >  1:3291341/1‑62 (MQ=255)

GGTGGCGATTCAGTGCATGAACATTGCTGCGAAACGCATCCATCAGCCGCCCGATGGTATCCAGCTCCCGTACGCCAGCGGTTTCCGGGAAAG  >  W3110S.gb/1055376‑1055468

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