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

GAATGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTCTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCTT  >  NC_000913/3477105‑3477193

gAATGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTa                                          <  1:1483021/49‑1 (MQ=255)

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

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

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

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

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

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

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

      ttGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGc                                     >  1:1149753/1‑48 (MQ=255)

        gTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCaa                                   <  1:2027917/48‑1 (MQ=255)

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

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

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

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

             gTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCAt                             <  1:370201/49‑1 (MQ=255)

                 gTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCt                          >  1:774972/1‑48 (MQ=255)

                  ttCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATcc                           <  1:451608/46‑1 (MQ=255)

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

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

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

                     gATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAaca                     >  1:1450750/1‑49 (MQ=255)

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

                       tctcTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAacac                    >  1:912427/1‑48 (MQ=255)

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

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

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

                             ggTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGc                <  1:1958862/46‑1 (MQ=255)

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

                              gTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGAt             >  1:1211921/1‑48 (MQ=255)

                              gTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGAt             <  1:1667509/48‑1 (MQ=255)

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

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

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

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

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

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

                                 ggAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGc         >  1:1742902/1‑49 (MQ=255)

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

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

                                      gTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATc    <  1:714040/49‑1 (MQ=255)

                                      gTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATc    >  1:1041499/1‑49 (MQ=255)

                                         tGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCtt  >  1:581549/1‑48 (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: `‑`