BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	2,093,973	A→G	18.8%	I293I (ATT→ATC)	tdcG ←	L‑serine dehydratase 3

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	2,093,973	0	A	G	18.8%	212.5 / 40.2	101	I293I (ATT→ATC)	tdcG	L‑serine dehydratase 3
Reads supporting (aligned to +/- strand): ref base A (46/36); new base G (14/5); total (60/41)
Fisher's exact test for biased strand distribution p-value = 2.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 108 total aligned reads displayed.

GAAAATAGCGGGCAATTGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAAC  >  Exported/2093903‑2094041

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

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

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

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

       gCGGGCAATTGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATg                                                                <  1:1789480/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)

                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)

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

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

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

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

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

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

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

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

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

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

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

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

                                ggACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGt                                      >  1:1516816/1‑71 (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: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:2066948/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:5502257/1‑66 (MQ=255)

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

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

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

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

                                              aTCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCgg                                            >  1:3465369/1‑51 (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)

                                               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)

                                                      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:2552742/71‑2 (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:182662/71‑1 (MQ=255)

                                                          cAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGtttt            >  1:1294979/1‑71 (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:3668176/1‑41 (MQ=255)

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

                                                               cTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGtt                                     >  1:3661615/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)

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

GAAAATAGCGGGCAATTGACCGCTCGTTTACCGGACGACGGAACTTATCGTAATAAGCCAGTACTGCCGGAATAATGCCGCACGCACCGTTAGTCGGTGCCGTTACCACGCGCCCGCCAGCTGCGTTTTCTTCACTAAC  >  Exported/2093903‑2094041

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