BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,093,712	G→C	8.1%	T380T (ACC→ACG)	tdcG ←	L‑serine dehydratase 3

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,093,712	0	G	C	8.1%	292.7 / 12.3	98	T380T (ACC→ACG)	tdcG	L‑serine dehydratase 3
Reads supporting (aligned to +/- strand): ref base G (52/38); new base C (5/3); total (57/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

Only 100 of 102 total aligned reads displayed.

CTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAG  >  minE/2093642‑2093780

cTGCTTTCAAGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAg                                                                      <  1:2499662/71‑1 (MQ=255)

  gCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGt                                                                    >  1:1794691/1‑71 (MQ=255)

  gCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGt                                                                    >  1:489675/1‑71 (MQ=255)

   cTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAgg                                                                     <  1:451817/69‑1 (MQ=255)

        aCGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGccc                                                              <  1:425164/71‑1 (MQ=255)

          ggCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGccc                                                              >  1:3905217/1‑69 (MQ=255)

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

           gCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAg                                                           >  1:1202718/1‑71 (MQ=255)

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

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

             aTTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGt                                                         >  1:4134366/1‑71 (MQ=255)

              ttAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGtt                                                        <  1:2191213/71‑1 (MQ=255)

              ttAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGt                                                         <  1:2911664/70‑1 (MQ=255)

                aaTGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTa                                                       <  1:3972130/70‑1 (MQ=255)

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

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

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

                 aTGGCGTTACGCTCGATGCATGGCACCTGTACCTGTCCGGCGACCGGGTCACACGTCAg                                                                 <  1:2715810/59‑1 (MQ=255)

                 aTGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATg                                                     >  1:4009992/1‑71 (MQ=255)

                  tGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAg                                                                 <  1:2023941/58‑1 (MQ=255)

                   ggCGTTACGCTCGATGCATGGCACCTGTACCTGTCCGGCGACCGGGTCACACGTCAgac                                                               >  1:640222/1‑57 (MQ=255)

                    gcgttACGCTCGATGCATGGCACCTGTACCTGTCCGGCGACCGGGTCACACGTCAGACCGAGGTTGTGCTc                                                  >  1:69646/4‑71 (MQ=255)

                    gCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTc                                                  <  1:2627347/71‑1 (MQ=255)

                     cATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTc                                                  >  1:1579481/1‑70 (MQ=255)

                       ttACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCt                                                   <  1:2258390/67‑1 (MQ=255)

                         aCGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAg                                                                      >  1:173999/1‑46 (MQ=255)

                         aCGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAg                                                                      >  1:57597/1‑46 (MQ=255)

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

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

                          cGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCaa                                                            >  1:3920832/1‑55 (MQ=255)

                          cGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATcg                                             >  1:1117471/1‑70 (MQ=255)

                            ttCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATg                                                     >  1:858876/1‑60 (MQ=255)

                            ttCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATg                                                     >  1:2011541/1‑60 (MQ=255)

                            ttCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGa                                          <  1:2503273/71‑1 (MQ=255)

                             tCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGAt                                         >  1:2310371/1‑71 (MQ=255)

                              cAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAg                                                           >  1:1413031/1‑52 (MQ=255)

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

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

                                 tGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTcc                                     >  1:2888998/1‑71 (MQ=255)

                                 tGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTcc                                     >  1:2745009/1‑71 (MQ=255)

                                 tGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTcc                                     <  1:1536656/71‑1 (MQ=255)

                                 tGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTc                                      >  1:1670312/1‑70 (MQ=255)

                                   cACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGt                                                         >  1:2691534/1‑49 (MQ=255)

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

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

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

                                   cACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTccg                                    <  1:397361/70‑1 (MQ=255)

                                    aCGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTccgcc                                  >  1:3321422/1‑71 (MQ=255)

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

                                     cGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTccgccg                                 >  1:3809978/1‑71 (MQ=255)

                                       ggATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCAt                                               >  1:2908358/1‑55 (MQ=255)

                                       ggATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCa                               <  1:859819/71‑1 (MQ=255)

                                       ggATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCa                               <  1:1040543/71‑1 (MQ=255)

                                       ggATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCa                               <  1:1789189/71‑1 (MQ=255)

                                       ggATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCa                               >  1:2310296/1‑71 (MQ=255)

                                        gATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCAt                              >  1:2266082/1‑71 (MQ=255)

                                         aTTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGt                                                         <  1:226628/43‑1 (MQ=255)

                                         aTTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGt                                                         <  1:3803092/43‑1 (MQ=255)

                                         aTTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCAtt                             <  1:3079567/71‑1 (MQ=255)

                                         aTTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCAtt                             <  1:2655062/71‑1 (MQ=255)

                                         aTTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCAtt                             <  1:392822/71‑1 (MQ=255)

                                          tttGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCAtt                             <  1:1633957/70‑1 (MQ=255)

                                           ctGTACCTGTCCGGCGACCGGGTCACACGTCAGACCGAGGTTGTGCTCCATGGCGATTTCcgccgc                                <  1:101106/65‑1 (MQ=12)

                                            tGTACCTGTCCGGCGACCGGGTCACACGTCAGACCGAGGTTGTGCTCCATGGCGATTTCcgccgc                                <  1:766565/65‑1 (MQ=12)

                                            tGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAg                                                           >  1:1869080/1‑38 (MQ=255)

                                            tGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTAt                                                      >  1:4073972/1‑43 (MQ=255)

                                            tGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTc                                      >  1:2712877/1‑59 (MQ=255)

                                            tGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCAtt                             >  1:2051032/1‑68 (MQ=255)

                                            tGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCa                          >  1:3514626/1‑71 (MQ=255)

                                             gTACCTGTCCGGCGACCGGGTCACACGTCAGACCGAGGTTGTGCTCCATGGCGATTTCCGCCGCGATGCa                          >  1:1487172/1‑70 (MQ=255)

                                             gTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCa                          >  1:3698576/1‑70 (MQ=255)

                                              tACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCcgccgc                                <  1:3578752/63‑1 (MQ=255)

                                              tACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGc                           <  1:1918185/68‑1 (MQ=255)

                                              tACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGc                           <  1:769223/68‑1 (MQ=255)

                                              tACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGc                           <  1:3268612/68‑1 (MQ=255)

                                              tACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGc                           <  1:3447487/68‑1 (MQ=255)

                                              tACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATa                        <  1:2384791/71‑1 (MQ=255)

                                               aCCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCAtt                             >  1:3439989/1‑65 (MQ=255)

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

                                                 cTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGt                                                         <  1:31938/35‑1 (MQ=255)

                                                   gTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTa                                                       >  1:53049/1‑35 (MQ=255)

                                                   gTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGc                           >  1:1330896/1‑63 (MQ=255)

                                                   gTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGc                           >  1:503312/1‑63 (MQ=255)

                                                    tCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTgcg                  >  1:1264869/1‑71 (MQ=255)

                                                       ggCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTc                                      <  1:225504/48‑1 (MQ=255)

                                                        gCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCa                               <  1:3575496/54‑1 (MQ=255)

                                                         cAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGAGATTTc                                      >  1:4093129/1‑46 (MQ=255)

                                                          aaCCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGAc            >  1:2078979/1‑71 (MQ=255)

                                                           aCCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCAt                         <  1:2454734/57‑1 (MQ=255)

                                                            ccGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATATCTGCGCCGGACTg          >  1:2669204/1‑71 (MQ=255)

                                                            ccGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCAGCCGCATTGCATACCTGCGCCGGACTg          >  1:2219078/1‑71 (MQ=255)

                                                             cGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCAtt                             <  1:4117087/51‑1 (MQ=255)

                                                             cGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGcc                <  1:2440865/64‑1 (MQ=255)

                                                             cGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTgc         >  1:335540/1‑71 (MQ=255)

                                                              ggATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCAtt                             <  1:2891229/50‑1 (MQ=255)

                                                                aTCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGc                           >  1:1735763/1‑50 (MQ=255)

                                                                  cGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCAt                              >  1:829668/1‑45 (MQ=255)

                                                                  cGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCAt                              >  1:2929045/1‑45 (MQ=255)

                                                                    cAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCag  >  1:3170170/1‑71 (MQ=255)

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

CTGCTTTCACGGCATTAATGGCATTACGTTCAATGCACGGGATTTGTACCTGTCCGGCAACCGGATCGCAGGTCAGCCCAAGGTTATGCTCCATCGCGATTTCCGCCGCATTGCATACCTGCGCCGGACTGCCGCCCAG  >  minE/2093642‑2093780

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