BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,366,640	C→T	43.3%	A271A (GCG→GCA)	pepQ ←	proline dipeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,366,640	0	C	T	43.3%	33.5 / 80.2	106	A271A (GCG→GCA)	pepQ	proline dipeptidase
Reads supporting (aligned to +/- strand): ref base C (36/24); new base T (18/28); total (54/52)
Fisher's exact test for biased strand distribution p-value = 4.95e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 108 total aligned reads displayed.

GCTGACGCCTGCTTTCATGGTGGCGATCAGCGCCAGTTGTTCATCATTAACGTCCTTCACCAGCTGCGCGTAGTCGTTGTCACTTTTTGCCGACCAGGTACGGGTCAGGTCAGCCGCATAGCCGTTATAT  >  minE/2366574‑2366703

gCTGACGCCTGCTTTCATGGTGGCGATCAGCGCCAGTTGTTCATCATTAACGTCCTTCACCAGCTgcgc                                                               >  1:1471043/1‑69 (MQ=255)

     cGCCTGCTTTCATGGTGGCGATCAGCGCCAGTTGTTCATCATTAACGTCCTTCACCAGCTGCGCGTAGt                                                          >  1:738628/1‑69 (MQ=255)

     cGCCTGCTTTCATGGTGGCGATCAGCGCCAGTTGTTCATCATTAACGTCCTTCACCAGCTGCGCGTAGt                                                          >  1:3289871/1‑69 (MQ=255)

     cGCCTGCTTTCATGGTGGCGATCAGCGCCAGTTGTTCATCATTAACGTCCTTCACCAGCTGCGCGTAGt                                                          >  1:1914606/1‑69 (MQ=255)

     cGCCTGCTTTCATGGTGGCGATCAGCGCCAGTTGTTCATCATTAACGTCCTTCACCAGCTGCGCGTAGt                                                          >  1:408790/1‑69 (MQ=255)

     cGCCTGCTTTCATGGTGGCGATCAGCGCCAGTTGTTCATCATTAACGTCCTTCACCAGCTGCGCGTAGt                                                          >  1:2223819/1‑69 (MQ=255)

     cGCCTGCTTTCATGGTGGCGATCAGCGCCAGTTGTTCATCATTAACGTCCTTCACCAGCTGCGCGTAGt                                                          >  1:548641/1‑69 (MQ=255)

        cTGCTTTCATGGTGGCGATCAGCGCCAGTTGCTCATCATTTACGTCTTTCACCAGCTGTGCGTAGTCGt                                                       >  1:2357669/1‑69 (MQ=255)

        cTGCTTTCATGGTGGCGATCAGCGCCAGTTGCTCATCATTTACGTCTTTCACCAGCTGTGCGTAGTCGt                                                       >  1:1781320/1‑69 (MQ=255)

        cTGCTTTCATGGTGGCGATCAGCGCCAGTTGCTCATCATTTACGTCTTTCACCAGCTGTGCGTAGTCGt                                                       >  1:741394/1‑69 (MQ=255)

        cTGCTTTCATGGTGGCGATCAGCGCCAGTTGCTCATCATTTACGTCTTTCACCAGCTGTGCGTAGTCGt                                                       >  1:1240817/1‑69 (MQ=255)

         tGCTTTCATGGTGGCGATCAGCGCCAGTTGTTCATCATTAACGTCCTTCACCAGCTGCGCGTAGTCGtt                                                      >  1:2388451/1‑69 (MQ=255)

         tGCTTTCATGGTGGCGATCAGCGCCAGTTGTTCATCATTAACGTCCTTCACCAGCTGCGCGTAGTCGtt                                                      >  1:97441/1‑69 (MQ=255)

         tGCTTTCATGGTGGCGATCAGCGCCAGTTGTTCATCATTAACGTCCTTCACCAGCTGCGCGTAGTCGtt                                                      >  1:2440508/1‑69 (MQ=255)

         tGCTTTCATGGTGGCGATCAGCGCCAGTTGTTCATCATTAACGTCCTTCACCAGCTGCGCGTAGTCGtt                                                      >  1:2486834/1‑69 (MQ=255)

         tGCTTTCATGGTGGCGATCAGCGCCAGTTGTTCATCATTAACGTCCTTCACCAGCTGCGCGTAGTCGtt                                                      >  1:2088719/1‑69 (MQ=255)

         tGCTTTCATGGTGGCGATCAGCGCCAGTTGTTCATCATTAACGTCCTTCACCAGCTGCGCGTAGTCGtt                                                      >  1:2867775/1‑69 (MQ=255)

         tGCTTTCATGGTGGCGATCAGCGCCAGTTGTTCATCATTAACGTCCTTCACCAGCTGCGCGTAGTCGtt                                                      >  1:2960638/1‑69 (MQ=255)

         tGCTTTCATGGTGGCGATCAGCGCCAGTTGTTCATCATTAACGTCCTTCACCAGCTGCGCGTAGTCGtt                                                      >  1:1844028/1‑69 (MQ=255)

         tGCTTTCATGGTGGCGATCAGCGCCAGTTGTTCATCATTAACGTCCTTCACCAGCTGCGCGTAGTCGtt                                                      >  1:1829327/1‑69 (MQ=255)

         tGCTTTCATGGTGGCGATCAGCGCCAGTTGTTCATCATTAACGTCCTTCACCAGCTGCGCGTAGTCGtt                                                      >  1:3033615/1‑69 (MQ=255)

         tGCTTTCATGGTGGCGATCAGCGCCAGTTGTTCATCATTAACGTCCTTCACCAGCTGCGCGTAGTCGtt                                                      >  1:3140172/1‑69 (MQ=255)

         tGCTTTCATGGTGGCGATCAGCGCCAGTTGTTCATCATTAACGTCCTTCACCAGCTGCGCGTAGTCGt                                                       >  1:341839/1‑68 (MQ=255)

              tCATGGTGGCGATCAGCGCCAGTTGCTCATCATTTACGTCTTTCACCAGCTGTGCGTAGTCGTTGTcgc                                                 >  1:1474393/1‑67 (MQ=255)

              tCATGGTGGCGATCAGCGCCAGTTGCTCATCATTTACGTCTTTCACCAGCTGTGCGTAGTCGTTGTcgc                                                 >  1:517002/1‑67 (MQ=255)

              tCATGGTGGCGATCAGCGCCAGTTGCTCATCATTTACGTCTTTCACCAGCTGTGCGTAGTCGTTGTcgc                                                 >  1:1355122/1‑67 (MQ=255)

              tCATGGTGGCGATCAGCGCCAGTTGCTCATCATTTACGTCTTTCACCAGCTGTGCGTAGTCGTTGTcgc                                                 >  1:1358553/1‑67 (MQ=255)

              tCATGGTGGCGATCAGCGCCAGTTGCTCATCATTTACGTCTTTCACCAGCTGTGCGTAGTCGTTGTcgc                                                 >  1:574774/1‑67 (MQ=255)

              tCATGGTGGCGATCAGCGCCAGTTGCTCATCATTTACGTCTTTCACCAGCTGTGCGTAGTCGTTGTcgc                                                 >  1:2885810/1‑67 (MQ=255)

              tCATGGTGGCGATCAGCGCCAGTTGCTCATCATTTACGTCTTTCACCAGCTGTGCGTAGTCGTTGTcgc                                                 >  1:885521/1‑67 (MQ=255)

              tCATGGTGGCGATCAGCGCCAGTTGCTCATCATTTACGTCTTTCACCAGCTGTGCGTAGTCGTTGTcgc                                                 >  1:3174928/1‑67 (MQ=255)

              tCATGGTGGCGATCAGCGCCAGTTGCTCATCATTTACGTCTTTCACCAGCTGTGCGTAGTCGTTGTcgc                                                 >  1:3094585/1‑67 (MQ=255)

              tCATGGTGGCGATCAGCGCCAGTTGCTCATCATTTACGTCTTTCACCAGCTGTGCGTAGTCGTTGTcgc                                                 >  1:711108/1‑67 (MQ=255)

              tCATGGTGGCGATCAGCGCCAGTTGCTCATCATTTACGTCTTTCACCAGCTGTGCGTAGTCGTTGTcgc                                                 >  1:72151/1‑67 (MQ=255)

              tCATGGTGGCGATCAGCGCCAGTTGCTCATCATTTACGTCTTTCACCAGCTGTGCGTAGTCGTTGTcgc                                                 >  1:1567342/1‑67 (MQ=255)

                   gtggCGATCAGCGCCAGTTGTTCATCATTAACGTCCTTCACCAGCTgcgc                                                               >  1:175399/1‑50 (MQ=255)

                   gtggCGATCAGCGCCAGTTGTTCATCATTAACGTCCTTCACCAGCTGCGCGTAGTCGTTGTCACttttt                                            >  1:91588/1‑69 (MQ=255)

                   gtggCGATCAGCGCCAGTTGTTCATCATTAACGTCCTTCACCAGCTGCGCGTAGTCGTTGTCACttttt                                            >  1:1889379/1‑69 (MQ=255)

                   gtggCGATCAGCGCCAGTTGTTCATCATTAACGTCCTTCACCAGCTGCGCGTAGTCGTTGTCACttttt                                            >  1:2771701/1‑69 (MQ=255)

                   gtggCGATCAGCGCCAGTTGTTCATCATTAACGTCCTTCACCAGCTGCGCGTAGTCGTTGTCACttttt                                            >  1:1173093/1‑69 (MQ=255)

                   gtggCGATCAGCGCCAGTTGTTCATCATTAACGTCCTTCACCAGCTGCGCGTAGTCGTTGTCACttttt                                            >  1:1437707/1‑69 (MQ=255)

                   gtggCGATCAGCGCCAGTTGTTCATCATTAACGTCCTTCACCAGCTGCGCGTAGTCGTTGTCACttttt                                            >  1:1246940/1‑69 (MQ=255)

                   gtggCGATCAGCGCCAGTTGTTCATCATTAACGTCCTTCACCAGCTGCGCGTAGTCGTTGTCACttttt                                            >  1:1393182/1‑69 (MQ=255)

                      gCGATCAGCGCCAGTTGTTCATCATTAACGTCCTTCACCAGCTGCGCGTAGTCGTTGTCACTTTTTGcc                                         >  1:792303/1‑69 (MQ=255)

                         aTCAGCGCCAGTTGTTCATCATTAACGTCCTTCACCAGCTGCGCGTAGTCGTTGTCACTTTTTGCCGAc                                      <  1:1607323/69‑1 (MQ=255)

                         aTCAGCGCCAGTTGTTCATCATTAACGTCCTTCACCAGCTGCGCGTAGTCGTTGTCACTTTTTGCCGAc                                      <  1:2362996/69‑1 (MQ=255)

                          tCAGCGCCAGTTGTTCATCATTAACGTCCTTCACCAGCTGCGCGTAGTCGTTGTCACTTTTTGCCGAc                                      <  1:1168595/68‑1 (MQ=255)

                           cAGCGCCAGTTGTTCATCATTAACGTCCTTCACCAGCTGCGCGTAGTCGTTGTCACTTTTTGCCGACCa                                    >  1:512433/1‑69 (MQ=255)

                           cAGCGCCAGTTGTTCATCATTAACGTCCTTCACCAGCTGCGCGTAGTCGTTGTCACTTTTTGCCGACCa                                    >  1:416825/1‑69 (MQ=255)

                           cAGCGCCAGTTGTTCATCATTAACGTCCTTCACCAGCTGCGCGTAGTCGTTGTCACTTTTTGCCGACCa                                    >  1:3050827/1‑69 (MQ=255)

                           cAGCGCCAGTTGTTCATCATTAACGTCCTTCACCAGCTGCGCGTAGTCGTTGTCACTTTTTGCCGACCa                                    >  1:1488143/1‑69 (MQ=255)

                           cAGCGCCAGTTGTTCATCATTAACGTCCTTCACCAGCTGCGCGTAGTCGTTGTCACTTTTTGCCGACCa                                    >  1:2310686/1‑69 (MQ=255)

                                      gttCATCATTAACGTCCTTCACCAGCTGCGCGTAGTCGTTGTCACTTTTTGCCGACCAGGTACGggtca                         >  1:1866249/1‑69 (MQ=255)

                                               tAACGTCCTTCACCAGCTGCGCGTAGTCGTTGTCACTTTTTGCCGACCAGGTACGGGTCAGGTCAGCCg                <  1:2381558/69‑1 (MQ=255)

                                               tAACGTCCTTCACCAGCTGCGCGTAGTCGTTGTCACTTTTTGCCGACCAGGTACGGGTCAGGTCAGCCg                <  1:3173568/69‑1 (MQ=255)

                                               tAACGTCCTTCACCAGCTGCGCGTAGTCGTTGTCACTTTTTGCCGACCAGGTACGGGTCAGGTCAGCCg                <  1:2198828/69‑1 (MQ=255)

                                                aaCGTCCTTCACCAGCTGCGCGTAGTCGTTGTCACTTTTTGCCGACCAGGTACGGGTCAGGTCAGCCGc               <  1:2786919/69‑1 (MQ=255)

                                                aaCGTCCTTCACCAGCTGCGCGTAGTCGTTGTCACTTTTTGCCGACCAGGTACGGGTCAGGTCAGCCGc               <  1:1505545/69‑1 (MQ=255)

                                                aaCGTCCTTCACCAGCTGCGCGTAGTCGTTGTCACTTTTTGCCGACCAGGTACGGGTCAGGTCAGCCGc               <  1:1489203/69‑1 (MQ=255)

                                                aaCGTCCTTCACCAGCTGCGCGTAGTCGTTGTCACTTTTTGCCGACCAGGTACGGGTCAGGTCAGCCGc               <  1:2908513/69‑1 (MQ=255)

                                                aaCGTCCTTCACCAGCTGCGCGTAGTCGTTGTCACTTTTTGCCGACCAGGTACGGGTCAGGTCAGCCGc               <  1:596717/69‑1 (MQ=255)

                                                aaCGTCCTTCACCAGCTGCGCGTAGTCGTTGTCACTTTTTGCCGACCAGGTACGGGTCAGGTCAGCCGc               <  1:1942566/69‑1 (MQ=255)

                                                aaCGTCCTTCACCAGCTGCGCGTAGTCGTTGTCACTTTTTGCCGACCAGGTACGGGTCAGGTCAGCCGc               <  1:2689964/69‑1 (MQ=255)

                                                aaCGTCCTTCACCAGCTGCGCGTAGTCGTTGTCACTTTTTGCCGACCAGGTACGGGTCAGGTCAGCCGc               <  1:203128/69‑1 (MQ=255)

                                                aaCGTCCTTCACCAGCTGCGCGTAGTCGTTGTCACTTTTTGCCGACCAGGTACGGGTCAGGTCAGCCGc               <  1:2575520/69‑1 (MQ=255)

                                                aaCGTCCTTCACCAGCTGCGCGTAGTCGTTGTCACTTTTTGCCGACCAGGTACGGGTCAGGTCAGCCGc               <  1:2290939/69‑1 (MQ=255)

                                                aaCGTCCTTCACCAGCTGCGCGTAGTCGTTGTCACTTTTTGCCGACCAGGTACGGGTCAGGTCAGCCGc               <  1:779336/69‑1 (MQ=255)

                                                   gTCCTTCACCAGCTGCGCGTAGTCGTTGTCACTTTTTGCCGACCAGGTACGGGTCAGGTCAGCCGCATa            <  1:379188/69‑1 (MQ=255)

                                                   gTCCTTCACCAGCTGCGCGTAGTCGTTGTCACTTTTTGCCGACCAGGTACGGGTCAGGTCAGCCGCATa            <  1:623417/69‑1 (MQ=255)

                                                    tCCTTCACCAGCTGCGCGTAGTCGTTGTCACTTTTTGCCGACCAGGTACGGGTCAGGTCAGCCGCATAg           >  1:863852/1‑69 (MQ=255)

                                                    tCCTTCACCAGCTGCGCGTAGTCGTTGTCACTTTTTGCCGACCAGGTACGGGTCAGGTCAGCCGCATAg           >  1:2415407/1‑69 (MQ=255)

                                                    tCCTTCACCAGCTGCGCGTAGTCGTTGTCACTTTTTGCCGACCAGGTACGGGTCAGGTCAGCCGCATAg           >  1:1428100/1‑69 (MQ=255)

                                                        tCACCAGCTGCGCGTAGTCGTTGTCACTTTTTGCCGACCAGGTACGGGTCAGGTCAGCCGCATAGCCGt       <  1:1810331/69‑1 (MQ=255)

                                                        tCACCAGCTGCGCGTAGTCGTTGTCACTTTTTGCCGACCAGGTACGGGTCAGGTCAGCCGCATAGCCGt       <  1:1318052/69‑1 (MQ=255)

                                                        tCACCAGCTGCGCGTAGTCGTTGTCACTTTTTGCCGACCAGCTACGGGTCAGGTCAGCCGCATAGCCGt       <  1:1021858/69‑1 (MQ=255)

                                                            cAGCTGCGCGTAGTCGTTGTCACTTTTTGCCGACCAGGTACGGGTCAGGTCAGCCGCATAGCCGTtata   <  1:2274765/69‑1 (MQ=255)

                                                             aGCTGTGCGTAGTCGTTGTCGCTTTTTGCCGACCAGGTACTGGTTAGATCAGCCGCATAGCCGTTatat  <  1:1000996/69‑1 (MQ=255)

                                                             aGCTGTGCGTAGTCGTTGTCGCTTTTTGCCGACCAGGTACGGGTTAGATCAGCCGCATAGCCTTTatat  <  1:7113/69‑1 (MQ=255)

                                                             aGCTGTGCGTAGTCGTTGTCGCTTTTTGCCGACCAGGTACGGGTTAGATCAGCCGCATAGCCGTtatgt  <  1:2366187/69‑3 (MQ=255)

                                                             aGCTGTGCGTAGTCGTTGTCGCTTTTTGCCGACCAGGTACGGGTTAGATCAGCCGCATAGCCGTTatat  <  1:636397/69‑1 (MQ=255)

                                                             aGCTGTGCGTAGTCGTTGTCGCTTTTTGCCGACCAGGTACGGGTTAGATCAGCCGCATAGCCGTTatat  <  1:1342632/69‑1 (MQ=255)

                                                             aGCTGTGCGTAGTCGTTGTCGCTTTTTGCCGACCAGGTACGGGTTAGATCAGCCGCATAGCCGTTatat  <  1:1124697/69‑1 (MQ=255)

                                                             aGCTGTGCGTAGTCGTTGTCGCTTTTTGCCGACCAGGTACGGGTTAGATCAGCCGCATAGCCGTTatat  <  1:864624/69‑1 (MQ=255)

                                                             aGCTGTGCGTAGTCGTTGTCGCTTTTTGCCGACCAGGTACGGGTTAGATCAGCCGCATAGCCGTTatat  <  1:1571262/69‑1 (MQ=255)

                                                             aGCTGTGCGTAGTCGTTGTCGCTTTTTGCCGACCAGGTACGGGTTAGATCAGCCGCATAGCCGTTatat  <  1:1040354/69‑1 (MQ=255)

                                                             aGCTGTGCGTAGTCGTTGTCGCTTTTTGCCGACCAGGTACGGGTTAGATCAGCCGCATAGCCGTTatat  <  1:737463/69‑1 (MQ=255)

                                                             aGCTGTGCGTAGTCGTTGTCGCTTTTTGCCGACCAGGTACGGGTTAGATCAGCCGCATAGCCGTTatat  <  1:1655130/69‑1 (MQ=255)

                                                             aGCTGTGCGTAGTCGTTGTCGCTTTTTGCCGACCAGGTACGGGTTAGATCAGCCGCATAGCCGTTatat  <  1:1015064/69‑1 (MQ=255)

                                                             aGCTGTGCGTAGTCGTTGTCGCTTTTTGCCGACCAGGTACGGGTTAGATCAGCCGCATAGCCGTTatat  <  1:1713928/69‑1 (MQ=255)

                                                             aGCTGTGCGTAGTCGTTGTCGCTTTTTGCCGACCAGGTACGGGTTAGATCAGCCGCATAGCCGTTatat  <  1:3102015/69‑1 (MQ=255)

                                                             aGCTGTGCGTAGTCGTTGTCGCTTTTTGCCGACCAGGTACGGGTTAGATCAGCCGCATAGCCGTTatat  <  1:1900851/69‑1 (MQ=255)

                                                             aGCTGTGCGTAGTCGTTGTCGCTTTTTGCCGACCAGGTACGGGTTAGATCAGCCGCATAGCCGTTatat  <  1:1940986/69‑1 (MQ=255)

                                                             aGCTGTGCGTAGTCGTTGTCGCTTTTTGCCGACCAGGTACGGGTTAGATCAGCCGCATAGCCGTTatat  <  1:2361918/69‑1 (MQ=255)

                                                             aGCTGTGCGTAGTCGTTGTCGCTTTTTGCCGACCAGGTACGGGTTAGATCAGCCGCATAGCCGTTatat  <  1:2659742/69‑1 (MQ=255)

                                                             aGCTGTGCGTAGTCGTTGTCGCTTTTTGCCGACCAGGTACGGGTTAGATCAGCCGCATAGCCGTTatat  <  1:2669223/69‑1 (MQ=255)

                                                             aGCTGTGCGTAGTCGTTGTCGCTTTTTGCCGACCAGGTACGGGTTAGATCAGCCGCATAGCCGTTatat  <  1:303057/69‑1 (MQ=255)

                                                             aGCTGTGCGTAGTCGTTGTCGCTTTTTGCCGACCAGGTACGGGTTAGATCAGCCGCATAGCCGTTatat  <  1:304539/69‑1 (MQ=255)

                                                             aGCTGTGCGTAGTCGTTGTCGCTTTTTGCCGACCAGGTACGGGTTAGATCAGCCGCATAGCCGTTatat  <  1:3081977/69‑1 (MQ=255)

                                                             aGCTGTGCGTAGTCGTTGTCGCTTTTTGCCGACCAGGTACGGGTCAGATCAGCCGCATAGCCGTTatat  <  1:730710/69‑1 (MQ=255)

                                                              gCTGTGCGTAGTCGTTGTCGCTTTTTGCCGACCAGGTACGGGTTAGATCAGCCGCATAGCCGTTatat  <  1:3226966/68‑1 (MQ=255)

GCTGACGCCTGCTTTCATGGTGGCGATCAGCGCCAGTTGTTCATCATTAACGTCCTTCACCAGCTGCGCGTAGTCGTTGTCACTTTTTGCCGACCAGGTACGGGTCAGGTCAGCCGCATAGCCGTTATAT  >  minE/2366574‑2366703

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