BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,093,223	A→G	V72V (GTT→GTC)	rhaM ←	L‑rhamnose mutarotase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,093,223	0	A	G	100.0%	258.5 / NA	65	V72V (GTT→GTC)	rhaM	L‑rhamnose mutarotase
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (34/31); total (34/31)

GGGTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAATCT  >  NC_000913/4093176‑4093269

gggTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGAc                                              <  1:51250/50‑1 (MQ=255)

gggTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGAc                                              <  1:1439735/50‑1 (MQ=255)

  gTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACAt                                            <  1:3536594/50‑1 (MQ=255)

  gTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACAt                                            >  1:912609/1‑50 (MQ=255)

  gTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACAt                                            >  1:811275/1‑50 (MQ=255)

  gTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACAt                                            <  1:770840/50‑1 (MQ=255)

   ttAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACAt                                            >  1:2691095/1‑49 (MQ=255)

    tAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACAt                                            >  1:896584/1‑48 (MQ=255)

      gCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAg                                         <  1:1379643/49‑1 (MQ=255)

      gCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAg                                         <  1:3502795/49‑1 (MQ=255)

      gCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGt                                        <  1:2045557/50‑1 (MQ=255)

        gggCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGc                                      >  1:3267986/1‑50 (MQ=255)

        gggCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGc                                      >  1:553152/1‑50 (MQ=255)

        gggCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGc                                      <  1:3126304/50‑1 (MQ=255)

         ggCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCt                                     <  1:142072/50‑1 (MQ=255)

           cATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTgg                                   >  1:3206414/1‑50 (MQ=255)

           cATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTgg                                   <  1:3707761/50‑1 (MQ=255)

           cATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTgg                                   >  1:2377771/1‑50 (MQ=255)

            aTAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTgg                                   >  1:1031057/1‑49 (MQ=255)

            aTAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGc                                  <  1:4337670/50‑1 (MQ=255)

             tAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGc                                  >  1:3600562/1‑49 (MQ=255)

              aaCATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCaa                                >  1:429018/1‑50 (MQ=255)

               aCATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCaa                                <  1:3876177/49‑1 (MQ=255)

               aCATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAAc                               >  1:1185963/1‑50 (MQ=255)

               aCATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAAc                               >  1:1135196/1‑50 (MQ=255)

                cATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCGAcc                              <  1:425886/50‑1 (MQ=255)

                cATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAAcc                              <  1:4439106/50‑1 (MQ=255)

                cATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAAcc                              <  1:3109228/50‑1 (MQ=255)

                 aTCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCg                             >  1:3394621/1‑50 (MQ=255)

                 aTCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCg                             >  1:1702124/1‑50 (MQ=255)

                 aTCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCg                             <  1:4227305/50‑1 (MQ=255)

                    ggTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCAt                          >  1:2439963/1‑50 (MQ=255)

                     gTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCAtt                         <  1:1933162/50‑1 (MQ=255)

                      tCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTc                        >  1:3772517/1‑50 (MQ=255)

                        atatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTcc                       <  1:2191852/49‑1 (MQ=255)

                        atatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTc                        >  1:3178979/1‑48 (MQ=255)

                        atatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCa                      <  1:1012869/50‑1 (MQ=255)

                           tatTTCCACCAACGTTGGCAGACCTCAGTGCTGGCAACCGCATTCCAGCg                   >  1:653164/1‑50 (MQ=39)

                           tatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCg                   >  1:2892601/1‑50 (MQ=255)

                           tatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCg                   <  1:2368166/50‑1 (MQ=255)

                             tttCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGtt                 <  1:399806/50‑1 (MQ=255)

                             tttCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGtt                 <  1:4247619/50‑1 (MQ=255)

                               tCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttct               >  1:3622300/1‑50 (MQ=255)

                               tCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttc                <  1:2038649/49‑1 (MQ=255)

                                ccaccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt              >  1:1485197/1‑50 (MQ=255)

                                 caccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt              >  1:3754043/1‑49 (MQ=255)

                                  accaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt              >  1:3204649/1‑48 (MQ=255)

                                  accaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt              >  1:3720378/1‑48 (MQ=255)

                                  accaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCa            <  1:3295043/50‑1 (MQ=255)

                                   ccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAg           <  1:2510795/50‑1 (MQ=255)

                                   ccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAg           <  1:2550961/50‑1 (MQ=255)

                                    caACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGa          <  1:627233/50‑1 (MQ=255)

                                     aaCGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAt         <  1:4077651/50‑1 (MQ=255)

                                     aaCGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAt         >  1:1073185/1‑50 (MQ=255)

                                       cGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAtt        >  1:3959022/1‑49 (MQ=255)

                                        gTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTc       >  1:3361203/1‑49 (MQ=255)

                                        gTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTc       >  1:3517915/1‑49 (MQ=255)

                                         ttGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCaa     >  1:2357394/1‑50 (MQ=255)

                                          tGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAt    >  1:2706367/1‑50 (MQ=255)

                                          tGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAt    >  1:94190/1‑50 (MQ=255)

                                           ggCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAtc   >  1:2313507/1‑50 (MQ=255)

                                           ggCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAt    <  1:4265959/49‑1 (MQ=255)

                                            gCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAtct  <  1:4099102/50‑1 (MQ=255)

                                            gCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAtct  >  1:3426946/1‑50 (MQ=255)

                                            gCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAtct  <  1:2843366/50‑1 (MQ=255)

GGGTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAATCT  >  NC_000913/4093176‑4093269

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