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%	169.9 / NA	43	K90K (AAG→AAA)	fkpA	peptidyl‑prolyl cis‑trans isomerase FkpA
Reads supporting (aligned to +/- strand): ref base C (0/0); new base T (10/33); total (10/33)

GAATGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTCTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCTT  >  NC_000913/3477105‑3477193

gAATGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTa                                          >  1:2327462/1‑49 (MQ=255)

gAATGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTa                                          >  1:174766/1‑49 (MQ=255)

 aaTGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTAt                                         <  1:2198175/49‑1 (MQ=255)

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

 aaTGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTAt                                         <  1:1371472/49‑1 (MQ=255)

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

   tGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTAt                                         >  1:1950262/1‑47 (MQ=255)

   tGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCa                                       <  1:1924290/49‑1 (MQ=255)

   tGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCa                                       <  1:2422722/49‑1 (MQ=255)

   tGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCa                                       <  1:1024483/49‑1 (MQ=255)

     cTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAg                                      >  1:575804/1‑48 (MQ=255)

     cTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAg                                      <  1:2485236/48‑1 (MQ=255)

     cTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAg                                      <  1:1010537/48‑1 (MQ=255)

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

        gTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCa                                    <  1:2082534/47‑1 (MQ=255)

        gTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCa                                    <  1:621104/47‑1 (MQ=255)

        gTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCa                                    <  1:1090352/47‑1 (MQ=255)

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

            agTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCa                              <  1:2374107/49‑1 (MQ=255)

              tCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATc                            <  1:80146/49‑1 (MQ=255)

              tCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATc                            <  1:2290635/49‑1 (MQ=255)

                   tCGATCTCTTGGCCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGaa                       >  1:1641609/1‑49 (MQ=39)

                    cGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGaa                       <  1:2107104/48‑1 (MQ=255)

                    cGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGaa                       <  1:638588/48‑1 (MQ=255)

                    cGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGaa                       <  1:1146714/48‑1 (MQ=255)

                    cGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGaa                       <  1:413982/48‑1 (MQ=255)

                     gATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAac                      <  1:1127773/48‑1 (MQ=255)

                      aTCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAacac                    >  1:217755/1‑49 (MQ=255)

                      aTCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAacac                    >  1:2387931/1‑49 (MQ=255)

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

                           ttGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAg                 <  1:2166340/47‑1 (MQ=255)

                           ttGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAg                 <  1:2808592/47‑1 (MQ=255)

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

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

                               tCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCa           <  1:2220197/49‑1 (MQ=255)

                                 ggAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGc         <  1:385397/49‑1 (MQ=255)

                                 ggAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGc         <  1:813075/49‑1 (MQ=255)

                                  gagaGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCt        >  1:547065/1‑49 (MQ=255)

                                  gagaGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCt        <  1:2593079/49‑1 (MQ=255)

                                  gagaGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCt        <  1:2538498/49‑1 (MQ=255)

                                  gagaGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCt        <  1:204767/49‑1 (MQ=255)

                                       tttGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCt   <  1:2732340/49‑1 (MQ=255)

                                         tGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCtt  <  1:2145034/48‑1 (MQ=255)

GAATGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTCTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCTT  >  NC_000913/3477105‑3477193

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