BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	504,214	A→G	S161S (AGT→AGC)	fsr ←	fosmidomycin efflux pump

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	504,214	0	A	G	100.0%	124.4 / NA	36	S161S (AGT→AGC)	fsr	fosmidomycin efflux pump
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (19/17); total (19/17)

TTCGGTGCTGTGCCGAGTACCAACGGCTGATTTGCGCCAACACCACGATCGCCAGCAGTGCCGCAAGCACAAACCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACTGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCACGCGAGAAGATTCCGGATGAAAGACCGATGAACCGGTACCGACCAGCGCCGCCGC  >  NC_000913/504078‑504342

ttcGGTGCTGTGCCGAGTACCAACGGCTGATTTGCGCCAACACCACGATCGCCAGCAGTGCCGCAAGCACAAACCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAGCTg                                                                                                                               >  2:29575/1‑140 (MQ=255)

 tcGGTGCTGTGCCGAGTACCAACGGCTGATTTGCGCCAACACCACGATCGCCAGCAGTGCCGCAAGCACAAACCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAGCTGc                                                                                                                              >  1:186200/1‑140 (MQ=255)

               agcaCCAACGGCTGATTTGCGCCAACACCACGATCGCCAGCAGTGCCGCAAGCACAAACCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAGCTGCCAAAGTTGccgccg                                                                                                                >  1:508942/4‑140 (MQ=255)

               aGTACCAACGGCTGATTTGCGCCAACACCACGATCGCCAGCAGTGCCGCAAGCACAAACCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAGCTGCCAAAGTTGccgccg                                                                                                                >  2:815011/1‑140 (MQ=255)

               aGTACCAACGGCTGATTTGCGCCAACACCACGATCGCCAGCAGTGCCGCAAGCACAAACCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAGCTGCCAAAGTTGccgccg                                                                                                                >  2:515637/1‑140 (MQ=255)

               aGTACCAACGGCTGATTTGCGCCAACACCACGATCGCCAGCAGTGCCGCAAGCACAAACCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAGCTGCCAAAGTTGccgccg                                                                                                                >  1:497696/1‑140 (MQ=255)

               aGTACCAACGGCTGATTTGCGCCAACACCACGATCGCCAGCAGTGCCGCAAGCACAAACCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAGCTGCCAAAGTTGccgccg                                                                                                                >  1:640120/1‑140 (MQ=255)

                 tACCAACGGCTGATTTGCGCCAACACCACGATCGCCAGCAGTGCCGCAAGCACAAACCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAGCTGCCAAAGTTGCCGCCGAc                                                                                                              <  2:6974/140‑1 (MQ=255)

                 tACCAACGGCTGATTTGCGCCAACACCACGATCGCCAGCAGTGCCGCAAGCACAAACCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAGCTGCCAAAGTTGCCGCCGAc                                                                                                              <  2:640120/140‑1 (MQ=255)

                 tACCAACGGCTGATTTGCGCCAACACCACGATCGCCAGCAGTGCCGCAAGCACAAACCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAGCTGCCAAAGTTGCCGCCGAc                                                                                                              <  1:515637/140‑1 (MQ=255)

                   ccAACGGCTGATTTGCGCCAACACCACGATCGCCAGCAGTGCCGCAAGCACAAACCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAGCTGCCAAAGTTGCCGCCGACCt                                                                                                            <  1:815011/140‑1 (MQ=255)

                    cAACGGCTGATTTGCGCCAACACCACGATCGCCAGCAGTGCCGCAAGCACAAACCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAGCTGCCAAAGTTGCCGCCGACCTg                                                                                                           >  1:13786/1‑140 (MQ=255)

                           tGATTTGCGCCAACACCACGATCGCCAGCAGTGCCGCAAGCACAAACCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAGCTGCCAAAGTTGCCGCCGACCTGAAaga                                                                                                      >  1:479884/1‑138 (MQ=255)

                           tGATTTGCGCCAACACCACGATCGCCAGCAGTGCCGCAAGCACAAACCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAGCTGCCAAAGTTGCCGCCGACCTGAAaga                                                                                                      <  2:479884/138‑1 (MQ=255)

                             aTTTGCGCCAACACCACGATCGCCAGCAGTGCCGCAAGCACAAACCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAGCTGCCAAAGTTGCCGCCGACCTGAAagataga                                                                                                  <  2:508942/140‑1 (MQ=255)

                              tttGCGCCAACACCACGATCGCCAGCAGTGCCGCAAGCACAAACCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAGCTGCCAAAGTTGCCGCCGACCTGAAagatagat                                                                                                 >  2:643911/1‑140 (MQ=255)

                                        caccacGATCGCCAGCAGTGCCGCAAGCACAAACCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAGCTGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGc                                                                                       >  1:443582/1‑140 (MQ=255)

                                                   ccAGCAGTGCCGCAAGCACAAACCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAGCTGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGccgcccg                                                                            >  1:567429/1‑140 (MQ=255)

                                                     agcagTGCCGCAAGCACAAACCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAGCTGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCcgcccgcc                                                                          <  2:497696/140‑1 (MQ=255)

                                                         gTGCCGCAAGCACAAACCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAGCTGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGaa                                                                      <  1:24997/140‑1 (MQ=255)

                                                          tGCCGCAAGCACAAACCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAGCTGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAg                                                                     >  2:862243/1‑140 (MQ=255)

                                                          tGCCGCAAGCACAAACCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAGCTGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAg                                                                     >  2:209535/1‑140 (MQ=255)

                                                            ccGCAAGCACAAACCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAGCTGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGcc                                                                   >  2:105670/1‑140 (MQ=255)

                                                               cAAGCACAAACCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAGCTGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATa                                                                <  2:443582/140‑1 (MQ=255)

                                                                   cacaAACCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAGCTGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACggg                                                            <  1:862243/140‑1 (MQ=255)

                                                                   cacaAACCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAGCTGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACggg                                                            <  1:209535/140‑1 (MQ=255)

                                                                        aCCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAGCTGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCAcg                                                       <  2:13786/140‑1 (MQ=255)

                                                                        aCCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAGCTGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCAcg                                                       <  1:105670/140‑1 (MQ=255)

                                                                            gggaaCGTCGCATTTGCCATAAGGCGCGATAATCCCCGCCGCCAGCAAGGGTCCCAGGGAGCTGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCACGCgag                                                   <  2:567429/137‑1 (MQ=255)

                                                                                  gTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAGCTGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCACGCGAGAAGAtt                                             >  2:310081/1‑140 (MQ=255)

                                                                                      ccTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAGCTGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCACGCGAGAAGATTCCgg                                         >  1:423952/1‑140 (MQ=255)

                                                                                           gCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAGCTGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCACGCGAGAAGATTCCGGATGaa                                    >  2:499417/1‑140 (MQ=255)

                                                                                                      cgATAATCACCGCCGCCAGCAAGGGTCCCAGGGAGCTGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCACGCGAGAAGATTCCGGATGAAAGACCGATGaa                         >  2:809791/1‑140 (MQ=255)

                                                                                                       gATAATCACCGCCGCCAGCAAGGGTCCCAGGGAGCTGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCACGCGAGAAGATTCCGGATGAAAGACCGATGAAc                        <  2:186200/140‑1 (MQ=255)

                                                                                                                           ggggTCCCAGGGAGCTGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCACGCGAGAAGATTCCGGATGAAAGACCGATGAACCGGTACCGACCAGcgccgcc    <  2:429231/139‑1 (MQ=255)

                                                                                                                             ggTCCCAGGGAGCTGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCACGCGAGAAGATTCCGGATGAAAGACCGATGAACCGGTACCGACCAGcgccgccgc  <  2:843087/140‑1 (MQ=255)

TTCGGTGCTGTGCCGAGTACCAACGGCTGATTTGCGCCAACACCACGATCGCCAGCAGTGCCGCAAGCACAAACCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACTGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCACGCGAGAAGATTCCGGATGAAAGACCGATGAACCGGTACCGACCAGCGCCGCCGC  >  NC_000913/504078‑504342

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