BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,093,718	C→A	15.0%	G378G (GGG→GGT)	tdcG ←	L‑serine dehydratase 3

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,093,718	0	C	A	15.0%	283.2 / 22.9	100	G378G (GGG→GGT)	tdcG	L‑serine dehydratase 3
Reads supporting (aligned to +/- strand): ref base C (37/47); new base A (12/3); total (50/50)
Fisher's exact test for biased strand distribution p-value = 1.22e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 103 total aligned reads displayed.

TCACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTT  >  minE/2093648‑2093785

tCACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGc                                                                     >  1:5304007/1‑71 (MQ=255)

    ggCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCaa                                                                 <  1:1316666/71‑1 (MQ=255)

    ggCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCaa                                                                 <  1:4238672/71‑1 (MQ=255)

      cATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAgg                                                               >  1:1290078/1‑71 (MQ=255)

       aTTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGt                                                              <  1:3518293/71‑1 (MQ=255)

          aaTGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTAt                                                           >  1:62536/1‑71 (MQ=255)

           aTGGCGTTACGCTCGATGCATGGCACCTGTACCTGTCCGGCGACCGGGTCACACGTCAGaccg                                                                  <  1:4700496/63‑2 (MQ=255)

            tGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCa                                                                  <  1:2850388/62‑1 (MQ=255)

            tGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCa                                                                  <  1:4082327/62‑1 (MQ=255)

             ggCGTTACGCTCGATGCATGGCACCTGTACCTGTCCGGCGACCGGGTCACACGTCAGACCGAGGTTGTGCt                                                        >  1:72115/1‑71 (MQ=255)

             ggCGTTACGCTCGATGCATGGCACCTGTACCTGTCCGGCGACCGGGTCACACGTCAGACCGAGGTTGTGCt                                                        >  1:5369122/1‑71 (MQ=255)

             ggCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCt                                                        >  1:2510150/1‑71 (MQ=255)

                 ttACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAgg                                                               <  1:5595479/60‑1 (MQ=255)

                 ttACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAgg                                                               <  1:1899907/60‑1 (MQ=255)

                 ttACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAgg                                                               <  1:1350301/60‑1 (MQ=255)

                 ttACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAgg                                                               <  1:3464394/60‑1 (MQ=255)

                 ttACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTAt                                                           <  1:4430465/64‑1 (MQ=255)

                 ttACGCTCGATGCATGGCACCTGTACCTGTCCGGCGACCGGGTCACACGTCAGACCGAgg                                                               <  1:5253238/60‑1 (MQ=255)

                 ttACGCTCGATGCATGGCACCTGTACCTGTCCGGCGACCGGGTCACACGTCAGACCGAGGTTGTGCTCCAt                                                    <  1:2779624/71‑1 (MQ=255)

                   aCGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGccc                                                                   >  1:4994782/1‑54 (MQ=255)

                   aCGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATcg                                                  >  1:5653179/1‑71 (MQ=255)

                   aCGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATcg                                                  >  1:114394/1‑71 (MQ=255)

                   aCGCTCGATGCATGGCACCTGTACCTGTCCGGCGACCGGGTCACACGTCAGACCGAGGtt                                                             >  1:605005/1‑60 (MQ=255)

                   aCGCTCGATGCATGGCACCTGTACCTGTCCGGCGACCGGGTCACACGTCAGACCGAGGTTGTGCt                                                        >  1:346/1‑65 (MQ=255)

                   aCGCTCGATGCATGGCACCTGTACCTGTCCGGCGACCGGGTCACACGTCAGACCGAGGTTGTGCTCCAtgg                                                  >  1:1677202/1‑69 (MQ=255)

                    cGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATcgc                                                 <  1:3410395/71‑1 (MQ=255)

                    cGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATcgc                                                 <  1:1106021/71‑1 (MQ=255)

                        cAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTa                                                            >  1:3693864/1‑56 (MQ=255)

                        cAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTa                                                            >  1:4494426/1‑56 (MQ=255)

                        cAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTa                                                            >  1:4516193/1‑56 (MQ=255)

                        cAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGAtt                                             >  1:3590660/1‑71 (MQ=255)

                        cAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGAtt                                             >  1:3584056/1‑71 (MQ=255)

                        cAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGAtt                                             >  1:1451747/1‑71 (MQ=255)

                        cAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGAt                                              <  1:4537438/70‑1 (MQ=255)

                         aaTGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGAttt                                            <  1:1529390/71‑1 (MQ=255)

                           tGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGt                                                              <  1:226556/51‑1 (MQ=255)

                           tGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGt                                                              <  1:4681619/51‑1 (MQ=255)

                           tGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGt                                                              <  1:4015058/51‑1 (MQ=255)

                            gCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTAt                                                           <  1:236864/53‑1 (MQ=255)

                             cACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTccgc                                        >  1:1744041/1‑71 (MQ=255)

                              aCGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGc                                                         >  1:216264/1‑53 (MQ=255)

                               cGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTccgccg                                      <  1:5343954/71‑1 (MQ=255)

                               cGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTccgccg                                      <  1:381/71‑1 (MQ=255)

                                gggATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTcc                                          >  1:1219149/1‑66 (MQ=255)

                                  gATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTAt                                                           <  1:4799612/47‑1 (MQ=255)

                                   aTTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTa                                                            <  1:686516/45‑1 (MQ=255)

                                   aTTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTa                                                            >  1:2327529/1‑45 (MQ=255)

                                    tttGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCGAAGGTTATGCTCCATCGCGATTTCCGCCGCATTg                                 >  1:1978175/1‑71 (MQ=255)

                                    tttGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTAt                                                           <  1:3246839/45‑1 (MQ=255)

                                    tttGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTAt                                                           <  1:3523781/45‑1 (MQ=255)

                                    tttGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTcc                                                      <  1:5240538/50‑1 (MQ=255)

                                    tttGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTcc                                                      <  1:4156908/50‑1 (MQ=255)

                                    tttGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTcc                                                      <  1:20459/50‑1 (MQ=255)

                                     ttGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTg                                 >  1:5319644/1‑70 (MQ=255)

                                      tGTACCTGTCCGGCGACCGGGTCACACGTCAGACCGAGGTTGTGCTCCa                                                     >  1:3967197/1‑49 (MQ=255)

                                      tGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCAt                                                    >  1:2946154/1‑50 (MQ=255)

                                       gTACCTGTCCGGCGACCGGGTCACACGTCAGACCGAGGTtgtg                                                          >  1:5605557/1‑43 (MQ=255)

                                       gTACCTGTCCGGCGACCGGGTCACACGTCAGACCGAGGTTGTGCTCCATGGCGATTTCCGCCGCGATGCa                               >  1:2198475/1‑70 (MQ=12)

                                       gTACCTGTCCGGCGACCGGGTCACACGTCAGACCGAGGTTGTGCTCCATGGCGATTTCCGCCGCGATGCa                               >  1:2405313/1‑70 (MQ=12)

                                       gTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCa                               <  1:1369682/70‑1 (MQ=255)

                                          ccTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATAcc                           >  1:885311/1‑71 (MQ=255)

                                           cTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGTGATTTCCGCCGCATTGCAt                              <  1:774486/67‑1 (MQ=255)

                                           cTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCAt                                   >  1:3664762/1‑62 (MQ=255)

                                           cTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCAt                                   >  1:1351605/1‑62 (MQ=255)

                                           cTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCAt                              <  1:2939965/67‑1 (MQ=255)

                                             gTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATAc                            >  1:810687/1‑67 (MQ=255)

                                              tCCGGCGACCGGGTCACACGTCAGACCGAGGTTGTGCTCCATGGCGATTTCCGCCGCGATGCACACCTgcg                       >  1:2494244/1‑71 (MQ=255)

                                              tCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGAttt                                            >  1:4332631/1‑50 (MQ=255)

                                              tCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGAttt                                            >  1:206288/1‑50 (MQ=255)

                                              tCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGAttt                                            >  1:2369391/1‑50 (MQ=255)

                                              tCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTg                                 <  1:1450476/61‑1 (MQ=255)

                                              tCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTg                                 <  1:3523794/61‑1 (MQ=255)

                                              tCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTg                                 <  1:1704363/61‑1 (MQ=255)

                                              tCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTg                                 <  1:5115985/61‑1 (MQ=255)

                                              tCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCAt                              <  1:1091913/64‑1 (MQ=255)

                                              tCCGCCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCAt                              <  1:5143837/64‑1 (MQ=255)

                                               ccGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTccgc                                        >  1:2417862/1‑53 (MQ=255)

                                               ccGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTccgc                                        >  1:3493960/1‑53 (MQ=255)

                                               ccGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTccgc                                        >  1:5102481/1‑53 (MQ=255)

                                               ccGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTccgc                                        >  1:5033261/1‑53 (MQ=255)

                                               ccGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTccgc                                        >  1:3931271/1‑53 (MQ=255)

                                                cGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGAttt                                            >  1:1296817/1‑48 (MQ=255)

                                                cGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCa                               >  1:4030355/1‑61 (MQ=255)

                                                cGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGcc                     >  1:5285012/1‑71 (MQ=255)

                                                cGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGcc                     >  1:1777827/1‑71 (MQ=255)

                                                  gcgaCCGGGTCACACGTCAGACCGAGGTTGTGCTCCATGGCGATTTCCGCCGCGATGCACACCTGCGCCgg                   >  1:3939444/4‑71 (MQ=9)

                                                  gcgaCCGGGTCACACGTCAGACCGAGGTTGTGCTCCATGGCGATTTCCGCCGCGATGCACACCTGCGCCgg                   <  1:5587325/68‑1 (MQ=9)

                                                        ggATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTg                         >  1:3450218/1‑59 (MQ=255)

                                                          aTCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTgccgc           <  1:3127635/71‑1 (MQ=255)

                                                             gCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTcc                                          >  1:2162419/1‑37 (MQ=255)

                                                             gCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCgg                   <  1:1952331/60‑1 (MQ=255)

                                                             gCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCgg                   <  1:2574925/60‑1 (MQ=255)

                                                             gCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCgg                   <  1:3908258/60‑1 (MQ=255)

                                                             gCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCgg                   <  1:4934977/60‑1 (MQ=255)

                                                             gCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCa        <  1:3822503/71‑1 (MQ=255)

                                                              cAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCag       <  1:3440753/71‑1 (MQ=255)

                                                                  tCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCagtagt   <  1:3609381/71‑1 (MQ=255)

                                                                  gcAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTg                                 <  1:5623993/40‑1 (MQ=255)

                                                                    aGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGtt  <  1:577449/70‑1 (MQ=255)

                                                                     gCCCAAGGTTATGCTCCATCGCGATTTCCGCCGGATAGCATAc                            >  1:322804/1‑43 (MQ=38)

TCACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAGTAGTT  >  minE/2093648‑2093785

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 14 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`