BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,563,246	A→G	45.2%	F277F (TTT→TTC)	lldP ←	L‑lactate permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,563,246	0	A	G	45.2%	6.7 / 50.2	62	F277F (TTT→TTC)	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 = 1.00e+00

ACACCTGACCCGCAGTGTAACCGGTGTGGGCCAGCGTCATATCAACCTGTGACGCCCCCAAATCACCAAAACGGAATACACGCACTGGCTGCCAGCGTTTGAGGAACAGCGTCAGGCAGAGCAGTGAT  >  minE/2563179‑2563306

acacCTGACCCGCAGTGTAACCGGTGTGGGCCAGCGTCATATCAACCTGTGACGCCCCCAAATCACCaa                                                             >  1:696899/1‑69 (MQ=255)

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

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

     tGACCCGCAGTGTAACCGGTGTGGGCCAGCGTCATATCAACCTGTGACGCCCCCAAATCACCAAAACgg                                                        <  1:1404918/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:1522772/64‑1 (MQ=255)

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

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

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

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

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

                        tgtgAGCCAGCGTCATATCAACCTGTGACGCCCCTAAGTCACCGAAGCGGAATACGCGCACTgg                                          <  1:181453/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:2995817/1‑68 (MQ=255)

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

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

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

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

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

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

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

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

                                                       ccctAAGTCACCGAAGCGGAATACGCGCACTGGTTGCCAGCGTTTGAGGAACAGTGTCAGACAGAGCa       >  1:295479/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:1599262/1‑68 (MQ=255)

                                                       ccctAAGTCACCGAAGCGGAATACGCGCACTGGTTGCCAGCGTTTGAGGAACAGTGTCAGACAGAGCa       >  1:2132249/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:1635502/69‑1 (MQ=255)

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

ACACCTGACCCGCAGTGTAACCGGTGTGGGCCAGCGTCATATCAACCTGTGACGCCCCCAAATCACCAAAACGGAATACACGCACTGGCTGCCAGCGTTTGAGGAACAGCGTCAGGCAGAGCAGTGAT  >  minE/2563179‑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: `‑`