BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	2,094,810	G→A	10.5%	P14P (CCC→CCT)	tdcG ←	L‑serine dehydratase 3

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	2,094,810	0	G	A	10.5%	313.7 / 19.2	114	P14P (CCC→CCT)	tdcG	L‑serine dehydratase 3
Reads supporting (aligned to +/- strand): ref base G (43/59); new base A (8/4); total (51/63)
Fisher's exact test for biased strand distribution p-value = 1.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 117 total aligned reads displayed.

CTTTCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCTCGGATTGCCGTTCAGTGA  >  Exported/2094748‑2094880

cTTTCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCaa                                                                 <  1:4652208/70‑1 (MQ=255)

cTTTCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAAt                                                                >  1:5313544/1‑71 (MQ=255)

  ttCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATcc                                                              >  1:5114875/1‑71 (MQ=255)

   tCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATcc                                                              <  1:1358814/70‑1 (MQ=255)

   tCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATcc                                                              >  1:5202564/1‑70 (MQ=255)

   tCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCACACGGTATGCGAACTGGAGGGACCAATcc                                                              >  1:491368/1‑70 (MQ=255)

            tCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGaa                                                    <  1:4121931/71‑1 (MQ=255)

              aaTAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGaaaa                                                  <  1:4592910/71‑1 (MQ=255)

                  aaaCTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATc                                              <  1:5076281/71‑1 (MQ=255)

                   aaCTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCg                                             <  1:3942990/71‑1 (MQ=255)

                   aaCTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCg                                             <  1:1765634/71‑1 (MQ=255)

                   aaCTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCg                                             <  1:1802592/71‑1 (MQ=255)

                    aCTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCg                                             <  1:4506995/70‑1 (MQ=255)

                    aCTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGa                                            <  1:218863/71‑1 (MQ=255)

                      tttttCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAAt                                                                >  1:4217547/1‑49 (MQ=255)

                        tttCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAAt                                                          >  1:3459377/1‑53 (MQ=255)

                        tttCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAAt                                                          >  1:5331385/1‑53 (MQ=255)

                            cGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACt                                     <  1:3990211/70‑1 (MQ=255)

                            cGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACt                                     <  1:4448373/70‑1 (MQ=255)

                             ggCATTCATTGGCCCCACGGTATGCGAACTGGAGGGAc                                                                    >  1:249946/1‑38 (MQ=255)

                             ggCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTaa                                   <  1:3967628/71‑1 (MQ=255)

                              gctttcattGGTCCAACGGTATGAGAACTGGAAGGGCCAATGCCGATTTTGAAAATATCGAATACGCTAAt                                  >  1:3441019/4‑71 (MQ=255)

                               ctttcattGGTCCAACGGTATGAGAACTGGAAGGGCCAATGCCGATTTTGAAAATATCTAAta                                         >  1:569464/3‑62 (MQ=255)

                               ctttcattGGTCCAACGGTATGAGAACTGGAAGGGCCAATGCCGATTTTGAAAATATCGAAta                                         >  1:4994076/3‑62 (MQ=255)

                                tttcattGGTCCAACGGTATGAGAACTGGAAGGGCCAATGCCGATTTTGAAAATATCGAAta                                         >  1:2722853/2‑61 (MQ=255)

                                attcattGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCa                                >  1:1953523/1‑71 (MQ=255)

                                 ttcattGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGaa                                                    <  1:3304151/50‑1 (MQ=255)

                                 ttcattGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGaa                                                    <  1:2567559/50‑1 (MQ=255)

                                 ttcattGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCg                                             <  1:4444206/57‑1 (MQ=255)

                                 ttcattGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCg                                             <  1:1480598/57‑1 (MQ=255)

                                 ttcattGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATg                                         <  1:2780132/61‑1 (MQ=255)

                                 ttcattGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATg                                         <  1:4967718/61‑1 (MQ=255)

                                 ttcattGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATg                                         <  1:4680636/61‑1 (MQ=255)

                                 ttcattGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATg                                         <  1:4988320/61‑1 (MQ=255)

                                 ttcattGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATg                                         <  1:5178284/61‑1 (MQ=255)

                                 ttcattGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATg                                         <  1:4937353/61‑1 (MQ=255)

                                 ttcattGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCat                               <  1:1549673/71‑1 (MQ=255)

                                  tcattGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCata                              >  1:185753/1‑71 (MQ=255)

                                  tcattGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCata                              >  1:4600281/1‑71 (MQ=255)

                                   cattGGTCCAACGGTATGAGAACTGGAAGGGCCAATGCCGATTTTGaaaa                                                  >  1:2021963/1‑50 (MQ=255)

                                    attGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCatat                             <  1:3673656/70‑1 (MQ=255)

                                    attGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATc                            <  1:2866449/71‑1 (MQ=255)

                                    attGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATc                            <  1:3152344/71‑1 (MQ=255)

                                     ttGGTCCAACGGTATGAGAACTGGAAGGGCCAATGCCGATTTTGAAAATATCGAATACGCTAATCatagga                           <  1:1323247/71‑4 (MQ=255)

                                      tGGTCCAACGGTATGAGAACTGGAAGGGCCAATGCCGATTTTGAAAATATCGAATACGCTAATCatagga                           <  1:3516884/70‑4 (MQ=255)

                                       ggCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAatat                                               >  1:335444/1‑49 (MQ=255)

                                       ggCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAat                                                 >  1:3137799/1‑47 (MQ=255)

                                       ggCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCca                          <  1:1537781/70‑1 (MQ=255)

                                           ccACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCCACAccc                     >  1:5661671/1‑71 (MQ=255)

                                           ccACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCCACAcc                      >  1:4446386/1‑70 (MQ=255)

                                            cACGGTATGCGAACTGGAGGGACCAATCCCAATTTTg                                                      <  1:5688508/37‑1 (MQ=255)

                                            cACGGTATGCGAACTGGAGGGACCAATCCCAATTTTg                                                      <  1:4752548/37‑1 (MQ=255)

                                            cACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAat                                                 >  1:5219107/1‑42 (MQ=255)

                                            cACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCt                    >  1:3152767/1‑71 (MQ=255)

                                             aCGGTATGCGAACTGGAGGGACCAATCCCAATTTTGaaaa                                                  >  1:745586/1‑40 (MQ=255)

                                             aCGGTATGCGAACTGGAGGGACCAATCCCAATTTTGaaaa                                                  >  1:770185/1‑40 (MQ=255)

                                             aCGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATg                                         >  1:4065286/1‑49 (MQ=255)

                                             aCGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCtg                   >  1:5020014/1‑70 (MQ=255)

                                              cGGTATGCGAACTGGAGGGACCAATCCCAATTTTGaa                                                    <  1:4042463/37‑1 (MQ=255)

                                              cGGTATGCGAACTGGAGGGACCAATCCCAATTTTGaa                                                    <  1:4805529/37‑1 (MQ=255)

                                              cGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAatat                                               >  1:2013054/1‑42 (MQ=255)

                                              cGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATc                            >  1:3286851/1‑61 (MQ=255)

                                              cGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCTCg                  <  1:2278327/71‑1 (MQ=255)

                                              cGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCTCg                  <  1:4705420/71‑1 (MQ=255)

                                              cGGTATGAGAACTGGAAGGGCCAATGCCGATTTTGAAAat                                                 <  1:5658574/40‑1 (MQ=16)

                                               ggTATGAGAACTGGAAGGGCCAATGCCGATTTTGaaaa                                                  >  1:3166002/1‑38 (MQ=16)

                                                gTATGCGAACTGGAGGGACCAATCCCAATTTTGaaaa                                                  >  1:663403/1‑37 (MQ=255)

                                                gTATGCGAACTGGAGGGACCAATCCCAATTTTGaaaa                                                  >  1:1508416/1‑37 (MQ=255)

                                                gTATGCGAACTGGAGGGACCAATCCCAATTTTGaaaa                                                  >  1:4926920/1‑37 (MQ=255)

                                                gTATGCGAACTGGAGGGACCAATCCCAATTTTGaaaa                                                  >  1:5399223/1‑37 (MQ=255)

                                                gTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAat                                                 >  1:1046720/1‑38 (MQ=255)

                                                gTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCTCGGa                <  1:2819036/71‑1 (MQ=255)

                                                gTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCTCGGa                <  1:4223379/71‑1 (MQ=255)

                                                 tATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCTCGGAt               >  1:2882038/1‑71 (MQ=255)

                                                   tGCGAACTGGAGGGACCAATCCCAATTTTGAAAat                                                 <  1:3104218/35‑1 (MQ=255)

                                                   tGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGaa                                           >  1:834518/1‑41 (MQ=255)

                                                    gCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATc                            <  1:4857633/55‑1 (MQ=255)

                                                     cGAACTGGAGGGACCAATCCCAATTTTGAAAATATCg                                             <  1:3115342/37‑1 (MQ=255)

                                                     cGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCTCGGATTGcc           >  1:2376258/1‑71 (MQ=255)

                                                      gAACTGGAGGGACCAATCCCAATTTTGAAAATATCGaa                                           >  1:927602/1‑38 (MQ=255)

                                                      gAACTGGAAGGGCCAATGCCGATTTTGAAAATATCGAATACGCTAATCa                                >  1:5605382/1‑49 (MQ=19)

                                                       aaCTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCTCg                  <  1:4262326/62‑1 (MQ=255)

                                                       aaCTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCTCg                  <  1:217615/62‑1 (MQ=255)

                                                       aaCTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCTCGGATTGCCGt         >  1:468979/1‑71 (MQ=255)

                                                       aaCTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCTCGGATTGCCGt         >  1:2702154/1‑71 (MQ=255)

                                                       aaCTGGAAGGGCCAATGCCGATTTTGAAAATATCGAATACGCTaa                                   <  1:4230516/45‑1 (MQ=16)

                                                        aCTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCTCGGATTGCCGtt        <  1:531360/71‑1 (MQ=255)

                                                         cTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAAt                                  <  1:1651229/44‑1 (MQ=255)

                                                         cTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAAt                                  <  1:236812/44‑1 (MQ=255)

                                                         cTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCTCGGATTGCCGtt        >  1:4034579/1‑70 (MQ=255)

                                                          tGGAGGGACCAATCCCAATTTTGAAAATATCGAATGc                                        <  1:1646663/37‑1 (MQ=255)

                                                          tGGAGGGACCAATCCCAATTTTGAAAATATCGAATGc                                        <  1:1482948/37‑1 (MQ=255)

                                                           ggAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCTCGGATTGCCGTTCAg     >  1:2716294/1‑71 (MQ=255)

                                                            gAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATc                                 >  1:1207226/1‑42 (MQ=255)

                                                            gAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCTCGGATTGCCGTTCAg     <  1:4689144/70‑1 (MQ=255)

                                                            gAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCTCGGATTGCCGTTCAg     <  1:180410/70‑1 (MQ=255)

                                                            gAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCTCGGATTGCCGTTCAGt    >  1:1683090/1‑71 (MQ=255)

                                                              gggACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCTCGGa                >  1:349613/1‑57 (MQ=255)

                                                              gggACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCTCGGATTGCCGTTCAGTg   >  1:961928/1‑70 (MQ=255)

                                                              gggACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCTCGGATTGCCGTTCAGTGa  <  1:2657072/71‑1 (MQ=255)

CTTTCCAGCCGATCAATAAAACTTTTTCCGGCATTCATTGGCCCCACGGTATGCGAACTGGAGGGACCAATCCCAATTTTGAAAATATCGAATGCACTAATCATATCCACACCCTCGGATTGCCGTTCAGTGA  >  Exported/2094748‑2094880

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