BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,093,958	A→C	16.7%	A298A (GCT→GCG)	tdcG ←	L‑serine dehydratase 3

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,093,958	0	A	C	16.7%	182.5 / 40.8	96	A298A (GCT→GCG)	tdcG	L‑serine dehydratase 3
Reads supporting (aligned to +/- strand): ref base A (38/41); new base C (12/4); total (50/46)
Fisher's exact test for biased strand distribution p-value = 5.88e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.96e-01

Only 100 of 102 total aligned reads displayed.

CCCGCGGCCAGAAAATAGCGGGCAATTGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGT  >  minE/2093893‑2094028

cccGCGGCCAGAAAATAGCGGGCAATTGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAg                                                                   <  1:5310707/71‑1 (MQ=255)

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

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

        cAGAAAATAGCGGGCAATTGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGc                                                                      <  1:868925/60‑1 (MQ=255)

        cAGAAAATAGCGGGCAATTGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGc                                                                      <  1:2698739/60‑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:3536218/71‑1 (MQ=255)

                        aTTGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCcgcacgc                                           <  1:3938239/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:2532291/71‑1 (MQ=255)

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

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

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

                               gCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCa                                                                    >  1:3419113/1‑39 (MQ=255)

                               gCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCa                                                                    >  1:563282/1‑39 (MQ=255)

                               gCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAg                                                                   <  1:243809/40‑1 (MQ=255)

                               gCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAg                                                                   <  1:1302585/40‑1 (MQ=255)

                               gCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAg                                                                   <  1:3204843/40‑1 (MQ=255)

                               gCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAg                                                                   <  1:2404873/40‑1 (MQ=255)

                               gCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAg                                                                   <  1:5287495/40‑1 (MQ=255)

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

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

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

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

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

                                 tCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGt                                  >  1:347994/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:3184121/42‑1 (MQ=255)

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

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

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

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

                                          ggACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGt                         <  1:351582/71‑1 (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:376702/1‑59 (MQ=255)

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

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

                                              gacgGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTAcc                     <  1:330598/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:2729027/52‑1 (MQ=25)

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

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

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

                                                    aCTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACgcgc               >  1:4843/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:2066948/4‑71 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

                                                        aTCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTa                       >  1:4688323/1‑59 (MQ=255)

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

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

                                                        aTCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGcgcccgc           >  1:4637198/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)

CCCGCGGCCAGAAAATAGCGGGCAATTGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGT  >  minE/2093893‑2094028

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