BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,382,555	C→T	S203S (AGC→AGT)	ompG →	outer membrane porin G

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,382,555	0	C	T	100.0%	166.3 / NA	42	S203S (AGC→AGT)	ompG	outer membrane porin G
Reads supporting (aligned to +/- strand): ref base C (0/0); new base T (16/26); total (16/26)

TTGCGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGCCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTATTTG  >  NC_000913/1382508‑1382600

ttGCGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCg                                             >  1:4935087/1‑50 (MQ=255)

   cGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCaa                                          >  1:1170697/1‑50 (MQ=255)

   cGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCaa                                          <  1:3133474/50‑1 (MQ=255)

   cGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCa                                           <  1:2214280/49‑1 (MQ=255)

    gAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAAt                                         <  1:4426218/50‑1 (MQ=255)

     aGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATa                                        <  1:3807479/50‑1 (MQ=255)

     aGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATa                                        <  1:3915827/50‑1 (MQ=255)

         aaCTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACgg                                    >  1:5012974/1‑50 (MQ=255)

          aCTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGt                                   <  1:3536320/50‑1 (MQ=255)

              ttatCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGt                               <  1:139844/50‑1 (MQ=255)

              ttatCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGt                               <  1:4229278/50‑1 (MQ=255)

              ttatCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGt                               <  1:261223/50‑1 (MQ=255)

               tatCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGtt                              >  1:3765489/1‑50 (MQ=255)

               tatCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGtt                              >  1:2700890/1‑50 (MQ=255)

                 tctcGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGttt                             <  1:2683991/49‑1 (MQ=255)

                 tctcGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGttt                             <  1:1909376/49‑1 (MQ=255)

                 tctcGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGttt                             <  1:5131448/49‑1 (MQ=255)

                     gAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCa                         <  1:3921702/49‑1 (MQ=255)

                     gAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCa                         <  1:1117409/49‑1 (MQ=255)

                      aGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACg                       >  1:2068876/1‑50 (MQ=255)

                      aGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACg                       >  1:5028689/1‑50 (MQ=255)

                         gcgGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCaa                    <  1:1439037/50‑1 (MQ=255)

                          cgGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGc                      >  1:3551480/1‑47 (MQ=255)

                           ggCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGa                  <  1:1517662/50‑1 (MQ=255)

                           ggCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGa                  <  1:1407951/50‑1 (MQ=255)

                            gCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGa                  >  1:213801/1‑49 (MQ=255)

                            gCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGa                  >  1:3739870/1‑49 (MQ=255)

                                cAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAAtt              <  1:3395977/49‑1 (MQ=255)

                                 aaTATGGACTACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTcg            >  1:719500/1‑50 (MQ=255)

                                 aaTATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTc             <  1:2185800/49‑1 (MQ=255)

                                  atatGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTcgc           <  1:1645182/50‑1 (MQ=255)

                                  atagGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTcgc           <  1:2930573/50‑1 (MQ=255)

                                   tatGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTcgcg          >  1:2005418/1‑50 (MQ=255)

                                   tatGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTcgcg          >  1:813652/1‑50 (MQ=255)

                                     tGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCgcgc         <  1:587477/49‑1 (MQ=255)

                                      ggACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGcc        >  1:2409426/1‑49 (MQ=255)

                                       gacgacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCt       <  1:5003727/49‑1 (MQ=255)

                                         cgacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTAtt    >  1:2360591/1‑50 (MQ=255)

                                          gacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTAttt   <  1:2881658/50‑1 (MQ=255)

                                          gacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTAttt   >  1:4328240/1‑50 (MQ=255)

                                           acAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTATTTg  <  1:3651833/50‑1 (MQ=255)

                                           acAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTATTTg  <  1:695206/50‑1 (MQ=255)

TTGCGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGCCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTATTTG  >  NC_000913/1382508‑1382600

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