BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,295,597	T→C	19.7%	T67T (ACA→ACG)	yijE ←	predicted permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,295,597	0	T	C	19.7%	144.6 / 25.5	101	T67T (ACA→ACG)	yijE	predicted permease
Reads supporting (aligned to +/- strand): ref base T (55/26); new base C (6/14); total (61/40)
Fisher's exact test for biased strand distribution p-value = 4.06e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

Only 100 of 101 total aligned reads displayed.

GCGCCAGACCAACCATCCCGCAGGTTTGTAACAGGGCAATGGCTAAGGTGTATTTAAACGGTGTCGGGCGCATTCCGCGACCACGTAATAAAAGGACGATGAATAAAACGAGAGCGCCGAAAA  >  minE/2295536‑2295658

gcgcCAGACCAACCATCCCGCAGGTTTGTAACAGGGCAATGGCTAAGGTGTATTTAAACGGTGTCggg                                                         >  1:2662528/1‑68 (MQ=255)

gcgcCAGACCAACCATCCCGCAGGTTTGTAACAGGGCAATGGCTAAGGTGTATTTAAACGGTGTCggg                                                         >  1:816700/1‑68 (MQ=255)

gcgcCAGACCAACCATCCCGCAGGTTTGTAACAGGGCAATGGCTAAGGTGTATTTAAACGGTGTCggg                                                         >  1:2070666/1‑68 (MQ=255)

gcgcCAGACCAACCATCCCGCAGGTTTGTAACAGGGCAATGGCTAAGGTGTATTTAAACGGTGTCggg                                                         >  1:2930047/1‑68 (MQ=255)

