BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,213,803	T→C	S253P (TCT→CCT)	rpoD →	RNA polymerase, sigma 70 (sigma D) factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,213,803	0	T	C	100.0%	535.9 / NA	130	S253P (TCT→CCT)	rpoD	RNA polymerase, sigma 70 (sigma D) factor
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (66/64); total (66/64)

Only 100 of 130 total aligned reads displayed.

CGAAAGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTA  >  NC_000913/3213756‑3213849

cGAAAGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCt                                              >  1:8601832/1‑50 (MQ=255)

cGAAAGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCt                                              >  1:8561183/1‑50 (MQ=255)

cGAAAGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCt                                              >  1:4204348/1‑50 (MQ=255)

 gAAAGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCt                                              >  1:462553/1‑49 (MQ=255)

 gAAAGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCt                                              >  1:1938315/1‑49 (MQ=255)

   aaGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGaa                                           <  1:7083691/50‑1 (MQ=255)

   aaGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGaa                                           <  1:3832783/50‑1 (MQ=255)

    aGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAg                                          <  1:1201542/50‑1 (MQ=255)

    aGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAg                                          <  1:5775854/50‑1 (MQ=255)

     ggTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAg                                          >  1:8225912/1‑49 (MQ=255)

     ggTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGt                                         <  1:7969413/50‑1 (MQ=255)

      gTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTa                                        <  1:4408550/50‑1 (MQ=255)

      gTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTa                                        <  1:5633669/50‑1 (MQ=255)

       tCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTAt                                       >  1:5188155/1‑50 (MQ=255)

       tCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTAt                                       <  1:2594815/50‑1 (MQ=255)

        cGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTAtt                                      >  1:1738170/1‑50 (MQ=255)

        cGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTAtt                                      >  1:5381373/1‑50 (MQ=255)

        cGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTAtt                                      >  1:7373541/1‑50 (MQ=255)

        cGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTAtt                                      >  1:7941314/1‑50 (MQ=255)

         gCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTAtt                                      <  1:2889451/49‑1 (MQ=255)

         gCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTAtt                                      >  1:4650230/1‑49 (MQ=255)

         gCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTAtt                                      <  1:1373625/49‑1 (MQ=255)

          cAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTc                                     >  1:6554242/1‑49 (MQ=255)

          cAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTc                                     >  1:6349339/1‑49 (MQ=255)

          cAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTc                                     >  1:4421660/1‑49 (MQ=255)

          cAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCa                                    <  1:6852459/50‑1 (MQ=255)

          cAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCa                                    <  1:7435663/50‑1 (MQ=255)

          cAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCa                                    <  1:1134800/50‑1 (MQ=255)

          cAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCa                                    <  1:3764158/50‑1 (MQ=255)

           aGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCaa                                   <  1:6134276/50‑1 (MQ=255)

           aGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCaa                                   <  1:4341878/50‑1 (MQ=255)

           aGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCa                                    >  1:6016029/1‑49 (MQ=255)

            gTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCaaa                                  <  1:5238164/50‑1 (MQ=255)

            gTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCaaa                                  <  1:3175873/50‑1 (MQ=255)

            gTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCaaa                                  >  1:6331848/1‑50 (MQ=255)

            gTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCaaa                                  <  1:4739669/50‑1 (MQ=255)

            gTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCaaa                                  <  1:6941502/50‑1 (MQ=255)

            gTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCaaa                                  <  1:823061/50‑1 (MQ=255)

             tCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCaaa                                  >  1:1778866/1‑49 (MQ=255)

             tCACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAAc                                 <  1:8417022/50‑1 (MQ=255)

              cACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCaaa                                  >  1:5841583/1‑48 (MQ=255)

              cACGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCaaa                                  >  1:5426203/1‑48 (MQ=255)

               aCGTTACCGCTCAGTAAGTGATCCTGAAACTGCCTGAAGTATTCAAACa                                <  1:8078543/49‑1 (MQ=37)

               aCGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAg                               >  1:2217525/1‑50 (MQ=255)

               aCGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAg                               <  1:1717835/50‑1 (MQ=255)

               aCGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAg                               >  1:5724398/1‑50 (MQ=255)

                cGCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGt                              >  1:6613144/1‑50 (MQ=255)

                 gCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGtt                             <  1:5094419/50‑1 (MQ=255)

                 gCTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGtt                             >  1:4074718/1‑50 (MQ=255)

                  cTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTc                            >  1:6698886/1‑50 (MQ=255)

                  cTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTc                            >  1:1592394/1‑50 (MQ=255)

                  cTACCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTc                            <  1:4561287/50‑1 (MQ=255)

                    aCCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCg                          >  1:7522261/1‑50 (MQ=255)

                    aCCGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCg                          >  1:6722059/1‑50 (MQ=255)

                     ccGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTcc                           <  1:8475391/48‑1 (MQ=255)

                     ccGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCg                          <  1:497973/49‑1 (MQ=255)

                     ccGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGc                         >  1:5062055/1‑50 (MQ=255)

                      cGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGcc                        <  1:1714908/50‑1 (MQ=255)

                      cGCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGcc                        <  1:3717226/50‑1 (MQ=255)

                       gCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCt                       >  1:5607629/1‑50 (MQ=255)

                       gCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCt                       >  1:5186001/1‑50 (MQ=255)

                       gCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCt                       <  1:250777/50‑1 (MQ=255)

                       gCTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCt                       >  1:6741956/1‑50 (MQ=255)

                        cTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTg                      >  1:7764747/1‑50 (MQ=255)

                        cTCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTg                      <  1:5073961/50‑1 (MQ=255)

                         tCAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTgg                     >  1:6052727/1‑50 (MQ=255)

                          cAGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGt                    >  1:7738562/1‑50 (MQ=255)

                           aGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTgg                     <  1:7789297/48‑1 (MQ=255)

                           aGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTg                   <  1:1243906/50‑1 (MQ=255)

                           aGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTg                   <  1:3241461/50‑1 (MQ=255)

                           aGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTg                   >  1:1979708/1‑50 (MQ=255)

                           aGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTg                   <  1:6471724/50‑1 (MQ=255)

                           aGGAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTg                   <  1:7571585/50‑1 (MQ=255)

                             gAAGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGcc                 <  1:8289545/50‑1 (MQ=255)

                              aaGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGcc                 >  1:6238541/1‑49 (MQ=255)

                              aaGAGATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCg                <  1:1820905/50‑1 (MQ=255)

                               agagATCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCg                >  1:3320178/1‑49 (MQ=255)

                                gagaTCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGaa              <  1:7376996/50‑1 (MQ=255)

                                 agaTCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAg             <  1:1859369/50‑1 (MQ=255)

                                  gaTCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAg             >  1:7632379/1‑49 (MQ=255)

                                  gaTCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAg             >  1:5845449/1‑49 (MQ=255)

                                  gaTCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGc            <  1:4687871/50‑1 (MQ=255)

                                  gaTCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGc            <  1:2526493/50‑1 (MQ=255)

                                  gaTCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGc            <  1:1472166/50‑1 (MQ=255)

                                   aTCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAg             >  1:8532857/1‑48 (MQ=255)

                                   aTCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAg             >  1:1215974/1‑48 (MQ=255)

                                    tCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAg          >  1:1745564/1‑50 (MQ=255)

                                    tCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAg          <  1:8426060/50‑1 (MQ=255)

                                    tCCTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAg          >  1:6197251/1‑50 (MQ=255)

                                     ccTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGt         >  1:1336425/1‑50 (MQ=255)

                                      cTGAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGtt        <  1:6481379/50‑1 (MQ=255)

                                        gAAACTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCTGAAGCAGttt       >  1:8227017/1‑49 (MQ=38)

                                         aaaCTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGttt       <  1:5046634/48‑1 (MQ=255)

                                         aaaCTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGa     <  1:4737007/50‑1 (MQ=255)

                                          aaCTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGAc    <  1:4170184/50‑1 (MQ=255)

                                          aaCTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGAc    >  1:3120711/1‑50 (MQ=255)

                                           aCTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGAc    >  1:3276031/1‑49 (MQ=255)

                                           aCTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACt   <  1:5216082/50‑1 (MQ=255)

                                           aCTGCCTGAAGTATTCAAACAGCTCCGCCTGGTGCCCAAGCAGTTTGACt   <  1:7638900/50‑1 (MQ=38)

                                            cTGCCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTa  >  1:5330732/1‑50 (MQ=255)

CGAAAGGTCGCAGTCACGCTACCGCTCAGGAAGAGATCCTGAAACTGTCTGAAGTATTCAAACAGTTCCGCCTGGTGCCGAAGCAGTTTGACTA  >  NC_000913/3213756‑3213849

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