BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,904,674	T→A	50.2%	A416A (GCA→GCT)	ptr ←	protease III

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,904,674	0	T	A	50.2%	‑4.8 / 215.1	168	A416A (GCA→GCT)	ptr	protease III
Reads supporting (aligned to +/- strand): ref base T (23/60); new base A (27/57); total (51/117)
Fisher's exact test for biased strand distribution p-value = 6.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.58e-01

Only 100 of 169 total aligned reads displayed.

GATCGGCAATATTGACTGCATCCAGCGTATGCTCAACAGGAACGCGAATCATGGTATCTGCCAGCCATTCGACGTAATCCATATCACGGGTGATCGACGGATAACGGAAGTCGATATCCAGCACA  >  minE/1904616‑1904740

gATCGGCAATATTGACTGCATCCAGTGTATGCTCAACAGGAACGCGAATCATGGTATCAGCCAGCCAtt                                                          <  1:2833737/69‑1 (MQ=255)

gATCGGCAATATTGACTGCATCCAGTGTATGCTCAACAGGAACGCGAATCATGGTATCAGCCAGCCAtt                                                          <  1:1634974/69‑1 (MQ=255)

gATCGGCAATATTGACTGCATCCAGTGTATGCTCAACAGGAACGCGAATCATGGTATCAGCCAGCCAtt                                                          <  1:179964/69‑1 (MQ=255)

gATCGGCAATATTGACTGCATCCAGTGTATGCTCAACAGGAACGCGAATCATGGTATCAGCCAGCCAtt                                                          <  1:1823043/69‑1 (MQ=255)

gATCGGCAATATTGACTGCATCCAGTGTATGCTCAACAGGAACGCGAATCATGGTATCAGCCAGCCAtt                                                          <  1:1871524/69‑1 (MQ=255)

gATCGGCAATATTGACTGCATCCAGTGTATGCTCAACAGGAACGCGAATCATGGTATCAGCCAGCCAtt                                                          <  1:197581/69‑1 (MQ=255)

gATCGGCAATATTGACTGCATCCAGTGTATGCTCAACAGGAACGCGAATCATGGTATCAGCCAGCCAtt                                                          <  1:2180554/69‑1 (MQ=255)

gATCGGCAATATTGACTGCATCCAGTGTATGCTCAACAGGAACGCGAATCATGGTATCAGCCAGCCAtt                                                          <  1:2305166/69‑1 (MQ=255)

gATCGGCAATATTGACTGCATCCAGTGTATGCTCAACAGGAACGCGAATCATGGTATCAGCCAGCCAtt                                                          <  1:2743343/69‑1 (MQ=255)

gATCGGCAATATTGACTGCATCCAGTGTATGCTCAACAGGAACGCGAATCATGGTATCAGCCAGCCAtt                                                          <  1:2557727/69‑1 (MQ=255)

gATCGGCAATATTGACTGCATCCAGTGTATGCTCAACAGGAACGCGAATCATGGTATCAGCCAGCCAtt                                                          <  1:2904718/69‑1 (MQ=255)

gATCGGCAATATTGACTGCATCCAGTGTATGCTCAACAGGAACGCGAATCATGGTATCAGCCAGCCAtt                                                          <  1:3097458/69‑1 (MQ=255)

gATCGGCAATATTGACTGCATCCAGTGTATGCTCAACAGGAACGCGAATCATGGTATCAGCCAGCCAtt                                                          <  1:410615/69‑1 (MQ=255)

gATCGGCAATATTGACTGCATCCAGTGTATGCTCAACAGGAACGCGAATCATGGTATCAGCCAGCCAtt                                                          <  1:786716/69‑1 (MQ=255)

gATCGGCAATATTGACTGCATCCAGTGTATGCTCAACAGGAACGCGAATCATGGTATCAGCCAGCCAtt                                                          <  1:821263/69‑1 (MQ=255)

gATCGGCAATATTGACTGCATCCAGTGTATGCTCAACAGGAACGCGAATCATGGTATCAGCCAGCCAtt                                                          <  1:919333/69‑1 (MQ=255)

