BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	1,787,235	G→A	S72S (AGC→AGT)	yfjG ←	conserved hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,787,235	0	G	A	100.0%	192.4 / NA	58	S72S (AGC→AGT)	yfjG	conserved hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (0/0); new base A (37/21); total (37/21)

CACCAGATTCATAAGAATACTTTGGTTACTGGTCAACTGGTTGCGGGTAGTAAACGTTTTGCTGATCCCAGCCTTAGAGACATCTACCGCCGCAGTCATCTGCCCAGGAGTGGACTCCAGAAT  >  minE/1787175‑1787297

caccaGATTCATAAGAATACTTTGGTTACTGGTCAACTGGTTGCGGGTAGTAAACGTTTTACTGATCCCAg                                                      >  1:1289707/1‑71 (MQ=255)

caccaGATTCATAAGAATACTTTGGTTACTGGTCAACTGGTTGCGGGTAGTAAACGTTTTACTGATCCCAg                                                      >  1:883381/1‑71 (MQ=255)

caccaGATTCATAAGAATACTTTGGTTACTGGTCAACTGGTTGCGGGTAGTAAACGTTTTACTGATCCCAg                                                      >  1:761762/1‑71 (MQ=255)

caccaGATTCATAAGAATACTTTGGTTACTGGTCAACTGGTTGCGGGTAGTAAACGTTTTACTGATCCCAg                                                      >  1:581275/1‑71 (MQ=255)

caccaGATTCATAAGAATACTTTGGTTACTGGTCAACTGGTTGCGGGTAGTAAACGTTTTACTGATCCCAg                                                      >  1:412165/1‑71 (MQ=255)

caccaGATTCATAAGAATACTTTGGTTACTGGTCAACTGGTTGCGGGTAGTAAACGTTTTACTGATCCCAg                                                      >  1:3533663/1‑71 (MQ=255)

caccaGATTCATAAGAATACTTTGGTTACTGGTCAACTGGTTGCGGGTAGTAAACGTTTTACTGATCCCAg                                                      >  1:2986067/1‑71 (MQ=255)

caccaGATTCATAAGAATACTTTGGTTACTGGTCAACTGGTTGCGGGTAGTAAACGTTTTACTGATCCCAg                                                      >  1:2532063/1‑71 (MQ=255)

caccaGATTCATAAGAATACTTTGGTTACTGGTCAACTGGTTGCGGGTAGTAAACGTTTTACTGATCCCAg                                                      >  1:149279/1‑71 (MQ=255)

