BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,970,359	G→A	51.0%	L145L (CTG→TTG)	speB ←	agmatinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,970,359	0	G	A	51.0%	‑3.7 / 94.5	102	L145L (CTG→TTG)	speB	agmatinase
Reads supporting (aligned to +/- strand): ref base G (31/19); new base A (7/45); total (38/64)
Fisher's exact test for biased strand distribution p-value = 4.58e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

Only 100 of 105 total aligned reads displayed.

GAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCATAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAGCAGCGGCAGCGTAACAAAGTGGTCACC  >  minE/1970291‑1970422

gAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCATAGGTATCGGTGTGGGCGTCAAAGTGTACCAg                                                                 >  1:245730/1‑69 (MQ=255)

gAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCATAGGTATCGGTGTGGGCGTCAAAGTGTACCAg                                                                 >  1:124877/1‑69 (MQ=255)

      aGTGCCGTGGTCAAATTCACAACCGTTCGCATAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCAt                                                           <  1:2451185/69‑1 (MQ=255)

      aGTGCCGTGGTCAAATTCACAACCGTTCGCATAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCAt                                                           <  1:1142694/69‑1 (MQ=255)

      aGTGCCGTGGTCAAATTCACAACCGTTCGCATAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCAt                                                           <  1:552676/69‑1 (MQ=255)

         gCCGTGGTCAAATTCACAACCGTTAGCGTAGGTATCGGTGTGGGCATCAAAGTGTACCAACGCCAtttt                                                        >  1:1753014/1‑69 (MQ=255)

         gCCGTGGTCAAATTCACAACCGTTAGCGTAGGTATCGGTGTGGGCATCAAAGTGTACCAACGCCAtttt                                                        >  1:1112333/1‑69 (MQ=255)

         gCCGTGGTCAAATTCACAACCGTTAGCGTAGGTATCGGTGTGGGCATCAAAGTGTACCAACGCCAtttt                                                        >  1:817956/1‑69 (MQ=255)

         gCCGTGGTCAAATTCACAACCGTTAGCGTAGGTATCGGTGTGGGCATCAAAGTGTACCAACGCCAtttt                                                        >  1:2514388/1‑69 (MQ=255)

         gCCGTGGTCAAATTCACAACCGTTAGCGTAGGTATCGGTGTGGGCATCAAAGTGTACCAACGCCAtttt                                                        >  1:1289773/1‑69 (MQ=255)

         gCCGTGGTCAAATTCACAACCGTTAGCGTAGGTATCGGTGTGGGCATCAAAGTGTACCAACGCCAtttt                                                        >  1:2928915/1‑69 (MQ=255)

         gCCGTGGTCAAATTCACAACCGTTAGCGTAGGTATCGGTGTGGGCATCAAAGTGTACCAACGCCAtttt                                                        >  1:305247/1‑69 (MQ=255)

             tGGTCAAATTCACAACCGTTAGCGTAGGTATCGGTGTGGGCATCAAAGTGTACCAACGCCATTTTGCCg                                                    <  1:32810/69‑1 (MQ=255)

             tGGTCAAATTCACAACCGTTAGCGTAGGTATCGGTGTGGGCATCAAAGTGTACCAACGCCATTTTGCCg                                                    <  1:1849621/69‑1 (MQ=255)

             tGGTCAAATTCACAACCGTTAGCGTAGGTATCGGTGTGGGCATCAAAGTGTACCAACGCCATTTTGCCg                                                    <  1:236146/69‑1 (MQ=255)

             tGGTCAAATTCACAACCGTTAGCGTAGGTATCGGTGTGGGCATCAAAGTGTACCAACGCCATTTTGCCg                                                    <  1:2709940/69‑1 (MQ=255)

             tGGTCAAATTCACAACCGTTAGCGTAGGTATCGGTGTGGGCATCAAAGTGTACCAACGCCATTTTGCCg                                                    <  1:3153102/69‑1 (MQ=255)

             tGGTCAAATTCACAACCGTTAGCGTAGGTATCGGTGTGGGCATCAAAGTGTACCAACGCCATTTTGCCg                                                    <  1:1052721/69‑1 (MQ=255)

             tGGTCAAATTCACAACCGTTAGCGTAGGTATCGGTGTGGGCATCAAAGTGTACCAACGCCATTTTGCCg                                                    <  1:40637/69‑1 (MQ=255)

             tGGTCAAATTCACAACCGTTAGCGTAGGTATCGGTGTGGGCATCAAAGTGTACCAACGCCATTTTGCCg                                                    <  1:7738/69‑1 (MQ=255)

             tGGTCAAATTCACAACCGTTAGCGTAGGTATCGGTGTGGGCATCAAAGTGTACCAACGCCATTTTGCCg                                                    <  1:830361/69‑1 (MQ=255)

             tGGTCAAATTCACAACCGTTAGCGTAGGTATCGGTGTGGGCATCAAAGTGTACCAACGCCATTTTGCCg                                                    <  1:843088/69‑1 (MQ=255)

             tGGTCAAATTCACAACCGTTAGCGTAGGTATCGGTGTGGGCATCAAAGTGTACCAACGCCATTTTGCCg                                                    <  1:936784/69‑1 (MQ=255)

             tGGTCAAATTCACAACCGTTAGCGTAGGTATCGGTGTGGGCATCAAAGTGTACCAACGCCATTTTGCCg                                                    <  1:941592/69‑1 (MQ=255)

             tGGTCAAATTCACAACCGTTAGCGTAGGTATCGGTGTGGGCATCAAAGTGTACCAACGCCATTTTGCCg                                                    <  1:1058394/69‑1 (MQ=255)

             tGGTCAAATTCACAACCGTTAGCGTAGGTATCGGTGTGGGCATCAAAGTGTACCAACGCCATTTTGCCg                                                    <  1:11975/69‑1 (MQ=255)

             tGGTCAAATTCACAACCGTTAGCGTAGGTATCGGTGTGGGCATCAAAGTGTACCAACGCCATTTTGCCg                                                    <  1:1583200/69‑1 (MQ=255)

             tGGTCAAATTCACAACCGTTAGCGTAGGTATCGGTGTGGGCATCAAAGTGTACCAACGCCATATTGCCg                                                    <  1:441462/69‑1 (MQ=255)

             tGGGCAAATTCACAACCGTTAGCGTAGGTATCGGTGTGGGCATCAAAGTGTACCAACGCCATTTTGCCg                                                    <  1:2211694/69‑1 (MQ=255)

              ggTCAAATTCACAACCGTTAGCGTAGGTATCGGTGTGGGCATCAAAGTGTACCAACGCCATTTTGCCg                                                    <  1:1698364/68‑1 (MQ=255)

                     ttCACAACCGTTCGCATAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCtt                                            >  1:1837137/1‑69 (MQ=255)

                     ttCACAACCGTTCGCATAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCtt                                            >  1:2755800/1‑69 (MQ=255)

                       cacaACCGTTCGCATAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTc                                           >  1:1231747/1‑68 (MQ=255)

                       cacaACCGTTCGCATAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCg                                          >  1:1004493/1‑69 (MQ=255)

                       cacaACCGTTCGCATAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCg                                          >  1:1956779/1‑69 (MQ=255)

                       cacaACAGTTCGCATAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCg                                          >  1:438920/1‑69 (MQ=255)

                         caACCGTTCGCATAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGc                                         <  1:2892205/68‑1 (MQ=255)

                         caACCGTTCGCATAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGc                                         <  1:1270855/68‑1 (MQ=255)

                         caACCGTTCGCATAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGc                                         <  1:1413585/68‑1 (MQ=255)

                         caACCGTTCGCATAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGc                                         <  1:2237393/68‑1 (MQ=255)

                             cGTTCGCATAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAg                                    <  1:884316/69‑1 (MQ=255)

                                             ggTGTGTGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAgcagcggcag                    >  1:436223/1‑69 (MQ=255)

                                             ggTGTGGGCGTCAAAGTGTACCAGCTCCATTTTGCCGAAATGCTTCGCATGAGCACGCAgcagctgcag                    >  1:1519271/1‑69 (MQ=255)

                                             ggTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAgcagcggcag                    >  1:3248715/1‑69 (MQ=255)

                                             ggTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAgcagcggcag                    >  1:3387983/1‑69 (MQ=255)

                                             ggTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAgcagcggcag                    >  1:1870991/1‑69 (MQ=255)

                                             ggTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAgcagcggcag                    >  1:2878211/1‑69 (MQ=255)

                                             ggTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAgcagcggcag                    >  1:557729/1‑69 (MQ=255)

                                             ggTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAgcagcggcag                    >  1:1596995/1‑69 (MQ=255)

                                             ggTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAgcagcggcag                    >  1:85395/1‑69 (MQ=255)

                                             ggTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAgcagcggcag                    >  1:247791/1‑69 (MQ=255)

                                             ggTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAgcagcggcag                    >  1:2451133/1‑69 (MQ=255)

                                             ggTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAgcagcggcag                    >  1:1668823/1‑69 (MQ=255)

                                             ggTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAgcagcggcag                    >  1:2346351/1‑69 (MQ=255)

                                             ggTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAgcagcggcag                    >  1:2458651/1‑69 (MQ=255)

                                             ggTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAgcagcggcag                    >  1:2425417/1‑69 (MQ=255)

                                                gtgGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAgcagcggcagcg                  >  1:1057252/1‑68 (MQ=255)

                                                gtgGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAGCAGCGGCAGCGt                 >  1:1636957/1‑69 (MQ=255)

                                                gtgGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAGCAGCGGCAGCGt                 >  1:343168/1‑69 (MQ=255)

                                                gtgGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAGCAGCGGCAGCGt                 >  1:2343811/1‑69 (MQ=255)

                                                gtgGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAGCAGCGGCAGCGt                 >  1:1339582/1‑69 (MQ=255)

                                                gtgGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAGCAGCGGCAGCGt                 >  1:2426057/1‑69 (MQ=255)

                                                gtgGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAACAGCGGAAGCGt                 >  1:1187205/1‑69 (MQ=255)

                                                     cGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAGCAGCGGCAGCGTAACaa            <  1:2925833/69‑1 (MQ=255)

                                                     cGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAGCAGCGGCAGCGTAACaa            <  1:3003341/69‑1 (MQ=255)

                                                         aaaGTGTACCAACGCCATTTTGCCGAAATGCTTCGCATGAGCTCGCAGCAGCGGCAGCGTAACAAAg          <  1:2871122/67‑1 (MQ=255)

                                                         aaaGTGTACCAACGCCATTTTGCCGAAATGCTTCGCATGAGCGCGCAGCAGCGGCAGCGTAACAAAg          <  1:603138/67‑1 (MQ=255)

                                                         aaaGTGTACCAACGCCATTTTGCCGAAATGCTTCGCATGAGCGCGCAGCAGCGGCAGCGTAACAAAg          <  1:1842109/67‑1 (MQ=255)

                                                         aaaGTGTACCAACGCCATTTTGCCGAAATGCTTCGCATGAGCGCGCAGCAGCGGCAGCGTAACAAAg          <  1:469727/67‑1 (MQ=255)

                                                         aaaGTGTACCAACGCCATTTTGCCGAAATGCTTCGCATGAGCGCGCAGCAGCGGCAGCGTAACAAAg          <  1:707788/67‑1 (MQ=255)

                                                         aaaGTGTACCAACGCCATTTTGCCGAAATGCTTCGCATGAGCGCGCAGCAGCGGCAGCGTAACAAAg          <  1:734482/67‑1 (MQ=255)

                                                         aaaGTGTACCAACGCCATTTTGCCGAAATGCTTCGCATGAGCGCGCAGCAGCGGCAGCGTAACAAAg          <  1:798696/67‑1 (MQ=255)

                                                         aaaGTGTACCAACGCCATTTTGCCGAAATGCTTCGCATGAGCGCGCAGCAGCGGCAGCGTAACAAAg          <  1:3307292/67‑1 (MQ=255)

                                                         aaaGTGTACCAACGCCATTTTGCCGAAATGCTTCGCATGAGCGCGCAGCAGCGGCAGCGTAACAAAg          <  1:2272679/67‑1 (MQ=255)

                                                         aaaGTGTACCAACGCCATTTTGCCGAAATGCTTCGCATGAGCGCGCAGCAGCGGCAGCGTAACAAAg          <  1:3243484/67‑1 (MQ=255)

                                                         aaaGTGTACCAACGCCATTTTGCCGAAATGCTTCGCATGAGCGCGCAGCAGCGGCAGCGTAACAAAg          <  1:3166264/67‑1 (MQ=255)

                                                         aaaGTGTACCAACGCCATTTTGCCGAAATGCTTCGCATGAGCGCGCAGCAGCGGCAGCGTAACAAAg          <  1:2208655/67‑1 (MQ=255)

                                                         aaaGTGTACCAACGCCATTTTGCCGAAATGCTTCGCATGAGCGCGCAGCAGCGGCAGCGTAACAAAg          <  1:251282/67‑1 (MQ=255)

                                                         aaaGTGTACCAACGCCATTTTGCCGAAATGCTTCGCATGAGCGCGCAGCAGCGGCAGCGTAACAAAg          <  1:1480436/67‑1 (MQ=255)

                                                         aaaGTGTACCAACGCCATTTTGCCGAAATGCTTCGCATGAGCGCGCAGCAGCGGCAGCGTAACAAAg          <  1:3098964/67‑1 (MQ=255)

                                                         aaaGTGTACCAACGCCATTTTGCCGAAATGCTTCGCATGAGCGCGCAGCAGCGGCAGCGTAACAAAg          <  1:2560700/67‑1 (MQ=255)

                                                         aaaGTGTACCAACGCCATTTTGCCGAAATGCTTCGCATGAGCGCGCAGCAGCGGCAGCGTAACAAAg          <  1:3047176/67‑1 (MQ=255)

                                                         aaaGTGTACCAACGCCATTTTGCCGAAATGCTTCGCATGAGCGCGCAGCAGCGGCAGCGTAACAAAg          <  1:2077667/67‑1 (MQ=255)

                                                         aaaGTGTACCAACGCCATTTTGCCGAAATGCTTCGCATGAGCGCGCAGCAGCGGCAGCGTAACAAAg          <  1:2180333/67‑1 (MQ=255)

                                                         aaaGTGTACCAACGCCATTTTGCCGAAATGCTTCGCATGAGCGCGCAGCAGCGGCAGCGTAACAAAg          <  1:2604381/67‑1 (MQ=255)

                                                         aaaGTGTACCAACGCCATTTTGCCGAAATGCTTCGCATGAGCGCGCAGCAGCGGCAGCGTAACAAAg          <  1:2627269/67‑1 (MQ=255)

                                                         aaaGTGTACCAACGCCATTTTGCCGAAATGCTTCGCATGAGCGCGCAGCAGCGGCAGCGTAACAAAg          <  1:1475515/67‑1 (MQ=255)

                                                         aaaGTGTACCAACGCCATTTTGCCGAAATGCTTCGCATGAGCGCGCAGCAGCGGCAGCGTAACAAAg          <  1:2701839/67‑1 (MQ=255)

                                                         aaaGTGTACCAACGCCATTTTGCCGAAATGCTTCGCATGAGCCCGCAGCAGCGTCAGCGTAACAAAg          <  1:3411204/67‑1 (MQ=255)

                                                             tgtACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAGCAGCGGCAGCGTAACAAAGTGGTc     <  1:584321/68‑1 (MQ=255)

                                                              gtACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAGCAGCGGCAGCGTAACAAAGTGGTcac   <  1:562538/69‑1 (MQ=255)

                                                              gtACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAGCAGCGGCAGCGTAACAAAGTGGTcac   <  1:1584356/69‑1 (MQ=255)

                                                              gtACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAGCAGCGGCAGCGTAACAAAGTGGTcac   <  1:1164103/69‑1 (MQ=255)

                                                              gtACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAGCAGCGGCAGCGTAACAAAGTGGTcac   <  1:1120726/69‑1 (MQ=255)

                                                              gtACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAGCAGCGGCAGCGTAACAAAGTGGTcac   <  1:1902453/69‑1 (MQ=255)

                                                              gtACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAGCAGCGGCAGCGTAACAAAGTGGTcac   <  1:1316948/69‑1 (MQ=255)

                                                              gtACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAGCAGCGGCAGCGTAACAAAGTGGTcac   <  1:234199/69‑1 (MQ=255)

                                                               tACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAGCAGCGGCAGCGTAACAAAGTGGTcacc  >  1:3121805/1‑69 (MQ=255)

                                                               tACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAGCAGCGGCAGCGTAACAAAGTGGTcacc  >  1:1584214/1‑69 (MQ=255)

                                                                aCCAACGCCATTTTGCCGAAATGCTTCGCATGAGCGCGCAGCAGCGGCAGCGTAACAAAg          <  1:1923844/60‑1 (MQ=255)

GAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCATAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAGCAGCGGCAGCGTAACAAAGTGGTCACC  >  minE/1970291‑1970422

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