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

TTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAATCT  >  NC_000913/4093179‑4093269

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

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

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

       gCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTg                                    >  1:3064794/1‑50 (MQ=255)

       gCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTg                                    >  1:622210/1‑50 (MQ=255)

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

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

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

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

               tCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCg                             >  1:453811/1‑49 (MQ=255)

                cGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCa                           >  1:935143/1‑50 (MQ=255)

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

                  gTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCAtt                         >  1:515421/1‑50 (MQ=255)

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

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

                    cATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTc                        >  1:2843802/1‑49 (MQ=255)

                    cATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTc                        >  1:1467166/1‑49 (MQ=255)

                     atatatTTCCCCCAACGTTGGCAGACATCAGTGCTGGCCACCGCATTCCa                      <  1:1503098/50‑1 (MQ=38)

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

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

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

                      tatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAg                     >  1:1992607/1‑50 (MQ=255)

                      tatatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAg                     <  1:900675/50‑1 (MQ=255)

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

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

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

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

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

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

                          tttCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGtt                 >  1:781180/1‑50 (MQ=255)

                          tttCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGt                  <  1:2301481/49‑1 (MQ=255)

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

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

                             ccaccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttct               <  1:2779626/49‑1 (MQ=255)

                              caccaACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttcttc             <  1:589468/50‑1 (MQ=255)

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

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

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

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

                                 caACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttctt              >  1:379248/1‑46 (MQ=255)

                                 caACGTTGGCAGAAATCAGTGCTGGCAACCGCATTCCAGCGttctt              >  1:2587051/1‑46 (MQ=38)

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

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

                                   aCGTTGGCAGACATCAGTGCTGGCTACCGCATTCCAGCGTTCTTCAGAtt        <  1:3377321/50‑1 (MQ=255)

                                   aCGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGAtt        <  1:1364254/50‑1 (MQ=255)

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

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

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

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

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

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

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

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

TTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAATCT  >  NC_000913/4093179‑4093269

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