BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,477,150	C→T	K90K (AAG→AAA)	fkpA ←	peptidyl‑prolyl cis‑trans isomerase FkpA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,477,150	0	C	T	100.0%	430.2 / NA	109	K90K (AAG→AAA)	fkpA	peptidyl‑prolyl cis‑trans isomerase FkpA
Reads supporting (aligned to +/- strand): ref base C (0/0); new base T (37/72); total (37/72)

Only 100 of 109 total aligned reads displayed.

GAATGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTCTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCTTTA  >  NC_000913/3477105‑3477195

gAATGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTAt                                           <  1:6202021/50‑1 (MQ=255)

gAATGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTAt                                           >  1:170381/1‑50 (MQ=255)

 aaTGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTAt                                           >  1:3461534/1‑49 (MQ=255)

 aaTGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTAt                                           >  1:5542409/1‑49 (MQ=255)

 aaTGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATc                                          <  1:1225091/50‑1 (MQ=255)

 aaTGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATc                                          <  1:2420609/50‑1 (MQ=255)

 aaTGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATc                                          <  1:6632152/50‑1 (MQ=255)

 aaTGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATc                                          <  1:2567185/50‑1 (MQ=255)

  aTGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTAt                                           >  1:710231/1‑48 (MQ=255)

  aTGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCa                                         <  1:788364/50‑1 (MQ=255)

  aTGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCa                                         <  1:5590283/50‑1 (MQ=255)

   tGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAg                                        <  1:1057424/50‑1 (MQ=255)

   tGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAg                                        <  1:1464425/50‑1 (MQ=255)

    gCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAg                                        >  1:2803740/1‑49 (MQ=255)

    gCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAg                                        >  1:5794450/1‑49 (MQ=255)

    gCTTGTAGAGTCTGCTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAg                                        >  1:416021/1‑49 (MQ=255)

     cTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGc                                       <  1:85861/49‑1 (MQ=255)

     cTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGc                                       <  1:4913162/49‑1 (MQ=255)

     cTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGc                                       >  1:1278999/1‑49 (MQ=255)

     cTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCa                                      <  1:4236469/50‑1 (MQ=255)

     cTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCa                                      <  1:2127373/50‑1 (MQ=255)

     cTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCa                                      <  1:807139/50‑1 (MQ=255)

        gTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCaaa                                    <  1:2523385/49‑1 (MQ=255)

        gTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAAt                                   >  1:2803377/1‑50 (MQ=255)

        gTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAAt                                   >  1:4506952/1‑50 (MQ=255)

         tAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAAt                                   <  1:3102616/49‑1 (MQ=255)

         tAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATg                                  <  1:6699494/50‑1 (MQ=255)

          agagTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATg                                  >  1:2501045/1‑49 (MQ=255)

           gagTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGc                                 <  1:834098/49‑1 (MQ=255)

            agTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCAt                               <  1:1687922/50‑1 (MQ=255)

            agTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCAt                               <  1:1657387/50‑1 (MQ=255)

            agTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCAt                               <  1:3996964/50‑1 (MQ=255)

            agTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCAt                               <  1:5308200/50‑1 (MQ=255)

             gTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATc                              >  1:3167759/1‑50 (MQ=255)

             gTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATc                              <  1:2224830/50‑1 (MQ=255)

             gTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATc                              <  1:5111050/50‑1 (MQ=255)

              tCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATcg                             <  1:4711874/50‑2 (MQ=255)

              tCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATcc                             >  1:6307790/1‑50 (MQ=255)

                 gTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTg                           >  1:6254424/1‑49 (MQ=255)

                  ttCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCt                            <  1:4461601/47‑1 (MQ=255)

                  ttCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTg                           <  1:1414663/48‑1 (MQ=255)

                  ttCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTg                           <  1:4796854/48‑1 (MQ=255)

                  ttCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGaa                         >  1:6498794/1‑50 (MQ=255)

                  ttCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGa                          <  1:5734066/49‑1 (MQ=255)

                  ttCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGa                          <  1:4038636/49‑1 (MQ=255)

                   tCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGaa                         >  1:5322787/1‑49 (MQ=255)

                   tCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGaa                         >  1:1807893/1‑49 (MQ=255)

                    cGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAaca                       >  1:255610/1‑50 (MQ=255)

                     gATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAacac                      >  1:1289475/1‑50 (MQ=255)

                     gATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAacac                      >  1:1266919/1‑50 (MQ=255)

                     gATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAacac                      >  1:3764588/1‑50 (MQ=255)

                     gATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAaca                       <  1:5903236/49‑1 (MQ=255)

                         tctTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCTTGAACACCAGc                  >  1:4916033/1‑50 (MQ=255)

                          ctTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGc                  >  1:2054838/1‑49 (MQ=255)

                           ttGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAacac                      <  1:3958595/44‑1 (MQ=255)

                           ttGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACAcc                     <  1:21185/45‑1 (MQ=255)

                           ttGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACAcc                     <  1:5872932/45‑1 (MQ=255)

                           ttGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCac                   <  1:2556370/47‑2 (MQ=255)

                           ttGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCa                    <  1:6205825/46‑1 (MQ=255)

                           ttGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCa                    <  1:6193163/46‑1 (MQ=255)

                           ttGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCa                    <  1:555714/46‑1 (MQ=255)

                           ttGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGc                  <  1:425375/48‑1 (MQ=255)

                           ttGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGc                  <  1:1658452/48‑1 (MQ=255)

                           ttGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGc                  <  1:1541161/48‑1 (MQ=255)

                            tGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGAt               >  1:4737829/1‑50 (MQ=255)

                             ggTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACGCCAGCGAt               >  1:1650829/1‑49 (MQ=39)

                             ggTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCg                 <  1:3113848/47‑1 (MQ=255)

                             ggTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCg                 <  1:3078377/47‑1 (MQ=255)

                             ggTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGa                <  1:5602394/48‑1 (MQ=255)

                             ggTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGAt               >  1:3416065/1‑49 (MQ=255)

                              gTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATc              <  1:283796/49‑1 (MQ=255)

                              gTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATc              <  1:2567285/49‑1 (MQ=255)

                              gTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCa             <  1:2226077/50‑1 (MQ=255)

                               tCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAg            <  1:6415860/50‑1 (MQ=255)

                               tCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAg            <  1:2808300/50‑1 (MQ=255)

                                cGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGc           >  1:4493994/1‑50 (MQ=255)

                                 ggAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCt          <  1:2393750/50‑1 (MQ=255)

                                 ggAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCt          <  1:4460643/50‑1 (MQ=255)

                                 ggAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCt          <  1:1800444/50‑1 (MQ=255)

                                 ggAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCt          <  1:1064815/50‑1 (MQ=255)

                                  gagaGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTg         <  1:26881/50‑1 (MQ=255)

                                   agagTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGa        >  1:3100848/1‑50 (MQ=255)

                                   agagTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGa        <  1:3814911/50‑1 (MQ=255)

                                    gagTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGAt       <  1:3618719/50‑1 (MQ=255)

                                    gagTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGAt       <  1:111886/50‑1 (MQ=255)

                                    gagTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGAt       <  1:44433/50‑1 (MQ=255)

                                    gagTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGAt       >  1:833817/1‑50 (MQ=255)

                                     agTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATc      >  1:4353020/1‑50 (MQ=255)

                                     agTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATc      <  1:254830/50‑1 (MQ=255)

                                     agTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATc      <  1:5039368/50‑1 (MQ=255)

                                      gTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCt     >  1:1207134/1‑50 (MQ=255)

                                      gTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCt     >  1:2840269/1‑50 (MQ=255)

                                      gTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCt     <  1:509416/50‑1 (MQ=255)

                                       tttGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCtt    <  1:701176/50‑1 (MQ=255)

                                       tttGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCtt    <  1:1113824/50‑1 (MQ=255)

                                       tttGCTTTTATCAGCAAATGCATACTGAACACCAGCGATCAGCTGATCtt    <  1:6209689/50‑1 (MQ=38)

                                         tGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCttt   <  1:5470670/49‑1 (MQ=255)

                                         tGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCTTTa  >  1:5602830/1‑50 (MQ=255)

                                          gCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCTTTa  <  1:5803889/49‑1 (MQ=255)

                                          gCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCTTTa  <  1:1327249/49‑1 (MQ=255)

GAATGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTCTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCTTTA  >  NC_000913/3477105‑3477195

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