BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	455,019	A→G	35.3%	Y22Y (TAT→TAC)	ybgH ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	455,019	0	A	G	35.3%	109.8 / 105.4	116	Y22Y (TAT→TAC)	ybgH	predicted transporter
Reads supporting (aligned to +/- strand): ref base A (37/38); new base G (15/26); total (52/64)
Fisher's exact test for biased strand distribution p-value = 2.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 133 total aligned reads displayed.

AGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAAT  >  Exported/454950‑455087

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

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

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

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

        tGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTcc                                                                >  1:776356/1‑68 (MQ=255)

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

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

                 gTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAg                                                     <  1:5079031/70‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

                      ttGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAAc                                               >  1:2081188/1‑71 (MQ=255)

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

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

                          tGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGAttt                                                         <  1:3099311/57‑1 (MQ=255)

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

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

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

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

                            agaTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGatat                                          <  1:3291639/70‑1 (MQ=255)

                            agaTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGatat                                          <  1:3771237/70‑1 (MQ=255)

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

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

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

                              aTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGAt                                                           <  1:2896687/51‑1 (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)

                               tAGAGAATCAGCAGTGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTa                                      <  1:1399617/71‑1 (MQ=255)

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

                                agagAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAg                                                             >  1:4100467/1‑47 (MQ=255)

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

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

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

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

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

                                 gagaATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTata                                    <  1:700401/71‑3 (MQ=255)

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

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

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

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

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

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

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

                                  agaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa                                    >  1:4696071/1‑70 (MQ=255)

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

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

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

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

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

                                         gcagGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa                                    >  1:494033/1‑63 (MQ=255)

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

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

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

                                           aaggCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCgg                          >  1:1049272/3‑71 (MQ=255)

                                           aaggCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCgg                          >  1:7948/3‑71 (MQ=255)

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

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

                                              gCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtg                       <  1:3150876/71‑1 (MQ=255)

                                               cACGCATGCCGTAAAAGCTAAAATATTTCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtgt                      <  1:1217570/71‑1 (MQ=255)

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

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

                                               cACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtgt                      <  1:4509645/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)

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

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

                                                   cATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAg                                        <  1:3572438/49‑1 (MQ=255)

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

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

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

                                                       ccGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCa                                                 >  1:4173295/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)

                                                            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)

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

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

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

                                                                   aaaTATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTa          >  1:5443618/1‑63 (MQ=255)

                                                                   aaaTATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCataat  <  1:4807160/71‑1 (MQ=255)

AGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAAT  >  Exported/454950‑455087

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