BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,969,173	A→G	19.7%	V89V (GTA→GTG)	yggG →	predicted peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,969,173	0	A	G	19.7%	204.1 / 36.3	136	V89V (GTA→GTG)	yggG	predicted peptidase
Reads supporting (aligned to +/- strand): ref base A (41/68); new base G (10/17); total (51/85)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.74e-01

Only 100 of 181 total aligned reads displayed.

AGCGAATACGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGTATATGGCGAAGGATGTGAACGCCTTTGCAATGGCTAACGGCTG  >  minE/1969105‑1969229

aGCGAATACGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTa                                                          >  1:2126462/1‑69 (MQ=255)

aGCGAATACGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTa                                                          >  1:1497716/1‑69 (MQ=255)

aGCGAATACGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTa                                                          >  1:2544153/1‑69 (MQ=255)

aGCGAATACGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTa                                                          >  1:1491754/1‑69 (MQ=255)

aGCGAATACGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTa                                                          >  1:1754266/1‑69 (MQ=255)

aGCGAATACGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTa                                                          >  1:527392/1‑69 (MQ=255)

aGCGAATACGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTa                                                          >  1:2981901/1‑69 (MQ=255)

aGCGAATACGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTa                                                          >  1:1781071/1‑69 (MQ=255)

 gCGAATACGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTaa                                                         >  1:3396234/1‑69 (MQ=255)

 gCGAATACGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTaa                                                         >  1:1050378/1‑69 (MQ=255)

 gCGAATACGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTaa                                                         >  1:3360512/1‑69 (MQ=255)

 gCGAATACGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTaa                                                         >  1:1774398/1‑69 (MQ=255)

 gCGAATACGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTaa                                                         >  1:1787644/1‑69 (MQ=255)

 gCGAATACGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTaa                                                         >  1:2839445/1‑69 (MQ=255)

 gCGAATACGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTaa                                                         >  1:2602849/1‑69 (MQ=255)

 gCGAATACGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTaa                                                         >  1:1974644/1‑69 (MQ=255)

 gCGAATACGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTaa                                                         >  1:2524347/1‑69 (MQ=255)

 gCGAATACGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTaa                                                         >  1:211353/1‑69 (MQ=255)

 gCGAATACGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTaa                                                         >  1:216664/1‑69 (MQ=255)

 gCGAATACGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTaa                                                         >  1:342766/1‑69 (MQ=255)

 gCGAATACGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTaa                                                         >  1:1097480/1‑69 (MQ=255)

 gCGAATACGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTaa                                                         >  1:663736/1‑69 (MQ=255)

 gCGAATACGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTaa                                                         >  1:1456848/1‑69 (MQ=255)

 gCGAATACGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTaa                                                         >  1:1427231/1‑69 (MQ=255)

 gCGAATACGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTaa                                                         >  1:1423008/1‑69 (MQ=255)

 gCGAATACGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTaa                                                         >  1:474498/1‑69 (MQ=255)

 gCGAATACGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTaa                                                         >  1:58896/1‑69 (MQ=255)

 gCGAATACGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTa                                                          >  1:289125/1‑68 (MQ=255)

        cGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACaa                                                  >  1:2230999/1‑69 (MQ=255)

        cGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACaa                                                  >  1:2479017/1‑69 (MQ=255)

        cGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACaa                                                  >  1:765551/1‑69 (MQ=255)

        cGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACaa                                                  >  1:1207571/1‑69 (MQ=255)

        cGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACaa                                                  >  1:2947207/1‑69 (MQ=255)

            aaaCGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAgtg                                              >  1:3017790/1‑69 (MQ=255)

            aaaCGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAgtg                                              >  1:1181644/1‑69 (MQ=255)

            aaaCGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAgt                                               >  1:227390/1‑68 (MQ=255)

               cGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAGAGTGta                                            <  1:442162/68‑1 (MQ=255)

               cGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTTta                                            <  1:256857/68‑1 (MQ=255)

               cGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGta                                            <  1:490834/68‑1 (MQ=255)

               cGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGta                                            <  1:3067074/68‑1 (MQ=255)

               cGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGta                                            <  1:3285746/68‑1 (MQ=255)

               cGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGta                                            <  1:377515/68‑1 (MQ=255)

               cGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGta                                            <  1:799215/68‑1 (MQ=255)

               cGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGta                                            <  1:2425290/68‑1 (MQ=255)

               cGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGta                                            <  1:3418199/68‑1 (MQ=255)

               cGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGta                                            <  1:1147433/68‑1 (MQ=255)

               cGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGta                                            <  1:1861564/68‑1 (MQ=255)

               cGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGta                                            <  1:1227148/68‑1 (MQ=255)

               cGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGta                                            <  1:1373361/68‑1 (MQ=255)

               cGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGta                                            <  1:1817856/68‑1 (MQ=255)

               cGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGta                                            <  1:1940331/68‑1 (MQ=255)

               cGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGta                                            <  1:1812657/68‑1 (MQ=255)

               cGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGta                                            <  1:1713977/68‑1 (MQ=255)

                     aCAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGta                                            >  1:116514/1‑62 (MQ=255)

                          tATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGTATATGGCGAAGGa                                <  1:1080538/69‑1 (MQ=255)

                          tATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGTATATGGCGAAGGa                                <  1:244996/69‑1 (MQ=255)

                           aTTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGTATATGGCGAAgg                                 >  1:1581051/1‑67 (MQ=255)

                           aTTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGTATATGGCGAAgg                                 >  1:3121703/1‑67 (MQ=255)

                           aTTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGTATATGGCGAAgg                                 >  1:3139599/1‑67 (MQ=255)

                           aTTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGTATATGGCGAAgg                                 >  1:2955102/1‑67 (MQ=255)

                           aTTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGTATATGGCGAAgg                                 >  1:3173949/1‑67 (MQ=255)

                           aTTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGTATATGGCGAAgg                                 >  1:1511545/1‑67 (MQ=255)

                           aTTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGTATATGGCGAAgg                                 >  1:1554696/1‑67 (MQ=255)

                           aTTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGTATATGGCGAAgg                                 >  1:2214732/1‑67 (MQ=255)

                              gCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGTATATGGCGAAGGAtgtt                            <  1:2377454/69‑2 (MQ=255)

                              gCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGTATATGGCGAAGGAtgtg                            <  1:1570734/69‑1 (MQ=255)

                               ccAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGTATATGGCGAAGGAtgtg                            <  1:722936/68‑1 (MQ=255)

                                   tGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAgtg                                              <  1:1928685/46‑1 (MQ=255)

                                   tGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAgtg                                              <  1:2597856/46‑1 (MQ=255)

                                   tGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAgtg                                              <  1:3035925/46‑1 (MQ=255)

                                   tGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAgtg                                              <  1:888070/46‑1 (MQ=255)

                                   tGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAgtg                                              <  1:2413664/46‑1 (MQ=255)

                                   tGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAgtg                                              <  1:2405285/46‑1 (MQ=255)

                                   tGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGTATATGGCGAAGGATGTGAACGc                       <  1:634806/69‑1 (MQ=255)

                                   tGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGTATATGGCGAAGGATGTGAACGc                       <  1:3228983/69‑1 (MQ=255)

                                   tGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGTATATGGCGAAGGATGTGAACGc                       <  1:1665909/69‑1 (MQ=255)

                                    gcgcTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGTATATGGCGAAGGATGTGAACGc                       <  1:91269/68‑1 (MQ=255)

                                           gCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGTATATGGCGAAGGATGTGAACGc                       <  1:359760/61‑1 (MQ=255)

                                                    tCAACGGTCAGCCGGTAAATTACAAAGTGTATATGGCGAAGGATGTGAACGCCTTTGCAATGGCTAACg      >  1:990657/1‑69 (MQ=255)

                                                     cAACGGTCAGCCGGTGAATTACAATGTGTATATGGCGAAGGATGTGAACGCCTTTGc                 <  1:2035212/57‑1 (MQ=255)

                                                     cAACGGTCAGCCGGTGAATTACAAAGTGTATATGGCGAAGGATGTGAACGCCTTTGc                 <  1:2358315/57‑1 (MQ=255)

                                                     cAACGGTCAGCCGGTGAATTACAAAGTGTATATGGCGAAGGATGTGAACGCCTTTGc                 <  1:968697/57‑1 (MQ=255)

                                                     cAACGGTCAGCCGGTGAATTACAAAGTGTATATGGCGAAGGATGTGAACGCCTTTGc                 <  1:3292898/57‑1 (MQ=255)

                                                     cAACGGTCAGCCGGTGAATTACAAAGTGTATATGGCGAAGGATGTGAACGCCTTTGc                 <  1:2641448/57‑1 (MQ=255)

                                                     cAACGGTCAGCCGGTGAATTACAAAGTGTATATGGCGAAGGATGTGAACGCCTTTGc                 <  1:1368667/57‑1 (MQ=255)

                                                     cAACGGTCAGCCGGTGAATTACAAAGTGTATATGGCGAAGGATGTGAACGCCTTTGc                 <  1:1616815/57‑1 (MQ=255)

                                                     cAACGGTCAGCCGGTGAATTACAAAGTGTATATGGCGAAGGATGTGAACGCCTTTGc                 <  1:1896953/57‑1 (MQ=255)

                                                      aaCGGTCAGCCGGTGATTTACAAAGTGTATATGGCg                                     >  1:63992/1‑36 (MQ=38)

                                                      aaCGGTCAGCCGGTGAATTACAAAGTGTATATGGCGAAGGa                                >  1:2063536/1‑41 (MQ=255)

                                                      aaCGGTCAGCCGGTGAATTACAAAGTGTATATGGCGAAGGATGTGAACGCCTTTGCaa               >  1:2731695/1‑58 (MQ=255)

                                                      aaCGGTCAGCCGGTGAATTACAAAGTGTATATGGCGAAGGATGTGAACGCCTTTGCaa               >  1:704921/1‑58 (MQ=255)

                                                      aaCGGTCAGCCGGTGAATTACAAAGTGTATATGGCGAAGGATGTGAACGCCTTTGCaa               >  1:464359/1‑58 (MQ=255)

                                                      aaCGGTCAGCCGGTGAATTACAAAGTGTATATGGCGAAGGATGTGAACGCCTTTGCaa               >  1:784011/1‑58 (MQ=255)

                                                      aaCGGTCAGCCGGTGAATTACAAAGTGTATATGGCGAAGGATGTGAACGCCTTTGCaa               >  1:1030291/1‑58 (MQ=255)

                                                      aaCGGTCAGCCGGTGAATTACAAAGTGTATATGGCGAAGGATGTGAACGCCTTTGCaa               >  1:3350527/1‑58 (MQ=255)

                                                        cGGTCAGCCGGTAAATTACAAAGTGTATATGGCGAAGGATGTGAACGCCTTTGCAATGGCTAACGGCTg  >  1:1325101/1‑69 (MQ=255)

                                                        cGGTCAGCCGGTAAATTACAAAGTGTATATGGCGAAGGATGTGAACGCCTTTGCAATGGCTAACGGCTg  >  1:1860512/1‑69 (MQ=255)

                                                        cGGTCAGCCGGTAAATTACAAAGTGTATATGGCGAAGGATGTGAACGCCTTTGCAATGGCTAACGGCTg  >  1:1978788/1‑69 (MQ=255)

                                                        cGGTCAGCCGGTAAATTACAAAGTGTATATGGCGAAGGATGTGAACGCCTTTGCAATGGCTAACGGCTg  >  1:2514772/1‑69 (MQ=255)

                                                                cGGTAAATTACAAAGTGTATATGGCGAAGGAtgtg                            <  1:1637590/35‑1 (MQ=255)

AGCGAATACGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGTATATGGCGAAGGATGTGAACGCCTTTGCAATGGCTAACGGCTG  >  minE/1969105‑1969229

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