BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,563,249	A→G	45.2%	R276R (CGT→CGC)	lldP ←	L‑lactate permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,563,249	0	A	G	45.2%	6.5 / 49.9	62	R276R (CGT→CGC)	lldP	L‑lactate permease
Reads supporting (aligned to +/- strand): ref base A (17/17); new base G (17/11); total (34/28)
Fisher's exact test for biased strand distribution p-value = 4.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.99e-01

TGACCCGCAGTGTAACCGGTGTGGGCCAGCGTCATATCAACCTGTGACGCCCCCAAATCACCAAAACGGAATACACGCACTGGCTGCCAGCGTTTGAGGAACAGCGTCAGGCAGAGCAGTGAT  >  minE/2563184‑2563306

tGACCCGCAGTGTAACCGGTGTGGGCCAGCGTCATATCAACCTGTGACGCCCCCAAATCACCAAAACgg                                                        <  1:2182380/69‑1 (MQ=255)

tGACCCGCAGTGTAACCGGTGTGGGCCAGCGTCATATCAACCTGTGACGCCCCCAAATCACCAAAACgg                                                        <  1:1404918/69‑1 (MQ=255)

tGACCCGCAGTGTAACCGGTGTGGGCCAGCGTCATATCAACCTGTGACGCCCCCAAATCACCAAAACgg                                                        <  1:591845/69‑1 (MQ=255)

    ccGCAGTGTAACCGGTGTGGGCCAGCGTCATATCAACCTGTGACGCCCCCAAATCACCAAAACGGaa                                                      <  1:1024073/67‑1 (MQ=255)

    ccGCAGTGTAACCGGTGTGGGCCAGCGTCATATCAACCTGTGACGCCCCCAAATCACCAAAACGGaa                                                      <  1:2901773/67‑1 (MQ=255)

    ccGCAGTGTAACCGGTGTGGGCCAGCGTCATATCAACCTGTGACGCCCCCAAATCACCAAAACGGaa                                                      <  1:3333508/67‑1 (MQ=255)

    ccGCAGTGTAACCGGTGTGGGCCAGCGTCATATCAACCTGTGACGCCCCCAAATCACCAAAACGGaa                                                      <  1:480476/67‑1 (MQ=255)

    ccGCAGTGTAACCGGTGTGGGCCAGCGTCATATCAACCTGTGACGCCCCCAAATCACCAAAACGGaa                                                      <  1:799091/67‑1 (MQ=255)

    ccGCAGTGTAACCGGTGTGGGCCAGCGTCATATCAACCTGTGACGCCCCAAAATCACCAAAACGGaa                                                      <  1:155850/67‑1 (MQ=255)

        aGTGTAACCGGTGTGGGCCAGCGTCATATCAACCTGTGACGCCCCCAAATCACCAAAACGGAATACACg                                                <  1:1594524/69‑1 (MQ=255)

        aGTGTAACCGGTGTGGGCCAGCGTCATATCAACCTGTGACGCCCCCAAATCACCAAAACGGAATACACg                                                <  1:1751823/69‑1 (MQ=255)

            tAACCGGTGTGGGCCAGCGTCATATCAACCTGTGACGCCCCCAAATCACCAAAACGGAATACACGCACt                                            >  1:3168196/1‑69 (MQ=255)

            tAACCGGTGTGGGCCAGCGTCATATCAACCTGTGACGCCCCCAAATCACCAAAACGGAATACACGCACt                                            >  1:293282/1‑69 (MQ=255)

            tAACCGGTGTGGGCCAGCGTCATATCAACCTGTGACGCCCCCAAATCACCAAAACGGAATACACGCACt                                            >  1:2679508/1‑69 (MQ=255)

            tAACCGGTGTGGGCCAGCGTCATATCAACCTGTGACGCCCCCAAATCACCAAAACGGAATACACGCACt                                            >  1:2390762/1‑69 (MQ=255)

            tAACCGGTGTGGGCCAGCGTCATATCAACCTGTGACGCCCCCAAATCACCAAAACGGAATACACGCACt                                            >  1:2670049/1‑69 (MQ=255)

            tAACCGGTGTGGGCCAGCGTCATATCAACCTGTGACGCCCCCAAATCACCAAAACGGAATACACGCACt                                            >  1:1678735/1‑69 (MQ=255)

            tAACCGGTGTGGGCCAGCGTCATATCAACCTGTGACGCCACCAAATCACCAAAACGGAATACACGCACt                                            >  1:1537212/1‑69 (MQ=255)

                   tgtgGGCCAGCGTCATATCAACCTGTGACGCCCCCAAATCACCAAAACGGAATACACGCACTGGCTGcc                                     >  1:768093/1‑69 (MQ=255)

                   tgtgGGCCAGCGTCATATCAACCTGTGACGCCCCCAAATCACCAAAACGGAATACACGCACTGGCTGc                                      >  1:2007057/1‑68 (MQ=255)

                   tgtgAGCCAGCGTCATATCAACCTGTGACGCCCCTAAGTCACCGAAGCGGAATACGCGCACTgg                                          <  1:1492246/64‑1 (MQ=255)

                   tgtgAGCCAGCGTCATATCAACCTGTGACGCCCCTAAGTCACCGAAGCGGAATACGCGCACTgg                                          <  1:1016775/64‑1 (MQ=255)

                   tgtgAGCCAGCGTCATATCAACCTGTGACGCCCCTAAGTCACCGAAGCGGAATACGCGCACTgg                                          <  1:1522772/64‑1 (MQ=255)

                   tgtgAGCCAGCGTCATATCAACCTGTGACGCCCCTAAGTCACCGAAGCGGAATACGCGCACTgg                                          <  1:640991/64‑1 (MQ=255)

                   tgtgAGCCAGCGTCATATCAACCTGTGACGCCCCTAAGTCACCGAAGCGGAATACGCGCACTgg                                          <  1:165382/64‑1 (MQ=255)

                   tgtgAGCCAGCGTCATATCAACCTGTGACGCCCCTAAGTCACCGAAGCGGAATACGCGCACTgg                                          <  1:1147497/64‑1 (MQ=255)

                   tgtgAGCCAGCGTCATATCAACCTGTGACGCCCCTAAGTCACCGAAGCGGAATACGCGCACTgg                                          <  1:181453/64‑1 (MQ=255)

                   tgtgAGCCAGCGTCATATCAACCTGTGACGCCCCTAAGTCACCGAAGCGGAATACGCGCACTgg                                          <  1:3368507/64‑1 (MQ=255)

                   tgtgAGCCAGCGTCATATCAACCTGTGACGCCCCTAAGTCACCGAAGCGGAATACGCGCACTgg                                          <  1:803413/64‑1 (MQ=255)

                   tgtgAGCCAGCGTCATATCAACCTGTGACGCCCCTAAGTCACCGAAGCGGAATACGCGCACTgg                                          <  1:96077/64‑1 (MQ=255)

                    gtgAGCCAGCGTCATATCAACCTGTGACGCCCCTAAGTCACCGAAGCGGAATACGCGCACTgg                                          <  1:2084432/63‑1 (MQ=255)

                              gTCATATCAACCTGTGACGCCCCCAAATCACCAAAACGGAATACACGCACTGGCTGCCAGCGTTTGa                            >  1:657541/1‑67 (MQ=255)

                              gTCATATCAACCTGTGACGCCCCCAAATCACCAAAACGGAATACACGCACTGGCTGCCAGCGTTTGAgg                          >  1:2922130/1‑69 (MQ=255)

                              gTCATATCAACCTGTGACGCCCCCAAATCACCAAAACGGAATACACGCACTGGCTGCCAGCGTTTGAgg                          <  1:2671139/69‑1 (MQ=255)

                                    tCAACCTGTGACGCCCCCAAATCACCAAAACGGAATACACGCACTGGCTGCCAGCGTTTGAGGAACAg                     >  1:385708/1‑68 (MQ=255)

                                    tCAACCTGTGACGCCCCCAAATCACCAAAACGGAATACACGCACTGGCTGCCAGCGTTTGAGGAACAGc                    >  1:1304369/1‑69 (MQ=255)

                                    tCAACCTGTGACGCCCCCAAATCACCAAAACGGAATACACGCACTGGCTGCCAGCGTTTGAGGAACAGc                    >  1:1240080/1‑69 (MQ=255)

                                      aaCCTGTGACGCCCCCAAATCACCAAAACGGAATACACGCACTGGCTGCCAGCGTTTGAGGAACAGc                    >  1:2264526/1‑67 (MQ=255)

                                      aaCCTGTGACGCCCCCAAATCACCAAAACGGAATACACGCACTGGCTGCCAGCGTTTGAGGAACAGc                    >  1:1932790/1‑67 (MQ=255)

                                              aCGCCCCCAAATCACCAAAACGGAATACACGCACTGGCTGCCAGCGTTTGAGGAACAGCGTCAGGCag           >  1:2265145/1‑68 (MQ=255)

                                                  ccctAAGTCACCGAAGCGGAATACGCGCACTGGTTGCCAGCGTTTGAGGAACAGTGTCAGACAGAGCa       >  1:750828/1‑68 (MQ=255)

                                                  ccctAAGTCACCGAAGCGGAATACGCGCACTGGTTGCCAGCGTTTGAGGAACAGTGTCAGACAGAGCa       >  1:731090/1‑68 (MQ=255)

                                                  ccctAAGTCACCGAAGCGGAATACGCGCACTGGTTGCCAGCGTTTGAGGAACAGTGTCAGACAGAGCa       >  1:297779/1‑68 (MQ=255)

                                                  ccctAAGTCACCGAAGCGGAATACGCGCACTGGTTGCCAGCGTTTGAGGAACAGTGTCAGACAGAGCa       >  1:754379/1‑68 (MQ=255)

                                                  ccctAAGTCACCGAAGCGGAATACGCGCACTGGTTGCCAGCGTTTGAGGAACAGTGTCAGACAGAGCa       >  1:511429/1‑68 (MQ=255)

                                                  ccctAAGTCACCGAAGCGGAATACGCGCACTGGTTGCCAGCGTTTGAGGAACAGTGTCAGACAGAGCa       >  1:496878/1‑68 (MQ=255)

                                                  ccctAAGTCACCGAAGCGGAATACGCGCACTGGTTGCCAGCGTTTGAGGAACAGTGTCAGACAGAGCa       >  1:1144460/1‑68 (MQ=255)

                                                  ccctAAGTCACCGAAGCGGAATACGCGCACTGGTTGCCAGCGTTTGAGGAACAGTGTCAGACAGAGCa       >  1:3354651/1‑68 (MQ=255)

                                                  ccctAAGTCACCGAAGCGGAATACGCGCACTGGTTGCCAGCGTTTGAGGAACAGTGTCAGACAGAGCa       >  1:2995817/1‑68 (MQ=255)

                                                  ccctAAGTCACCGAAGCGGAATACGCGCACTGGTTGCCAGCGTTTGAGGAACAGTGTCAGACAGAGCa       >  1:295479/1‑68 (MQ=255)

                                                  ccctAAGTCACCGAAGCGGAATACGCGCACTGGTTGCCAGCGTTTGAGGAACAGTGTCAGACAGAGCa       >  1:131391/1‑68 (MQ=255)

                                                  ccctAAGTCACCGAAGCGGAATACGCGCACTGGTTGCCAGCGTTTGAGGAACAGTGTCAGACAGAGCa       >  1:137612/1‑68 (MQ=255)

                                                  ccctAAGTCACCGAAGCGGAATACGCGCACTGGTTGCCAGCGTTTGAGGAACAGTGTCAGACAGAGCa       >  1:2875769/1‑68 (MQ=255)

                                                  ccctAAGTCACCGAAGCGGAATACGCGCACTGGTTGCCAGCGTTTGAGGAACAGTGTCAGACAGAGCa       >  1:2595515/1‑68 (MQ=255)

                                                  ccctAAGTCACCGAAGCGGAATACGCGCACTGGTTGCCAGCGTTTGAGGAACAGTGTCAGACAGAGCa       >  1:2132249/1‑68 (MQ=255)

                                                  ccctAAGTCACCGAAGCGGAATACGCGCACTGGTTGCCAGCGTTTGAGGAACAGTGTCAGACAGAGCa       >  1:1599262/1‑68 (MQ=255)

                                                  ccctAAGTCACCGAAGCGGAATACGCGCACTGGTTGCCAGCGTTTGAGGAACAGTGTCAGACAGAGCa       >  1:1715631/1‑68 (MQ=255)

                                                      aaaTCACCAAAACGGAATACACGCACTGGCTGCCAGCGTTTGAGGAACAGCGTCAGGCAGAGCAGTGAt  <  1:1260014/69‑1 (MQ=255)

                                                      aaaTCACCAAAACGGAATACACGCACTGGCTGCCAGCGTTTGAGGAACAGCGTCAGGCAGAGCAGTGAt  <  1:567478/69‑1 (MQ=255)

                                                      aaaTCACCAAAACGGAATACACGCACTGGCTGCCAGCGTTTGAGGAACAGCGTCAGGCAGAGCAGTGAt  <  1:630068/69‑1 (MQ=255)

                                                      aaaTCACCAAAACGGAATACACGCACTGGCTGCCAGCGTTTGAGGAACAGCGTCAGGCAGAGCAGTGAt  <  1:855016/69‑1 (MQ=255)

                                                      aaaTCACCAAAACGGAATACACGCACTGGCTGCCAGCGTTTGAGGAACAGCGTCAGGCAGAGCAGTGAt  <  1:1635502/69‑1 (MQ=255)

TGACCCGCAGTGTAACCGGTGTGGGCCAGCGTCATATCAACCTGTGACGCCCCCAAATCACCAAAACGGAATACACGCACTGGCTGCCAGCGTTTGAGGAACAGCGTCAGGCAGAGCAGTGAT  >  minE/2563184‑2563306

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 33 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`