BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	455,016	A→G	36.6%	F23F (TTT→TTC)	ybgH ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	455,016	0	A	G	36.6%	104.0 / 118.6	123	F23F (TTT→TTC)	ybgH	predicted transporter
Reads supporting (aligned to +/- strand): ref base A (37/41); new base G (19/26); total (56/67)
Fisher's exact test for biased strand distribution p-value = 7.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.66e-01

Only 100 of 143 total aligned reads displayed.

GTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTAT  >  Exported/454947‑455080

gTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTaa                                                                 <  1:3644714/71‑1 (MQ=255)

gTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTaa                                                                 <  1:4791161/71‑1 (MQ=255)

 tGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTaaa                                                                <  1:1362548/71‑1 (MQ=255)

   aGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAat                                                              <  1:3626119/71‑1 (MQ=255)

     tATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAatat                                                            <  1:3961095/71‑1 (MQ=255)

      aTCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTaa                                                                 <  1:3515809/65‑1 (MQ=255)

       tCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTaa                                                                 <  1:46847/64‑1 (MQ=255)

        cGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTaa                                                                 <  1:4189223/63‑1 (MQ=255)

         gTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTcc                                                         >  1:2946076/1‑70 (MQ=255)

         gTTGTATTTTAGATGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATAtt                                                           >  1:5247165/1‑68 (MQ=255)

           tGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCa                                                       <  1:5350036/70‑1 (MQ=255)

                tttAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTcc                                                         <  1:1029401/63‑1 (MQ=255)

                 ttAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTg                                                 <  1:414974/70‑1 (MQ=255)

                     ttGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAgcg                                            >  1:2372739/1‑71 (MQ=255)

                     ttGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAg                                              <  1:2972474/69‑1 (MQ=255)

                      tGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAgcg                                            >  1:1587338/1‑70 (MQ=255)

                       gATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTaaaa                                                               >  1:5095262/1‑50 (MQ=255)

                       gATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTaaaa                                                               >  1:2646619/1‑50 (MQ=255)

                       gATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAatat                                                            <  1:782649/53‑1 (MQ=255)

                       gATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAgcgc                                           >  1:3082149/1‑70 (MQ=255)

                         ttGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTaaaa                                                               <  1:4013472/48‑1 (MQ=255)

                         ttGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTaaaa                                                               <  1:1316045/48‑1 (MQ=255)

                         ttGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTaaaa                                                               <  1:686172/48‑1 (MQ=255)

                         tgggTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCAc                                        <  1:3862299/69‑1 (MQ=255)

                         tgggTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCAc                                        <  1:3353582/69‑1 (MQ=255)

                         tgggTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCAc                                        <  1:3263095/69‑1 (MQ=255)

                            gTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGAt                                                    >  1:4363715/1‑56 (MQ=255)

                             tGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACa                                    <  1:1790594/71‑1 (MQ=255)

                              gagaTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAgta                                                             >  1:507627/1‑45 (MQ=255)

                              gagaTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAgta                                                             >  1:5553697/1‑45 (MQ=255)

                              gagaTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAgta                                                             >  1:5420728/1‑45 (MQ=255)

                              gagaTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAgta                                                             >  1:5342075/1‑45 (MQ=255)

                              gagaTAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTaaa                                                                <  1:1035054/42‑1 (MQ=255)

                                gaTAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAg                                 <  1:265496/71‑1 (MQ=255)

                                gaTAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAg                                 >  1:3082007/1‑71 (MQ=255)

                                 aTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGAt                                                    <  1:2419711/51‑1 (MQ=255)

                                 aTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGAtata                                  <  1:2225637/69‑1 (MQ=255)

                                 aTAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAg                                 >  1:2010975/1‑70 (MQ=255)

                                  tAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGAtata                                  <  1:5071689/68‑1 (MQ=255)

                                  tAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGt                                <  1:4606702/70‑1 (MQ=255)

                                  tAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTaa                                                                 <  1:2885980/37‑1 (MQ=255)

                                  tAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCTCAACGACATAGTa                               <  1:1553546/71‑1 (MQ=255)

                                  tAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTACCAGATTTGCAGCGCAACGACATAGTa                               <  1:4030605/71‑1 (MQ=255)

                                   agagAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTa                               >  1:136720/1‑70 (MQ=255)

                                   agagAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTa                               <  1:4914839/70‑1 (MQ=255)

                                    gagaATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGata                                    >  1:2970162/1‑62 (MQ=255)

                                    gagaATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGata                                    >  1:3129678/1‑62 (MQ=255)

                                    gagaATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGt                                >  1:4002479/1‑68 (MQ=255)

                                    gagaATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTat                              <  1:700044/70‑2 (MQ=255)

                                    gagaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa                             <  1:3352644/71‑1 (MQ=255)

                                    gagaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa                             >  1:328063/1‑71 (MQ=255)

                                    gagaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa                             <  1:762769/71‑1 (MQ=255)

                                    gagaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa                             <  1:2087903/71‑1 (MQ=255)

                                    gagaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa                             <  1:3796277/71‑1 (MQ=255)

                                    gagaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa                             <  1:4845549/71‑1 (MQ=255)

                                      gaATCAGCAGGGCACGCATGCCGTAAAAGCTaaaa                                                               <  1:1031592/35‑1 (MQ=255)

                                       aaTCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATa                                  >  1:1265823/1‑63 (MQ=255)

                                       aaTCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATa                                  >  1:4607506/1‑63 (MQ=255)

                                        aTCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAg                          <  1:1708427/70‑1 (MQ=255)

                                          cagcagGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAgcgc                       >  1:2903474/1‑71 (MQ=255)

                                          cagcagGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAgcgc                       >  1:4318729/1‑71 (MQ=255)

                                           agcagGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAgcgcg                      >  1:1768082/1‑71 (MQ=255)

                                           agcagGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAgcgcg                      >  1:626790/1‑71 (MQ=255)

                                             cagGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGcgcgcg                    >  1:925877/1‑71 (MQ=255)

                                              agGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACa                                    <  1:283032/54‑1 (MQ=255)

                                              agGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCgg                   >  1:2968215/1‑71 (MQ=255)

                                              agGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCgg                   >  1:1929373/1‑71 (MQ=255)

                                               aggCACGCATGCCGTAAAAACTGAAGTACTCCCAGAttt                                                  >  1:501213/2‑39 (MQ=21)

                                                  cACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtgt               <  1:4503315/71‑1 (MQ=255)

                                                  cACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtgt               <  1:2786885/71‑1 (MQ=255)

                                                  cACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtg                <  1:1933962/70‑1 (MQ=255)

                                                  cACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCtgt               <  1:3778408/71‑1 (MQ=255)

                                                   aCGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAg                                              >  1:1352737/1‑39 (MQ=255)

                                                   aCGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAg                                              >  1:5562898/1‑39 (MQ=255)

                                                    cGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGa             <  1:70855/71‑1 (MQ=255)

                                                     gCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGag            >  1:2119683/1‑70 (MQ=255)

                                                      cATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAAc                                        >  1:3208739/1‑42 (MQ=255)

                                                      cATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGTTGTGATg           >  1:3128724/1‑71 (MQ=255)

                                                      cATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGagg           <  1:121246/71‑3 (MQ=255)

                                                        tGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTa                               <  1:4810159/49‑1 (MQ=255)

                                                        tGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTatat                            <  1:4975880/52‑1 (MQ=255)

                                                        tGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGagg           <  1:1998365/69‑3 (MQ=255)

                                                          ccGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCa                                          >  1:4484290/1‑36 (MQ=255)

                                                          ccGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCa                                          >  1:4221607/1‑36 (MQ=255)

                                                          ccGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCgtg       >  1:442478/1‑71 (MQ=255)

                                                             tAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATa                                  <  1:4168569/41‑1 (MQ=255)

                                                             tAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGt                                <  1:941405/43‑1 (MQ=255)

                                                             tAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGt                                <  1:1091625/43‑1 (MQ=255)

                                                             tAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa                             <  1:4982778/46‑1 (MQ=255)

                                                             tAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCtgt               <  1:5459242/60‑1 (MQ=255)

                                                               aaaaCTGAAGTACTCCCAGATTTGTATCGCCACGATATAGTATATCGCGCGCGGCTGTGag            >  1:3834926/1‑60 (MQ=255)

                                                               aaaaCTGAAGTACTCCCAGATTTGGATCGCCACGATATAg                                 <  1:236756/40‑1 (MQ=255)

                                                               aaaaCTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGA‑GGGTGttttc  <  1:1798406/70‑3 (MQ=255)

                                                               aaaaCTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGA‑GGGTGttttc  <  1:4056179/70‑3 (MQ=255)

                                                               aaaaCTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGA‑GGGTGttttc  <  1:3652356/70‑3 (MQ=255)

                                                               aaaGCTAAAATATTCCCAGATTTGCAGCGCAACGACAt                                   >  1:776803/1‑38 (MQ=255)

                                                               aaaGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAg                          >  1:1839905/1‑47 (MQ=255)

                                                               aaaGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTAt  >  1:3473302/1‑71 (MQ=255)

                                                                  gCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGc                  <  1:2532398/52‑1 (MQ=255)

                                                                    tAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtgt               <  1:515501/53‑1 (MQ=255)

GTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTAT  >  Exported/454947‑455080

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