BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	257,794	T→C	intergenic (+23/‑35)	frsA → / → crl	fermentation‑respiration switch protein/RNA polymerase holoenzyme assembly factor Crl

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	257,794	0	T	C	100.0%	164.8 / NA	41	intergenic (+23/‑35)	frsA/crl	fermentation‑respiration switch protein/RNA polymerase holoenzyme assembly factor Crl
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (23/18); total (23/18)

TGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACC  >  NC_000913/257748‑257839

tGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTg                                             >  1:1482806/1‑49 (MQ=255)

  gATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTg                                             <  1:1276377/47‑1 (MQ=255)

  gATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTg                                             <  1:527643/47‑1 (MQ=255)

  gATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGt                                           <  1:1382265/49‑1 (MQ=255)

   aTCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGt                                           >  1:234818/1‑48 (MQ=255)

    tCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTa                                          >  1:641127/1‑48 (MQ=255)

     cGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTaa                                         >  1:2192092/1‑48 (MQ=255)

      gAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTaaaa                                       >  1:2613326/1‑49 (MQ=255)

      gAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTaaa                                        <  1:877964/48‑1 (MQ=255)

      gAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTaaa                                        <  1:536174/48‑1 (MQ=255)

      gAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTaaa                                        <  1:1175609/48‑1 (MQ=255)

       aaaaaCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTaaaa                                       >  1:2298886/1‑48 (MQ=255)

         aaaCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAg                                    <  1:499501/49‑1 (MQ=255)

         aaaCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAg                                    <  1:76005/49‑1 (MQ=255)

         aaaCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAg                                    >  1:874156/1‑49 (MQ=255)

          aaCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGt                                   >  1:176912/1‑49 (MQ=255)

          aaCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGt                                   <  1:1697293/49‑1 (MQ=255)

           aCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGtt                                  >  1:2605663/1‑49 (MQ=255)

             gCTTGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGc                                <  1:354887/49‑1 (MQ=255)

               ttGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCAt                              >  1:1987557/1‑49 (MQ=255)

               ttGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCAt                              >  1:2421603/1‑49 (MQ=255)

               ttGTGTTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCAt                              >  1:2364152/1‑49 (MQ=255)

                  tgtTAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATca                            <  1:1281336/48‑1 (MQ=255)

                    ttAAAAATTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCacaa                         >  1:793241/1‑49 (MQ=255)

                      aaaaaTTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCacaac                        <  1:1915923/48‑1 (MQ=255)

                        aaaTTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAg                      >  1:401969/1‑48 (MQ=255)

                        aaaTTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAg                      >  1:403000/1‑48 (MQ=255)

                         aaTTTGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAg                      >  1:2121803/1‑47 (MQ=255)

                             tGCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGgaga                  >  1:1839087/1‑47 (MQ=255)

                              gCTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATa                <  1:1217226/48‑1 (MQ=255)

                               cTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAg               <  1:810073/48‑1 (MQ=255)

                               cTAAATTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAg               <  1:255874/48‑1 (MQ=255)

                                 aaaTTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCaa            >  1:2833443/1‑49 (MQ=255)

                                   aTTTTGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAAt           >  1:1889310/1‑48 (MQ=255)

                                      ttGCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGAc        <  1:1258415/48‑1 (MQ=255)

                                        gCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGt      >  1:1294438/1‑48 (MQ=255)

                                        gCCAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGt      <  1:203644/48‑1 (MQ=255)

                                         ccAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTa    >  1:1905958/1‑49 (MQ=255)

                                          cAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTAc   <  1:465147/49‑1 (MQ=255)

                                          cAATCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTAc   >  1:1742125/1‑49 (MQ=255)

                                           aaTCTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTAcc  >  1:1381152/1‑49 (MQ=255)

TGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACC  >  NC_000913/257748‑257839

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: `‑`