BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,470,435	C→T	5.3%	Q429Q (CAG→CAA)	yfaL ←	adhesin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,470,435	0	C	T	5.3%	418.9 / 7.4	132	Q429Q (CAG→CAA)	yfaL	adhesin
Reads supporting (aligned to +/- strand): ref base C (58/67); new base T (4/3); total (62/70)
Fisher's exact test for biased strand distribution p-value = 7.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 143 total aligned reads displayed.

CCAGCACGTAGCTGCCGTTTTGCGCAATGGTGAGCTGACCAGCACCTTCGATGATGCCAGCAAAACTGCCCTGATTAACAGTAACGTTGCCACCAGCATCGATATTTAAAGTGCCATTCTGAAAGCCCGTCAATGCGTG  >  minE/1470365‑1470503

 cAGCACGTAGCTGCCGTTTTGCGCAATGGTGAGCTGACCAGCACCTTCGATGATGCCAGCAAAACTGCCCt                                                                     <  1:3359808/71‑1 (MQ=255)

 cAGCACGTAGCTGCCGTTTTGCGCAATGGTGAGCTGACCAGCACCTTCGATGATGCCAGCAAAACTGCCCt                                                                     <  1:2259726/71‑1 (MQ=255)

     aCGTAGCTGCCGTTTTGCGCAATGGTGAGCTGACCAGCACCTTCGATGATTCCAGCAAAACTGCCCTGAtt                                                                 <  1:3137634/71‑1 (MQ=255)

       gTAGCTGCCGTTTTGCGCAATGGTGAGCTGACCAGCACCTTCGATGATGCCAGCAAAACTGCCCTGATTaa                                                               <  1:951750/71‑1 (MQ=255)

           cTGCCGTTTTGCGCAATGGTGAGCTGACCAGCACCTTCGATGATGCCAGCAAAACTGCCTTGATTAACAGt                                                           <  1:1897969/71‑1 (MQ=255)

              ccGTTTTGCGCAATGGTGAGCTGACCAGCACCTTCGATGATGCCAGCAAAACTGCCCTGATTAACAGTaa                                                         <  1:112533/70‑1 (MQ=255)

              ccGTTTTGCGCAATGGTGAGCTGACCAGCACCTTCGATGATGCCAGCAAAACTGCCCTGATTAACAGTa                                                          <  1:1024228/69‑1 (MQ=255)

              ccGTTTTGCGCAATGGTGAGCTGACCAGCACCTTCGATGATGCCAGCAAAACTGCCCTGATTAACAGTa                                                          <  1:3095340/69‑1 (MQ=255)

              ccGTTTTGCGCAATGGTGAGCTGACCAGCACCTTCGATGATGCCAGCAAAACTGCCCTGATTAACAGTAAc                                                        >  1:3154705/1‑71 (MQ=255)

              ccGTTTTGCGCAATGGTGAGCTGACCAGCACCTTCGATGATGCCAGCAAAACTGCCCTGATTAACAGTAAc                                                        >  1:1394434/1‑71 (MQ=255)

               cGTTTTGCGCAATGGTGAGCTGACCAGCACCTTCGATGATGCCAGCAAAACTGCCCTGATTAACAGTAACg                                                       >  1:1748408/1‑71 (MQ=255)

                gTTTTGCGCAATGGTGAGCTGACCAGCACCTTCGATGATGCCAGCAAAACTGCCCTGATTAACAGTAACGt                                                      >  1:383925/1‑71 (MQ=255)

                 ttttGCGCAATGGTGAGCTGACCAGCACCTTCGATGATGCCAGCAAAACTGCCCTGATTAACAGt                                                           <  1:2430016/65‑1 (MQ=255)

                  tttGCGCAATGGTGAGCTGACCAGCACCTTCGATGATGCCAGCAAAACTGCCTTGATTAACAGTAACGTTg                                                    >  1:1717988/1‑71 (MQ=255)

                  tttGCGCAATGGTGAGCTGACCAGCACCTTCGATGATGCCAGCAAAACTGCCCTGATTAACAGTAACGTTg                                                    >  1:2820316/1‑71 (MQ=255)

                  tttGCGCAATGGTGAGCTGACCAGCACCTTCGATGATGCCAGCAAAACTGCCCTGATTAACAGTAACGTTg                                                    >  1:505484/1‑71 (MQ=255)

                   ttGCGCAATGGTGATCTGACCAGCACCTTCGATGATGCCAGCAAAACTGCCCTGATTAACAGTAACGTTGc                                                   >  1:1698197/1‑71 (MQ=255)

                      cgcAATGGTGAGCTGACCAGCACCTTCGATGATGCCAGCAAAACTGCCCTGATTAACAGTAACGTTGccac                                                <  1:2605455/71‑1 (MQ=255)

                      cgcAATGGTGAGCTGACCAGCACCTTCGATGATGCCAGCAAAACTGCCCTGATTAACAGTAACGTTGccac                                                <  1:3782947/71‑1 (MQ=255)

                      cgcAATGGTGAGCTGACCAGCACCTTCGATGATGCCAGCAAAACTGCCCTGATTAACAGTAACGTTGccac                                                <  1:1846413/71‑1 (MQ=255)

                         aaTGGTGAGCTGACCAGCACCTTCGATGATGCCAGCAAAACTGCCCTGATTAACAGTAACGTTGCCACCAg                                             <  1:80203/71‑1 (MQ=255)

                         aaTGGTGAGCTGACCAGCACCTTCGATGATGCCAGCAAAACTGCCCTGATTAACAGTAACGTTGCCACCAg                                             <  1:3505817/71‑1 (MQ=255)

                           tGGTGAGCTGACCAGCACCTTCGATGATGCCAGCAAAACTGCCCTGAt                                                                  <  1:269950/48‑1 (MQ=255)

                           tGGTGAGCTGACCAGCACCTTCGATGATGCCAGCAAAACTGCCCTGATTAACAGt                                                           <  1:981754/55‑1 (MQ=255)

                           tGGTGAGCTGACCAGCACCTTCGATGATGCCAGCAAAACTGCCCTGATTAACAGTAACGTTGCCACCAGCa                                           >  1:3414080/1‑71 (MQ=255)

                            ggTGAGCTGACCAGCACCTTCGATGATGCCAGCAAAACTGCCCTGATTAACAGTAACGTTGCCACCAGCAt                                          >  1:3384869/1‑71 (MQ=255)

                            ggTGAGCTGACCAGCACCTTCGATGATGCCAGCAAAACTGCCCTGATTAACAGTAACGTTGCCACCAGCAt                                          >  1:2070498/1‑71 (MQ=255)

                             gTGAGCTGACCAGCACCTTCGATGATGCCAGCAAAACTGCCCTGATTAACAGTAACGTTGCCACCAGCATc                                         <  1:1579006/71‑1 (MQ=255)

                              tGAGCTGACCAGCACCTTCGATGATGCCAGCAAAACTGCCCTGATTaa                                                               <  1:4105944/48‑1 (MQ=255)

                              tGAGCTGACCAGCACCTTCGATGATGCCAGCAAAACTGCCCTGATTAACAGTAACGTTGCCACCAGCATCg                                        <  1:4468249/71‑1 (MQ=255)

                               gAGCTGACCAGCACCTTCGATGATGCCAGCAAAACTGCCCTGATTAACAGTAACGTTGCCACCAGCATCGa                                       >  1:1032705/1‑71 (MQ=255)

                                aGCTGACCAGCACCTTCGATGATGCCAGCAAAACTGCCCTGATTAACAGTAACGTTGCCACCAGCATCGat                                      <  1:683805/71‑1 (MQ=255)

                                aGCTGACCAGCACCTTCGATGATGCCAGCAAAACTGCCCTGATTAACAGTAACGTTGCCACCAGCATCGat                                      <  1:970895/71‑1 (MQ=255)

                                aGCTGACCAGCACCTTCGATGATGCCAGCAAAACTGCCCTGATTAACAGTAACGTTGCCACCAGCATCGa                                       <  1:3915982/70‑1 (MQ=255)

                                  cTGACCAGCACCTTCGATGATGCCAGCAAAACTGCCCTGATTAACAGTAACGTTGCCACCAGCATCGatat                                    >  1:1246733/1‑71 (MQ=255)

                                   tGACCAGCACCTTCGATGATGCCAGCAAAACTGCCCTGATTAACAGTAACGTTGCCACCAGCATc                                         <  1:2924418/65‑1 (MQ=255)

                                     aCCAGCACCTTCGATGATGCCAGCAAAACTGCCCTGATTAACAGTAACGt                                                      >  1:4380400/1‑50 (MQ=255)

                                     aCCAGCACCTTCGATGATGCCAGCAAAACTGCCCTGATTAACAGTAACGt                                                      >  1:247087/1‑50 (MQ=255)

                                     aCCAGCACCTTCGATGATGCCAGCAAAACTGCCCTGATTAACAGTAACGTTGCCACCAGCATc                                         <  1:3699316/63‑1 (MQ=255)

                                     aCCAGCACCTTCGATGATGCCAGCAAAACTGCCCTGATTAACAGTAACGTTGCCACCAGCATCGATAttt                                  <  1:110187/70‑1 (MQ=255)

                                     aCCAGCACCTTCGATGATGCCAGCAAAACTGCCCTGATTAACAGTAACGTTGCCACCAGCATCGATATTTa                                 <  1:3296996/71‑1 (MQ=255)

                                     aCCAGCACCTTCGATGATGCCAGCAAAACTGCCCTGATTAACAGTAACGTTGCCACCAGCATCGATATTTa                                 <  1:1689223/71‑1 (MQ=255)

                                     aCCAGCACCTTCGATGATGCCAGCAAAACTGCCCTGATTAACAGT‑ACGTTGCCACCAGCATCGATATTTa                                 <  1:61155/70‑1 (MQ=255)

                                     aCCAGAACCTTCGATGATGCCAGCAAAACTGCCCTGATTAACAGTAACGTTGCCACCAGCATCGATATTTa                                 <  1:1552923/71‑1 (MQ=255)

                                      ccAGCACCTTCGATGATGCCATCAAAACTGCCCTGATTAACAGTAACGTTGCCACCAGCATCGATATTTa                                 >  1:3533606/1‑70 (MQ=255)

                                      ccAGCACCTTCGATGATGCCAGCAAAACTGCCTTGATTAACAGTAACGTTGCCACCAGCATCGATATTTaa                                <  1:578021/71‑1 (MQ=255)

                                      ccAGCACCTTCGATGATGCCAGCAAAACTGCCTTGATTAACAGTAACGTTGCCACCAGCATCGATATTTa                                 >  1:1351055/1‑70 (MQ=255)

                                      ccAGCACCTTCGATGATGCCAGCAAAACTGCCCTGATTAACAGTAACGTTGCCACCAGCATCGATATTTaa                                <  1:2057067/71‑1 (MQ=255)

                                          cACCTTCGATGATGCCAGCAAAACTGCCCTGATTAACAGTaa                                                         >  1:597902/1‑42 (MQ=255)

                                          cACCTTCGATGATGCCAGCAAAACTGCCCTGATTAACAGTAACGTTGCCACCAGCATCGATATTTaa                                <  1:1645314/67‑1 (MQ=255)

                                           aCCTTCGATGATGCCAGCAAAACTGCCCTGATTAACAGTa                                                          <  1:3262486/40‑1 (MQ=255)

                                           aCCTTCGATGATGCCAGCAAAACTGCCCTGATTAACAGTAA‑GTTGCCACCAGCATCGATATTTAAAGTGc                           >  1:3895360/1‑70 (MQ=255)

                                            ccTTCGATGATGCCAGCAAAACTGCCCTGATTAACAGTAACGt                                                      >  1:3321007/1‑43 (MQ=255)

                                            ccTTCGATGATGCCAGCAAAACTGCCCTGATTAACAGTAACGTTGCCACCAGCATCg                                        <  1:902910/57‑1 (MQ=255)

                                             cTTCGATGATGCCAGCAAAACTGCCCTGATTAACAGTAACGTTGcca                                                 <  1:3967230/47‑1 (MQ=255)

                                             cTTCGATGATGCCAGCAAAACTGCCCTGATTAACAGTAACGTTGCCACCAGCATCGatat                                    <  1:94022/60‑1 (MQ=255)

                                             cTTCGATGATGCCAGCAAAACTGCCCTGATTAACAGTAACGTTGCCACCAGCATCGATATTTAAAGTg                            <  1:4232519/68‑1 (MQ=255)

                                              ttCGATGATGCCAGCAAAACTGCCCTGATTAACAGTAACGTTGCCACCAGCATCg                                        <  1:1802614/55‑1 (MQ=255)

                                              ttCGATGATGCCAGCAAAACTGCCCTGATTAACAGTAACGTTGCCACCAGCATCGATATTTAAAGTGCCAt                        >  1:2439514/1‑71 (MQ=255)

                                              ttCGATGATGCCAGCAAAACTGCCCTGATTAACAGTAACGTTGCCACCAGCATCGATATTTAAAGTGCCAt                        >  1:185963/1‑71 (MQ=255)

                                               tCGATGATGCCAGCAAAACTGCCTTGATTAACAGTaa                                                         >  1:4287154/1‑37 (MQ=255)

                                               tCGATGATGCCAGCAAAACTGCCCTGATTAACAGTAACGtt                                                     >  1:1281577/1‑41 (MQ=255)

                                               tCGATGATGCCAGCAAAACTGCCCTGATTAACAGTAACGTTGCCACCAGCATCGa                                       >  1:434568/1‑55 (MQ=255)

                                               tCGATGATGCCAGCAAAACTGCCCTGATTAACAGTAACGTTGCCACCAGCATCGATATTTAAAGTGCCAtt                       >  1:4107687/1‑71 (MQ=255)

                                               tCGATGATGCCAGCAAAACTGCCCTGATTAACAGTAACGTTGCCACCAGCATCGATATTTAAAGTGCCAtt                       >  1:909022/1‑71 (MQ=255)

                                               tCGATGATGCCAGCAAAACTGCCCTGATTAACAGTAACGTTGCCACCAGCATCGATATTTAAAGTGCCAtt                       >  1:58434/1‑71 (MQ=255)

                                               tCGATGATGCCAGCAAAACTGCCCTGATTAACAGTAACGTTGCCACCAGCATCGATATTTAAAGTGCCAt                        >  1:1528331/1‑70 (MQ=255)

                                               tCGATGATGCCAGCAAAACTGCCCTGATTAACAGTAACGTTGCCACCAGCATCGATATTTAAAGTGACAtt                       >  1:2981650/1‑71 (MQ=255)

                                                cGATGATGCCAGCAAAACTGCCCTGATTAACAGTAACGt                                                      >  1:2992851/1‑39 (MQ=255)

                                                cGATGATGCCAGCAAAACTGCCCTGATTAACAGTAACGTTGCCACCAGCATCGATATTTAAAGTGCCAtt                       <  1:3970514/70‑1 (MQ=255)

                                                 gatgatGCCAGCAAAACTGCCCTGATTAACAGTAACGTTGCCACCAGCATCGATATTTAAAGTGCCATTc                      <  1:247732/70‑1 (MQ=255)

                                                  atgatgCCAGCAAAACTGCCCTGATTAACAGTAACGTTGCCACCAGCATCGATATTTAAAGTGCCATTCTg                    <  1:2309130/71‑1 (MQ=255)

                                                   tgatgCCAGCAAAACTGCCCTGATTAACAGTAACGTTGCCACCAGCATCGATATTTAAAGTGCCATTCTGa                   <  1:1415565/71‑1 (MQ=255)

                                                      tgCCAGCAAAACTGCCCTGATTAACAGTAACGTTGCCACCAGCATCGATATTTAAAGTGCCATTCTGAAAg                >  1:561512/1‑71 (MQ=255)

                                                      tgCCAGCAAAACTGCCCTGATTAACAGTAACGTTGCCACCAGCATCGATATTTAAAGTGCCATTCTGAAAg                <  1:503650/71‑1 (MQ=255)

                                                       gCCAGCAAAACTGCCCTGATTAACAGTAACGTTGCCACCAGCa                                           >  1:2528618/1‑43 (MQ=255)

                                                       gCCAGCAAAACTGCCCTGATTAACAGTAACGTTGCCACCAGCATCGATATTTAAAGTGCCATTCTGAAAGc               >  1:3110663/1‑71 (MQ=255)

                                                         cAGCAAAACTGCCCTGATTAACAGTAACGTTGCCACCAGCa                                           <  1:38927/41‑1 (MQ=255)

                                                           gCAAAGCTGCCCTGATTAACAGTAACGTTGCCACCAGCATCGATATTTaa                                <  1:2904124/50‑1 (MQ=255)

                                                           gCAAAACTGCCCTGATTAACAGTAACGTTGCCACCAGCATCGATATTTAAAGTGCCATTCTGAAAGCCCGt           >  1:1600984/1‑71 (MQ=255)

                                                            cAAAACTGCCCTGATTAACAGTAACGTTGCCACCAGCATCGa                                       >  1:2013522/1‑42 (MQ=255)

                                                            cAAAACTGCCCTGATTAACAGTAACGTTGCCACCAGCATCGATATTTAAAGTGCCATTCTGAAAGCCCGt           >  1:2492155/1‑70 (MQ=255)

                                                            cAAAACTGCCCTGATTAACAGTAACGTTGCCACCAGCATCGATATTTAAAGTGCCATTCTGAAAGCCCGt           >  1:1664002/1‑70 (MQ=255)

                                                             aaaaCTGCCCTGATTAACAGTAACGTTGCCACCAGCATCGATATTTAAAGTGCCATTCTGAAAGCCCGTCa         <  1:1020345/71‑1 (MQ=255)

                                                             aaaaCTGCCCTGATTAACAGTAACGTTGCCACCAGCATCGATATTTAAAGTGCCATTCTGAAAGACCGTc          >  1:1825243/1‑70 (MQ=255)

                                                              aaaCTGCCCTGATTAACAGTAACGTTGCCACCAGCATCGATATTTAAAg                              >  1:2726940/1‑49 (MQ=255)

                                                              aaaCTGCCCTGATTAACAGTAACGTTGCCACCAGCATCGATATTTAAAg                              >  1:3428182/1‑49 (MQ=255)

                                                               aaCTGCCCTGATTAACAGTAACGTTGCCACCAGCATCg                                        <  1:1614550/38‑1 (MQ=255)

                                                               aaCTGCCCTGATTAACAGTAACGTTGCCACCAGCATCg                                        <  1:3914060/38‑1 (MQ=255)

                                                                aCTGCCCTGATTAACAGTAACGTTGCCACCAGCATCGATATTTAAAGTGCCATTCTGAAAGCCCGTCAATg      <  1:881442/71‑1 (MQ=255)

                                                                 cTGCCTTGATTAACAGTAACGTTGCCACCAGCATCGATATTTAAAGTGCCATTCTGAAAGCCCGTCAATGc     <  1:2104221/71‑1 (MQ=255)

                                                                 cTGCCCTGATTAACAGTAACGTTGCCACCAGCATCGATATTTAAAGTGCCATTCTGAAAGCCCGTCAATGc     >  1:449344/1‑71 (MQ=255)

                                                                 cTGCCCTGATTAACAGTAACGTTGCCACCAGCATCGATATTTAAAGTGCCATTCTGAAAGCCCGTCAATGc     >  1:3794277/1‑71 (MQ=255)

                                                                   gCCCTGATTAACAGTAACGTTGCCACCAGCATCGATATTTAAAGTGCCATTCTGAAAGCCCGTCAATGCGt   >  1:4009408/1‑71 (MQ=255)

                                                                    cccTGATTAACAGTAACGTTGCCACCAGCATCGATATTTAAAGTGCCATTCTGAAAGCCCGTc          <  1:3357767/63‑1 (MQ=255)

                                                                    cccTGATTAACAGTAACGTTGCCACCAGCATCGATATTTAAAGTGCCATTCTGAAAGCCCGTCAATGCg    <  1:2184059/69‑1 (MQ=255)

                                                                    cccTGATTAACAGTAACGTTGCCACCAGCATCGATATTTAAAGTGCCATTCTGAAAGCCCGTCAATGCGTg  <  1:86400/71‑1 (MQ=255)

                                                                      cTGATTAACAGTAACGTTGCCACCAGCATCGATATTTaa                                <  1:3660409/39‑1 (MQ=255)

                                                                      cTGATTAACAGTAACGTTGCCACCAGCATCGATATTTAAAGTGCCATTCTGAAAGCCCGTc          <  1:4151199/61‑1 (MQ=255)

                                                                      cTGATTAACAGTAACGTTGCCACCAGCATCGATATTTAAAGTGCCATTCTGAAAGCCCGTc          <  1:701408/61‑1 (MQ=255)

CCAGCACGTAGCTGCCGTTTTGCGCAATGGTGAGCTGACCAGCACCTTCGATGATGCCAGCAAAACTGCCCTGATTAACAGTAACGTTGCCACCAGCATCGATATTTAAAGTGCCATTCTGAAAGCCCGTCAATGCGTG  >  minE/1470365‑1470503

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