BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,091,305	A→G	S58G (AGC→GGC)	rsmE →	16S rRNA m(3)U1498 methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,091,305	0	A	G	100.0%	163.2 / NA	42	S58G (AGC→GGC)	rsmE	16S rRNA m(3)U1498 methyltransferase
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (23/19); total (23/19)

GTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACG  >  NC_000913/3091259‑3091350

gTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGc                                             >  1:2607619/1‑49 (MQ=255)

 tttGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCa                                            >  1:3917549/1‑49 (MQ=255)

   tGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaa                                          >  1:1688261/1‑49 (MQ=255)

    gACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaa                                         <  1:1163635/49‑1 (MQ=255)

     aCGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaa                                        >  1:1373249/1‑49 (MQ=255)

      cGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaaaa                                      <  1:3827974/50‑1 (MQ=255)

      cGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaa                                        >  1:3313831/1‑48 (MQ=255)

      cGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaa                                        >  1:1066816/1‑48 (MQ=255)

       ggTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAg                                     <  1:131760/50‑1 (MQ=255)

        gTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGc                                    >  1:3543904/1‑50 (MQ=255)

        gTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGc                                    >  1:2405318/1‑50 (MQ=255)

            cAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTg                                 >  1:1175098/1‑49 (MQ=255)

              aCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTgg                                >  1:3933575/1‑48 (MQ=255)

              aCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGaa                              <  1:2183605/50‑1 (MQ=255)

              aCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGaa                              <  1:520031/50‑1 (MQ=255)

               ccAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAg                             <  1:3653648/50‑1 (MQ=255)

                cAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGt                            >  1:319805/1‑50 (MQ=255)

                 aGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGt                            >  1:3537711/1‑49 (MQ=255)

                   gTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTg                           >  1:1484648/1‑48 (MQ=255)

                   gTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGaa                         <  1:2983783/50‑1 (MQ=255)

                    tCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAg                        <  1:312992/50‑1 (MQ=255)

                    tCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAg                        >  1:1204907/1‑50 (MQ=255)

                     cTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAgg                       >  1:2073879/1‑50 (MQ=255)

                      tttGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGt                      >  1:3625153/1‑50 (MQ=255)

                       ttGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTg                     >  1:670185/1‑50 (MQ=255)

                        tGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAgg                       <  1:894121/47‑1 (MQ=255)

                         gACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCt                   >  1:3557769/1‑50 (MQ=255)

                            gCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTgg                 >  1:2008829/1‑49 (MQ=255)

                            gCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTgg                 >  1:3801106/1‑49 (MQ=255)

                             ccGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGaa               <  1:729966/50‑1 (MQ=255)

                               gAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAgg             <  1:3498267/50‑1 (MQ=255)

                               gAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAgg             <  1:3374617/50‑1 (MQ=255)

                                  aTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc            >  1:1330898/1‑48 (MQ=255)

                                  aTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCa          <  1:2774082/50‑1 (MQ=255)

                                      ccAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGAt       <  1:3223049/49‑1 (MQ=255)

                                      ccAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGAt       <  1:652947/49‑1 (MQ=255)

                                      ccAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATc      <  1:2638008/50‑1 (MQ=255)

                                       cAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcg     <  1:3378122/50‑1 (MQ=255)

                                        aGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcga    >  1:3929620/1‑50 (MQ=255)

                                         gcgcCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcgac   >  1:2905057/1‑50 (MQ=255)

                                          cgcCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcgacg  <  1:2824549/50‑1 (MQ=255)

                                          cgcCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcgacg  <  1:996970/50‑1 (MQ=255)

GTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACG  >  NC_000913/3091259‑3091350

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