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	93.0%	348.1 / 11.1	100	S203S (AGC→AGT)	ompG	outer membrane porin G
Reads supporting (aligned to +/- strand): ref base C (4/3); new base T (49/44); total (53/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

Only 100 of 101 total aligned reads displayed.

CCTTGCGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGCCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTATTTGC  >  NC_000913/1382506‑1382601

ccTTGCGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGc                                                >  1:2993902/1‑50 (MQ=255)

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

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

    gCGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCg                                              <  1:4008499/48‑1 (MQ=255)

    gCGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCg                                              <  2:2776691/48‑1 (MQ=255)

    gCGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGc                                             <  1:3533782/49‑1 (MQ=255)

    gCGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGc                                             >  1:4206882/1‑49 (MQ=255)

    gCGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGc                                             <  2:4347570/49‑1 (MQ=255)

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

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

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

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

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

      gAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGCCGCAAt                                          <  1:3363050/50‑1 (MQ=255)

        gTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATaa                                        >  1:1983474/1‑50 (MQ=255)

        gTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATaa                                        <  2:4013263/50‑1 (MQ=255)

        gTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATaa                                        >  1:3578054/1‑50 (MQ=255)

         tGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAAc                                       <  1:931206/50‑1 (MQ=255)

         tGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAAc                                       <  2:5202820/50‑1 (MQ=255)

         tGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAAc                                       <  1:3723454/50‑1 (MQ=255)

         tGAACTATTATCGCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAAc                                       <  2:1690849/50‑1 (MQ=255)

          gAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACg                                      >  2:614649/1‑50 (MQ=255)

           aaCTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACg                                      <  1:4447965/49‑1 (MQ=255)

           aaCTATTATCTCGAGCGCGGCTTCAATATGGACGACAGCCGCAATAACgg                                     <  1:3689495/50‑1 (MQ=255)

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

            aCTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGt                                    >  1:3520613/1‑50 (MQ=255)

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

             cTATTATCTCGAGCGCGGCTTCAATATGGACGACAGCCGCAATAACGGTg                                   <  1:3917382/50‑1 (MQ=255)

              tattatCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGt                                    >  2:2167413/1‑48 (MQ=255)

              tattatCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGt                                    >  2:140775/1‑48 (MQ=255)

              tattatCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGt                                    >  1:5320884/1‑48 (MQ=255)

              tattatCTCGAGCGCGGCTTCAATATGGACGACAGCCGCAATAACGGt                                    >  2:764004/1‑48 (MQ=255)

               attatCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAg                                 >  1:3450602/1‑50 (MQ=255)

                ttatCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGt                                >  2:1549883/1‑50 (MQ=255)

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

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

                  atCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGtt                               <  2:2736562/49‑1 (MQ=255)

                   tctcGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGtttt                             >  2:697611/1‑50 (MQ=255)

                   tctcGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGtttt                             >  2:1175016/1‑50 (MQ=255)

                      cGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCa                          >  2:3237061/1‑50 (MQ=255)

                       gAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTcc                           <  2:3461601/48‑1 (MQ=255)

                       gAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTcc                           <  2:3854657/48‑1 (MQ=255)

                       gAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTc                            <  1:3475135/47‑1 (MQ=255)

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

                       gAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCAc                         >  1:916542/1‑50 (MQ=255)

                       gAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCAc                         >  1:1287126/1‑50 (MQ=255)

                       gAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCAc                         >  1:2431682/1‑50 (MQ=255)

                       gAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCAc                         <  1:866002/50‑1 (MQ=255)

                       gAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCAc                         <  1:3552994/50‑1 (MQ=255)

                       gAGCGCGGCTTCAATATGGACGACAGTCACAATAAGGGTGAGtttt                             <  2:425922/46‑1 (MQ=255)

                        aGCGCGGCTTCAATATGGACGACAGTCGGAATAACGGTGAGTTTTCCACg                        >  2:3171305/1‑50 (MQ=255)

                        aGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACg                        <  2:1330446/50‑1 (MQ=255)

                         gcgcGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGc                       <  1:4637207/50‑1 (MQ=255)

                           gcgGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGc                       >  2:22309/1‑48 (MQ=255)

                           gcgGCTCCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGc                       >  2:2750818/1‑48 (MQ=255)

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

                            cgGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAg                    <  1:169624/50‑1 (MQ=255)

                             ggCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAg                    >  2:3443814/1‑49 (MQ=255)

                             ggCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAg                    >  1:4487414/1‑49 (MQ=255)

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

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

                               cTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGaaa                 <  1:1267206/50‑1 (MQ=255)

                               cTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGaaa                 >  2:5212276/1‑50 (MQ=255)

                               cTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGaaa                 <  1:2765060/50‑1 (MQ=255)

                                ttCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAAt                >  1:1913988/1‑50 (MQ=255)

                                ttCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAAt                >  2:5028467/1‑50 (MQ=255)

                                 tCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAAtt               >  2:2441715/1‑50 (MQ=255)

                                 tCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAAtt               >  2:5313483/1‑50 (MQ=255)

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

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

                                  cAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTc              >  1:3044937/1‑50 (MQ=255)

                                   aaTATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTcg             >  2:690863/1‑50 (MQ=255)

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

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

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

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

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

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

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

                                     tatGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTcgcg           <  1:563934/50‑1 (MQ=255)

                                      atGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCgcgc          >  1:2001772/1‑50 (MQ=255)

                                       tGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGcc         >  1:1483604/1‑50 (MQ=255)

                                       tGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGcc         <  2:4543980/50‑1 (MQ=255)

                                       tGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGcc         >  2:4690079/1‑50 (MQ=255)

                                       tGGACGACAGCCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGcc         >  2:5043509/1‑50 (MQ=255)

                                       tGGACGACAGCCGCAATAACGGTGAGTTTTCAACGCAAGAAATTCGCGcc         >  2:4899529/1‑50 (MQ=255)

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

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

                                         gacgacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTa       <  1:5165843/50‑1 (MQ=255)

                                         gacgacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTa       >  2:1495550/1‑50 (MQ=255)

                                          acgacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTAt      >  2:1489691/1‑50 (MQ=255)

                                           cgacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTAtt     <  2:3975604/50‑1 (MQ=255)

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

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

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

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

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

                                             acAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTATTTg   >  1:3192547/1‑50 (MQ=255)

                                              cAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTATTTGc  >  2:3216509/1‑50 (MQ=255)

                                              cAGCCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTATTTGc  >  2:5202782/1‑50 (MQ=255)

CCTTGCGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGCCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTATTTGC  >  NC_000913/1382506‑1382601

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