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%	257.7 / NA	62	S203S (AGC→AGT)	ompG	outer membrane porin G
Reads supporting (aligned to +/- strand): ref base C (0/0); new base T (22/40); total (22/40)

ATCTCGAGCGCGGCTTCAATATGGACGACAGCCGCAATAACGGTGAGTTTTCCACGCAAGAA > NC_000913/1382524‑1382585

atCTCGAGCGCGGCTTCAATATGGACGACAGTCGc                             >  1:1490918/1‑35 (MQ=255)

 tctcGAGCGCGGCTTCAATATGGACGACAGTCGCa                            >  1:3133269/1‑35 (MQ=255)

  ctcGAGCGCGGCTTCAATATGGACGACAGTCGCaa                           <  1:916699/35‑1 (MQ=39)

  ctcGAGCGCGGCTTCAATATGGACGACAGTCGCaa                           <  1:4439673/35‑1 (MQ=255)

  ctcGAGCGCGGCTTCAATATGGACGACAGTCGCaa                           <  1:5019102/35‑1 (MQ=255)

  ctcGAGCGCGGCTTCAATATGGACGACAGTCGCaa                           <  1:6535647/35‑1 (MQ=255)

  ctcGAGCGCGGCTTCAATATGGACGACAGTCGCaa                           <  1:7620063/35‑1 (MQ=255)

   tcGAGCGCGGCTTCAATATGGACGACAGTCGCAAt                          <  1:6454630/35‑1 (MQ=39)

   tcGAGCGCGGCTTCAATATGGACGACAGTCGCAAt                          >  1:7889267/1‑35 (MQ=39)

   tcGAGCGCGGCTTCAATATGGACGACAGTCGCAAt                          >  1:2414863/1‑35 (MQ=255)

   tcGAGCGCGGCTTCAATATGGACGACAGTCGCAAt                          <  1:2447194/35‑1 (MQ=255)

   tcGAGCGCGGCTTCAATATGGACGACAGTCGCAAt                          <  1:1096216/35‑1 (MQ=39)

   tcGAGCGCGGCTTCAATATGGACGACAGTCGCAAt                          <  1:3590792/35‑1 (MQ=255)

   tcGAGCGCGGCTTCAATATGGACGACAGTCGCAAt                          <  1:3242732/35‑1 (MQ=255)

   tcGAGCGCGGCTTCAATATGGACGACAGTCGCAAt                          >  1:3173067/1‑35 (MQ=255)

    cGAGCGCGGCTTCAATATGGACGACAGTCGCAATa                         <  1:2931971/35‑1 (MQ=255)

     gAGCGCGGCTTCAATATGGACGACAGTCGCAATaa                        >  1:6293086/1‑35 (MQ=39)

     gAGCGCGGCTTCAATATGGACGACAGTCGCAATaa                        <  1:5236240/35‑1 (MQ=39)

     gAGCGCGGCTTCAATATGGACGACAGTCGCAATaa                        <  1:2226314/35‑1 (MQ=255)

      aGCGCGGCTTCAATATGGACGACAGTCGCAATAAc                       >  1:3662592/1‑35 (MQ=255)

      aGCGCGGCTTCAATATGGACGACAGTCGCAATAAc                       <  1:1025915/35‑1 (MQ=255)

        cgcgGCTTCAATATGGACGACAGTCGCAATAACgg                     <  1:2534663/35‑1 (MQ=255)

        cgcgGCTTCAATATGGACGACAGTCGCAATAACgg                     <  1:6324054/35‑1 (MQ=39)

        cgcgGCTTCAATATGGACGACAGTCGCAATAACgg                     <  1:5501806/35‑1 (MQ=39)

        cgcgGCTTCAATATGGACGACAGTCGCAATAACgg                     <  1:4283584/35‑1 (MQ=39)

         gcgGCTTCAATATGGACGACAGTCGCAATAACGGt                    <  1:658484/35‑1 (MQ=255)

         gcgGCTTCAATATGGACGACAGTCGCAATAACGGt                    <  1:3599787/35‑1 (MQ=255)

         gcgGCTTCAATATGGACGACAGTCGCAATAACGGt                    >  1:4201172/1‑35 (MQ=39)

         gcgGCTTCAATATGGACGACAGTCGCAATAACGGt                    <  1:4268741/35‑1 (MQ=39)

         gcgGCTTCAATATGGACGACAGTCGCAATAACGGt                    <  1:326919/35‑1 (MQ=255)

            gCTTCAATATGGACGACAGTCGCAATAACGGTGAg                 <  1:5837328/35‑1 (MQ=39)

             cTTCAATATGGACGACAGTCGCAATAACGGTGAGt                >  1:3670888/1‑35 (MQ=39)

             cTTCAATATGGACGACAGTCGCAATAACGGTGAGt                <  1:5234347/35‑1 (MQ=255)

              ttCAATATGGACGACAGTCGCAATAACGGTGAGtt               <  1:4090204/35‑1 (MQ=39)

                cAATATGGACGACAGTCGCAATAACGGTGAGtttt             <  1:5985134/35‑1 (MQ=39)

                cAATATGGACGACAGTCGCAATAACGGTGAGtttt             <  1:5970750/35‑1 (MQ=255)

                 aaTATGGACGACAGTCGCAATAACGGTGAGTTTTc            >  1:2662692/1‑35 (MQ=255)

                 aaTATGGACGACAGTCGCAATAACGGTGAGTTTTc            <  1:6070506/35‑1 (MQ=255)

                 aaTATGGACGACAGTCGCAATAACGGTGAGTTTTc            <  1:7677813/35‑1 (MQ=255)

                  atatGGACGACAGTCGCAATAACGGTGAGTTTTcc           >  1:5207367/1‑35 (MQ=39)

                   tatGGACGACAGTCGCAATAACGGTGAGTTTTCCa          >  1:2729494/1‑35 (MQ=38)

                   tatGGACGACAGTCGCAATAACGGTGAGTTTTCCa          <  1:1192870/35‑1 (MQ=39)

                   tatGGACGACAGTCGCAATAACGGTGAGTTTTCCa          >  1:4142289/1‑35 (MQ=39)

                   tatGGACGACAGTCGCAATAACGGTGAGTTTTCCa          >  1:3846148/1‑35 (MQ=39)

                    atGGACGACAGTCGCAATAACGGTGAGTTTTCCAc         >  1:1010680/1‑35 (MQ=39)

                    atGGACGACAGTCGCAATAACGGTGAGTTTTCCAc         >  1:1804774/1‑35 (MQ=255)

                     tGGACGACAGTCGCAATAACGGTGAGTTTTCCACg        >  1:1338485/1‑35 (MQ=255)

                     tGGACGACAGTCGCAATAACGGTGAGTTTTCCACg        >  1:1514780/1‑35 (MQ=39)

                     tGGACGACAGTCGCAATAACGGTGAGTTTTCCACg        >  1:2354799/1‑35 (MQ=255)

                      ggACGACAGTCGCAATAACGGTGAGTTTTCCACGc       <  1:3804630/35‑1 (MQ=255)

                      ggACGACAGTCGCAATAACGGTGAGTTTTCCACGc       <  1:7434774/35‑1 (MQ=39)

                         cgacAGTCGCAATAACGGTGAGTTTTCCACGCAAg    >  1:762990/1‑35 (MQ=255)

                         cgacAGTCGCAATAACGGTGAGTTTTCCACGCAAg    <  1:2713178/35‑1 (MQ=255)

                         cgacAGTCGCAATAACGGTGAGTTTTCCACGCAAg    <  1:1719770/35‑1 (MQ=255)

                         cgacAGTCGCAATAACGGTGAGTTTTCCACGCAAg    >  1:2696487/1‑35 (MQ=39)

                          gacAGTCGCAATAACGGTGAGTTTTCCACGCAAGa   <  1:6707319/35‑1 (MQ=39)

                          gacAGTCGCAATAACGGTGAGTTTTCCACGCAAGa   <  1:7089277/35‑1 (MQ=39)

                          gacAGTCGCAATAACGGTGAGTTTTCCACGCAAGa   <  1:4874895/35‑1 (MQ=39)

                           acAGTCGCAATAACGGTGAGTTTTCCACGCAAGaa  <  1:2212468/35‑1 (MQ=38)

                           acAGTCGCAATAACGGTGAGTTTTCCACGCAAGaa  <  1:518194/35‑1 (MQ=255)

                           acAGTCGCAATAACGGTGAGTTTTCCACGCAAGaa  <  1:7771385/35‑1 (MQ=255)

                           acAGTCGCAATAACGCTGAGTTTTCCACGCAAGaa  >  1:6636327/1‑35 (MQ=37)

ATCTCGAGCGCGGCTTCAATATGGACGACAGCCGCAATAACGGTGAGTTTTCCACGCAAGAA > NC_000913/1382524‑1382585

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