caccaGATTCATAAGAATACTTTGGTTACTGGTCAACTGGTTGCGGGTAGTAAACGATTTACTGATCCCAg                                                      >  1:1796780/1‑71 (MQ=255)

   caGATTCATAAGAATACTTTGGTTACTGGTCAACTGGTTGCGGGTAGTAAACGTTTTACTGATCCCAGCCt                                                   <  1:3248438/71‑1 (MQ=255)

   caGATTCATAAGAATACTTTGGTTACTGGTCAACTGGTTGCGGGTAGTAAACGTTTTACTGATCCCAGCCt                                                   <  1:811460/71‑1 (MQ=255)

   caGATTCATAAGAATACTTTGGTTACTGGTCAACTGGTTGCGGGTAGTAAACGTTTTACTGATCCCAGCCt                                                   <  1:139883/71‑1 (MQ=255)

       ttCATAAGAATACTTTGGTTACTGGTCAACTGGTTGCGGGTAGTAAACGTTTTACTGATCCCAGCCtt                                                  >  1:1252954/1‑68 (MQ=255)

       ttCATAAGAATACTTTGGTTACTGGTCAACTGGTTGCGGGTAGTAAACGTTTTACTGATCCCAGCCTTa                                                 >  1:586194/1‑69 (MQ=255)

       ttCATAAGAATACTTTGGTTACTGGTCAACTGGTTGCGGGTAGTAAACGTTTTACTGATCCCAGCCTTa                                                 >  1:3905138/1‑69 (MQ=255)

                  aCTTTGGTTACTGGTCAACTGGTTGCGGGTAGTAAACGTTTTACTGATCCCAGCCTTAGAGACATCTAccg                                    >  1:1996742/1‑71 (MQ=255)

                  aCTTTGGTTACTGGTCAACTGGTTGCGGGTAGTAAACGTTTTACTGATCCCAGCCTTAGAGACATCTAccg                                    >  1:2903038/1‑71 (MQ=255)

                  aCTTTGGTTACTGGTCAACTGGTTGCGGGTAGTAAACGTTTTACTGATCCCAGCCTTAGAGACATCTAccg                                    >  1:3787412/1‑71 (MQ=255)

                  aCTTTGGTTACTGGTCAACTGGTTGCGGGTAGTAAACGTTTTACTGATCCCAGCCTTAGAGACATCTAccg                                    >  1:1456360/1‑71 (MQ=255)

                  aCTTTGGTTACTGGTCAACTGGTTGCGGGTAGTAAACGTTTTACTGATCCCAGCCTTAGAGACATCTAcc                                     >  1:1598381/1‑70 (MQ=255)

                      tGGTTACTGGTCAACTGGTTGCGGGTAGTAAATGTTTTACTGATCCCAGCCt                                                   >  1:636590/1‑52 (MQ=255)

                      tGGTTACTGGTCAACTGGTTGCGGGTAGTAAACGTTTTACTGATCCCAGCCt                                                   >  1:2509779/1‑52 (MQ=255)

                      tGGTTACTGGTCAACTGGTTGCGGGTAGTAAACGTTTTACTGATCCCAGCCt                                                   >  1:4278111/1‑52 (MQ=255)

                      tGGTTACTGGTCAACTGGTTGCGGGTAGTAAACGTTTTACTGATCCCAGCCt                                                   >  1:453968/1‑52 (MQ=255)

                      tGGTTACTGGTCAACTGGTTGCGGGTAGTAAACGTTTTACTGATCCCAGCCt                                                   >  1:3689037/1‑52 (MQ=255)

                      tGGTTACTGGTCAACTGGTTGCGGGTAGTAAACGTTTTACTGATAACAGCCt                                                   >  1:4195199/1‑52 (MQ=255)

                               gTCAACTGGTTTCGGGTAGTAAACGTTTTACTGATCCCAGCCTTAGAGACATCTACCGCCGCAGTCATCTg                       <  1:2558420/71‑1 (MQ=255)

                               gTCAACTGGTTGCGGGTAGTAAACGTTTTACTGATCCCAGCCTTAGAGACATCTACCGCCGCAGTCATCTg                       <  1:4341998/71‑1 (MQ=255)

                               gTCAACTGGTTGCGGGTAGTAAACGTTTTACTGATCCCAGCCTTAGAGACATCTACCGCCGCAGTCATCTg                       <  1:568942/71‑1 (MQ=255)

                               gTCAACTGGTTGCGGGTAGTAAACGTTTTACTGATCCCAGCCTTAGAGACATCTACCGCCGCAGTCATCTg                       <  1:3703655/71‑1 (MQ=255)

                               gTCAACTGGTTGCGGGTAGTAAACGTTTTACTGATCCCAGCCTTAGAGACATCTACCGCCGCAGTCATCTg                       <  1:146544/71‑1 (MQ=255)

                               gTCAACTGGTTGCGGGTAGTAAACGTTTTACTGATCCCAGCCTTAGAGACATCTACCGCCGCAGTCATCTg                       <  1:2077355/71‑1 (MQ=255)

                               gTCAACTGGTTGCGGGTAGTAAACGTTTTACTGATCCCAGCCTTAGAGACATCTACCGCCGCAGTCATCTg                       <  1:1613199/71‑1 (MQ=255)

                               gTCAACTGGTTGCGGGTAGTAAACGTTTTACTGATCCCAGCCTTAGAGACATCTACCGCCGCAGTCATCTg                       <  1:3127508/71‑1 (MQ=255)

                               gTCAACTGGTTGCGGGTAGTAAACGTTTTACTGATCCCAGCCTTAGAGACATCTACCGCCGCAGTCATCTg                       <  1:299125/71‑1 (MQ=255)

                               gTCAACTGGTTGCGGGTAGTAAACGTTTTACTGATCCCAGCCTTAGAGACATCTACCGCCGCAGTCATCTg                       <  1:338846/71‑1 (MQ=255)

                               gTCAACTGGTTGCGGGTAGTAAACGTTTTACTGATCCCAGCCTTAGAGACATCTACCGCCGCAGTCATCTg                       <  1:151581/71‑1 (MQ=255)

                                       gTTGCGGGTAGTAAACGTTTTACTGATCCCAGCCTTAGAGACATCTACCGCCGCAGTCATCTGCCCAgg                 <  1:4134041/69‑1 (MQ=255)

                                       gTTGCGGGTAGTAAACGTTTTACTGATCCCAGCCTTAGAGACATCTACCGCCGCAGTCATCTGCCCAgg                 <  1:2229362/69‑1 (MQ=255)

                                       gTTGCGGGTAGTAAACGTTTTACTGATCCCAGCCTTAGAGACATCTACCGCCGCAGTCATCTGCCCAgg                 <  1:1651471/69‑1 (MQ=255)

                                       gTTGCGGGTAGTAAACGTTTTACTGATCCCAGCCTTAGAGACATCTACCGCCGCAGTCATCTGCCCAgg                 <  1:3961317/69‑1 (MQ=255)

                                       gTTGCGGGTAGTAAACGTTTTACTGATCCCAGCCTTAGAGACATCTACCGCCGCAGTCATCTGCCCAgg                 <  1:273156/69‑1 (MQ=255)

                                       gTTGCGGGTAGTAAACGTTTTACTGATCCCAGCCTTAGAGACATCTACCGCCGCAGTCATCTGCCCAgg                 <  1:1564532/69‑1 (MQ=255)

                                                   aaaCGTTTTACTGATCCCAGCCTTAGAGACATCTACCGCCGCAGTCATCTg                       <  1:722817/51‑1 (MQ=255)

                                                    aaCGTTTTACTGATCCCAGCCTTAGAGACATCTACCGCCGCAGTCATCTGCCCAGGAGTGGACTCCAGAAt  >  1:452648/1‑71 (MQ=255)

                                                      cGTTTTACTGATCCCAGCCTTAGAGACATCTACCGCCGCAGTCATCTg                       >  1:2140585/1‑48 (MQ=255)

                                                      cGTTTTACTGATCCCAGCCTTAGAGACATCTACCGCCGCAGTCATCTGc                      >  1:2156192/1‑49 (MQ=255)

                                                      cGTTTTACTGATCCCAGCCTTAGAGACATCTACCGCCGCAGTCATCTGc                      >  1:1326619/1‑49 (MQ=255)

                                                      cGTTTTACTGATCCCAGCCTTAGAGACATCTACCGCCGCAGTCATCTGc                      >  1:1656309/1‑49 (MQ=255)

                                                      cGTTTTACTGATCCCAGCCTTAGAGACATCTACCGCCGCAGTCATCTGc                      >  1:1726852/1‑49 (MQ=255)

                                                      cGTTTTACTGATCCCAGCCTTAGAGACATCTACCGCCGCAGTCATCTGc                      >  1:510312/1‑49 (MQ=255)

                                                      cGTTTTACTGATCCCAGCCTTAGAGACATCTACCGCCGCAGTCATCTGc                      >  1:498490/1‑49 (MQ=255)

                                                      cGTTTTACTGATCCCAGCCTTAGAGACATCTACCGCCGCAGTCATCTGc                      >  1:2532050/1‑49 (MQ=255)

                                                      cGTTTTACTGATCCCAGCCTTAGAGACATCTACCGCCGCAGTCATCTGc                      >  1:4178381/1‑49 (MQ=255)

                                                      cGTTTTACTGATCCCAGCCTTAGAGACATCTACCGCCGCAGTCATCTGc                      >  1:2878173/1‑49 (MQ=255)

                                                      cGTTTTACTGATCCCAGCCTTAGAGACATCTACCGCCGCAGTCATCTGc                      >  1:3584176/1‑49 (MQ=255)

                                                      cGTTTTACTGATCCCAGCCTTAGAGACATCTACCGCCGCAGTCATCTGc                      >  1:1555041/1‑49 (MQ=255)

CACCAGATTCATAAGAATACTTTGGTTACTGGTCAACTGGTTGCGGGTAGTAAACGTTTTGCTGATCCCAGCCTTAGAGACATCTACCGCCGCAGTCATCTGCCCAGGAGTGGACTCCAGAAT  >  minE/1787175‑1787297

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