BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,705,682	A→T	8.1%	D167V (GAT→GTT)	dam →	DNA adenine methylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,705,682	0	A	T	8.1%	483.8 / 28.2	173	D167V (GAT→GTT)	dam	DNA adenine methylase
Reads supporting (aligned to +/- strand): ref base A (90/69); new base T (7/7); total (97/76)
Fisher's exact test for biased strand distribution p-value = 7.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

Only 100 of 183 total aligned reads displayed.

CCGGAAGCAGAGTTGTATCACTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCGCCTTATGCACCGCTGTCTGCGACC  >  minE/2705612‑2705752

ccGGAAGCAGAGTTGTATCACTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGa                                                                        <  1:3101540/71‑1 (MQ=255)

   gAAGCAGAGTTGTATCACTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAg                                                                     <  1:1336613/71‑1 (MQ=255)

    aaGCAGAGTTGTATCACTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAg                                                                     >  1:2970905/1‑70 (MQ=255)

    aaGCAGAGTTGTATCACTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAg                                                                     >  1:2741238/1‑70 (MQ=255)

    aaGCAGAGTTGTATCACTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGc                                                                    <  1:2385351/71‑1 (MQ=255)

     aGCAGAGTTGTATCACTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGTTAGCa                                                                   <  1:3914419/71‑1 (MQ=255)

        agagTTGTATCACTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGTTAGCATgg                                                                >  1:3861250/1‑71 (MQ=255)

          agTTGTATCACTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGgcg                                                              >  1:585531/1‑71 (MQ=255)

              gTATCACTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGAt                                                                       >  1:2195869/1‑58 (MQ=255)

                aTCACTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGgcgcgcg                                                          >  1:1109384/1‑69 (MQ=255)

                 tCACTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGgcgcgcg                                                          >  1:2115528/1‑68 (MQ=255)

                 tCACTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGgcgcgcg                                                          >  1:2598282/1‑68 (MQ=255)

                 tCACTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCa                                                        <  1:799093/70‑1 (MQ=255)

                 tCACTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCa                                                        <  1:3642654/70‑1 (MQ=255)

                 tCACTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCa                                                        <  1:286163/70‑1 (MQ=255)

                 tCACTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAg                                                       <  1:2419617/71‑1 (MQ=255)

                 tCACTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAg                                                       <  1:3361924/71‑1 (MQ=255)

                 tCACTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAg                                                       <  1:66578/71‑1 (MQ=255)

                  cACTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATa                                                                      <  1:2742687/55‑1 (MQ=255)

                   aCTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGTTAGCATGGCGCGCGCAgat                                                     >  1:3263304/1‑71 (MQ=255)

                   aCTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGTTAGCATGGCGCGCGCAgat                                                     <  1:1517382/71‑1 (MQ=255)

                   aCTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAgat                                                     <  1:522962/71‑1 (MQ=255)

                   aCTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAgat                                                     <  1:3745474/71‑1 (MQ=255)

                    cTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATATCATGGCGCGCGCAgatg                                                    <  1:3044385/71‑1 (MQ=255)

                    cTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAgatg                                                    <  1:2837789/71‑1 (MQ=255)

                    cTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAgat                                                     >  1:3593212/1‑70 (MQ=255)

                     ttCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAgatga                                                   >  1:3356786/1‑71 (MQ=255)

                     ttCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAgatga                                                   >  1:25705/1‑71 (MQ=255)

                     ttCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAgatga                                                   >  1:2444609/1‑71 (MQ=255)

                     ttCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAgatga                                                   >  1:3053020/1‑71 (MQ=255)

                      tCGCTGAAAAAGCGCAGAATGCCTTTTTCTAATGTGAGTCTTACGCCGATa                                                                      >  1:2036968/1‑51 (MQ=255)

                       cGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAgatga                                                   <  1:1842792/69‑1 (MQ=255)

                       cGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAgatga                                                   <  1:3394512/69‑1 (MQ=255)

                       cGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGatgatg                                                 >  1:3217947/1‑71 (MQ=255)

                       cGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGatgatg                                                 >  1:2306242/1‑71 (MQ=255)

                           gAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAgatg                                                    <  1:1989743/64‑1 (MQ=255)

                           gAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGc                                                >  1:2526501/1‑68 (MQ=255)

                            aaaaaGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGTTAGCATGGCGCGCGCAGATGATGCATcc                                            >  1:846590/1‑71 (MQ=255)

                            aaaaaGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCa                                                                   >  1:3300286/1‑48 (MQ=255)

                            aaaaaGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATcc                                            <  1:1877256/71‑1 (MQ=255)

                             aaaaGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGTTAGCATGGCGCGCGCAGATGATGCATCcg                                           <  1:286228/71‑1 (MQ=255)

                              aaaGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATcc                                            >  1:805702/1‑69 (MQ=255)

                                   gcAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAg                                                       >  1:2519234/1‑53 (MQ=255)

                                   gcAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCCgtcgtc                                      >  1:462295/1‑70 (MQ=255)

                                   gcAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCCgtcgtc                                      >  1:1669903/1‑70 (MQ=255)

                                   gcAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCCgtcgtc                                      <  1:3172311/70‑1 (MQ=255)

                                    cAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCt                                     <  1:2696498/70‑1 (MQ=255)

                                     aGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAg                                                                     <  1:1347180/37‑1 (MQ=255)

                                     aGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCa                                                                   <  1:3702152/39‑1 (MQ=255)

                                     aGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTAt                                   >  1:1164796/1‑71 (MQ=255)

                                       aaTGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCcgt                                          <  1:1391705/62‑1 (MQ=255)

                                       aaTGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTg                                 >  1:3031259/1‑71 (MQ=255)

                                       aaTGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCACAGATGATGCATCCGTCGTCTATTg                                 <  1:1286527/71‑1 (MQ=255)

                                          gCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGa                              >  1:2629883/1‑71 (MQ=255)

                                           ccTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGa                              <  1:3158732/70‑1 (MQ=255)

                                             tttttCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATcc                           >  1:1219966/1‑71 (MQ=255)

                                              ttttCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGc                                <  1:382199/65‑1 (MQ=255)

                                              ttttCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCg                          >  1:641812/1‑71 (MQ=255)

                                               tttCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCGc                         >  1:3946022/1‑71 (MQ=255)

                                                ttCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCcgt                                          >  1:1998428/1‑53 (MQ=255)

                                                ttCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGAt                             >  1:2697633/1‑66 (MQ=255)

                                                ttCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCGc                         >  1:3559899/1‑70 (MQ=255)

                                                 tCTATTGTGAGTCTTACGCCGTTAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCGcc                        <  1:3775952/70‑1 (MQ=255)

                                                 tCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCGCCt                       >  1:1797774/1‑71 (MQ=255)

                                                 tCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCGCCt                       >  1:1203300/1‑71 (MQ=255)

                                                  cTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCGcc                        <  1:3318720/69‑1 (MQ=255)

                                                  cTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCGCCt                       >  1:3719060/1‑70 (MQ=255)

                                                   tATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGAt                             <  1:1710585/63‑1 (MQ=255)

                                                     ttGTGAGTCTTACGCCGATAGCATGGCGCGCGCAgatga                                                   >  1:220967/1‑39 (MQ=255)

                                                     ttGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCGCCtt                      <  1:570287/68‑1 (MQ=255)

                                                        tgAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCGCCtt                      <  1:2799373/65‑1 (MQ=255)

                                                          aGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCGCCTTATGCa                 >  1:1478400/1‑68 (MQ=255)

                                                          aGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCGCCTTATGCACCg              >  1:1053277/1‑71 (MQ=255)

                                                              ttACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGc                                >  1:3118011/1‑49 (MQ=255)

                                                              ttACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGa                              <  1:2039998/51‑1 (MQ=255)

                                                              ttACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCGCCTTATGCACCGCTg           >  1:1201211/1‑70 (MQ=255)

                                                              ttACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCGCCTTATGCACCGCTGt          >  1:3378991/1‑71 (MQ=255)

                                                              ttACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCGCCTTATGCACCGCTGt          >  1:385792/1‑71 (MQ=255)

                                                               tACGCCGATAGCATGGCGCTCGCAGATGATGCATCCGTCGTCTATTGCGATCCGcc                        <  1:2927822/56‑1 (MQ=255)

                                                                aCGCCGTTAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCGCCTTATGCACCGCTGTCt        >  1:2309433/1‑71 (MQ=255)

                                                                aCGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGa                              >  1:3754371/1‑49 (MQ=255)

                                                                aCGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCGCGTTATGCACCGCTGTCt        >  1:3890118/1‑71 (MQ=255)

                                                                aCGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCGCCTTATGCACCGCTGt          <  1:288732/69‑1 (MQ=255)

                                                                aCGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCGCCTTATGCACCGCTGTc         <  1:3461571/70‑1 (MQ=255)

                                                                aCGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCGCCTTATGCACCGCTGTc         >  1:76027/1‑70 (MQ=255)

                                                                aCGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCGCCTTATGCACCGCTGTCt        >  1:952657/1‑71 (MQ=255)

                                                                aCGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCGCCTTATGCACCGCTGTCt        >  1:489556/1‑71 (MQ=255)

                                                                 cGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCGCCTTATGCACCGCTGTCTg       >  1:1005761/1‑71 (MQ=255)

                                                                  gCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCGCCTTATGCACCGCTGTCTGc      >  1:2064480/1‑71 (MQ=255)

                                                                   ccGTTAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCGCCTTATGCACCGCTGTCTGCg     <  1:1920420/71‑1 (MQ=255)

                                                                    cGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCt                                     >  1:1441503/1‑38 (MQ=255)

                                                                    cGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCGCCtt                      >  1:91286/1‑53 (MQ=255)

                                                                    cGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCGCCTTATGCACCGCTGt          <  1:309022/65‑1 (MQ=255)

                                                                    cGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCGCCTTATGCACCGCTGTCTGCGa    >  1:2502450/1‑71 (MQ=255)

                                                                    cGATAGCATGGCGCGCGCAGATGATGCACCCGTCGTCTATTGCGATCCGCCTTATGCACCGCTGTCTGCGa    >  1:3710197/1‑71 (MQ=255)

                                                                     gATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCGCCTTATGCACCGCTGTCTGCGAc   <  1:3184069/71‑1 (MQ=255)

                                                                      aTAGCATGGCGCGCGCAGATGATGCATCCgtcgtc                                      <  1:2409655/35‑1 (MQ=255)

                                                                      aTAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCGCCt                       <  1:403128/50‑1 (MQ=255)

                                                                      aTAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCGCCTTa                     >  1:1441295/1‑52 (MQ=255)

                                                                      aTAGCATGGCGCGCGCAGATAATGCATCCGTCGTCTATTGCGATCCGCCTTATGCACCGCTGTCTGCGAc   <  1:879226/70‑1 (MQ=255)

CCGGAAGCAGAGTTGTATCACTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCGCCTTATGCACCGCTGTCTGCGACC  >  minE/2705612‑2705752

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