BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,505,809	C→A	12.7%	G510G (GGG→GGT)	yfbS ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,505,809	0	C	A	12.7%	136.0 / 6.5	62	G510G (GGG→GGT)	yfbS	predicted transporter
Reads supporting (aligned to +/- strand): ref base C (28/26); new base A (5/3); total (33/29)
Fisher's exact test for biased strand distribution p-value = 7.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

AACACCGCCGTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCATCAGCCCTTTCACCGCCAGCGC  >  minE/1505747‑1505871

aaCACCGCCGTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGc                                                          >  1:985875/1‑69 (MQ=255)

aaCACCGCCGTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGc                                                          >  1:2361729/1‑69 (MQ=255)

aaCACCGCCGTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGc                                                          >  1:3166209/1‑69 (MQ=255)

aaCACCGCCGTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGc                                                          >  1:403168/1‑69 (MQ=255)

    ccgccgTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGcatca                                                      >  1:896376/1‑69 (MQ=255)

    ccgccgTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGcatca                                                      >  1:1496723/1‑69 (MQ=255)

    ccgccgTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGcatca                                                      <  1:592760/69‑1 (MQ=255)

    ccgccgTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGcatca                                                      >  1:3315647/1‑69 (MQ=255)

    ccgccgTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGcatca                                                      <  1:2889013/69‑1 (MQ=255)

    ccgccgTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGcatca                                                      <  1:2610844/69‑1 (MQ=255)

      gccgTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCata                                                    <  1:858879/69‑1 (MQ=255)

      gccgTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCata                                                    <  1:1407476/69‑1 (MQ=255)

      gccgTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCata                                                    <  1:1722437/69‑1 (MQ=255)

      gccgTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGAAAAGACAAACAAACACCCCAGCATCata                                                    <  1:1507605/69‑1 (MQ=255)

       ccgTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCata                                                    <  1:941142/68‑1 (MQ=255)

                tATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGt                                          <  1:1416547/69‑1 (MQ=255)

                tATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGt                                          <  1:2435173/69‑1 (MQ=255)

                 aTTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTa                                         <  1:153997/69‑1 (MQ=255)

                 aTTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTa                                         <  1:1891557/69‑1 (MQ=255)

                  ttAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTa                                         <  1:1102302/68‑1 (MQ=255)

                    agagATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGc                                              >  1:1085710/1‑61 (MQ=255)

                    agagATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAacc                                      >  1:2598292/1‑69 (MQ=255)

                    agagATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAacc                                      >  1:1505080/1‑69 (MQ=255)

                       gaTAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCata                                                    <  1:1015040/52‑1 (MQ=255)

                          aaaTAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGcc                                <  1:1738058/69‑1 (MQ=255)

                          aaaTAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGcc                                <  1:2917957/69‑1 (MQ=255)

                           aaTAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGcc                                <  1:1571196/68‑1 (MQ=255)

                             tAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAAt                             >  1:2951204/1‑69 (MQ=255)

                              aGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATg                            <  1:1335371/69‑1 (MQ=255)

                                    aTAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTa                                         <  1:1393920/50‑1 (MQ=255)

                                    aTAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTa                                         <  1:3197563/50‑1 (MQ=255)

                                    aTAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTa                                         <  1:2720259/50‑1 (MQ=255)

                                    aTAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTa                                         <  1:1131323/50‑1 (MQ=255)

                                      aaCCGCCGACAAGACAAACAAACACCCCAGCATCata                                                    <  1:3273754/37‑1 (MQ=255)

                                      aaCCGCCGACAAGACAAACAAACACCCCAGCATCata                                                    <  1:1908795/37‑1 (MQ=255)

                                      aaCCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCa                        >  1:2069350/1‑65 (MQ=255)

                                      aaCCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCATCAg                    >  1:648902/1‑69 (MQ=255)

                                      aaCCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCATCAg                    >  1:504616/1‑69 (MQ=255)

                                      aaCCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCATCAg                    >  1:421036/1‑69 (MQ=255)

                                      aaCCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCATCAg                    >  1:713820/1‑69 (MQ=255)

                                      aaCCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCATCAg                    >  1:1030110/1‑69 (MQ=255)

                                      aaCCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCATCAg                    >  1:323409/1‑69 (MQ=255)

                                      aaCCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCATCAg                    >  1:781215/1‑69 (MQ=255)

                                      aaCCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCATCAg                    >  1:2050538/1‑69 (MQ=255)

                                      aaCCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCATCAg                    >  1:2642875/1‑69 (MQ=255)

                                      aaCCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCATCAg                    >  1:196487/1‑69 (MQ=255)

                                      aaCCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCATCAg                    >  1:2240270/1‑69 (MQ=255)

                                      aaCCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCATCAg                    >  1:1955936/1‑69 (MQ=255)

                                      aaCCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCATCAg                    >  1:2521224/1‑69 (MQ=255)

                                      aaCCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCATCAg                    >  1:2672126/1‑69 (MQ=255)

                                      aaCCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCATCAg                    >  1:1815605/1‑69 (MQ=255)

                                      aaCCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCATAAg                    >  1:880651/1‑69 (MQ=255)

                                              acaagacaaACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCATCAGCCCTTTc             <  1:2837577/68‑1 (MQ=255)

                                                  gacaaACAAACAACCCAGCATCATATGTGGCCCGTAACCGCCGCCAATGTCCATCAGCCCTTTCACCGc        <  1:1551177/69‑1 (MQ=255)

                                                  gacaaACAAACAACCCAGCATCATATGTGGCCCGTAACCGCCGCCAATGTCCATCAGCCCTTTCACCGc        <  1:1829722/69‑1 (MQ=255)

                                                  gacaaACAAACAACCCAGCATCATATGTGGCCCGTAACCGCCGCCAATGTCCATCAGCCCTTTCACCGc        <  1:2158003/69‑1 (MQ=255)

                                                   acaaacaaacaaCCCAGCATCATATGTGGCCCGTAACCGCCGCCAATGTCCATCAGCCCTTTCACCGc        <  1:468041/68‑1 (MQ=255)

                                                    caaacaaaCAACCCAGCATCATATGTGGCCCGTAACCGCCGCCAATGTCCATCAGCCCTTTCACCGCCa      >  1:2798308/1‑69 (MQ=255)

                                                    caaacaaaCAACCCAGCATCATATGTGGCCCGTAACCGCCGCCAATGTCCATCAGCCCTTTCACCGCCa      >  1:695049/1‑69 (MQ=255)

                                                    caaacaaaCAACCCAGCATCATATGTGGCCCGTAACCGCCGCCAATGTCCATCAGCCCTTTCACCGCCa      >  1:2944992/1‑69 (MQ=255)

                                                    caaacaaaCAACCCAGCATCATATGTGGCCCGTAACCGCCGCCAATGTCCATCAGCCCTTTCACCGCCa      >  1:2139193/1‑69 (MQ=255)

                                                    caaacaaaCAACCCAGCATCATATGTGGCCCGTAACCGCCGCCAATGTCCATCAGCCCTTTCACCGCCa      >  1:2982651/1‑69 (MQ=255)

                                                        caaacaCCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCATCAGCCCTTTCACCGCCAgcgc  <  1:3138759/69‑1 (MQ=255)

AACACCGCCGTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCATCAGCCCTTTCACCGCCAGCGC  >  minE/1505747‑1505871

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