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

Only 100 of 104 total aligned reads displayed.

GGGTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAATCT  >  NC_000913/4093176‑4093269

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

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

 ggTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACa                                             <  1:779374/50‑1 (MQ=255)

 ggTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACa                                             <  1:6027322/50‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

       cGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTg                                       <  1:4400726/50‑1 (MQ=255)

       cGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTg                                       <  1:7748853/50‑1 (MQ=255)

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

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

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

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

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

          gCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTg                                    <  1:3706123/50‑1 (MQ=255)

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

          gCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTg                                    <  1:1081084/50‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

               aCATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTGCTGGCAAc                               <  1:8053548/50‑1 (MQ=255)

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

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

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

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

                 aTCGGTCATATATTTCCACCAACGTTGGCAGACATCAGGGCTGGCAACCg                             >  1:4905539/1‑50 (MQ=39)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                           tatTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGc                    <  1:5432858/49‑1 (MQ=255)

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

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

                            atTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGt                  >  1:5409424/1‑50 (MQ=255)

                            atTTCCACCAACGTTGGCAGACATCAGTGCTGGCAACCGCATTCCAGCGt                  >  1:1407472/1‑50 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                            gCAGACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAAtc   <  1:8157294/49‑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: `‑`