BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	2,093,979	G→C	19.4%	G291G (GGC→GGG)	tdcG ←	L‑serine dehydratase 3

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	2,093,979	0	G	C	19.4%	219.2 / 44.2	103	G291G (GGC→GGG)	tdcG	L‑serine dehydratase 3
Reads supporting (aligned to +/- strand): ref base G (47/36); new base C (14/6); total (61/42)
Fisher's exact test for biased strand distribution p-value = 3.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 107 total aligned reads displayed.

GCGGGCAATTGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGC  >  Exported/2093910‑2094044

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

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

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

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

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

       aTTGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCcgcacg                                                            >  1:2563832/1‑70 (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:5650663/71‑1 (MQ=255)

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

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

              gCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGt                                                      >  1:221553/1‑69 (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)

                         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: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: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)

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

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

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

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

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

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

                                      tATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTa                                                    >  1:136755/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:3890246/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)

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

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

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

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

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

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

                                        tCGTAGTACGCCAGAACTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTcaccac                                   >  1:1527158/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)

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

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

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

                                                 gCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCg                   <  1:3037432/69‑1 (MQ=255)

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

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

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

                                                  ccAGAACTGCCGGGATAATCCCGCACGCACCGTTAGTCGGCGCAGTCACCACGCGACCGCCAGCAGCGttc                >  1:2620699/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:3581412/1‑71 (MQ=255)

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

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

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

                                                        cTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGtt                                        >  1:3664791/1‑41 (MQ=255)

                                                        cTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGtt                                        >  1:3661615/1‑41 (MQ=255)

                                                        cTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGtt                                        >  1:3668176/1‑41 (MQ=255)

                                                        cTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTttctt            >  1:1208784/1‑69 (MQ=255)

                                                         tGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTa                                       <  1:556242/41‑1 (MQ=255)

                                                         tGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTa                                       <  1:3043430/41‑1 (MQ=255)

                                                          gCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGtt                 <  1:1173057/62‑1 (MQ=255)

                                                          gCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCAc         >  1:3455443/1‑70 (MQ=255)

                                                           ccGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTa                                       <  1:2388891/39‑1 (MQ=255)

                                                            cGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTaa      <  1:384785/71‑1 (MQ=255)

                                                             ggAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAAc     >  1:3191954/1‑71 (MQ=255)

                                                              ggataatCCCGCACGCACCGTTAGTCGGCGCAGTCACCACGCGACCGCCAGCAGCGtt                 <  1:5387372/56‑1 (MQ=255)

                                                              gAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTaa      >  1:4677484/1‑69 (MQ=255)

                                                                ataatGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGc  >  1:4308065/1‑71 (MQ=255)

                                                                ataatGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGc  >  1:5546509/1‑71 (MQ=255)

                                                                  aatGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTttct             <  1:604785/58‑1 (MQ=255)

                                                                  aatGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTttct             <  1:4361282/58‑1 (MQ=255)

                                                                  aatGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTttct             <  1:2004588/58‑1 (MQ=255)

GCGGGCAATTGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAACCGC  >  Exported/2093910‑2094044

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