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%	160.6 / NA	40	V72V (GTT→GTC)	rhaM	L‑rhamnose mutarotase
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (16/24); total (16/24)

GGTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAATCT  >  NC_000913/4093177‑4093269

ggTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACa                                             >  1:2725909/1‑50 (MQ=255)

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

  ttAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATc                                           <  1:3471498/50‑1 (MQ=255)

    aGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCa                                          <  1:2654424/49‑1 (MQ=255)

    aGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAg                                         >  1:936199/1‑50 (MQ=255)

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

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

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

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

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

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

        ggCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCt                                     >  1:1671348/1‑50 (MQ=255)

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

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

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

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

            tAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCa                                 <  1:89919/50‑1 (MQ=255)

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

                   ggTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCAt                          <  1:1125934/50‑1 (MQ=255)

                   ggTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCAt                          <  1:4110880/50‑1 (MQ=255)

                         atatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGc                    >  1:2181920/1‑50 (MQ=255)

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

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

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

                                caccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt              <  1:407386/49‑1 (MQ=255)

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

                                caccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt              <  1:515887/49‑1 (MQ=255)

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

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

                                  ccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt              >  1:207781/1‑47 (MQ=255)

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

                                   caACGTTGGCAGACATCAGTGCTGGCGACCGCATTCCAGCGTTCTTCAGa          <  1:2628336/50‑1 (MQ=255)

                                     aCGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAtt        >  1:4237403/1‑50 (MQ=255)

                                      cGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTc       <  1:4290269/50‑1 (MQ=255)

                                      cGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTc       <  1:1949679/50‑1 (MQ=255)

                                       gTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCa      <  1:454778/50‑1 (MQ=255)

                                        ttGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCaa     <  1:2770436/50‑1 (MQ=255)

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

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

                                           gCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAtc   <  1:1024757/49‑1 (MQ=255)

GGTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAATCT  >  NC_000913/4093177‑4093269

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