BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	454,992	C→T	16.2%	L31L (CTG→CTA) ‡	ybgH ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	454,992	0	C	T	16.2%	162.2 / 21.6	68	L31L (CTG→CTA) ‡	ybgH	predicted transporter
Reads supporting (aligned to +/- strand): ref base C (26/31); new base T (3/8); total (29/39)
Fisher's exact test for biased strand distribution p-value = 3.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCG  >  minE/454923‑455060

gTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATcagc                                                                      <  1:1765773/70‑1 (MQ=255)

      gcTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCAc                                                               <  1:2538141/71‑1 (MQ=255)

       cTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACg                                                              <  1:2607543/71‑1 (MQ=255)

          aaTAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCAt                                                           <  1:139096/71‑1 (MQ=255)

           aTAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCAt                                                           <  1:1014072/70‑1 (MQ=255)

           aTAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATg                                                          <  1:2499243/71‑1 (MQ=255)

            tAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGc                                                         >  1:2646086/1‑71 (MQ=255)

            tAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGc                                                         <  1:754617/71‑1 (MQ=255)

              aCTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCg                                                       <  1:2713603/71‑1 (MQ=255)

               cTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCg                                                       >  1:3646856/1‑70 (MQ=255)

                  gTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTaa                                                    <  1:1299840/70‑1 (MQ=255)

                  gTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTaa                                                    >  1:598975/1‑70 (MQ=255)

                  gTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTaa                                                    >  1:2412742/1‑70 (MQ=255)

                  gTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTaa                                                    >  1:2358149/1‑70 (MQ=255)

                  gTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTaa                                                    >  1:594306/1‑70 (MQ=255)

                  gTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTaa                                                    <  1:2143446/70‑1 (MQ=255)

                  gTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTaa                                                    >  1:444932/1‑70 (MQ=255)

                  gTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTaa                                                    >  1:459620/1‑70 (MQ=255)

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

                        gTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTa                                              >  1:521910/1‑70 (MQ=255)

                         tGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTaa                                             >  1:872901/1‑70 (MQ=255)

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

                            gTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGt                                                      <  1:3472690/58‑1 (MQ=255)

                            gTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAata                                         <  1:3885647/71‑1 (MQ=255)

                            gTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAata                                         >  1:1974156/1‑71 (MQ=255)

                             tATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGt                                                      >  1:3491445/1‑57 (MQ=255)

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

                             tATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAatat                                        >  1:808394/1‑71 (MQ=255)

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

                               tCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCAt                                                           >  1:3520769/1‑50 (MQ=255)

                                cGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTc                                      >  1:3907370/1‑70 (MQ=255)

                                cGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTc                                      >  1:3902917/1‑70 (MQ=255)

                                       ttttAGTTGATTGGTGAGATAGAGAATCAGCAggg                                                                  <  1:2219114/35‑1 (MQ=255)

                                       ttttAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGAttt                              >  1:3281715/1‑71 (MQ=255)

                                       ttttAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGAttt                              <  1:3918520/71‑1 (MQ=255)

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

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

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

                                          tAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCa                                                            >  1:2754271/1‑38 (MQ=255)

                                            gTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTaaaa                                           >  1:1556433/1‑53 (MQ=255)

                                                atgggTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGAtcgc                       <  1:2110966/66‑1 (MQ=255)

                                                atgggTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGAtcgc                       <  1:594407/66‑1 (MQ=255)

                                                 ttGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCaa                     >  1:3571166/1‑70 (MQ=255)

                                                  tGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCaa                     >  1:2071687/1‑69 (MQ=255)

                                                  gggTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTg                             >  1:3637716/2‑61 (MQ=255)

                                                   ggTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACg                   <  1:3413146/70‑1 (MQ=255)

                                                   ggTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACg                   <  1:1588269/70‑1 (MQ=255)

                                                     tGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGata                <  1:479848/71‑3 (MQ=255)

                                                      gagaTAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACa                >  1:2183757/1‑70 (MQ=255)

                                                       agaTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGAtata              <  1:1031011/71‑1 (MQ=255)

                                                         aTAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGc                            >  1:2157446/1‑55 (MQ=255)

                                                         aTAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGt            <  1:2289002/71‑1 (MQ=255)

                                                          tAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTa           >  1:1461850/1‑71 (MQ=255)

                                                          tAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAata                                         <  1:2343091/41‑1 (MQ=255)

                                                          tAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAata                                         <  1:47862/41‑1 (MQ=255)

                                                          tAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAata                                         <  1:43173/41‑1 (MQ=255)

                                                          tAGAGAATCAGCAGGACACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTa           <  1:220749/71‑1 (MQ=255)

                                                           agagAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAgcgc                       <  1:2383786/58‑1 (MQ=255)

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

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

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

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

                                                             agaATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGAtata              <  1:3048127/65‑1 (MQ=255)

                                                             agaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATAtt                                       >  1:3915264/1‑40 (MQ=255)

                                                             agaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATAtt                                       >  1:2882151/1‑40 (MQ=255)

                                                             agaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATAtt                                       >  1:1788054/1‑40 (MQ=255)

                                                              gaATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTatatc       <  1:3408650/71‑2 (MQ=255)

                                                              gaATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTatatc       >  1:665482/1‑70 (MQ=255)

                                                              gaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAAt        <  1:1390922/70‑1 (MQ=255)

                                                                   agcagGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAg                          >  1:2611027/1‑47 (MQ=255)

                                                                   agcagGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAgcgcg  <  1:954568/71‑1 (MQ=255)

GTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCG  >  minE/454923‑455060

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