BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,614,238	G→T	L133I (CTA→ATA)	yjjI ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,614,238	0	G	T	100.0%	609.9 / NA	182	L133I (CTA→ATA)	yjjI	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (0/0); new base T (107/75); total (107/75)

Only 100 of 182 total aligned reads displayed.

AGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAGAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTAC  >  NC_000913/4614166‑4614310

aGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAAt                                                                       >  1:3554834/1‑76 (MQ=255)

 ggCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAAtt                                                                      >  2:1270564/1‑76 (MQ=255)

 ggCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAAtt                                                                      <  2:5288021/76‑1 (MQ=255)

  gCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTc                                                                     <  1:3029807/76‑1 (MQ=255)

    gTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTa                                                                   >  2:3716209/1‑76 (MQ=255)

    gTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTa                                                                   >  1:3372084/1‑76 (MQ=255)

    gTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTa                                                                   >  2:1748097/1‑76 (MQ=255)

    gTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTa                                                                   >  2:3019893/1‑76 (MQ=255)

    gTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTa                                                                   <  2:2659486/76‑1 (MQ=255)

     tCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTaa                                                                  <  2:2274489/76‑1 (MQ=255)

     tCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTaa                                                                  <  2:5330695/76‑1 (MQ=255)

      cTGGCAGTGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAAc                                                                 >  1:1054753/1‑76 (MQ=255)

        ggCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACAt                                                               >  1:3025280/1‑76 (MQ=255)

         gCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATa                                                              <  2:1411177/76‑1 (MQ=255)

         gCAGGGGTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATa                                                              >  1:1723885/1‑76 (MQ=255)

          cAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATAc                                                             <  1:4120005/76‑1 (MQ=255)

           aGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACg                                                            >  1:5474068/1‑76 (MQ=255)

            gggTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACgg                                                           <  1:2607369/76‑1 (MQ=255)

            gggTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACgg                                                           >  2:1653914/1‑76 (MQ=255)

             ggTTCTGTCGAAGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGt                                                          >  2:5574318/1‑76 (MQ=255)

              gTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGtt                                                         <  2:3045848/76‑1 (MQ=255)

               ttCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTg                                                        <  2:4361913/76‑1 (MQ=255)

                tCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGc                                                       <  1:1510481/76‑1 (MQ=255)

                tCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGc                                                       <  2:2495020/76‑1 (MQ=255)

                tCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGc                                                       <  1:4468806/76‑1 (MQ=255)

                 ttGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGca                                                      <  2:4001989/75‑1 (MQ=255)

                 cTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGca                                                      >  2:4928067/1‑76 (MQ=255)

                 cTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGca                                                      <  2:691998/76‑1 (MQ=255)

                 cTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGca                                                      >  1:5380796/1‑76 (MQ=255)

                  tGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGcaa                                                     >  2:2326605/1‑76 (MQ=255)

                    tCGAGGTAACGCCAGAAACGTTTTATTTGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGt                                                          <  1:4251925/69‑1 (MQ=255)

                    tCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGcaaca                                                   >  2:25263/1‑76 (MQ=255)

                    tCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGcaaca                                                   >  2:3330981/1‑76 (MQ=255)

                     cGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGcaacaa                                                  >  1:5619706/1‑76 (MQ=255)

                      gatgTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCaacaac                                                 <  1:4245824/73‑1 (MQ=255)

                      gAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCaacaac                                                 >  2:1655362/1‑76 (MQ=255)

                      gAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCaacaac                                                 >  1:3953706/1‑76 (MQ=255)

                      gAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCaacaac                                                 <  1:852454/76‑1 (MQ=255)

                       aGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACg                                                >  1:4714224/1‑76 (MQ=255)

                        ggTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGc                                               <  2:3237603/76‑1 (MQ=255)

                        ggTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGc                                               >  2:2950805/1‑76 (MQ=255)

                         gTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCa                                              >  2:4859990/1‑76 (MQ=255)

                          tAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCAt                                             <  2:3522997/76‑1 (MQ=255)

                            aCGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATcc                                           >  1:909269/1‑76 (MQ=255)

                             cGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCa                                          <  2:5415913/76‑1 (MQ=255)

                              gCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAg                                         <  2:1052229/76‑1 (MQ=255)

                               ccAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGt                                        >  1:2245352/1‑76 (MQ=255)

                                cAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGtt                                       >  2:5050806/1‑76 (MQ=255)

                                cAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGtt                                       >  2:2908909/1‑76 (MQ=255)

                                cAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGtt                                       >  1:4929263/1‑76 (MQ=255)

                                cAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGtt                                       <  2:2721947/76‑1 (MQ=255)

                                 aGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTg                                      <  2:440068/76‑1 (MQ=255)

                                   aaaCGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGcc                                    >  1:3105363/1‑76 (MQ=255)

                                    aaCGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGccc                                   <  2:2203584/76‑1 (MQ=255)

                                    aaCGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGccc                                   >  2:2224870/1‑76 (MQ=255)

                                       gTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAg                                >  1:2146733/1‑76 (MQ=255)

                                       gTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAg                                >  2:4143346/1‑76 (MQ=255)

                                       gTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAg                                >  2:3418585/1‑76 (MQ=255)

                                       gTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAg                                >  2:2628491/1‑76 (MQ=255)

                                       gTTTTAGTAGAACATTGAGCTCGTCGTGTGCTATAAGTCTAACACACGGATGCAACAACGCATACAGGTGCGCCAg                                >  1:543328/1‑76 (MQ=255)

                                         tttATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCATTTGCCCCAGGt                              <  1:2637093/76‑1 (MQ=255)

                                           tATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAg                            <  1:2950805/76‑1 (MQ=255)

                                            aTTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGa                           >  2:3921843/1‑76 (MQ=255)

                                             ttCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGAc                          <  1:1794367/76‑1 (MQ=255)

                                             ttCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGAc                          >  1:3983821/1‑76 (MQ=255)

                                              tCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGAcc                         >  2:2857614/1‑76 (MQ=255)

                                              tCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGAcc                         <  2:4091446/76‑1 (MQ=255)

                                              tCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGAcc                         <  2:1297044/76‑1 (MQ=255)

                                               cGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCg                        >  1:230893/1‑76 (MQ=255)

                                               cGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCg                        >  2:2081619/1‑76 (MQ=255)

                                                 aaCATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGCCCGGc                      <  1:3773072/76‑1 (MQ=255)

                                                 aaCATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGc                      >  1:1144244/1‑76 (MQ=255)

                                                 aaCATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGc                      <  2:3554834/76‑1 (MQ=255)

                                                  aCATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCa                     >  1:5174526/1‑76 (MQ=255)

                                                   cATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGCCCGGCAt                    <  1:25263/76‑1 (MQ=255)

                                                   cATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCAt                    <  1:2988482/76‑1 (MQ=255)

                                                    aTCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATc                   <  2:3372084/76‑1 (MQ=255)

                                                    aTCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATc                   >  1:4600749/1‑76 (MQ=255)

                                                    aTCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATc                   <  2:2056379/76‑1 (MQ=255)

                                                     tCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCg                  >  1:5015621/1‑76 (MQ=255)

                                                      cGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGa                 >  1:5735000/1‑76 (MQ=255)

                                                       gATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGAt                >  2:710916/1‑76 (MQ=255)

                                                       gATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGAt                <  2:3264847/76‑1 (MQ=255)

                                                       gATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGAt                >  2:4478698/1‑76 (MQ=255)

                                                        aTCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGAtg               >  2:4723324/1‑76 (MQ=255)

                                                         tctcGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGAtgt              <  1:3914116/76‑1 (MQ=255)

                                                         tctcGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGAtgt              >  2:4676047/1‑76 (MQ=255)

                                                           tcGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGa            >  1:5200270/1‑76 (MQ=255)

                                                            cGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGAc           >  1:4337647/1‑76 (MQ=255)

                                                             gTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATTTGAcc          >  2:79243/1‑76 (MQ=255)

                                                             gTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGAcc          >  2:2822867/1‑76 (MQ=255)

                                                             gTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGAcc          <  1:3257845/76‑1 (MQ=255)

                                                             gTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCACTTGCCCCAGGTAGACCGGCATCGATGTGAcc          >  2:3086442/1‑76 (MQ=255)

                                                              tCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCg         >  2:3777912/1‑76 (MQ=255)

                                                                ttGTGTTATAATTCTAACATCCGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGAc       <  2:3959490/76‑1 (MQ=255)

                                                                ttGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGAc       >  1:1158694/1‑76 (MQ=255)

                                                                 tgtgTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACg      <  1:5556141/76‑1 (MQ=255)

                                                                   tgtTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGt    <  1:710916/76‑1 (MQ=255)

                                                                    gtTATAATTCTAACATAGTGTTTCAACAACGCATCCAGTTGCCCCAGGGAGACCGGCATCGATGTGACCGACGGTa   <  1:2764938/76‑1 (MQ=255)

                                                                     ttATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTAc  <  1:1435705/76‑1 (MQ=255)

AGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAGAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTAC  >  NC_000913/4614166‑4614310

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 12 ≤` `ATCG/ATCG` `< 29 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`