BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	826,746	T→A	11.0%	S381S (TCT→TCA)	ychP →	predicted invasin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	826,746	0	T	A	11.0%	285.3 / 23.1	100	S381S (TCT→TCA)	ychP	predicted invasin
Reads supporting (aligned to +/- strand): ref base T (53/36); new base A (4/7); total (57/43)
Fisher's exact test for biased strand distribution p-value = 1.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

Only 100 of 112 total aligned reads displayed.

CGGAGGGCTGGACGCTGATCATGCCTGACTGGCAGAACGGGGAAGGGGCAAGCAATCACTGGCGATTGTCTGTGGTGGTGGAAGATAACCAGGGGCAGCGTGTCTCCTCCAATGAGATCACGCTAACGCTTGTCGA  >  Exported/826676‑826811

cgGAGGGCTGGACGCTGATCATGCCTGACTGGCAGAACGGGGAAGGGGCAAGCAATCACTGGCGATTGTCt                                                                   <  1:1882100/71‑1 (MQ=255)

  gAGGGCTGGACGCTGATCATGCCTGACTGGCAGAACGGGGAAGGGGCAAGCAATCACTGGCGATTGTCtgt                                                                 >  1:4718191/1‑71 (MQ=255)

  gAGGGCTGGACGCTGATCATGCCTGACTGGCAGAACGGGGAAGGGGCAAGCAATCACTGGCGATTGTCtg                                                                  <  1:1744409/70‑1 (MQ=255)

     ggCTGGACGCTGATCATGCCTGACTGGCAGAACGGGGAAGGGGCAAGCAATCACTGGCGATTGTCTgtggt                                                              >  1:4297582/1‑71 (MQ=255)

     ggCTGGACGCTGATCATGCCTGACTGGCAGAACGGGGAAGGGGCAAGCAATCACTGGCGATTGTCTgtggt                                                              <  1:2195724/71‑1 (MQ=255)

      gCTGGACGCTGATCATGCCTGACTGGCAGAACGGGGAAGGGGCAAGCAATCACTGGCGATTGTCTgtggtg                                                             <  1:2204267/71‑1 (MQ=255)

           aCGCTGATCATGCCTGACTGGCAGAACGGGGAAGGGGCAAGCAATCACTGGCGATTGTCTGtggtggtgg                                                         >  1:4449204/1‑70 (MQ=255)

           aCGCTGATCATGCCTGACTGGCAGAACGGGGAAGGGGCAAGCAATCACTGGCGATTGTCTGtggtggtgg                                                         >  1:5165053/1‑70 (MQ=255)

           aCGCTGATCATGCCTGACTGGCAGAACGGGGAAGGGGCAAGCAATCACTGGCGATTGTCTGtggtggtgg                                                         >  1:721280/1‑70 (MQ=255)

           aCGCTGATCATGCCTGACTGGCAGAACGGGGAAGGGGCAAGCAATCACTGGCGATTGTCTGtggtggtgg                                                         >  1:226944/1‑70 (MQ=255)

           aCGCTGATCATGCCTGACTGGCAGAACGGGGAAGGGGCAAGCAATCACTGGCGATTGTCTGtggtggtgg                                                         >  1:2166535/1‑70 (MQ=255)

           aCGCTGATCATGCCTGACTGGCAGAACGGGGAAGGGGCAAGCAATCACTGGCGATTGTCTGtggtggtgg                                                         >  1:3464438/1‑70 (MQ=255)

           aCGCTGATCATGCCTGACTGGCAGAACGGGGAAGGGGCAAGCAATCACTGGCGATTGTCTGtggtggtgg                                                         >  1:3570425/1‑70 (MQ=255)

            cGCTGATCATGCCTGACTGGCAGAACGGGGAAGGGGCAAGCAATCACTGGCGATTGTCTGTGGTGGTGGa                                                        <  1:1643884/70‑1 (MQ=255)

            cGCTGATCATGCCTGACTGGCAGAACGGGGAAGGGGCAAGCAATCACTGGCGATTGTCTGTGGTGGTGGa                                                        <  1:2824707/70‑1 (MQ=255)

            cGCTGATCATGCCTGACTGGCAGAACGGGGAAGGGGCAAGCAATCACTGGCGATTGTCTGTGGTGGTGGa                                                        <  1:4048300/70‑1 (MQ=255)

             gCTGATCATGCCTGACTGGCAGAACGGGGAAGGGGCAAGCAATCACTGGCGATTGTCTGtggtggtgg                                                         >  1:670893/1‑68 (MQ=255)

             gCTGATCATGCCTGACTGGCAGAACGGGGAAGGGGCAAGCAATCACTGGCGATTGTCTGtggtggtgg                                                         >  1:1289157/1‑68 (MQ=255)

                   cATGCCTGACTGGCAGAACGGGGAAGGGGCAAGCAATCACTGGCGATTGTCTGTGGTGGTGGa                                                        <  1:4396616/63‑1 (MQ=255)

                    aTGCCTGACTGGCAGAACGGGGAAGGGGCAAGCAATCACTGGCGATTGTCTGTGGTGGTGGAAGATAACCa                                               <  1:1333607/71‑1 (MQ=255)

                       ccTGACTGGCAGCACGGGGAAGGGGCAAGCAATCACTGGCGATTGTCTGTGGTGGTGGAAGATAACCAggg                                            >  1:466942/1‑71 (MQ=255)

                       ccTGACTGGCAGAACGGGGAAGGGGCAAGCAATCACTGGCGATTGTCTGTGGTGGTGGAAGATAACCAggg                                            >  1:2245823/1‑71 (MQ=255)

                       ccTGACTGGCAGAACGGGGAAGGGGCAAGCAATCACTGGCGATTGTCTGTGGTGGTGGAAGATAACCAggg                                            >  1:878393/1‑71 (MQ=255)

                       ccTGACTGGCAGAACGGGGAAGGGGCAAGCAATCACTGGCGATTGTCTGTGGTGGTGGAAGATAACCAgg                                             >  1:1840714/1‑70 (MQ=255)

                       ccTGACTGGCAGAACGGGGAAGGGGCAAGCAATCACTGGCGATTGTCTGTGGTGGTGGAAGATAACCAgg                                             <  1:4713006/70‑1 (MQ=255)

                       ccTGACTGGCAGAACGGGGAAGGGGCAAGCAATCACTGGCGATTGTCTGTGGTGGTGGAAGATAACCAgg                                             <  1:3318184/70‑1 (MQ=255)

                         tGACTGGCAAAACGGGGAAGGCGCAAGCAATCACTGGCGATTGTCAGTGGTGGTGGAAGATAACCAgggg                                           >  1:2618852/1‑70 (MQ=255)

                         tGACTGGCAAAACGGGGAAGGCGCAAGCAATCACTGGCGATTGTCAGTGGTGGTGGAAGATAACCAgggg                                           >  1:2962748/1‑70 (MQ=255)

                         tGACTGGCAAAACGGGGAAGGCGCAAGCAATCACTGGCGATTGTCAGTGGTGGTGGAAGATAACCAgggg                                           >  1:2612645/1‑70 (MQ=255)

                         tGACTGGCAAAACGGGGAAGGCGCAAGCAATCACTGGCGATTGTCAGTGGTGGTGGAAGATAACCAggg                                            >  1:1315284/1‑69 (MQ=255)

                           aCTGGCAGAACGGGGAAGGGGCAAGCAATCACTGGCGATTGTCTGTGGTGGTGGAAGATAACCAgggg                                           >  1:3681987/1‑68 (MQ=255)

                           aCTGGCAGAACGGGGAAGGGGCAAGCAATCACTGGCGATTGTCTGTGGTGGTGGAAGATAACCAgggg                                           >  1:1078315/1‑68 (MQ=255)

                           aCTGGCAGAACGGGGAAGGGGCAAGCAATCACTGGCGATTGTCTGTGGTGGTGGAAGATAACCAgggg                                           >  1:81885/1‑68 (MQ=255)

                           aCTGGCAGAACGGGGAAGGGGCAAGCAATCACTGGCGATTGTCTGTGGTGGTGGAAGATAACCAGGGGCAg                                        >  1:395934/1‑71 (MQ=255)

                           aCTGGCAGAACGGGGAAGGGGCAAGCAATCACTGGCGATTGTCTGTGGTGGTGGAAGATAACCAGGGGCAg                                        >  1:46532/1‑71 (MQ=255)

                           aCTGGCAGAACGGGGAAGGGGCAAGCAATCACTGGCGATTGTCTGTGGTGGTGGAAGATAACCAGGGGCAg                                        >  1:674736/1‑71 (MQ=255)

                            cTGGCAGAACGGGGAAGGGGCAAGCAATCACTGGCGATTGTCTGTGGTGGTGGAAGATAACCAGGGGCAGc                                       >  1:2338321/1‑71 (MQ=255)

                             tGGCAGAACGGGGAAGGGGCAAGCAATCACTGGCGATTGTCtgtg                                                                >  1:4572935/1‑45 (MQ=255)

                             tGGCAGAACGGGGAAGGGGCAAGCAATCACTGGCGATTGTCtgtg                                                                >  1:3889551/1‑45 (MQ=255)

                             tGGCAGAACGGGGAAGGGGCAAGCAATCACTGGCGATTGTCtgtg                                                                >  1:3516559/1‑45 (MQ=255)

                             tGGCAGAACGGGGAAGGGGCAAGCAATCACTGGCGATTGTCtgtg                                                                >  1:1660646/1‑45 (MQ=255)

                             tGGCAGAACGGGGAAGGGGCAAGCAATCACTGGCGATTGTCtgtg                                                                >  1:2711388/1‑45 (MQ=255)

                             tGGCAGAACGGGGAAGGGGCAAGCAATCACTGGCGATTGTCtgtg                                                                >  1:2025382/1‑45 (MQ=255)

                             tGGCAGAACGGGGAAGGGGCAAGCAATCACTGGCGATTGTCtgtg                                                                >  1:54018/1‑45 (MQ=255)

                             tGGCAGAACGGGGAAGGGGCAAGCAATCACTGGCGATTGTCtgtg                                                                >  1:4932930/1‑45 (MQ=255)

                               gCAGAACGGGGAAGGGGCAAGCAATCACTGGCGATTGTCTGTGGTGGTGGAAGa                                                     >  1:2663189/1‑54 (MQ=255)

                               gCAGAACGGGGAAGGGGCAAGCAATCACTGGCGATTGTCTGTGGTGGTGGAAGa                                                     >  1:3100554/1‑54 (MQ=255)

                               gCAGAACGGGGAAGGGGCAAGCAATCACTGGCGATTGTCTGTGGTGGTGGAAGa                                                     >  1:2102987/1‑54 (MQ=255)

                               gCAGAACGGGGAAGGGGCAAGCAATCACTGGCGATTGTCTGTGGTGGTGGAAGa                                                     >  1:4900077/1‑54 (MQ=255)

                               gCAGAACGGGGAAGGGGCAAGCAATCACTGGCGATTGTCTGTGGTGGTGGAAGa                                                     >  1:4344736/1‑54 (MQ=255)

                                  gAACGGGGAAGGGGCAAGCAATCACTGGCGATTGTCTGTGGTGGTGGAAGATAACCAGGGGCAGCGTGtct                                 <  1:4485138/71‑1 (MQ=255)

                                  gAACGGGGAAGGGGCAAGCAATCACTGGCGATTGTCTGTGGTGGTGGAAGATAACCAGGGGCAGCGTGtct                                 <  1:1160472/71‑1 (MQ=255)

                                    aCGGGGAAGGGGCAAGCAATCACTGGCGATTGTCTGTGGTGGTGGAAGATAACCAGGGGCAGCGTGTctcc                               >  1:1771822/1‑71 (MQ=255)

                                    aCGGGGAAGGGGCAAGCAATCACTGGCGATTGTCTGTGGTGGTGGAAGATAACCAGGGGCAGCGTGTctcc                               >  1:2338931/1‑71 (MQ=255)

                                    aCGGGGAAGGGGCAAGCAATCACTGGCGATTGTCTGTGGTGGTGGAAGATAACCAGGGGCAGCGTGTctcc                               >  1:2719806/1‑71 (MQ=255)

                                     cGGGGAAGGGGCAAGCAATCACTGGCGATTGTCTGTGGTGGTGGAAGATa                                                   <  1:3987089/50‑1 (MQ=255)

                                     cGGGGAAGGGGCAAGCAATCACTGGCGATTGTCTGTGGTGGTGGAAGATa                                                   <  1:4749063/50‑1 (MQ=255)

                                     cGGGGAAGGGGCAAGCAATCACTGGCGATTGTCTGTGGTGGTGGAAGATAACCAGGGGCAGCGTGTctcct                              <  1:3226958/71‑1 (MQ=255)

                                     cGGGGAAGGGGCAAGCAATCACTGGCGATTGTCTGTGGTGGTGGAAGATAACCAGGGGCAGCGTGTctcct                              <  1:3913078/71‑1 (MQ=255)

                                     cGGGGAAGGGGCAAGCAATCACTGGCGATTGTCTGTGGTGGTGGAAGATAACCAGGGGCAGCGTGTctcct                              <  1:3544216/71‑1 (MQ=255)

                                     cGGGGAAGGGGCAAGCAATCACTGGCGATTGTCTGTGGTGGTGGAAGATAACCAGGGGCAGCGTGTctcct                              <  1:1329415/71‑1 (MQ=255)

                                          aaGGGGCAAGCAATCACTGGCGATTGTCTGTGGTGGTGGAAGATAACCAGGGGCAGCGTGTCTCCTCCAAt                         <  1:2659336/71‑1 (MQ=255)

                                                    caaTCACTGGCGATTGTCTGTGGTGGTGGAAGATAACCAGGGGCAGCGTGTCTCCTCCAATGAGATCACGc               >  1:4919781/1‑71 (MQ=255)

                                                     aaTCACTGGCGATTGTCTGTGGTGGTGGAAGATAACCAGGGGCAGCGTGTCTCCTCCAATGAGATCACGCt              <  1:2617305/71‑1 (MQ=255)

                                                         aCTGGCGATTGTCTGTGGTGGTGGAAGATAACCAGGGGCAGCgt                                     >  1:4240547/1‑44 (MQ=255)

                                                         aCTGGCGATTGTCTGTGGTGGTGGAAGATAACCAGGGGCAGCGTGTCTCCTCCAATgag                      >  1:3818958/1‑59 (MQ=255)

                                                         aCTGGCGATTGTCTGTGGTGGTGGAAGATAACCAGGGGCAGCGTGTCTCCTCCAATGAGATCACGCTAAc           >  1:4872028/1‑70 (MQ=255)

                                                         aCTGGCGATTGTCTGTGGTGGTGGAAGATAACCAGGGGCAGCGTGTCTCCTCCAATGAGATCACGCTAACg          >  1:4381833/1‑71 (MQ=255)

                                                         aCTGGCGATTGTCTGTGGTGGTGGAAGATAACCAGGGGCAGCGTGTCTCCTCCAATGAGATCACGCTAACg          >  1:5078025/1‑71 (MQ=255)

                                                         aCTGGCGATTGTCAGTGGTGGTGGAAGATAACCAGGGGCAGCGTGTCTCCTCCAATGAGATCACGCTAAc           <  1:5064215/70‑1 (MQ=255)

                                                         aCTGGCGATTGTCAGTGGTGGTGGAAGATAACCAGGGGCAGCGTGTCTCCTCCAATGAGATCACGCTAAc           <  1:867287/70‑1 (MQ=255)

                                                         aCTGGCGATTGTCAGTGGTGGTGGAAGATAACCAGGGGCAGCGTGTCTCCTCCAATGAGATCACGCTAAc           <  1:3068329/70‑1 (MQ=255)

                                                         aCTGGCGATTGTCAGTGGTGGTGGAAGATAACCAGGGGCAGCGTGTCTCCTCCAATGAGATCACGCTAAc           <  1:2944399/70‑1 (MQ=255)

                                                         aCTGGCGATTGTCAGTGGTGGTGGAAGATAACCAGGGGCAGCGTGTCTCCTCCAATGAGATCACGCTAAc           <  1:2662471/70‑1 (MQ=255)

                                                         aCTGGCGATTGTCAGTGGTGGTGGAAGATAACCAGGGGCAGCGTGTCTCCTCCAATGAGATCACGCTAAc           <  1:2085507/70‑1 (MQ=255)

                                                         aCTCGCGATTGTCTGTGGTGGTGGAAGATAACCAGGGGCAGCGTGTCTCCTCCAATGAGAt                    >  1:1135561/1‑61 (MQ=255)

                                                          cTGGCGATTGTCAGTGGTGGTGGAAGATAACCAGGGGCAGCGTGTCTCCTCCAATGAGATCACGCTAAc           <  1:2938176/69‑1 (MQ=255)

                                                           tGGCGATTGTCTGTGGTGGTGGAAGATAACCAGGGGCAGCGTGtct                                 <  1:3463103/46‑1 (MQ=255)

                                                           tGGCGATTGTCTGTGGTGGTGGAAGATAACCAGGGGCAGCGTGtc                                  <  1:1500205/45‑1 (MQ=255)

                                                           tGGCGATTGTCTGTGGTGGTGGAAGATAACCAGGGGCAGCGTGtc                                  <  1:1706789/45‑1 (MQ=255)

                                                           tGGCGATTGTCTGTGGTGGTGGAAGATAACCAGGGGCAGCGTGtc                                  <  1:4779093/45‑1 (MQ=255)

                                                           tGGCGATTGTCTGTGGTGGTGGAAGATAACCAGGGGCAGCGTGtc                                  <  1:2007190/45‑1 (MQ=255)

                                                           tGGCGATTGTCTGTGGTGGTGGAAGATAACCAGGGGCAGCGTGtc                                  <  1:687610/45‑1 (MQ=255)

                                                           tGGCGATTGTCTGTGGTGGTGGAAGATAACCAGGGGCAGCGTGtc                                  <  1:4693352/45‑1 (MQ=255)

                                                           tGGCGATTGTCTGTGGTGGTGGAAGATAACCAGGGGCAGCGTGtc                                  <  1:2325954/45‑1 (MQ=255)

                                                           tGGCGATTGTCTGTGGTGGTGGAAGATAACCAGGGGCAGCGTGtc                                  <  1:368028/45‑1 (MQ=255)

                                                           tGGCGATTGTCTGTGGTGGTGGAAGATAACCAGGGGCAGCGTGTCTCCTCCAATGAGATCACGCTAACGCt        >  1:4011542/1‑71 (MQ=255)

                                                            ggCGATTGTCTGTGGTGGTGGAAGATAACCAGGGGCAGCgtgt                                   >  1:4070974/1‑43 (MQ=255)

                                                             gCGATTGTCTGTGGTGGTGGAAGATAACCAGGGGCAGCGTGTctcc                               >  1:718889/1‑46 (MQ=255)

                                                             gCGATTGTCTGTGGTGGTGGAAGATAACCAGGGGCAGCGTGTctcc                               >  1:3316315/1‑46 (MQ=255)

                                                             gCGATTGTCTGTGGTGGTGGAAGATAACCAGGGGCAGCGTGTctcc                               >  1:4765244/1‑46 (MQ=255)

                                                             gCGATTGTCTGTGGTGGTGGAAGATAACCAGGGGCAGCGTGTctcc                               >  1:1979439/1‑46 (MQ=255)

                                                             gCGATTGTCTGTGGTGGTGGAAGATAACCAGGGGCAGCGTGTctcc                               >  1:1572262/1‑46 (MQ=255)

                                                             gCGATTGTCTGTGGTGGTGGAAGATAACCAGGGGCAGCGTGTCTCCTCCAATgag                      >  1:2251597/1‑55 (MQ=255)

                                                             gCGATTGTCTGTGGTGGTGGAAGATAACCAGGGGCAGCGTGTCTCCTCCAATgag                      >  1:1587934/1‑55 (MQ=255)

                                                             gCGATTGTCTGTGGTGGTGGAAGATAACCAGGGGCAGCGTGTCTCCTCCAATgag                      >  1:1069458/1‑55 (MQ=255)

                                                                aTTGTCTGTTGTGGTGGAAGATAACCAGGGGCAGCGTGTCTCCTCCAATGAGATCACGCTAACGCTTGTCg   <  1:2807632/71‑1 (MQ=255)

                                                                aTTGTCTGTGGTGGTGGAAGATAACCAGGGGCAGCGTGTCTCCTCCAATGAGATCACGCTAACGCTTGTCg   <  1:1006645/71‑1 (MQ=255)

                                                                aTTGTCTGTGGTGGTGGAAGATAACCAGGGGCAGCGTGTCTCCTCCAATGAGATCACGCTAACGCTTGTCg   <  1:1850337/71‑1 (MQ=255)

                                                                 ttGTCTGTGGTGGTGGAAGATAACCAGGGGCAGCGTGTCTCCTCCAATGAGATCACGCTAACGCTTGTCGa  <  1:4121292/71‑1 (MQ=255)

CGGAGGGCTGGACGCTGATCATGCCTGACTGGCAGAACGGGGAAGGGGCAAGCAATCACTGGCGATTGTCTGTGGTGGTGGAAGATAACCAGGGGCAGCGTGTCTCCTCCAATGAGATCACGCTAACGCTTGTCGA  >  Exported/826676‑826811

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