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	99.3%	587.0 / NA	146	V72V (GTT→GTC)	rhaM	L‑rhamnose mutarotase
Reads supporting (aligned to +/- strand): ref base A (0/0); major base G (55/90); minor base T (0/1); total (55/91)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 146 total aligned reads displayed.

GGGTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAATCT  >  NC_000913/4093176‑4093269

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

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

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

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

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

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

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

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

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

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

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

       cGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTg                                       >  2:3028744/1‑50 (MQ=255)

       cGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAGACATCAGTg                                       >  2:1435718/1‑50 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                  tCGGTCAGATATTTTCACCAACGTTGGCAGAAATCAGTGCTGGCAACCGc                            <  2:3984101/50‑1 (MQ=37)

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

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

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

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

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

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

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

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

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

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

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

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

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

                             tttCCACCAACGTTGGCAGACATCAGTGCTGGCCACCGCATTCCAGCGt                  <  1:3105307/49‑1 (MQ=255)

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

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

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

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

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

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

                                ccaccaACGTTGGCAGACATCAGTGCTGGCCACCGCATTCCAGCGttct               <  1:319994/49‑1 (MQ=255)

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

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

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

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

                                  accaACGTTGGCAGACATCAGTGCTGGCGACCGCATTCCAGCGTTCTTCa            <  1:1074566/50‑1 (MQ=255)

                                  accaACGTTGGCAGACATCAGTGCTGGCCGCCGCATTCCAGCGTTCTTCa            <  2:803764/50‑1 (MQ=38)

                                  accaACGTTGGCAGACATCAGTGCTGGCAGCCGCATTCCAGCGTTCTTCa            <  2:332388/50‑1 (MQ=255)

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

                                  accaACGTTGGCAGACATCAGTGCTGGCAACCGCATGCCAGCGTTCTTCa            <  2:2767929/50‑1 (MQ=38)

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

                                    caACGTTGGCAGACATCAGTGCTGGCTACCGCATTCCAGCGTTCTTCAGa          <  2:3164204/50‑1 (MQ=255)

                                     aaCGTTGGCAGACATCAGTGCTGGCATCCGCATTCCAGCGTTCTTCAGAt         <  2:1020173/50‑1 (MQ=255)

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

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

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

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

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

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

                                        gTTGGCAGACATCAGTGCTGGCTACCGCATTCCAGCGTTCTTCAGATTCa      <  2:3205401/50‑1 (MQ=255)

                                        gTTGGCAGACATCAGTGCTGGCAGCCGCATTCCAGCGTTCTTCAGATTCa      <  2:115800/50‑1 (MQ=255)

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

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

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

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

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

                                           ggCATACATCGGTGCTGGCAGCCGCATTCCAGCGTTCTTCAGATTCAAt    <  1:3894317/49‑1 (MQ=38)

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

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

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

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

                                            gCAGACATCAGTGCTGGCGACCGCATTCCAGCGTTCTTCAGATTCAAtc   <  2:2905360/49‑1 (MQ=39)

                                            gCAGACATCAGTGCTGGCAGCCGCATTCCAGCGTTCTTCAGATTCAAtct  <  2:2243814/50‑1 (MQ=255)

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

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

GGGTTAGCGGGCATAACATCGGTCATATATTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGCATTCCAGCGTTCTTCAGATTCAATCT  >  NC_000913/4093176‑4093269

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