BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,904,710	C→T	40.6%	S404S (TCG→TCA)	ptr ←	protease III

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,904,710	0	C	T	40.6%	50.2 / 79.8	113	S404S (TCG→TCA)	ptr	protease III
Reads supporting (aligned to +/- strand): ref base C (29/38); new base T (26/20); total (55/58)
Fisher's exact test for biased strand distribution p-value = 1.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.38e-01

Only 100 of 113 total aligned reads displayed.

CTCAACAGGAACGCGAATCATGGTATCTGCCAGCCATTCGACGTAATCCATATCACGGGTGATCGACGGATAACGGAAGTCGATATCCAGCACATTCGCCAGTTCATCGAAGTATTGTTTATC  >  minE/1904647‑1904769

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

              gAATCATGGTATCGGCCAGCCATTCGACGTAATCCATATCACGGGTGATTGACGGATAACGGAAGTCGa                                          >  1:2382311/1‑69 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                              cagccaTTCGACGTAATCCATATCACGGGTGATCGACGGATAACGGAAGTCGATATCCAGcaca                               <  1:2446597/64‑1 (MQ=255)

                                  caTTCGACGTAATCCATATCACGGGTGATCGACGGATAACGGAAGTCGATATCCAGCACATTCGCCAg                       >  1:218208/1‑68 (MQ=255)

                                  caTTCGACGTAATCCATATCACGGGTGATCGACGGATAACGGAAGTCGATATCCAGCACATTCGCCAg                       >  1:2210520/1‑68 (MQ=255)

                                  caTTCGACGTAATCCATATCACGGGTGATCGACGGATAACGGAAGTCGATATCCAGCACATTCGCCAg                       >  1:1250885/1‑68 (MQ=255)

                                  caTTCGACGTAATCCATATCACGGGTGATCGACGGATAACGGAAGTCGATATCCAGCACATTCGCCAg                       >  1:717215/1‑68 (MQ=255)

                                  caTTCGACGTAATCCATATCACGGGTGATCGACGGATAACGGAAGTCGATATCCAGCACATTCGCCAg                       >  1:1773789/1‑68 (MQ=255)

                                      cGACGTAATCCATATCACGGGTGATCGACGGATAACGGAAGTCGATATCCAGCACATTCGCCAGTTCAt                  <  1:1636880/69‑1 (MQ=255)

                                          gTAATCCATATCACGGGTGATCGACGGATAACGGAAGTCGATATCCAGCACATTCGCCAGTTCATCGa               >  1:1203676/1‑68 (MQ=255)

                                           tAATCCATATCACGGGTGATCGACGGATAACGGAAGTCGATATCCAGCACATTCGCCAGTTCATCGaa              >  1:2790726/1‑68 (MQ=255)

                                                cATATCACGGGTGATCGACGGATAACGGAAGTCGATATCCAGCACATTCGCCAGTTCATCGAAGTAtt         >  1:2158341/1‑68 (MQ=255)

                                                cATATCACGGGTGATCGACGGATAACGGAAGTCGATATCCAGCACATTCGCCAGTTCATCGAAGTATTg        >  1:701555/1‑69 (MQ=255)

                                                cATATCACGGGTGATCGACGGATAACGGAAGTCGATATCCAGCACATTCGCCAGTTCATCGAAGTATTg        >  1:1550340/1‑69 (MQ=255)

                                                cATATCACGGGTGATCGACGGATAACGGAAGTCGATATCCAGCACATTCGCCAGTTCATCGAAGTATTg        >  1:104140/1‑69 (MQ=255)

                                                cATATCACGGGTGATCGACGGATAACGGAAGTCGATATCCAGCACATTCGCCAGTTCATCGAAGTATTg        >  1:1684376/1‑69 (MQ=255)

                                                cATATCACGGGTGATCGACGGATAACGGAAGTCGATATCCAGCACATTCGCCAGTTCATCGAAGTATTg        >  1:2525746/1‑69 (MQ=255)

                                                cATATCACGGGTGATCGACGGATAACGGAAGTCGATATCCAGCACATTCGCCAGTTCATCGAAGTATTg        >  1:3112941/1‑69 (MQ=255)

                                                cATATCACGGGTGATCGACGGATAACGGAAGTCGATATCCAGCACATTCGCCAGTTCATCGAAGTATTg        >  1:1900574/1‑69 (MQ=255)

                                                cATATCACGGGTGATCGACGGATAACGGAAGTCGATATCCAGCACATTCGCCAGTTCATCGAAGTATTg        >  1:1816385/1‑69 (MQ=255)

                                                cATATCACGGGTGATCGACGGATAACGGAAGTCGATATCCAGCACATTCGCCAGTTCATCGAAGTATTg        >  1:2988568/1‑69 (MQ=255)

                                                cATATCACGGGTGATCGACGGATAACGGAAGTCGATATCCAGCACATTCGCCAGTTCATCGAAGTATTg        >  1:2969686/1‑69 (MQ=255)

                                                    tCACGGGTGATCGACGGATAACGGAAGTCGATATCCAGCACATTCGCCAGTTCATCGAAGTATTGttt     >  1:1770437/1‑68 (MQ=255)

                                                    tCACGGGTGATCGACGGATAACGGAAGTCGATATCCAGCACATTCGCCAGTTCATCGAAGTATTGttt     >  1:777499/1‑68 (MQ=255)

                                                    tCACGGGTGATCGACGGATAACGGAAGTCGATATCCAGCACATTCGCCAGTTCATCGAAGTATTGttt     >  1:775133/1‑68 (MQ=255)

                                                      aCGTGTGATCGACGGATAACGGAAGTCGATATCCAGCACATTCGCCAGTTCATCGAAGTATTGTTTATc  <  1:788712/69‑1 (MQ=255)

                                                      aCGGGTGATCGACGGATAACGTAAGTCGATATCCAGCACATTCGCCAGTTCATCGAAGTATTGTTTATc  <  1:2589031/69‑1 (MQ=255)

                                                      aCGGGTGATCGACGGATAACGGAAGTCGATATCCAGCACATTCGCCAGTTCATCGAAGTATTGTTTATc  <  1:574287/69‑1 (MQ=255)

                                                      aCGGGTGATCGACGGATAACGGAAGTCGATATCCAGCACATTCGCCAGTTCATCGAAGTATTGTTTATc  <  1:1215243/69‑1 (MQ=255)

                                                      aCGGGTGATCGACGGATAACGGAAGTCGATATCCAGCACATTCGCCAGTTCATCGAAGTATTGTTTATc  <  1:2116996/69‑1 (MQ=255)

                                                      aCGGGTGATCGACGGATAACGGAAGTCGATATCCAGCACATTCGCCAGTTCATCGAAGTATTGTTTATc  <  1:346302/69‑1 (MQ=255)

                                                      aCGGGTGATCGACGGATAACGGAAGTCGATATCCAGCACATTCGCCAGTTCATCGAAGTATTGTTTATc  <  1:3217105/69‑1 (MQ=255)

                                                      aCGGGTGATCGACGGATAACGGAAGTCGATATCCAGCACATTCGCCAGTTCATCGAAGTATTGTTTATc  <  1:3399108/69‑1 (MQ=255)

                                                       cGGGTGATCGACGGATAACGGAAGTCGATATCCAGCACATTCGCCAGTTCATCGAAGTATTGTTTATc  <  1:805488/68‑1 (MQ=255)

                                                       cGGGTGATCGACGGATAACGGAAGTCGATATCCAGCACATTCGCCAGTTCATCGAAGTATTGTTTATc  <  1:2483578/68‑1 (MQ=255)

                                                             aTCGACGGATAACGGAAGTCGATATCCAGCACATTCGCCAGTTCATCGAAGTATTGTTTATc  <  1:229897/62‑1 (MQ=255)

CTCAACAGGAACGCGAATCATGGTATCTGCCAGCCATTCGACGTAATCCATATCACGGGTGATCGACGGATAACGGAAGTCGATATCCAGCACATTCGCCAGTTCATCGAAGTATTGTTTATC  >  minE/1904647‑1904769

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