gATCGGCAATATTGACTGCATCCAGTGTATGCTCAACAGGAACGCGAATCATGGTATCAGCCAGCCAtt                                                          <  1:992186/69‑1 (MQ=255)

gATCGGCAATATTGACTGCATCCAGTGTATGCTCAACAGGAACGCGAATCATGGTATCAGCCAGCCAtt                                                          <  1:1044678/69‑1 (MQ=255)

gATCGGCAATATTGACTGCATCCAGTGTATGCTCAACAGGAACGCGAATCATGGTATCAGCCAGCCAtt                                                          <  1:1054408/69‑1 (MQ=255)

gATCGGCAATATTGACTGCATCCAGTGTATGCTCAACAGGAACGCGAATCATGCTATCAGCCAGCCAtt                                                          <  1:681060/69‑1 (MQ=255)

gATCGGCAATATTGACTGCATCCAGCGTATGCTCAACAGGAACGCGAATCATGGTATCTGCCAGCCAt                                                           <  1:2213167/68‑1 (MQ=255)

 aTCGGCAATATTGACTGCATCCAGTGTATGCTCAACAGGAACGCGAATCATGGTATCAGCCAGCCAtt                                                          <  1:1213428/68‑1 (MQ=255)

          atTGACTGCATCCAGCGTATGCTCAACAGGAACGCGAATCATGGTATCTGCCAGCCATTCGACGTAATc                                                >  1:549806/1‑69 (MQ=255)

           ttGACTGCATCCAGCGTATGCTCAACAGGAACGCGAATCATGGTATCTGCCAGCCATTCGACGTaa                                                  <  1:2765305/66‑1 (MQ=255)

           ttGACTGCATCCAGCGTATGCTCAACAGGAACGCGAATCATGGTATCTGCCAGCCATTCGACGTaa                                                  <  1:2584234/66‑1 (MQ=255)

           ttGACTGCATCCAGCGTATGCTCAACAGGAACGCGAATCATGGTATCTGCCAGCCATTCGACGTaa                                                  <  1:593997/66‑1 (MQ=255)

           ttGACTGCATCCAGCGTATGCTCAACAGGAACGCGAATCATGGTATCTGCCAGCCATTCGACGTAATc                                                >  1:1756682/1‑68 (MQ=255)

           ttGACTGCATCCAGCGTATGCTCAACAGGAACGCGAATCATGGTATCTGCCAGCCATTCGACGTAATc                                                >  1:2506325/1‑68 (MQ=255)

             gACTGCATCCAGCGTATGCTCAACAGGAACGCGAATCATGGTATCTGCCAGCCATTCGACGTAATCCa                                              <  1:3074/68‑1 (MQ=255)

                        gCGTATGCTCAACAGGAACGCGAATCATGGTATCTGCCAGCCATTCGACGTAATCCATATCACGGGTg                                   <  1:2277793/68‑1 (MQ=255)

                        gCGTATGCTCAACAGGAACGCGAATCATGGTATCTGCCAGCCATTCGACGTAATCCATATCACGGGTg                                   <  1:425883/68‑1 (MQ=255)

                        gCGTATGCTCAACAGGAACGCGAATCATGGTATCTGCCAGCCATTCGACGTAATCCATATCACGGGTg                                   <  1:3033390/68‑1 (MQ=255)

                        gCGTATGCTCAACAGGAACGCGAATCATGGTATCTGCCAGCCATTCGACGTAATCCATATCACGGGTg                                   <  1:2860875/68‑1 (MQ=255)

                        gCGTATGCTCAACAGGAACGCGAATCATGGTATCTGCCAGCCATTCGACGTAATCCATATCACGGGTg                                   <  1:862579/68‑1 (MQ=255)

                          gTATGCTCAACAGGAACGCGAATCATGGTATCTGCCAGCCATTCGACGTAATCCATATCACGtgtg                                   <  1:2352173/66‑1 (MQ=255)

                             tGCTCAACAGGAACGCGAATCATGGTATCTGCCAGCCATTCGACGTAATc                                                <  1:2432293/50‑1 (MQ=255)

                             tGCTCAACAGGAACGCGAATCATGGTATCTGCCAGCCATTCGACGTAATc                                                <  1:3000249/50‑1 (MQ=255)

                             tGCTCAACAGGAACGCGAATCATGGTATCTGCCAGCCATTCGACGTAATc                                                <  1:1098675/50‑1 (MQ=255)

                             tGCTCAACAGGAACGCGAATCATGGTATCTGCCAGCCATTCGACGTAATc                                                <  1:2642305/50‑1 (MQ=255)

                             tGCTCAACAGGAACGCGAATCATGGTATCTGCCAGCCATTCGACGTAATc                                                <  1:2712042/50‑1 (MQ=255)

                             tGCTCAACAGGAACGCGAATCATGGTATCTGCCAGCCATTCGACGTAATc                                                <  1:2668912/50‑1 (MQ=255)

                             tGCTCAACAGGAACGCGAATCATGGTATCTGCCAGCCATTCGACGTAATc                                                <  1:2099581/50‑1 (MQ=255)

                             tGCTCAACAGGAACGCGAATCATGGTATCTGCCAGCCATTCGACGTAATc                                                <  1:706205/50‑1 (MQ=255)

                             tGCTCAACAGGAACGCGAATCATGGTATCTGCCAGCCATTCGACGTAATc                                                <  1:757232/50‑1 (MQ=255)

                             tGCTCAACAGGAACGCGAATCATGGTATCTGCCAGCCATTCGACGTAATc                                                <  1:1702709/50‑1 (MQ=255)

                               cTCAACAGGAACGCGAATCATGGTATCTGCCAGCCATTCGACGTAATCCATATCACGGGTGATCGAc                             >  1:2610758/1‑67 (MQ=255)

                                     aGGAACGCGAATCATGGTATCTGCCAGCCATTCGACGTaa                                                  <  1:505680/40‑1 (MQ=255)

                                     aGGAACGCGAATCATGGTATCTGCCAGCCATTCGACGTa                                                   <  1:3162543/39‑1 (MQ=255)

                                     aGGAACGCGAATCATGGTATCTGCCAGCCATTCGACGTa                                                   <  1:2448912/39‑1 (MQ=255)

                                     aGGAACGCGAATCATGGTATCTGCCAGCCATTCGACGTa                                                   <  1:1478005/39‑1 (MQ=255)

                                     aGGAACGCGAATCATGGTATCTGCCAGCCATTCGACGTa                                                   <  1:1563033/39‑1 (MQ=255)

                                     aGGAACGCGAATCAAGGTATCTGCCAGCCATTCGACGTa                                                   <  1:979804/39‑1 (MQ=255)

                                      ggAACGCGAATCATGGTATCAGCCAGCCATTCGACGTAATCCATATCACGGGTGATTGACGGATAACgg                    <  1:3208059/69‑1 (MQ=255)

                                      ggAACGCGAATCATGGTATCAGCCAGCCATTCGACGTAATCCATATCACGGGTGATTGACGGATAACgg                    <  1:3407954/69‑1 (MQ=255)

                                      ggAACGCGAATCATGGTATCAGCCAGCCATTCGACGTAATCCATATCACGGGTGATTGACGGATAACgg                    <  1:1068391/69‑1 (MQ=255)

                                      ggAACGCGAATCATGGTATCAGCCAGCCATTCGACGTAATCCATATCACGGGTGATTGACGGATAACgg                    <  1:1062732/69‑1 (MQ=255)

                                      ggAACGCGAATCATGGTATCAGCCAGCCATTCGACGTAATCCATATCACGGGTGATTGACGGATAACgg                    <  1:1535658/69‑1 (MQ=255)

                                      ggAACGCGAATCATGGTATCAGCCAGCCATTCGACGTAATCCATATCACGGGTGATTGACGGATAACgg                    <  1:2167000/69‑1 (MQ=255)

                                      ggAACGCGAATCATGGTATCAGCCAGCCATTCGACGTAATCCATATCACGGGTGATTGACGGATAACgg                    <  1:23035/69‑1 (MQ=255)

                                      ggAACGCGAATCATGGTATCAGCCAGCCATTCGACGTAATCCATATCACGGGTGATTGACGGATAACgg                    <  1:1261370/69‑1 (MQ=255)

                                      ggAACGCGAATCATGGTATCAGCCAGCCATTCGACGTAATCCATATCACGGGTGATTGACGGATAACgg                    <  1:1203383/69‑1 (MQ=255)

                                      ggAACGCGAATCATGGTATCAGCCAGCCATTCGACGTAATCCATATCACGGGTGATTGACGGATAACgg                    <  1:2357567/69‑1 (MQ=255)

                                      ggAACGCGAATCATGGTATCAGCCAGCCATTCGACGTAATCCATATCACGGGTGATTGACGGATAACgg                    <  1:139876/69‑1 (MQ=255)

                                      ggAACGCGAATCATGGTATCAGCCAGCCATTCGACGTAATCCATATCACGGGTGATTGACGGATAACgg                    <  1:1729678/69‑1 (MQ=255)

                                      ggAACGCGAATCATGGTATCAGCCAGCCATTCGACGTAATCCATATCACGGGTGATTGACGGATAACgg                    <  1:2655725/69‑1 (MQ=255)

                                        aaCGCGAATCATGGTATCTGCCAGCCATTCGACGTAATCCATATCACGGGt                                    >  1:1071724/1‑51 (MQ=255)

                                        aaCGCGAATCATGGTATCTGCCAGCCATTCGACGTAATCCATATCACGGGTg                                   >  1:2932668/1‑52 (MQ=255)

                                        aaCGCGAATCATGGTATCTGCCAGCCATTCGACGTAATCCATATCACGGGTg                                   >  1:964405/1‑52 (MQ=255)

                                        aaCGCGAATCATGGTATCTGCCAGCCATTCGACGTAATCCATATCACGGGTg                                   >  1:1683705/1‑52 (MQ=255)

                                        aaCGCGAATCATGGTATCTGCCAGCCATTCGACGTAATCCATATCACGGGTg                                   >  1:258657/1‑52 (MQ=255)

                                        aaCGCGAATCATGGTATCTGCCAGCCATTCGACGTAATCCATATCACGGGTg                                   >  1:663831/1‑52 (MQ=255)

                                           gcgAATCATGGTATCTGCCAGCCATTCGACGTAATCCATATCACGGGTGATCGACGGATAACGGAAGTc               <  1:192342/69‑1 (MQ=255)

                                             gAATCATGGTATCAGCCAGCCATTCGACGTAATCCatat                                           >  1:695293/1‑39 (MQ=255)

                                             gAATCATGGTATCAGCCAGCCATTCGACGTAATCCATATCACGGGTGAt                                 >  1:1242069/1‑49 (MQ=255)

                                             gAATCATGGTATCAGCCAGCCATTCGACGTAATCCATATCACGGGTGATTGACGGATAACGGAAGTc               >  1:1926716/1‑67 (MQ=255)

                                             gAATCATGGTATCAGCCAGCCATTCGACGTAATCCATATCACGGGTGATTGACGGATAACGGAAGTGGa             >  1:1156915/1‑69 (MQ=255)

                                             gAATCATGGTATCAGCCAGCCATTCGACGTAATCCATATCACGGGTGATTGACGGATAACGGAAGTCg              >  1:1090589/1‑68 (MQ=255)

                                             gAATCATGGTATCAGCCAGCCATTCGACGTAATCCATATCACGGGTGATTGACGGATAACGGAAGTCGa             >  1:1904680/1‑69 (MQ=255)

                                             gAATCATGGTATCAGCCAGCCATTCGACGTAATCCATATCACGGGTGATTGACGGATAACGGAAGTCGa             >  1:632927/1‑69 (MQ=255)

                                             gAATCATGGTATCAGCCAGCCATTCGACGTAATCCATATCACGGGTGATTGACGGATAACGGAAGTCGa             >  1:1421481/1‑69 (MQ=255)

                                             gAATCATGGTATCAGCCAGCCATTCGACGTAATCCATATCACGGGTGATTGACGGATAACGGAAGTCGa             >  1:149296/1‑69 (MQ=255)

                                             gAATCATGGTATCAGCCAGCCATTCGACGTAATCCATATCACGGGTGATTGACGGATAACGGAAGTCGa             >  1:49264/1‑69 (MQ=255)

                                             gAATCATGGTATCAGCCAGCCATTCGACGTAATCCATATCACGGGTGATTGACGGATAACGGAAGTCGa             >  1:2876992/1‑69 (MQ=255)

                                             gAATCATGGTATCAGCCAGCCATTCGACGTAATCCATATCACGGGTGATTGACGGATAACGGAAGTCGa             >  1:2951984/1‑69 (MQ=255)

                                             gAATCATGGTATCAGCCAGCCATTCGACGTAATCCATATCACGGGTGATTGACGGATAACGGAAGTCGa             >  1:2226764/1‑69 (MQ=255)

                                             gAATCATGGTATCAGCCAGCCATTCGACGTAATCCATATCACGGGTGATTGACGGATAACGGAAGTCGa             >  1:3229253/1‑69 (MQ=255)

                                             gAATCATGGTATCAGCCAGCCATTCGACGTAATCCATATCACGGGTGATTGACGGATAACGGAAGTCGa             >  1:3289456/1‑69 (MQ=255)

                                             gAATCATGGTATCAGCCAGCCATTCGACGTAATCCATATCACGGGTGATTGACGGATAACGGAAGTCGa             >  1:3276178/1‑69 (MQ=255)

                                               aTCATGGTATCTGCCAGCCATTCGACGTAATCCATATCACGGGTGATCGACGGATAACGGAAGTCGata           <  1:1145960/69‑1 (MQ=255)

                                               aTCATGGTATCTGCCAGCCATTCGACGTAATCCATATCACGGGTGATCGACGGATAACGGAAGTCGata           <  1:2411936/69‑1 (MQ=255)

                                               aTCATGGTATCTGCCAGCCATTCGACGTAATCCATATCACGGGTGATCGACGGATAACGGAAGTCGata           <  1:3389931/69‑1 (MQ=255)

                                               aTCATGGTATCTGCCAGCCATTCGACGTAATCCATATCACGGGTGATCGACGGATAACGGAAGTCGata           <  1:218739/69‑1 (MQ=255)

                                               aTCATGGTATCTGCCAGCCATTCGACGTAATCCATATCACGGGTGATCGACGGATAACGGAAGTCGata           <  1:560917/69‑1 (MQ=255)

                                               aTCATGGTATCTGCCAGCCATTCGACGTAATCCATATCACGGGTGATCGACGGATAACGGAAGTCGata           <  1:2120959/69‑1 (MQ=255)

                                               aTCATGGTATCTGCCAGCCATTCGACGTAACCCATATCACGGGTGATCGACGGATAACGGAAGTCGata           <  1:249410/69‑1 (MQ=255)

                                                  aTGGTATCTGCCAGCCATTCGACGTAATCCATATCACGGGTGATCGACGGATAACGGAAGTCGATATcc        >  1:322543/1‑69 (MQ=255)

                                                        tCTGCCAGCCATTCGACGTAATCCATATCACGGGTGATCGACGGATAACGGAAGTCGATATCCAGcaca  <  1:657211/69‑1 (MQ=255)

                                                        tCTGCCAGCCATTCGACGTAATCCATATCACGGGTGATCGACGGATAACGGAAGTCGATATCCAGcaca  <  1:3134992/69‑1 (MQ=255)

                                                        tCTGCCAGCCATTCGACGTAATCCATATCACGGGTGATCGACGGATAACGGAAGTCGATATCCAGcaca  <  1:205885/69‑1 (MQ=255)

                                                        tCTGCCAGCCATTCGACGTAATCCATATCACGGGTGATCGACGGATAACGGAAGTCGATATCCAGcaca  <  1:1000354/69‑1 (MQ=255)

GATCGGCAATATTGACTGCATCCAGCGTATGCTCAACAGGAACGCGAATCATGGTATCTGCCAGCCATTCGACGTAATCCATATCACGGGTGATCGACGGATAACGGAAGTCGATATCCAGCACA  >  minE/1904616‑1904740

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