BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,052,352	C→A	W389C (TGG→TGT)	ubiH ←	2‑octaprenyl‑6‑methoxyphenol hydroxylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,052,352	0	C	A	100.0%	151.6 / NA	39	W389C (TGG→TGT)	ubiH	2‑octaprenyl‑6‑methoxyphenol hydroxylase
Reads supporting (aligned to +/- strand): ref base C (0/0); new base A (27/12); total (27/12)

AACACTTTGCATTGTTTATTCCTTAAAACCGCCTTCAACGCGCCACCCAACCGAGGGTGCGCTGCGCCAGCACATCGCGTGCCGGGGTGAAT  >  NC_000913/3052306‑3052397

aaCACTTTGCATTGTTTATTCCTTAAAACCGCCTTCAACGCGCcacaca                                             <  1:1855169/49‑1 (MQ=255)

 acacTTTGCATTGTTTATTCCTTAAAACCGCCTTCAACGCGCCACACaa                                            <  1:230596/49‑1 (MQ=255)

 acacTTTGCATTGTTTATTCCTTAAAACCGCCTTCAACGCGCCACACaa                                            >  1:182034/1‑49 (MQ=255)

  cacTTTGCATTGTTTATTCCTTAAAACCGCCTTCAACGCGCCACACAAc                                           >  1:423685/1‑49 (MQ=255)

  cacTTTGCATTGTTTATTCCTTAAAACCGCCTTCAACGCGCCACACAAc                                           >  1:2641482/1‑49 (MQ=255)

  cacTTTGCATTGTTTATTCCTTAAAACCGCCTTCAACGCGCCACACAAc                                           >  1:1449487/1‑49 (MQ=255)

  cacTTTGCATTGTTTATTCCTTAAAACCGCCTTCAACGCGCCACACAAc                                           >  1:250548/1‑49 (MQ=255)

   acTTTGCATTGTTTATTCCTTAAAACCGCCTTCAACGCGCCACACAAcc                                          >  1:2841762/1‑49 (MQ=255)

    cTTTGCATTGTTTATTCCTTAAAACCGCCTTCAACGCGCCACACAAcc                                          >  1:2522438/1‑48 (MQ=255)

      ttGCATTGTTTATTCCTTAAAACCGCCTTCAACGCGCCACACAACCGAg                                       >  1:1053076/1‑49 (MQ=255)

      ttGCATTGTTTATTCCTTAAAACCGCCTTCAACGCGCCACACAACCGAg                                       >  1:690256/1‑49 (MQ=255)

       tGCATTGTTTATTCCTTAAAACCGCCTTCAACGCGCCACACAACCGAgg                                      >  1:2338722/1‑49 (MQ=255)

       tGCATTGTTTATTCCTTAAAACCGCCTTCAACGCGCCACACAACCGAgg                                      >  1:2584518/1‑49 (MQ=255)

       tGCATTGTTTATTCCTTAAAACCGCCTTCAACGCGCCACACAACCGAg                                       <  1:885496/48‑1 (MQ=255)

        gCATTGTTTATTCCTTAAAACCGCCTTCAACGCGCCACACAACCGAgg                                      <  1:758021/48‑1 (MQ=255)

           ttGTTTATTCCTTAAAACCGCCTTCAACGCGCCACACAACCGAGGGtgc                                  <  1:1693179/49‑1 (MQ=255)

             gTTTATTCCTTAAAACCGCCTTCAACGCGCCACACAACCGAGGGtgcgc                                >  1:697890/1‑49 (MQ=255)

              tttATTCCTTAAAACCGCCTTCAACGCGCCACACAACCGAGGGtgcgct                               <  1:307840/49‑1 (MQ=255)

              tttATTCCTTAAAACCGCCTTCAACGCGCCACACAACCGAGGGtgcgct                               <  1:2795012/49‑1 (MQ=255)

              tttATTCCTTAAAACCGCCTTCAACGCGCCACACAACCGAGGGtgcgct                               <  1:1414608/49‑1 (MQ=255)

              tttATTCCTTAAAACCGCCTTCAACGCGCCACACAACCGAGGGtgcgct                               >  1:2568913/1‑49 (MQ=255)

                       tAAAACCGCCTTCAACGCGCCACACAACCGAGGGTGCGCTGCGCCAGca                      >  1:1610616/1‑49 (MQ=255)

                        aaaaCCGCCTTCAACGCGCCACACAACCGAGGGTGCGCTGCGCCAGca                      >  1:1008012/1‑48 (MQ=255)

                         aaaCCGCCTTCAACGCGCCACACAACCGAGGGTGCGCTGCGCCAGcaca                    >  1:2137390/1‑49 (MQ=255)

                          aaCCGCCTTCAACGCGCCACACAACCGAGGGTGCGCTGCGCCAGCACAt                   >  1:2407935/1‑49 (MQ=255)

                          aaCCGCCTTCAACGCGCCACACAACCGAGGGTGCGCTGCGCCAGCACAt                   <  1:1190787/49‑1 (MQ=255)

                            ccGCCTTCAACGCGCCACACAACCGAGGGTGCGCTGCGCCAGCACATcg                 <  1:1657764/49‑1 (MQ=255)

                            ccGCCTTCAACGCGCCACACAACCGAGGGTGCGCTGCGCCAGCACATcg                 <  1:2553611/49‑1 (MQ=255)

                             cGCCTTCAACGCGCCACACAACCGAGGGTGCGCTGCGCCAGCACATcgc                >  1:2686857/1‑49 (MQ=255)

                             cGCCTTCAACGCGCCACACAACCGAGGGTGCGCTGCGCCAGCACATcgc                >  1:2539184/1‑49 (MQ=255)

                             cGCCTTCAACGCGCCACACAACCGAGGGTGCGCTGCGCCAGCACATcgc                <  1:360008/49‑1 (MQ=255)

                                cTTCAACGCGCCACACAACCGAGGGTGCGCTGCGCCAGCACATCGCGTg             >  1:1526924/1‑49 (MQ=255)

                                 ttCAACGCGCCACACAACCGAGGGTGCGCTGCGCCAGCACATCGCGTGc            >  1:2583202/1‑49 (MQ=255)

                                  tCAACGCGCCACACAACCGAGGGTGCGCTGCGCCAGCACATCGCGTGcc           >  1:850888/1‑49 (MQ=255)

                                   cAACGCGCCACACAACCGAGGGTGCGCTGCGCCAGCACATCGCGTGCCg          >  1:1475915/1‑49 (MQ=255)

                                   cAACGCGCCACACAACCGAGGGTGCGCTGCGCCAGCACATCGCGTGCCg          >  1:2065016/1‑49 (MQ=255)

                                     aCGCGCCACACAACCGAGGGTGCGCTGCGCCAGCACATCGCGTGCCggg        >  1:2692197/1‑49 (MQ=255)

                                     aCGCGCCACACAACCGAGGGTGCGCTGCGCCAGCACATCGCGTGCCggg        >  1:1593226/1‑49 (MQ=255)

                                     aCGCGCCACACAACCGAGGGTGCGCTGCGCCAGCACATCGCGTGCCggg        <  1:183069/49‑1 (MQ=255)

                                         gcCACACAACCGAGGGTGCGCTGCGCCAGCACATCGCGTGCCGGGGTg     >  1:1315421/1‑48 (MQ=255)

                                           cacacaACCGAGGGTGCGCTGCGCCAGCACATCGCGTGCCGGGGTGaat  <  1:2115248/49‑1 (MQ=255)

                                           cacacaACCGAGGGTGCGCTGCGCCAGCACATCGCGTGCCGGGGTGaa   >  1:317427/1‑48 (MQ=255)

AACACTTTGCATTGTTTATTCCTTAAAACCGCCTTCAACGCGCCACCCAACCGAGGGTGCGCTGCGCCAGCACATCGCGTGCCGGGGTGAAT  >  NC_000913/3052306‑3052397

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