BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	281,478	T→C	25.8%	L179L (TTA→TTG)	ampG ←	muropeptide transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	281,478	0	T	C	25.8%	58.9 / 18.5	54	L179L (TTA→TTG)	ampG	muropeptide transporter
Reads supporting (aligned to +/- strand): ref base T (19/21); new base C (12/2); total (31/23)
Fisher's exact test for biased strand distribution p-value = 1.46e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GGAATGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAACCAGTACATGCCCTGCCAGCCCAGCCATTTATCTGCCAGCCACAGGGCTAGGCCGCC  >  minE/281410‑281537

ggAATGGTGTCGGTTGGTTCTGGTGCGAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCAtt                                                             >  1:2327439/1‑69 (MQ=255)

ggAATGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACGGTGCCGCCAtt                                                             >  1:3076646/1‑69 (MQ=255)

ggAATGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCAtt                                                             >  1:3193327/1‑69 (MQ=255)

ggAATGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCAtt                                                             >  1:2966599/1‑69 (MQ=255)

ggAATGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCAtt                                                             >  1:2881519/1‑69 (MQ=255)

ggAATGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCAtt                                                             >  1:478768/1‑69 (MQ=255)

ggAATGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCAtt                                                             >  1:3311158/1‑69 (MQ=255)

ggAATGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCAtt                                                             >  1:2600785/1‑69 (MQ=255)

ggAATGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCAtt                                                             >  1:2702835/1‑69 (MQ=255)

   aTGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTa                                                            >  1:2052063/1‑67 (MQ=255)

   aTGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTa                                                            >  1:2024702/1‑67 (MQ=255)

   aTGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTa                                                            >  1:1958790/1‑67 (MQ=255)

   aTGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTa                                                            >  1:2187735/1‑67 (MQ=255)

   aTGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTa                                                            >  1:2551033/1‑67 (MQ=255)

   aTGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTa                                                            >  1:2769641/1‑67 (MQ=255)

   aTGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTa                                                            >  1:881674/1‑67 (MQ=255)

   aTGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTa                                                            >  1:1356333/1‑67 (MQ=255)

   aTGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTa                                                            >  1:1224604/1‑67 (MQ=255)

   aTGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTa                                                            >  1:3041941/1‑67 (MQ=255)

   aTGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTa                                                            >  1:3106592/1‑67 (MQ=255)

      gtgtCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAACCAg                                                       <  1:3254765/69‑1 (MQ=255)

      gtgtCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAACCAg                                                       <  1:2286368/69‑1 (MQ=255)

      gtgtCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAACCAg                                                       <  1:1792466/69‑1 (MQ=255)

        gtCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAACCAGTa                                                     <  1:1875300/69‑1 (MQ=255)

        gtCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAACCAGTa                                                     <  1:882946/69‑1 (MQ=255)

        gtCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAACCAGTa                                                     <  1:115679/69‑1 (MQ=255)

        gtCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAACCAGTa                                                     <  1:3050349/69‑1 (MQ=255)

        gtCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAACCAGTa                                                     <  1:1506766/69‑1 (MQ=255)

        gtCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAACCAGTa                                                     <  1:412492/69‑1 (MQ=255)

        gtCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAACCAGTa                                                     <  1:1687028/69‑1 (MQ=255)

        gtCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAACCAGTa                                                     <  1:186186/69‑1 (MQ=255)

           ggttggttCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAAc                                                          <  1:181923/61‑1 (MQ=255)

           ggttggttCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAAc                                                          <  1:582603/61‑1 (MQ=255)

           ggttggttCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAAc                                                          <  1:3409504/61‑1 (MQ=255)

           ggttggttCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAAc                                                          <  1:1098489/61‑1 (MQ=255)

           ggttggttCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAAc                                                          <  1:1266137/61‑1 (MQ=255)

           ggttggttCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAACCAGTACAt                                                  >  1:2261835/1‑69 (MQ=255)

                   cTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGCGCCGCCATCAACCAGTACATGCCCTGcc                                          >  1:249397/1‑69 (MQ=255)

                   cTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGCGCCGCCATCAACCAGTACATGCCCTGcc                                          >  1:814367/1‑69 (MQ=255)

                   cTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGCGCCGCCATCAACCAGTACATGCCCTGcc                                          >  1:1084346/1‑69 (MQ=255)

                   cTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGCGCCGCCATCAACCAGTACATGCCCTGcc                                          >  1:3137295/1‑69 (MQ=255)

                   cTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGCGCCGCCATCAACCAGTACATGCCCTGcc                                          >  1:3105012/1‑69 (MQ=255)

                   cTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGCGCCGCCATCAACCAGTACATGCCCTGcc                                          >  1:1170801/1‑69 (MQ=255)

                   cTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGCGCCGCCATCAACCAGTACATGCCCTGcc                                          >  1:1415965/1‑69 (MQ=255)

                   cTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGCGCCGCCATCAACCAGTACATGCCCTGcc                                          >  1:335578/1‑69 (MQ=255)

                   cTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGCGCCGCCATCAACCAGTACATGCCCTGcc                                          >  1:1515012/1‑69 (MQ=255)

                      gTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGCGCCGCCATCAACCAGTACATGCCCTGccag                                        >  1:753285/1‑68 (MQ=255)

                      gTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGCGCCGCCATCAACCAGTACATGCCCTGccag                                        >  1:2824454/1‑68 (MQ=255)

                      gTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGCGCCGCCATCAACCAGTACATGCCCTGccag                                        >  1:1654413/1‑68 (MQ=255)

                         caagcaaCGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAACCAGTACATGCCCTGccagccca                                    >  1:1045652/1‑69 (MQ=255)

                         caagcaaCGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAACCAGTACATGCCCTGccagccca                                    >  1:361394/1‑69 (MQ=255)

                         caagcaaCGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAACCAGTACATGCCCTGccagccca                                    >  1:2491317/1‑69 (MQ=255)

                         caagcaaCGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAACCAGTACATGCCCTGccagccca                                    >  1:2108664/1‑69 (MQ=255)

                                 gTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAACCAGTACATGCCCTGCCAGCCCAGCCATTTa                            <  1:3219127/69‑1 (MQ=255)

                                 gTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAACCAGTACATGCCCTGCCAGCCCAGCCATTTa                            <  1:1488934/69‑1 (MQ=255)

                                 gTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAACCAGTACATGCCCTGCCAGCCCAGCCATTTa                            <  1:1013808/69‑1 (MQ=255)

                                          ataCAGGGGATCAACAGTGCCGCCATTAACCAGTACATGCCCTGCCAGCCCAGCCATTTATCTgccagc                   <  1:2824793/69‑1 (MQ=255)

                                          ataCAGGGGATCAACAGTGCCGCCATTAACCAGTACATGCCCTGCCAGCCCAGCCATTTATCTgccagc                   >  1:2615649/1‑69 (MQ=255)

                                            aCAGGGGATCAACAGTGCCGCCATTAACCAGTACATGCCCTGCCAGCCCAGCCATTTATCTgccagc                   >  1:2838540/1‑67 (MQ=255)

                                             cAGGGGATCAACAGTGTCGCCATTAACCAGTACATGCCCTGCCAGCCCAGCCATTTATCTGCCAGCcac                <  1:2651231/69‑1 (MQ=255)

                                                     cAACAGCGCCGCCATCAACCAGTACATGCCCTGCCAGCCCAGCCATTTATCTGCCAGCCACAGGGCTaa        <  1:2616209/69‑2 (MQ=255)

                                                      aaCAGCGCCGCCATCAACCAGTACATGCCCTGCCAGCCCAGCCATTTATCTGCCAGCCACAGGGCTaa        <  1:2592896/68‑2 (MQ=255)

                                                            gccgccATTAACCAGTACATGCCCTGCCAGCCCAGCCATTTATCTGCCAGCCACAGGGCTAGgccgcc  >  1:229712/1‑68 (MQ=255)

GGAATGGTGTCGGTTGGTTCTGGTGCAAGCAACGTCGCGATAATACAGGGGATCAACAGTGCCGCCATTAACCAGTACATGCCCTGCCAGCCCAGCCATTTATCTGCCAGCCACAGGGCTAGGCCGCC  >  minE/281410‑281537

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