BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,330,617	G→A	21.2%	A323A (GCC→GCT)	cpsG ←	phosphomannomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,330,617	0	G	A	21.2%	211.9 / 41.5	145	A323A (GCC→GCT)	cpsG	phosphomannomutase
Reads supporting (aligned to +/- strand): ref base G (50/64); new base A (18/13); total (68/77)
Fisher's exact test for biased strand distribution p-value = 2.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 145 total aligned reads displayed.

GTAAGCGAAATCACGGAAGTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGCGTGCC  >  minE/1330566‑1330682

gTAAGCGAAATCACGGAAGTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCAt                                                  >  1:3113859/1‑69 (MQ=255)

gTAAGCGAAATCACGGAAGTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCAt                                                  >  1:412660/1‑69 (MQ=255)

gTAAGCGAAATCACGGAAGTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCAt                                                  >  1:429457/1‑69 (MQ=255)

gTAAGCGAAATCACGGAAGTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCAt                                                  >  1:2616485/1‑69 (MQ=255)

gTAAGCGAAATCACGGAAGTAATGGTGGCCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCAt                                                  >  1:1390572/1‑69 (MQ=255)

                    aaTGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCt                                           <  1:492084/56‑1 (MQ=255)

                    aaTGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCt                                           <  1:2443155/56‑1 (MQ=255)

                    aaTGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCt                                           >  1:1468083/1‑56 (MQ=255)

                    aaTGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCt                                           <  1:3189745/56‑1 (MQ=255)

                    aaTGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCt                                           >  1:3005646/1‑56 (MQ=255)

                     aTGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCt                                           <  1:645253/55‑1 (MQ=255)

                        gtggGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTcc                          >  1:2832139/1‑69 (MQ=255)

                        gtggGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTcc                          >  1:27379/1‑69 (MQ=255)

                        gtggGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTcc                          >  1:2498442/1‑69 (MQ=255)

                        gtggGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTcc                          >  1:2222746/1‑69 (MQ=255)

                        gtggGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTcc                          >  1:1473748/1‑69 (MQ=255)

                        gtggGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTcc                          >  1:2981893/1‑69 (MQ=255)

                        gtggGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTcc                          >  1:3115308/1‑69 (MQ=255)

                        gtggGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTc                           >  1:1764713/1‑68 (MQ=255)

                                tCATTTCGCCACCAT‑‑ATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTTCGGTTTTCGAc                <  1:244677/69‑1 (MQ=255)

                                 cATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCg                  >  1:929843/1‑68 (MQ=255)

                                 cATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGa                 <  1:2365334/69‑1 (MQ=255)

                                 cATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGa                 <  1:324172/69‑1 (MQ=255)

                                 cATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGa                 >  1:2293168/1‑69 (MQ=255)

                                 cATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGa                 <  1:3246428/69‑1 (MQ=255)

                                 cATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGa                 <  1:3267664/69‑1 (MQ=255)

                                 cATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGa                 <  1:2458292/69‑1 (MQ=255)

                                 cATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGa                 <  1:2632478/69‑1 (MQ=255)

                                 cATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGa                 <  1:2639810/69‑1 (MQ=255)

                                 cATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGa                 <  1:1035737/69‑1 (MQ=255)

                                 cATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGa                 <  1:3207122/69‑1 (MQ=255)

                                 cATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGa                 <  1:294780/69‑1 (MQ=255)

                                 cATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGa                 <  1:2920587/69‑1 (MQ=255)

                                 cATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGa                 <  1:2795751/69‑1 (MQ=255)

                                 cATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGa                 <  1:1798146/69‑1 (MQ=255)

                                 cATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGa                 <  1:1185038/69‑1 (MQ=255)

                                 cATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGa                 <  1:600266/69‑1 (MQ=255)

                                 cATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGa                 <  1:823559/69‑1 (MQ=255)

                                 cATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGa                 <  1:90702/69‑1 (MQ=255)

                                 cATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGa                 >  1:165521/1‑69 (MQ=255)

                                 cATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGa                 <  1:1745072/69‑1 (MQ=255)

                                 cATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGa                 <  1:114423/69‑1 (MQ=255)

                                 cATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGa                 <  1:1819249/69‑1 (MQ=255)

                                 cATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGa                 <  1:1103832/69‑1 (MQ=255)

                                 cATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGa                 <  1:1895164/69‑1 (MQ=255)

                                 cATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGa                 <  1:1923026/69‑1 (MQ=255)

                                 cATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGa                 <  1:3271967/69‑1 (MQ=255)

                                  aTTTCGCCACCTTAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTGTCGa                 <  1:3311737/68‑1 (MQ=255)

                                  aTTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGa                 <  1:197060/68‑1 (MQ=255)

                                    ttCGCCACCGTAGATAGCGTCTTCCTTGCGCATACGTTCttt                                         >  1:595270/1‑42 (MQ=38)

                                    ttCGCCACCGTAGATAGCGTCTTCCTTGCGCATACGTTCTTtaat                                      >  1:164412/1‑45 (MQ=255)

                                    ttCGCCACCGTAGATAGCGTCTTCCTTGCGCATACGTTCTTtaat                                      >  1:3063721/1‑45 (MQ=255)

                                    ttCGCCACCGTAGATAGCGTCTTCCTTGCGCATACGTTCTTtaat                                      >  1:1176114/1‑45 (MQ=255)

                                    ttCGCCACCGTAGATAGCGTCTTCCTTGCGCATACGTTCTTtaat                                      >  1:1012224/1‑45 (MQ=255)

                                    ttCGCCACCGTAGATAGCGTCTTCCTTGCGCATACGTTCTTtaat                                      >  1:507133/1‑45 (MQ=255)

                                    ttCGCCACCGTAGATAGCGTCTTCCTTGCGCATACGTTCTTtaat                                      >  1:2701262/1‑45 (MQ=255)

                                    ttCGCCACCGTAGATAGCGTCTTCCTTGCGCATACGTTCTTtaat                                      >  1:2145298/1‑45 (MQ=255)

                                    ttCGCCACCGTAGATAGCGTCTTCCTTGCGCATACGTTCTTtaat                                      >  1:2610870/1‑45 (MQ=255)

                                    ttCGCCACCGTAGATAGCGTCTTCCTTGCGCATACGTTCTTtaat                                      >  1:328369/1‑45 (MQ=255)

                                    ttCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGAc                <  1:3184265/67‑1 (MQ=255)

                                    ttCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGAc                <  1:342751/67‑1 (MQ=255)

                                    ttCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGAc                <  1:3421052/67‑1 (MQ=255)

                                    ttCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGAc                <  1:961302/67‑1 (MQ=255)

                                    ttCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGAc                <  1:1752381/67‑1 (MQ=255)

                                    ttCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGAc                <  1:1049239/67‑1 (MQ=255)

                                    ttCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGAc                <  1:1164308/67‑1 (MQ=255)

                                    ttCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGAc                <  1:1606605/67‑1 (MQ=255)

                                    ttCGCCACCATAGATGGCGCCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGAc                <  1:890639/67‑1 (MQ=255)

                                     tCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACAtt             >  1:1620985/1‑69 (MQ=255)

                                     tCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACAtt             >  1:1387190/1‑69 (MQ=255)

                                     tCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACAtt             >  1:1900984/1‑69 (MQ=255)

                                          accaTAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGttt                     <  1:52958/56‑1 (MQ=255)

                                          accaTAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGttt                     <  1:1773336/56‑1 (MQ=255)

                                          accaTAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGttt                     <  1:1164997/56‑1 (MQ=255)

                                          accaTAGATGGCGGCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGttt                     <  1:1969828/56‑1 (MQ=255)

                                             aTAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGAc                <  1:2767097/58‑1 (MQ=255)

                                             aTAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGCGt     >  1:2058057/1‑69 (MQ=255)

                                             aTAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGCGt     >  1:1865226/1‑69 (MQ=255)

                                             aTAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGCGt     >  1:644441/1‑69 (MQ=255)

                                             aTAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGCGt     >  1:343835/1‑69 (MQ=255)

                                               aGATAGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGGGTGc   <  1:695880/69‑1 (MQ=255)

                                               aGATAGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGGGTGc   <  1:766822/69‑1 (MQ=255)

                                               aGATAGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGGGTGc   <  1:772426/69‑1 (MQ=255)

                                               aGATAGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGGGTGc   <  1:920561/69‑1 (MQ=255)

                                               aGATAGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGGGTGc   <  1:2724464/69‑1 (MQ=255)

                                               aGATAGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGGGTGc   <  1:307092/69‑1 (MQ=255)

                                               aGATAGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGGGTGc   <  1:3435025/69‑1 (MQ=255)

                                               aGATAGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGGGTGc   <  1:205735/69‑1 (MQ=255)

                                               aGATAGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGGGTGc   <  1:2306310/69‑1 (MQ=255)

                                               aGATAGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGGGTGc   <  1:3200649/69‑1 (MQ=255)

                                                gATGGCGTCTTCCTTGGGCATACGTTCTTTAATAAAgg                                 >  1:2170510/1‑38 (MQ=255)

                                                gATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCg                               >  1:306757/1‑40 (MQ=255)

                                                gATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCTGTTTTCGa                 >  1:1147570/1‑54 (MQ=255)

                                                gATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGa                 >  1:2183660/1‑54 (MQ=255)

                                                gATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGa                 >  1:801744/1‑54 (MQ=255)

                                                gATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGCGTGcc  >  1:1204799/1‑69 (MQ=255)

                                                gATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGCGTGcc  >  1:2874256/1‑69 (MQ=255)

                                                gATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGCGTGcc  >  1:3232131/1‑69 (MQ=255)

                                                gATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGCGTGcc  >  1:1093804/1‑69 (MQ=255)

                                                gATAGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGGGTGc   <  1:802731/68‑1 (MQ=255)

GTAAGCGAAATCACGGAAGTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGCGTGCC  >  minE/1330566‑1330682

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