BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,093,964	T→A	18.7%	V296V (GTA→GTT)	tdcG ←	L‑serine dehydratase 3

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,093,964	0	T	A	18.7%	243.0 / 40.4	101	V296V (GTA→GTT)	tdcG	L‑serine dehydratase 3
Reads supporting (aligned to +/- strand): ref base T (42/40); new base A (14/5); total (56/45)
Fisher's exact test for biased strand distribution p-value = 1.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 105 total aligned reads displayed.

GGCCAGAAAATAGCGGGCAATTGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTT  >  minE/2093898‑2094034

ggCCAGAAAATAGCGGGCAATTGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGt                                                                        <  1:4442776/67‑1 (MQ=255)

 gCCAGAAAATAGCGGGCAATTGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGc                                                                   <  1:2753528/71‑1 (MQ=255)

   cAGAAAATAGCGGGCAATTGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCg                                                                 <  1:4122573/71‑1 (MQ=255)

   cAGAAAATAGCGGGCAATTGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCg                                                                 <  1:4313325/71‑1 (MQ=255)

    aGAAAATAGCGGGCAATTGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCgg                                                                <  1:4537720/71‑1 (MQ=255)

     gAAAATAGCGGGCAATTGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCATTACTGCCGGa                                                               <  1:117037/71‑1 (MQ=255)

     gAAAATAGCGGGCAATTGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCgg                                                                >  1:624616/1‑70 (MQ=255)

       aaaTAGCGGGCAATTGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGaat                                                             <  1:427164/71‑1 (MQ=255)

       aaaTAGCGGGCAATTGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGaa                                                              >  1:1877559/1‑70 (MQ=255)

           aGCGGGCAATTGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGt                                                                        <  1:4436167/56‑1 (MQ=255)

            gCGGGCAATTGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATga                                                        >  1:3436402/1‑70 (MQ=255)

            gCGGGCAATTGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATg                                                         <  1:1789480/70‑1 (MQ=255)

            gCGGGCAATTGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATg                                                         <  1:3251234/70‑1 (MQ=255)

                 cAATTGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCcgcac                                                   >  1:1937486/1‑71 (MQ=255)

                  aaTTGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGaat                                                             <  1:3958981/60‑1 (MQ=255)

                   aTTGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCcgcacgc                                                 <  1:3938239/71‑1 (MQ=255)

                   aTTGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCcgcacgc                                                 <  1:3536218/71‑1 (MQ=255)

                   aTTGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCcgcacg                                                  >  1:2563832/1‑70 (MQ=255)

                     tGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTa                                                                       >  1:2855360/1‑47 (MQ=255)

                     tGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCgcacgcac                                               >  1:4523193/1‑71 (MQ=255)

                      gACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCcgcacgc                                                 <  1:2546472/68‑1 (MQ=255)

                        ccGCTCGTTTACCGGACGACGGAACTTATCGTAATATGCCAGTACTGCCGGAATAATGCCGCACGCACCGt                                            <  1:1435165/71‑1 (MQ=255)

                        ccGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGt                                            <  1:5650663/71‑1 (MQ=255)

                        ccGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGt                                            <  1:3285351/71‑1 (MQ=255)

                        ccGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGt                                            <  1:2532291/71‑1 (MQ=255)

                         cGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGt                                            >  1:2341850/1‑70 (MQ=255)

                          gCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGt                                            >  1:221553/1‑69 (MQ=255)

                           cTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGaat                                                             <  1:3132883/51‑1 (MQ=255)

                            tCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTa                                          >  1:4310112/1‑69 (MQ=255)

                            tCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTa                                          >  1:5664956/1‑69 (MQ=255)

                            tCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGt                                        >  1:347994/1‑71 (MQ=255)

                            tCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGt                                        >  1:4883251/1‑71 (MQ=255)

                              gTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGt                                        <  1:2980555/69‑1 (MQ=255)

                                 tACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCg                                      >  1:2082869/1‑68 (MQ=255)

                                   ccGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGcc                                 <  1:5308122/71‑1 (MQ=255)

                                    cGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCggt                                                               <  1:1635642/40‑2 (MQ=255)

                                    cGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGaat                                                             <  1:2474628/42‑1 (MQ=255)

                                    cGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGaat                                                             <  1:3184121/42‑1 (MQ=255)

                                     ggACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGt                               <  1:351582/71‑1 (MQ=255)

                                     ggACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGt                               >  1:815046/1‑71 (MQ=255)

                                     ggACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGt                               >  1:936761/1‑71 (MQ=255)

                                     ggACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGt                               >  1:1516816/1‑71 (MQ=255)

                                      gacgacGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGtt                              <  1:3766914/71‑1 (MQ=255)

                                      gacgacGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGt                               >  1:1967858/1‑70 (MQ=255)

                                       acgacgGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTa                             <  1:4912529/71‑1 (MQ=255)

                                         gacgGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTc                                       >  1:4076938/1‑59 (MQ=255)

                                         gacgGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTAcc                           <  1:330598/71‑1 (MQ=255)

                                         gacgGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTAcc                           <  1:5573537/71‑1 (MQ=255)

                                           ataaaCTTGTCGTAGTACGCCAGAACTGCCGGGATAATCCCGCACGCACCGTTAGTCGgcgc                                  <  1:4072408/59‑1 (MQ=255)

                                            taaaCTTGTCGTAGTACGCCAGAACTGCCGGGATAATCCCGCACGCACCGTTAGTCGgcgc                                  >  1:2086508/3‑61 (MQ=255)

                                            taaaCTTGTCGTAGTACGCCAGAACTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTCAc                            >  1:5127424/3‑67 (MQ=255)

                                            ggAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCAcg                        >  1:164002/1‑71 (MQ=255)

                                              aaCTTGTCGTAGTACGCCAGAACTGCCGGGATAATCCCGCACGCACCGTTAg                                         <  1:753465/52‑1 (MQ=25)

                                              aaCTTGTCGTAGTACGCCAGAACTGCCGGGATAATCCCGCACGCACCGTTAg                                         <  1:2729027/52‑1 (MQ=25)

                                              aaCTTATCGTAATAAGCCAGTACTGCCGGAATAATg                                                         <  1:4587790/36‑1 (MQ=255)

                                               aCTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACgcgc                     >  1:4843/1‑71 (MQ=255)

                                               aCTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACgcgc                     >  1:3708721/1‑71 (MQ=255)

                                                cTTGTCGTAGTACGCCAGAACTGCCGGGATAATCCCGCACGCACCGTTAGTc                                       >  1:1343018/1‑52 (MQ=25)

                                                cTTGTCGTAGTACGCCAGAACTGCCGGGATAATCCCGCACGCACCGATAgac                                       >  1:1826969/1‑50 (MQ=16)

                                                cTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAg                                         <  1:1925721/50‑1 (MQ=255)

                                                cTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAg                                         <  1:4392894/50‑1 (MQ=255)

                                                 ttgtCGTAGTACGCCAGAACTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTCAc                            <  1:3702030/59‑1 (MQ=255)

                                                 ttgtCGTAGTACGCCAGAACTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTCACCACGcgacc                   >  1:5468454/4‑71 (MQ=255)

                                                 ttATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTc                                       <  1:3882923/51‑1 (MQ=255)

                                                 ttATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTc                                       <  1:1120625/51‑1 (MQ=255)

                                                  tATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTa                                          >  1:5032431/1‑47 (MQ=255)

                                                  tATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTa                                          >  1:1510676/1‑47 (MQ=255)

                                                  tATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTa                                          >  1:136755/1‑47 (MQ=255)

                                                  tATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTa                                          >  1:125019/1‑47 (MQ=255)

                                                  tATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCAcgc                       >  1:3890246/1‑66 (MQ=255)

                                                  tATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCAcgc                       >  1:217407/1‑66 (MQ=255)

                                                  tATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCAcgc                       >  1:5502257/1‑66 (MQ=255)

                                                   gtCGTAGTACGCCAGAACTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTCACCa                          >  1:3169996/2‑62 (MQ=255)

                                                   gtCGTAGTACGCCAGAACTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTCACCa                          >  1:3176502/2‑62 (MQ=255)

                                                   aTCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCgg                                     >  1:3465369/1‑51 (MQ=255)

                                                   aTCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGcgcccgc                 >  1:4637198/1‑71 (MQ=255)

                                                   aTCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGcgcccgc                 >  1:2676583/1‑71 (MQ=255)

                                                   aTCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGcgcccgc                 >  1:3695590/1‑71 (MQ=255)

                                                   aTCGTAATAAGACAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTa                             >  1:168245/1‑59 (MQ=255)

                                                    tCGTAGTACGCCAGAACTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTcaccac                         >  1:1527158/1‑62 (MQ=255)

                                                    tCGTAGTACGCCAGAACTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTcaccac                         >  1:3411282/1‑62 (MQ=255)

                                                      gtagtaCGCCAGAACTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTCACCACGCGACCGCCAg              >  1:806680/1‑71 (MQ=255)

                                                         ataaTCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGt                                        <  1:790571/42‑1 (MQ=255)

                                                         ataaGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGt                                        <  1:2796035/42‑1 (MQ=255)

                                                         ataaGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTg           >  1:4459363/1‑71 (MQ=255)

                                                           aaGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCa                          >  1:3661131/1‑54 (MQ=255)

                                                            aGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCg         <  1:659736/70‑1 (MQ=255)

                                                            aGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGt        >  1:2956915/1‑71 (MQ=255)

                                                             gCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGtt       <  1:5166382/71‑1 (MQ=255)

                                                             gCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGtt       <  1:465832/71‑1 (MQ=255)

                                                              ccAGAACTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTCACCACGCGACCGCCAGCAGCGttc      <  1:2552742/71‑2 (MQ=255)

                                                              ccAGAACTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTCACCACGCGACCGCCAGCAGCGttc      >  1:2620699/1‑70 (MQ=255)

                                                              ccAGAACTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTCACCACGCGACCGCCAGCAGCGttc      >  1:2181308/1‑70 (MQ=255)

                                                               cAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGtttt     <  1:2560447/71‑1 (MQ=255)

                                                               cAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGtttt     >  1:1294979/1‑71 (MQ=255)

                                                               cAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGtttt     <  1:182662/71‑1 (MQ=255)

                                                                aGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTttc    >  1:58555/1‑71 (MQ=255)

                                                                aGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTttc    >  1:3581412/1‑71 (MQ=255)

                                                                  tACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTttctt  <  1:4390054/71‑1 (MQ=255)

                                                                  tACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTttctt  <  1:3112079/71‑1 (MQ=255)

GGCCAGAAAATAGCGGGCAATTGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTT  >  minE/2093898‑2094034

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