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

CGAATGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTCTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCTTTA  >  NC_000913/3477104‑3477195

cGAATGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCtttt                                             <  1:2741332/49‑1 (MQ=255)

cGAATGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCtttt                                             <  1:667149/49‑1 (MQ=255)

cGAATGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTa                                            >  1:1559460/1‑50 (MQ=255)

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

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

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

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

      cTTGTAGAGTCTGTTCGCTCTCTTGGTCGGAGAGTTTGCTTTTATCAGc                                       >  1:1865286/1‑49 (MQ=38)

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

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

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

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

          tAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATg                                  >  1:2565230/1‑50 (MQ=255)

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

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

            gagTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCa                                <  1:2537217/50‑1 (MQ=255)

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

               tCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATcc                             <  1:1568402/50‑1 (MQ=255)

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

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

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

                     cGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAac                        <  1:1642838/49‑1 (MQ=255)

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

                          tctTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGc                  >  1:847430/1‑50 (MQ=255)

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

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

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

                            ttGGTCGGAGAGGTTGCTTTTATCAGCAAATGCATCCTGAACACCa                    <  1:1823959/46‑1 (MQ=39)

                            atGGTCGGAGGGTTTGCTTGTATCAGCAAATGCATCCTGAACAcc                     <  1:974476/44‑1 (MQ=37)

                              ggTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGAt               <  1:995820/49‑1 (MQ=255)

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

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

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

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

                                tCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAg            >  1:1989707/1‑50 (MQ=255)

                                tCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAg            >  1:1055481/1‑50 (MQ=255)

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

                                  ggAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCt          >  1:159276/1‑50 (MQ=255)

                                  ggAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCt          >  1:1335990/1‑50 (MQ=255)

                                  ggAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCt          >  1:1156582/1‑50 (MQ=255)

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

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

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

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

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

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

CGAATGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTCTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCTTTA  >  NC_000913/3477104‑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: `‑`