BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,876,236	A→G	62.2%	S396S (TCT→TCC)	eno ←	enolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,876,236	0	A	G	62.2%	36.3 / 40.0	66	S396S (TCT→TCC)	eno	enolase
Reads supporting (aligned to +/- strand): ref base A (14/11); new base G (28/13); total (42/24)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.44e-01

CTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCTCATAGAACCAGTTTTGATCTGGCCTGCAGCAGTACCAACAGCCAGGTCAGCGATGGTAGCGTCTTCAGT  >  minE/1876181‑1876301

cttcttCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCTCATAGAACCAGTTTTg                                                       >  1:3095835/1‑68 (MQ=255)

   cttcGATACGAATCATCTGGTTGTATTTAGCAACACGGTCAGAACGGCTCATAGAACCAGTTTTg                                                       >  1:1933303/1‑65 (MQ=255)

   cttcGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGATCGGCTCATAGAACCAGTTTTg                                                       >  1:553402/1‑65 (MQ=255)

   cttcGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCTCATAGAACCAGtttt                                                        >  1:567636/1‑64 (MQ=255)

   cttcGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCTCATAGAACCAGTTTTg                                                       >  1:3380868/1‑65 (MQ=255)

   cttcGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCTCATAGAACCAGTTTTg                                                       >  1:2005887/1‑65 (MQ=255)

   cttcGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCTCATAGAACCAGTTTTg                                                       >  1:2899749/1‑65 (MQ=255)

   cttcGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCTCATAGAACCAGTTTTg                                                       >  1:149654/1‑65 (MQ=255)

   cttcGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCTCATAGAACCAGTTTTg                                                       >  1:3219364/1‑65 (MQ=255)

   cttcGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCTCATAGAACCAGTTTTg                                                       >  1:3247069/1‑65 (MQ=255)

   cttcGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCTCATAGAACCAGTTTTg                                                       >  1:3411246/1‑65 (MQ=255)

   cttcGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCTCATAGAACCAGTTTTg                                                       >  1:402732/1‑65 (MQ=255)

   cttcGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCTCATAGAACCAGTTTTg                                                       >  1:809769/1‑65 (MQ=255)

   cttcGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCTCATAGAACCAGTTTTg                                                       >  1:1205942/1‑65 (MQ=255)

      cGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCTCATGGAACCCGttt                                                         >  1:1658322/1‑60 (MQ=255)

      cGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCTCATGGAACCCGTTTTGATCTGGc                                                >  1:507072/1‑69 (MQ=255)

      cGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCTCATGGAACCCGTTTTGATCTGGc                                                >  1:559531/1‑69 (MQ=255)

      cGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCTCATGGAACCCGTTTTGATCTGGc                                                >  1:2430942/1‑69 (MQ=255)

      cGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCTCATGGAACCCGTTTTGATCTGGc                                                >  1:2095088/1‑69 (MQ=255)

                        tGTATTTAGCAACACGGTCAGAACGTCTCATGGAACCCGTTTTGATCTGGCCTgcagc                                         <  1:497751/58‑1 (MQ=255)

                        tGTATTTAGCAACACGGTCAGAACGGCTCATGGAACCCGTTTTGATCTGGCCTgcagc                                         <  1:939837/58‑1 (MQ=255)

                        tGTATTTAGCAACACGGTCAGAACGGCTCATGGAACCCGTTTTGATCTGGCCTgcagc                                         <  1:484693/58‑1 (MQ=255)

                        tGTATTTAGCAACACGGTCAGAACGGCTCATGGAACCCGTTTTGATCTGGCCTgcagc                                         <  1:783177/58‑1 (MQ=255)

                        tGTATTTAGCAACACGGTCAGAACGGCTCATGGAACCCGTTTTGATCTGGCCTgcagc                                         <  1:2507026/58‑1 (MQ=255)

                        tGTATTTAGCAACACGGTCAGAACGGCTCATGGAACCCGTTTTGATCTGGCCTgcagc                                         <  1:2199027/58‑1 (MQ=255)

                        tGTATTTAGCAACACGGTCAGAACGGCTCATGGAACCCGTTTTGATCTGGCCTgcagc                                         <  1:1992930/58‑1 (MQ=255)

                        tGTATTTAGCAACACGGTCAGAACGGCTCATGGAACCCGTTTTGATCTGGCCTgcagc                                         <  1:121707/58‑1 (MQ=255)

                        tGTATTTAGCAACACGGTCAGAACGGCTCATGGAACCCGTTTTGATCTGGCCTgcagc                                         <  1:1331460/58‑1 (MQ=255)

                        tGTATTTAGCAACACGGTCAGAACGGCTCATGGAACCCGTTTTGATCTGGCCTgcagc                                         <  1:1444155/58‑1 (MQ=255)

                        tGTATTTAGCAACACGGTCAGAACGGCTCATGGAACCCGTTTTGATCTGGCCTgcagc                                         <  1:1460480/58‑1 (MQ=255)

                        tGTATTTAGCAACACGGTCAGAACGGCTCATGGAACCCGTTTTGATCTGGCCTgcagc                                         <  1:1798048/58‑1 (MQ=255)

                         gTATTTAGCAACACGGTCAGAACGGCTCATGGAACCCGTTTTGATCTGGCCTgcagc                                         <  1:3073974/57‑1 (MQ=255)

                                  aaCACGGTCAGAACGGCTCATAGAACCAGTTTTGATCTGGCCTGCAGCAGTACCAACAGCCAGGTCAGc                    <  1:4022/69‑1 (MQ=255)

                                   acacGGTCAGAACGGCTCATGGACCCCGTTTTGATCTGGCCTGCAGCAGTACCAACAGCCAGGTCAGCg                   >  1:2737977/1‑69 (MQ=255)

                                   acacGGTCAGAACGGCTCATGGAACCCGTTTTGATCt                                                   >  1:1856774/1‑37 (MQ=38)

                                   acacGGTCAGAACGGCTCATGGAACCCGTTTTGATCTTGCCTgcagc                                         >  1:3315811/1‑47 (MQ=255)

                                   acacGGTCAGAACGGCTCATGGAACCCGTTTTGATCTGGCCTGCAGCAGTACCAACAg                              >  1:2682299/1‑58 (MQ=255)

                                   acacGGTCAGAACGGCTCATGGAACCCGTTTTGATCTGGCCTGCAGCAGTACCAACAg                              >  1:1323875/1‑58 (MQ=255)

                                   acacGGTCAGAACGGCTCATGGAACCCGTTTTGATCTGGCCTGCAGCAGTACCAACAGc                             >  1:180291/1‑59 (MQ=255)

                                   acacGGTCAGAACGGCTCATGGAACCCGTTTTGATCTGGCCTGCAGCAGTACCAACAGCCAGGTCAgat                   >  1:2381152/1‑67 (MQ=255)

                                   acacGGTCAGAACGGCTCATGGAACCCGTTTTGATCTGGCCTGCAGCAGTACCAACAGCCAGGTCAGCg                   >  1:1130045/1‑69 (MQ=255)

                                   acacGGTCAGAACGGCTCATGGAACCCGTTTTGATCTGGCCTGCAGCAGTACCAACAGCCAGGTCAGCg                   >  1:921216/1‑69 (MQ=255)

                                   acacGGTCAGAACGGCTCATGGAACCCGTTTTGATCTGGCCTGCAGCAGTACCAACAGCCAGGTCAGCg                   >  1:816908/1‑69 (MQ=255)

                                   acacGGTCAGAACGGCTCATGGAACCCGTTTTGATCTGGCCTGCAGCAGTACCAACAGCCAGGTCAGCg                   >  1:1243008/1‑69 (MQ=255)

                                   acacGGTCAGAACGGCTCATGGAACCCGTTTTGATCTGGCCTGCAGCAGTACCAACAGCCAGGTCAGCg                   >  1:614402/1‑69 (MQ=255)

                                   acacGGTCAGAACGGCTCATGGAACCCGTTTTGATCTGGCCTGCAGCAGTACCAACAGCCAGGTCAGCg                   >  1:2281895/1‑69 (MQ=255)

                                   acacGGTCAGAACGGCTCATGGAACCCGTTTTGATCTGGCCTGCAGCAGTACCAACAGCCAGGTCAGCg                   >  1:136685/1‑69 (MQ=255)

                                   acacGGTCAGAACGGCTCATGGAACCCGTTTTGATCTGGCCTGCAGCAGTACCAACAGCCAGGTCAGCg                   >  1:137127/1‑69 (MQ=255)

                                   acacGGTCAGAACGGCTCATGGAACCCGTTTTGATCTGGCCTGCAGCAGTACCAACAGCCAGGTCAGCg                   >  1:354483/1‑69 (MQ=255)

                                   acacGGTCAGAACGGCTCATGGAACCCGTTTTGATCTGGCCTGCAGCAGTACCAACAGCCAGGTCAGCg                   >  1:2657600/1‑69 (MQ=255)

                                   acacGGTCAGAACGGCTCATGGAACCCGTTTTGATCTGGCCTGCAGCAGTACCAACAGCCAGGTCAGCg                   >  1:2619201/1‑69 (MQ=255)

                                   acacGGTCAGAACGGCTCATGGAACCCGTTTTGATCTGGCCTGCAGCAGTACCAACAGCCAGGTCAGCg                   >  1:1822658/1‑69 (MQ=255)

                                   acacGGTCAGAACGGCTCATGGAACCCGTTTTGATCTGGCCTGCAGCAGTACCAACAGCCAGGTCAGCg                   >  1:2003816/1‑69 (MQ=255)

                                   acacGGTCAGAACGGCTCATGGAACCCGTTTTGATCTGGCCTGCAGCAGTACCAACAGCCAGGTCAGCg                   >  1:32952/1‑69 (MQ=255)

                                   acacGGTCAGAACGGCTCATGGAACCCGTTTTGATCTGGCCTGCAGCAGTACCAACAGCCAGGTCAGCg                   >  1:1486137/1‑69 (MQ=255)

                                   acacGGTCAGAACGGCTCATGGAACCCGTTTTGATCTGGCCTGCAGCAGTACCAACAGCCAGGTCAGCg                   >  1:3125726/1‑69 (MQ=255)

                                             aaCGGCTCATAGAACCAGTTTTGATCTGGCCTGCAGCAGTACCAACAGCCAGGTCAGCGATGGTAGCg          <  1:2384703/68‑1 (MQ=255)

                                             aaCGGCTCATAGAACCAGTTTTGATCTGGCCTGCAGCAGTACCAACAGCCAGGTCAGCGATGGTAGCg          <  1:1819945/68‑1 (MQ=255)

                                             aaCGGCTCATAGAACCAGTTTTGATCTGGCCTGCAGCAGTACCAACAGCCAGGTCAGCGATGGTAGCg          <  1:348800/68‑1 (MQ=255)

                                                    cATAGAACCAGTTTTGATCTGGCCTGCAGCAGTACCAACAGCCAGGTCAGCGATGGTAGCGTCTTCAGt  <  1:1492829/69‑1 (MQ=255)

                                                    cATAGAACCAGTTTTGATCTGGCCTGCAGCAGTACCAACAGCCAGGTCAGCGATGGTAGCGTCTTCAGt  <  1:74183/69‑1 (MQ=255)

                                                    cATAGAACCAGTTTTGATCTGGCCTGCAGCAGTACCAACAGCCAGGTCAGCGATGGTAGCGTCTTCAGt  <  1:77906/69‑1 (MQ=255)

                                                    cATAGAACCAGTTTTGATCTGGCCTGCAGCAGTACCAACAGCCAGGTCAGCGATGGTAGCGTCTTCAGt  <  1:1318754/69‑1 (MQ=255)

                                                    cATAGAACCAGTTTTGATCTGGCCTGCAGCAGTACCAACAGCCAGGTCAGCGATGGTAGCGTCTTCAGt  <  1:1954126/69‑1 (MQ=255)

                                                    cATAGAACCAGTTTTGATCTGGCCTGCAGCAGTACCAACAGCCAGGTCAGCGATGGTAGCGTCTTCAGt  <  1:1443420/69‑1 (MQ=255)

                                                    cATAGAACCAGTTTTGATCTGGCCTGCAGCAGTACCAACAGCCAGGTCAGCGATGGTAGCGTCTTCAGt  <  1:1446800/69‑1 (MQ=255)

CTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCTCATAGAACCAGTTTTGATCTGGCCTGCAGCAGTACCAACAGCCAGGTCAGCGATGGTAGCGTCTTCAGT  >  minE/1876181‑1876301

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