BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	AP009048	4,125,002	A→T	6.7%	D167V (GAT→GTT)	dam →	DNA adenine methylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AP009048	4,125,002	0	A	T	6.7%	336.8 / 13.8	120	D167V (GAT→GTT)	dam	DNA adenine methylase
Reads supporting (aligned to +/- strand): ref base A (59/53); new base T (4/4); total (63/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.37e-01

Only 100 of 125 total aligned reads displayed.

CGGAAGCAGAGTTGTATCACTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCGCCTTATGCACCGCTGTCTGCGACC  >  AP009048/4124933‑4125072

cGGAAGCAGAGTTGTATCACTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGAt                                                                       <  1:454830/71‑1 (MQ=255)

 ggAAGCAGAGTTGTATCACTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATa                                                                      <  1:2587636/71‑1 (MQ=255)

 ggAAGCAGAGTTGTATAACTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATa                                                                      <  1:445550/71‑1 (MQ=255)

  gAAGCAGAGTTGTATCACTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAg                                                                     >  1:2423854/1‑71 (MQ=255)

    aGCAGAGTTGTATCACTTCGCTGAAAAAGCGCAGAATTCCTTTTTCTATTGTGAGTCTTACGCCGATAGCa                                                                   <  1:2981070/71‑1 (MQ=255)

      cAGAGTTGTATCACTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGc                                                                    >  1:2386293/1‑68 (MQ=255)

      cAGAGTTGTATCACTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATg                                                                 >  1:4079955/1‑71 (MQ=255)

       agagTTGTATCACTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATa                                                                      >  1:2760762/1‑65 (MQ=255)

       agagTTGTATCACTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATgg                                                                >  1:2718685/1‑71 (MQ=255)

       agagTTGTATCACTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATgg                                                                <  1:160302/71‑1 (MQ=255)

        gagTTGTATCACTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCa                                                                   <  1:1932356/67‑1 (MQ=255)

        gagTTGTATCACTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGgc                                                               >  1:3433089/1‑71 (MQ=255)

          gTTGTATCACTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGgcgc                                                             >  1:3609449/1‑71 (MQ=255)

          gTTGTATCACTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGgcgc                                                             >  1:2142801/1‑71 (MQ=255)

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

                tCACTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGTTAGCATGGCGCGCGCAg                                                       <  1:2600749/71‑1 (MQ=255)

                tCACTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGgc                                                               <  1:2800996/63‑1 (MQ=255)

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

                  aCTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGTTAGCATGGCGCGCGCAga                                                      <  1:2667597/70‑1 (MQ=255)

                  aCTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCTCGCAgat                                                     <  1:2599749/71‑1 (MQ=255)

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

                  aCTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAga                                                      <  1:3113011/70‑1 (MQ=255)

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

                    ttCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCa                                                                   <  1:2518607/55‑1 (MQ=255)

                    ttCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCa                                                        >  1:4177672/1‑66 (MQ=255)

                    ttCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAgatag                                                   >  1:3713988/1‑69 (MQ=255)

                      agCTGAAAAAGCGCAGAATGCCTTTTTCTATTGAGAGTCTTACGCCGTTAGc                                                                    >  1:4235403/2‑52 (MQ=255)

                       gCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGatgatg                                                 >  1:4016756/1‑70 (MQ=255)

                       gCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGatgatg                                                 <  1:1521884/70‑1 (MQ=255)

                        cTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCa                                               <  1:4487191/71‑1 (MQ=255)

                        cTGAAAAAGAGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCa                                               >  1:3879876/1‑71 (MQ=255)

                         tGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGTTAGCATGGCGCGCGCAGATGATGCAt                                              >  1:3872771/1‑71 (MQ=255)

                         tGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAgat                                                     <  1:3837490/64‑1 (MQ=255)

                          gAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGgc                                                               >  1:43423/1‑53 (MQ=255)

                           aaaaaGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCAt                                              <  1:3865628/69‑1 (MQ=255)

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

                           aaaaaGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATcc                                            >  1:4339066/1‑71 (MQ=255)

                            aaaaGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCcg                                           >  1:1444743/1‑71 (MQ=255)

                             aaaGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCcgt                                          <  1:3346920/71‑1 (MQ=255)

                             aaaGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCcg                                           <  1:846163/70‑1 (MQ=255)

                             aaaGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCcg                                           >  1:4446495/1‑70 (MQ=255)

                              aaGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCcgtc                                         >  1:3882932/1‑71 (MQ=255)

                               aGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCcgtcg                                        >  1:838884/1‑71 (MQ=255)

                                    aGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAga                                                      >  1:4076056/1‑52 (MQ=255)

                                    aGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAga                                                      >  1:1837815/1‑52 (MQ=255)

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

                                      aaTGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCa                                                        <  1:784060/48‑1 (MQ=255)

                                      aaTGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAgat                                                     <  1:2261333/51‑1 (MQ=255)

                                       aTGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAgatgat                                                  <  1:2948245/53‑1 (MQ=255)

                                        tGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAgatga                                                   >  1:3888527/1‑51 (MQ=255)

                                        tGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCg                               <  1:1051375/71‑1 (MQ=255)

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

                                          ccTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCg                               >  1:1703207/1‑69 (MQ=255)

                                          ccTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGAt                             >  1:2110404/1‑71 (MQ=255)

                                          ccTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGAt                             >  1:1222140/1‑71 (MQ=255)

                                           cTTTTTCTATTGTGAGTCTTACGCCGATAGCATGgcgc                                                             <  1:2150838/38‑1 (MQ=255)

                                            tttttCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAgatgat                                                  >  1:1505356/1‑48 (MQ=255)

                                              tttCTATTGTGAGTCTTACGCCGTTAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTAtt                                  >  1:2079571/1‑62 (MQ=255)

                                              tttCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATcc                                            >  1:1411354/1‑52 (MQ=255)

                                              tttCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCcgt                                          >  1:110002/1‑54 (MQ=255)

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

                                              tttCTATTGTGAGTCTTACGCCGATAGCATGGCGCACGCAGATGATGCATCCGTCGTCTATTGCGATCCGc                         >  1:216977/1‑71 (MQ=255)

                                               ttCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCGcc                        >  1:930722/1‑71 (MQ=255)

                                                tCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCa                                               >  1:2137762/1‑47 (MQ=255)

                                                 cTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCGCCt                       <  1:4111930/70‑1 (MQ=255)

                                                  tATTGTGAGTCTTACGCCGTTAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCg                          <  1:300450/66‑1 (MQ=255)

                                                  tATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCGCCTTa                     >  1:3831822/1‑71 (MQ=255)

                                                    ttGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGa                              <  1:510642/60‑1 (MQ=255)

                                                    ttGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCGCCTTATg                   >  1:4392504/1‑71 (MQ=255)

                                                     tgtgAGTCTTACGCCGATAGCATGGCGCGCGCAga                                                      <  1:536819/35‑1 (MQ=255)

                                                     tgtgAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCGCCTTATg                   >  1:1081462/1‑70 (MQ=255)

                                                      gtgACTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCGCCTTATGCa                 >  1:3301426/1‑71 (MQ=255)

                                                        gAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCGCCTTATGCAcc               >  1:3095841/1‑71 (MQ=255)

                                                         aGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCGCCt                       <  1:3263113/62‑1 (MQ=255)

                                                          gTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCcgtcgt                                       >  1:970530/1‑45 (MQ=255)

                                                          gTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCGCCTTATGCACCg              <  1:2270229/70‑1 (MQ=255)

                                                           tCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCGCCTTATGCACCGCt            <  1:2308921/71‑1 (MQ=255)

                                                           tCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCGCCTTATGCACCGCt            <  1:943775/71‑1 (MQ=255)

                                                            cTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTa                                    >  1:2678972/1‑46 (MQ=255)

                                                             ttACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGc                                <  1:1070497/49‑1 (MQ=255)

                                                             ttACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATcc                           >  1:2895199/1‑54 (MQ=255)

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

                                                              tACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCg                          >  1:1009913/1‑54 (MQ=255)

                                                               aCGCCGTTAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCGCCTTATGCACCGCTGt          <  1:1698558/69‑1 (MQ=255)

                                                               aCGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCGCCTTAt                    >  1:620917/1‑59 (MQ=255)

                                                               aCGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCGCCTTATGCACCg              >  1:2230111/1‑65 (MQ=255)

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

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

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

                                                                cGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCGCCTTATGCACCGCTGTCTg       <  1:4164873/71‑1 (MQ=255)

                                                                cGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCGCCTTAAGCACCGCTGTCt        <  1:1498130/70‑1 (MQ=255)

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

                                                                   cGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCGCCTTATGCACCGCTGTCTGCGa    <  1:1357893/71‑1 (MQ=255)

                                                                    gATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCGCCTTATGCACCGCTGTCTGCGAc   >  1:1234020/1‑71 (MQ=255)

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

                                                                     aTAGCATGGCGCGCGCATATGATGCATCCGTCGTCTATTg                                 <  1:1972921/40‑1 (MQ=255)

                                                                     aTAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGa                              >  1:3601070/1‑43 (MQ=255)

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

                                                                     aTAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCGCCTTATGCACCGCTGTCTGCGAcc  >  1:1638682/1‑71 (MQ=255)

                                                                     aTAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCGCCTTATGCACCGCTGTCTGCGAcc  >  1:2105523/1‑71 (MQ=255)

CGGAAGCAGAGTTGTATCACTTCGCTGAAAAAGCGCAGAATGCCTTTTTCTATTGTGAGTCTTACGCCGATAGCATGGCGCGCGCAGATGATGCATCCGTCGTCTATTGCGATCCGCCTTATGCACCGCTGTCTGCGACC  >  AP009048/4124933‑4125072

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