BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,513,003	A→G	D212D (GAT→GAC)	trpS ←	tryptophan‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,513,003	0	A	G	100.0%	129.7 / NA	36	D212D (GAT→GAC)	trpS	tryptophan‑‑tRNA ligase
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (17/19); total (17/19)

TCGGAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGATCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACTTG  >  NC_000913/3512956‑3513048

tcgGAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTc                                             <  1:1192020/50‑1 (MQ=255)

     gTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTc                                          >  1:1012344/1‑48 (MQ=255)

     gTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCca                                        <  1:1748549/50‑1 (MQ=255)

        aGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCcagca                                     <  1:324612/50‑1 (MQ=255)

           gACCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGc                                  >  1:1309397/1‑50 (MQ=255)

            aCCGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGcc                                 <  1:2536695/50‑1 (MQ=255)

             ccGCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCg                                <  1:1855609/50‑1 (MQ=255)

               gCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGa                               >  1:2304587/1‑49 (MQ=255)

               gCACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGAt                              <  1:2551297/50‑1 (MQ=255)

                cACGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGAt                              >  1:2475240/1‑49 (MQ=255)

                 aCGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATaa                            <  1:1179433/50‑1 (MQ=255)

                  cGTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAAc                           >  1:869098/1‑50 (MQ=255)

                   gTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACg                          >  1:317959/1‑50 (MQ=255)

                   gTTTGATTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACg                          <  1:2383419/50‑1 (MQ=255)

                        aTTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTAtt                     <  1:2242815/50‑1 (MQ=255)

                        aTTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTAt                      <  1:834717/49‑1 (MQ=255)

                        aTTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTAt                      <  1:656892/49‑1 (MQ=255)

                        aTTTTCTTCACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTAt                      <  1:999062/49‑1 (MQ=255)

                            tcttcACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCg                  <  1:236415/49‑1 (MQ=255)

                             cttcACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGAt                >  1:2003743/1‑50 (MQ=255)

                              ttcACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGAt                >  1:1589097/1‑49 (MQ=255)

                               tcACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTa              >  1:329773/1‑50 (MQ=255)

                                cACTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGAtt               <  1:1363482/48‑1 (MQ=255)

                                 aCTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATc            >  1:899368/1‑50 (MQ=255)

                                 aCTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATc            >  1:276011/1‑50 (MQ=255)

                                 aCTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATc            <  1:1635620/50‑1 (MQ=255)

                                  cTACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCg           >  1:263259/1‑50 (MQ=255)

                                   tACCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGt          >  1:1987110/1‑50 (MQ=255)

                                    aCCGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGt          >  1:1814236/1‑49 (MQ=255)

                                     ccGATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTc         <  1:1059905/49‑1 (MQ=255)

                                       gATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGa      >  1:111789/1‑50 (MQ=255)

                                       gATTTCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGa      <  1:2010857/50‑1 (MQ=255)

                                           tCGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACTTg  >  1:2187213/1‑50 (MQ=255)

                                            cGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACTTg  <  1:264885/49‑1 (MQ=255)

                                            cGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACTTg  >  1:687424/1‑49 (MQ=255)

                                            cGGGTCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACTTg  <  1:1880529/49‑1 (MQ=255)

TCGGAGTCAGTGACCGCACGTTTGATTTTCTTCACTACCGATTTCGGATCTTCCAGCAGGCCGATAACGTTATTGCGATTATCGTCAGACTTG  >  NC_000913/3512956‑3513048

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