BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	4,563,589	T→A	100%	K207M (AAG→ATG)	iadA ←	isoaspartyl dipeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	4,563,589	0	T	A	100.0%	135.3 / NA	44	K207M (AAG→ATG)	iadA	isoaspartyl dipeptidase
Reads supporting (aligned to +/- strand): ref base T (0/0); new base A (35/9); total (35/9)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTT  >  W3110S.gb/4563522‑4563628

ggTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCattt                                      <  1:1609654/71‑1 (MQ=255)

ggTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCattt                                      <  1:2765269/71‑1 (MQ=255)

    ggCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTAc                                   >  1:238665/1‑70 (MQ=255)

    ggCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACt                                  >  1:961593/1‑71 (MQ=255)

    ggCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACt                                  >  1:821571/1‑71 (MQ=255)

    ggCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACt                                  >  1:726865/1‑71 (MQ=255)

    ggCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACt                                  >  1:1603999/1‑71 (MQ=255)

    ggCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACt                                  >  1:189933/1‑71 (MQ=255)

                     cGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGGTGTAACGCCATTTTACTGTCGCCCATGTGGAaca                 >  1:2687999/1‑71 (MQ=255)

                     cGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACt                                  >  1:400383/1‑54 (MQ=255)

                     cGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAgca                 >  1:502714/1‑71 (MQ=255)

                     cGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAaca                 >  1:1665215/1‑71 (MQ=255)

                     cGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAaca                 >  1:2396693/1‑71 (MQ=255)

                     cGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAaca                 >  1:2397402/1‑71 (MQ=255)

                     cGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAaca                 >  1:2975299/1‑71 (MQ=255)

                     cGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAaca                 >  1:2513577/1‑71 (MQ=255)

                     cGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAaca                 >  1:2544498/1‑71 (MQ=255)

                     cGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAaca                 >  1:131382/1‑71 (MQ=255)

                     cGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAaca                 >  1:84422/1‑71 (MQ=255)

                     cGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAaca                 >  1:2658409/1‑71 (MQ=255)

                     cGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTCTCGCCCATGTGGAaca                 >  1:2157853/1‑71 (MQ=255)

                         acaTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGt             <  1:3292481/71‑1 (MQ=255)

                         acaTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGt             <  1:2946047/71‑1 (MQ=255)

                         acaTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGt             <  1:992697/71‑1 (MQ=255)

                         acaTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGt             <  1:3317893/71‑1 (MQ=255)

                         acaTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGt             <  1:1042704/71‑1 (MQ=255)

                         acaTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGt             <  1:2547277/71‑1 (MQ=255)

                         acaTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGt             <  1:2053510/71‑1 (MQ=255)

                         acaTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGt             <  1:2606148/71‑1 (MQ=255)

                           aTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGt             <  1:3307913/69‑1 (MQ=255)

                              gCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGc        >  1:219227/1‑71 (MQ=255)

                                    ttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCAtgtg                      >  1:1071429/1‑51 (MQ=255)

                                    ttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGtt  >  1:1051243/1‑71 (MQ=255)

                                    ttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGtt  >  1:1260147/1‑71 (MQ=255)

                                    ttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGtt  >  1:1303243/1‑71 (MQ=255)

                                    ttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGtt  >  1:537839/1‑71 (MQ=255)

                                    ttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGtt  >  1:1828975/1‑71 (MQ=255)

                                    ttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGtt  >  1:1832775/1‑71 (MQ=255)

                                    ttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGtt  >  1:1954938/1‑71 (MQ=255)

                                    ttCCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGtt  >  1:2888618/1‑71 (MQ=255)

                                                   aaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAgg    >  1:466302/1‑54 (MQ=255)

                                                   aaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGt   >  1:1624375/1‑55 (MQ=255)

                                                   aaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGt   >  1:1602821/1‑55 (MQ=255)

                                                   aaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGt   >  1:159927/1‑55 (MQ=255)

                                                   aaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGt   >  1:3272674/1‑55 (MQ=255)

                                                   aaTAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGt   >  1:1911188/1‑55 (MQ=255)

GGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTT  >  W3110S.gb/4563522‑4563628

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