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

GGGTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAATC  >  NC_000913/4093176‑4093268

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

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

  gTTAGCGAGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACAt                                           <  1:246836/50‑1 (MQ=39)

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

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

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

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

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

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

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

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

                cATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAAcc                             >  1:2079070/1‑50 (MQ=255)

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

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

                  tCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGc                           <  1:2440350/50‑1 (MQ=255)

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

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

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

                      tCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTc                       <  1:2570910/50‑1 (MQ=255)

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

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

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

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

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

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

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

                              ttCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGttc               >  1:1657064/1‑50 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

GGGTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAATC  >  NC_000913/4093176‑4093268

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: `‑`