BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,907,399	C→T	E13K (GAG→AAG)	cspC ←	stress protein, member of the CspA family

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,907,399	0	C	T	100.0%	188.5 / NA	47	E13K (GAG→AAG)	cspC	stress protein, member of the CspA family
Reads supporting (aligned to +/- strand): ref base C (0/0); new base T (18/29); total (18/29)

TGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATT  >  NC_000913/1907357‑1907445

tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGtt                                             <  1:2421208/46‑1 (MQ=255)

tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGaa                                          <  1:1589780/49‑1 (MQ=255)

tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGaa                                          <  1:2600135/49‑1 (MQ=255)

tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGa                                           <  1:806500/48‑1 (MQ=255)

 gctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAAc                                         >  1:2433922/1‑49 (MQ=255)

  ctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAAcc                                        <  1:2490604/49‑1 (MQ=255)

    gcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACt                                     >  1:896932/1‑50 (MQ=255)

     ccATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACtt                                    <  1:1622629/50‑1 (MQ=255)

      cATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTa                                   <  1:2721196/50‑1 (MQ=255)

       aTCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTaa                                  <  1:2507264/50‑1 (MQ=255)

       aTCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTaa                                  <  1:2571080/50‑1 (MQ=255)

          aGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCt                               <  1:2242907/50‑1 (MQ=255)

          aGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCt                               <  1:2475586/50‑1 (MQ=255)

          aGCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCt                               >  1:1859269/1‑50 (MQ=255)

           gCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTg                              <  1:2466294/50‑1 (MQ=255)

           gCCGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTg                              <  1:1293606/50‑1 (MQ=255)

            ccGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGa                             >  1:2200124/1‑50 (MQ=255)

             cGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAc                            >  1:2622117/1‑50 (MQ=255)

             cGGAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGAc                            >  1:806359/1‑50 (MQ=255)

               gAGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCt                          <  1:283372/50‑1 (MQ=255)

                aGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCtt                         <  1:1667677/50‑1 (MQ=255)

                aGTAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCtt                         >  1:2394336/1‑50 (MQ=255)

                  tAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCttt                        >  1:1166928/1‑49 (MQ=255)

                  tAATGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTa                       <  1:940269/50‑1 (MQ=255)

                   aaTGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTaa                      >  1:842811/1‑50 (MQ=255)

                    aTGAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTaa                      <  1:1021210/49‑1 (MQ=255)

                      gAAGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCt                   >  1:1379942/1‑50 (MQ=255)

                       aaGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCtt                  <  1:2284190/50‑1 (MQ=255)

                       aaGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCt                   <  1:2329492/49‑1 (MQ=255)

                       aaGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCt                   <  1:2627869/49‑1 (MQ=255)

                        aGCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCttt                 <  1:1976272/50‑1 (MQ=255)

                         gCCAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTg                >  1:1997006/1‑50 (MQ=255)

                           cAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGcc              <  1:2167363/50‑1 (MQ=255)

                           cAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGcc              <  1:681672/50‑1 (MQ=255)

                           cAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGc               >  1:2664337/1‑49 (MQ=255)

                           cAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGc               >  1:2741055/1‑49 (MQ=255)

                           cAAAACCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGc               >  1:1215691/1‑49 (MQ=255)

                             aaaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAt            <  1:897465/50‑1 (MQ=255)

                             aaaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAt            >  1:397559/1‑50 (MQ=255)

                              aaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAtt           <  1:651601/50‑1 (MQ=255)

                              aaCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAtt           >  1:818926/1‑50 (MQ=255)

                               aCCTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAttt          <  1:1924138/50‑1 (MQ=255)

                                 cTTTAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGa        <  1:1273142/50‑1 (MQ=255)

                                  tttAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGa        >  1:2360462/1‑49 (MQ=255)

                                    tAGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaa      <  1:2221360/49‑1 (MQ=255)

                                     aGACTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGaaaa     <  1:229266/49‑1 (MQ=255)

                                       aCTTGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAAtt  >  1:1202785/1‑50 (MQ=255)

TGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTTGAAAAATT  >  NC_000913/1907357‑1907445

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