gcgcCAGACCAACCATCCCGCAGGTTTGTAACAGGGCAATGGCTAAGGTGTATTTAAACGGTGTCggg                                                         >  1:2956325/1‑68 (MQ=255)

   ccAGACCAACCATCCCGCAGGTTTGTAACAGGGCAATGGCTAAGGTGTATTTAAACGGTGTCGGgcg                                                       <  1:330824/67‑1 (MQ=255)

    cAGACCAACCATCCCGCAGGTTTGTAACAGGGCAATGGCTAAGGTGTATTTAAACGGTGTCGGgcg                                                       <  1:714280/66‑1 (MQ=255)

         caaccaTCCCGCAGGTTTGTAACAGGGCAATGGCTAAGGTGTATTTAAACGGCGTCGGGCGCATTCcgc                                               <  1:3023742/69‑1 (MQ=255)

         caaccaTCCCGCAGGTTTGTAACAGGGCAATGGCTAAGGTGTATTTAAACGGCGTCGGGCGCATTCcgc                                               <  1:2053059/69‑1 (MQ=255)

         caaccaTCCCGCAGGTTTGTAACAGGGCAATGGCTAAGGTGTATTTAAACGGCGTCGGGCGCATTCcgc                                               <  1:95021/69‑1 (MQ=255)

         caaccaTCCCGCAGGTTTGTAACAGGGCAATGGCTAAGGTGTATTTAAACGGCGTCGGGCGCATTCcgc                                               <  1:3176894/69‑1 (MQ=255)

         caaccaTCCCGCAGGTTTGTAACAGGGCAATGGCTAAGGTGTATTTAAACGGCGTCGGGCGCATTCcgc                                               <  1:330438/69‑1 (MQ=255)

         caaccaTCCCGCAGGTTTGTAACAGGGCAATGGCTAAGGTGTATTTAAACGGCGTCGGGCGCATTCcgc                                               <  1:2152079/69‑1 (MQ=255)

         caaccaTCCCGCAGGTTTGTAACAGGGCAATGGCTAAGGTGTATTTAAACGGCGTCGGGCGCATTCcgc                                               <  1:1660187/69‑1 (MQ=255)

         caaccaTCCCGCAGGTTTGTAACAGGGCAATGGCTAAGGTGTATTTAAACGGCGTCGGGCGCATTCcgc                                               <  1:2881906/69‑1 (MQ=255)

         caaccaTCCCGCAGGTTTGTAACAGGGCAATGGCTAAGGTGTATTTAAACGGCGTCGGGCGCATTCcgc                                               <  1:245985/69‑1 (MQ=255)

         caaccaTCCCGCAGGTTTGTAACAGGGCAATGGCTAAGGTGTATTTAAACGGCGTCGGGCGCATTCcgc                                               <  1:2464098/69‑1 (MQ=255)

         caaccaTCCCGCAGGTTTGTAACAGGGCAATGGCTAAGGTGTATTTAAACGGCGTCGGGCGCATTCcgc                                               <  1:468189/69‑1 (MQ=255)

         caaccaTCCCGCAGGTTTGTAACAGGGCAATGGCTAAGGTGTATTTAAACGGCGTCGGGCGCATTCcgc                                               <  1:2614378/69‑1 (MQ=255)

         caaccaTCCCGCAGGTTTGTAACAGGGCAATGGCTAAGGGGTATTTAAACGGCGTCGGGCGCATTCcgc                                               <  1:767301/69‑1 (MQ=255)

          aaccaTCCCGCAGGTTTGTAACAGGGCAATGGCTAAGGTGTATTTAAACGGCGTCGGGCGCATTCcgc                                               <  1:865477/68‑1 (MQ=255)

           accaTCCCGCAGGTTTGTAACAGGGCAATGGCTAAGGTGTATTTAAACGGTGTCGGGCGCATTCcgcg                                              >  1:1082589/1‑68 (MQ=255)

           accaTCCCGCAGGTTTGTAACAGGGCAATGGCTAAGGTGTATTTAAACGGTGTCGGGCGCATTCcgcg                                              >  1:3050720/1‑68 (MQ=255)

           accaTCCCGCAGGTTTGTAACAGGGCAATGGCTAAGGTGTATTTAAACGGTGTCGGGCGCATTCcgcg                                              >  1:3340314/1‑68 (MQ=255)

                   gCAGGTTTGTAACAGGGCAATGGCTAAGGTGTATTTAAACGGTGTCGGGCGCATTCCGTGACCACGta                                      <  1:601321/68‑1 (MQ=255)

                   gCAGGTTTGTAACAGGGCAATGGCTAAGGTGTATTTAAACGGTGTCGGGCGCATTCCGCGACCACGta                                      <  1:78496/68‑1 (MQ=255)

                   gCAGGTTTGTAACAGGGCAATGGCTAAGGTGTATTTAAACGGTGTCGGGCGCATTCCGCGACCACGta                                      <  1:1345070/68‑1 (MQ=255)

                   gCAGGTTTGTAACAGGGCAATGGCTAAGGTGTATTTAAACGGTGTCGGGCGCATTCCGCGACCACGta                                      <  1:2087750/68‑1 (MQ=255)

                   gCAGGTTTGTAACAGGGCAATGGCTAAGGTGTATTTAAACGGTGTCGGGCGCATTCCGCGACCACGta                                      <  1:2081238/68‑1 (MQ=255)

                    cAGGTTTGTAACAGGGCAATGGCTAAGGTGTATTTAAACGGTGTCGGGCGCATTCCGCGACCACGtaat                                    <  1:1723636/69‑1 (MQ=255)

                    cAGGTTTGTAACAGGGCAATGGCTAAGGTGTATTTAAACGGTGTCGGGCGCATTCCGCGACCACGtaat                                    <  1:2168677/69‑1 (MQ=255)

                      ggTTTGTAACAGGGCAATGGCTAAGGTGTATTTAAACGGTGTCGGGCGCATTCCGCGACCACGtaataa                                  <  1:1540196/69‑1 (MQ=255)

                       gTTTGTAACAGGGCAATGGCTAAGGTGTATTTAAACGGTGTCGGGCGCATTCCGCGACCACGtaataa                                  <  1:2813393/68‑1 (MQ=255)

                        tttGTAAGAGGGCAATGGCTAAGGTGTATTTTAACGGTGTCGGGCGCATTCCGCGACCACGtaat                                    <  1:308452/65‑1 (MQ=255)

                        tttGTAACAGGGCAATGGCTAAGGTGTATTTAAACGGTGTCGGGCGCATTCCGCGACCACGTAATaaa                                 >  1:2484279/1‑68 (MQ=255)

                         ttGTAACAGGGCAATGGCTAAGGTGTATTTAAACGGTGTCGGGCGCATTCCGCGACCACGTAATAAAAg                               <  1:1965117/69‑1 (MQ=255)

                                aGGGCAATGGCTAAGGTGTATTTAAACGGTGTCGGGCGCATTCCGCGACCACGTAATAAAAGGACGATg                        <  1:3168660/69‑1 (MQ=255)

                                aGGGCAATGGCTAAGGTGTATTTAAACGGTGTCGGGCGCATTCCGCGACCACGTAATAAAAGGACGATg                        <  1:1495773/69‑1 (MQ=255)

                                aGGGCAATGGCTAAGGTGTATTTAAACGGTGTCGGGCGCATTCCGCGACCACGTAATAAAAGGACGATg                        <  1:2520509/69‑1 (MQ=255)

                                        ggCTAAGGTGTATTTAAACGGTGTCGGGCTCATTCCGCGACCACGTAATAAAAGGACGATGAATaaa                  >  1:1446920/1‑67 (MQ=255)

                                        ggCTAAGGTGTATTTAAACGGTGTCGGGCGCATTCCGCGACCACGTAATAAAAGGACGATGAATaaa                  >  1:2865670/1‑67 (MQ=255)

                                        ggCTAAGGTGTATTTAAACGGTGTCGGGCGCATTCCGCGACCACGTAATAAAAGGACGATGAATaaa                  >  1:679481/1‑67 (MQ=255)

                                        ggCTAAGGTGTATTTAAACGGTGTCGGGCGCATTCCGCGACCACGTAATAAAAGGACGATGAATaaa                  >  1:3321191/1‑67 (MQ=255)

                                        ggCTAAGGTGTATTTAAACGGTGTCGGGCGCATTCCGCGACCACGTAATAAAAGGACGATGAATaaa                  >  1:440704/1‑67 (MQ=255)

                                        ggCTAAGGTGTATTTAAACGGTGTCGGGCGCATTCCGCGACCACGTAATAAAAGGACGATGAATaaa                  >  1:3400614/1‑67 (MQ=255)

                                        ggCTAAGGTGTATTTAAACGGTGTCGGGCGCATTCCGCGACCACGTAATAAAAGGACGATGAATaaa                  >  1:3400917/1‑67 (MQ=255)

                                        ggCTAAGGTGTATTTAAACGGTGTCGGGCGCATTCCGCGACCACGTAATAAAAGGACGATGAATaaa                  >  1:376833/1‑67 (MQ=255)

                                        ggCTAAGGTGTATTTAAACGGTGTCGGGCGCATTCCGCGACCACGTAATAAAAGGACGATGAATaaa                  >  1:3083678/1‑67 (MQ=255)

                                        ggCTAAGGTGTATTTAAACGGTGTCGGGCGCATTCCGCGACCACGTAATAAAAGGACGATGAATaaa                  >  1:2179631/1‑67 (MQ=255)

                                        ggCTAAGGTGTATTTAAACGGTGTCGGGCGCATTCCGCGACCACGTAATAAAAGGACGATGAATaaa                  >  1:1929803/1‑67 (MQ=255)

                                        ggCTAAGGTGTATTTAAACGGTGTCGGGCGCATTCCGCGACCACGTAATAAAAGGACGATGAATaaa                  >  1:1779815/1‑67 (MQ=255)

                                        ggCTAAGGTGTATTTAAACGGTGTCGGGCGCATTCCGCGACCACGTAATAAAAGGACGATGAATaaa                  >  1:1758785/1‑67 (MQ=255)

                                        ggCTAAGGTGTATTTAAACGGTGTCGGGCGCATTCCGCGACCACGTAATAAAAGGACGATGAATaaa                  >  1:1998987/1‑67 (MQ=255)

                                        ggCTAAGGTGTATTTAAACGGTGTCGGGCGCATTCCGCGACCACGTAATAAAAGGACGATGAATaaa                  >  1:2039870/1‑67 (MQ=255)

                                        ggCTAAGGTGTATTTAAACGGTGTCGGGCGCATTCCGCGACCACGTAATAAAAGGACGATGAATaaa                  >  1:1712834/1‑67 (MQ=255)

                                        ggCTAAGGTGTATTTAAACGGTGTCGGGCGCATTCCGCGACCACGTAATAAAAGGACGATGAATaaa                  >  1:216811/1‑67 (MQ=255)

                                        ggCTAAGGTGTATTTAAACGGTGTCGGGCGCATTCCGCGACCACGTAATAAAAGGACGATGAATaaa                  >  1:2282206/1‑67 (MQ=255)

                                        ggCTAAGGTGTATTTAAACGGTGTCGGGCGCATTCCGCGACCACGTAATAAAAGGACGATGAATaaa                  >  1:913184/1‑67 (MQ=255)

                                        ggCTAAGGTGTATTTAAACGGTGTCGGGCGCATTCCGCGACCACGTAATAAAAGGACGATGAATaaa                  >  1:1072219/1‑67 (MQ=255)

                                        ggCTAAGGTGTATTTAAACGGTGTCGGGCGCATTCCGCGACCACGTAATAAAAGGACGATGAATaaa                  >  1:2552131/1‑67 (MQ=255)

                                        ggCTAAGGTGTATTTAAACGGTGTCGGGCGCATTCCGCGACCACGTAATAAAAGGACGATGAATaaa                  >  1:1119319/1‑67 (MQ=255)

                                        ggCTAAGGTGTATTTAAACGGTGTCGGGCGCATTCCGCGACCACGTAATAAAAGGACGATGAATaa                   >  1:3325725/1‑66 (MQ=255)

                                         gCTAAGGTGTATTTAAACGGTGTCGGGCGCATTCCGCGACCACGTAATAAAAGGACGATGAATaaa                  >  1:2555669/1‑66 (MQ=255)

                                         gCTAAGGTGTATTTAAACGGTGTCGGGCGCATTCCGCGACCACGAAATAAAAGGACGATGAATaaa                  >  1:1291923/1‑66 (MQ=255)

                                             aGGTGTATTTAAACGGCGTCGGGCGCATTCCGCGACCACGTAATAAAAGGACGATGAATAAAACGagag           >  1:695083/1‑69 (MQ=255)

                                             aGGTGTATTTAAACGGCGTCGGGCGCATTCCGCGACCACGTAATAAAAGGACGATGAATAAAACGagag           >  1:1468018/1‑69 (MQ=255)

                                             aGGTGTATTTAAACGGCGTCGGGCGCATTCCGCGACCACGTAATAAAAGGACGATGAATAAAACGagag           >  1:1716805/1‑69 (MQ=255)

                                             aGGTGTATTTAAACGGCGTCGGGCGCATTCCGCGACCACGTAATAAAAGGACGATGAATAAAACGagag           >  1:1823328/1‑69 (MQ=255)

                                             aGGTGTATTTAAACGGCGTCGGGCGCATTCCGCGACCACGTAATAAAAGGACGATGAATAAAACGagag           >  1:2084100/1‑69 (MQ=255)

                                             aGGTGTATTTAAACGGCGTCGGGCGCATTCCGCGACCACGTAATAAAAGGACGATGAATAAAACGagag           >  1:2914135/1‑69 (MQ=255)

                                               gtgtATTTAAACGGTGTCGTGCGCATTCCGCGACCACGTAATAAAAGGACGATGAATAAAACGAGAgc          <  1:1313823/68‑1 (MQ=255)

                                               gtgtATTTAAACGGTGTCGGGCGCATTCCGCGACCACGTAATAAAAGGACGATGAATAAAACGAGAgcg         >  1:3300617/1‑69 (MQ=255)

                                               gtgtATTTAAACGGTGTCGGGCGCATTCCGCGACCACGTAATAAAAGGACGATGAATAAAACGAGAgcg         >  1:1983773/1‑69 (MQ=255)

                                               gtgtATTTAAACGGTGTCGGGCGCATTCCGCGACCACGTAATAAAAGGACGATGAATAAAACGAGAgcg         >  1:225847/1‑69 (MQ=255)

                                               gtgtATTTAAACGGTGTCGGGCGCATTCCGCGACCACGTAATAAAAGGACGATGAATAAAACGAGAgcg         >  1:869122/1‑69 (MQ=255)

                                               gtgtATTTAAACGGTGTCGGGCGCATTCCGCGACCACGTAATAAAAGGACGATGAATAAAACGAGAgcg         >  1:781384/1‑69 (MQ=255)

                                               gtgtATTTAAACGGTGTCGGGCGCATTCCGCGACCACGTAATAAAAGGACGATGAATAAAACGAGAgc          <  1:3260157/68‑1 (MQ=255)

                                                    tttAAACGGTGTCGGGCGCATTCTGCGACCACGTAATAAAAGGACGATGAATAAAACGAGAGCGCCGaa    >  1:258622/1‑69 (MQ=255)

                                                    tttAAACGGTGTCGGGCGCATTCCGCGACCACGTAATAAAAGGACGATGAATAAAACGAGAGCGCCGaa    >  1:930203/1‑69 (MQ=255)

                                                    tttAAACGGTGTCGGGCGCATTCCGCGACCACGTAATAAAAGGACGATGAATAAAACGAGAGCGCCGaa    >  1:768543/1‑69 (MQ=255)

                                                    tttAAACGGTGTCGGGCGCATTCCGCGACCACGTAATAAAAGGACGATGAATAAAACGAGAGCGCCGaa    >  1:952534/1‑69 (MQ=255)

                                                    tttAAACGGTGTCGGGCGCATTCCGCGACCACGTAATAAAAGGACGATGAATAAAACGAGAGCGCCGaa    >  1:952853/1‑69 (MQ=255)

                                                    tttAAACGGTGTCGGGCGCATTCCGCGACCACGTAATAAAAGGACGATGAATAAAACGAGAGCGCCGaa    <  1:3349093/69‑1 (MQ=255)

                                                    tttAAACGGTGTCGGGCGCATTCCGCGACCACGTAATAAAAGGACGATGAATAAAACGAGAGCGCCGaa    <  1:1319354/69‑1 (MQ=255)

                                                    tttAAACGGTGTCGGGCGCATTCCGCGACCACGTAATAAAAGGACGATGAATAAAACGAGAGCGCCGaa    >  1:2059592/1‑69 (MQ=255)

                                                    tttAAACGGTGTCGGGCGCATTCCGCGACCACGTAATAAAAGGACGATGAATAAAACGAGAGCGCCGaa    <  1:2092174/69‑1 (MQ=255)

                                                    tttAAACGGTGTCGGGCGCATTCCGCGACCACGTAATAAAAGGACGATGAATAAAACGAGAGCGCCGaa    >  1:2451445/1‑69 (MQ=255)

                                                    tttAAACGGTGTCGGGCGCATTCCGCGACCACGTAATAAAAGGACGATGAATAAAACGAGAGCGCCGaa    >  1:2492140/1‑69 (MQ=255)

                                                    tttAAACGGTGTCGGGCGCATTCCGCGACCACGTAATAAAAGGACGATGAATAAAACGAGAGCGCCGaa    <  1:2610106/69‑1 (MQ=255)

                                                    tttAAACGGTGTCGGGCGCATTCCGCGACCACGTAATAAAAGGACGATGAATAAAACGAGAGCGCCGaa    <  1:3052406/69‑1 (MQ=255)

                                                    tttAAACGGTGTCGGGCGCATTCCGCGACCACGTAATAAAAGGACGATGAATAAAACGAGAGCGCCGaa    <  1:3065449/69‑1 (MQ=255)

                                                    tttAAACGGTGTCGGGCGCATTCCGCGACCACGTAATAAAAGGACGATGAATAAAACGAGAGCGCCGaa    >  1:3147425/1‑69 (MQ=255)

                                                    tttAAACGGTGTCGGGCGCATTCCGCGACCACGTAATAAAAGGACGATGAATAAAACGAGAGCGCCGaa    <  1:3156186/69‑1 (MQ=255)

                                                    tttAAACGGTGTCGGGCGCATTCCGCGACCACGTAATAAAAGGACGATGAATAAAACGAGAGCGCCGaa    >  1:332303/1‑69 (MQ=255)

                                                    tttAAACGGTGTCGGGCGCATTCCGCGACCACGTAATAAAAGGACGATGAATAAAACGAGAGCGCCGaa    >  1:336337/1‑69 (MQ=255)

                                                      tAAACGGTGTCGGGCGCATTCCGCGCCCACGTAATAAAAGGACGATGAATAAAACGAGAGCGCCGaaaa  >  1:1229687/1‑69 (MQ=255)

                                                      tAAACGGTGTCGGGCGCATTCCGCGACCACGTAATAAAAGGACGATGAATAAAACGAGAGCGCCGaaaa  >  1:2959582/1‑69 (MQ=255)

                                                      tAAACGGTGTCGGGCGCATTCCGCGACCACGTAATAAAAGGACGATGAATAAAACGAGAGCGCCGaaaa  >  1:2654889/1‑69 (MQ=255)

                                                      tAAACGGTGTCGGGCGCATTCCGCGACCACGTAATAAAAGGACGATGAATAAAACGAGAGCGCCGaaaa  >  1:2531853/1‑69 (MQ=255)

                                                          cGGTGTCGGGCGCATTCCGCGACCACGTAATAAAAGGACGATg                        <  1:2495454/43‑1 (MQ=255)

GCGCCAGACCAACCATCCCGCAGGTTTGTAACAGGGCAATGGCTAAGGTGTATTTAAACGGTGTCGGGCGCATTCCGCGACCACGTAATAAAAGGACGATGAATAAAACGAGAGCGCCGAAAA  >  minE/2295536‑2295658

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