BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,665,980	T→C	42.4%	Q259Q (CAA→CAG)	guaB ←	IMP dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,665,980	0	T	C	42.4%	13.2 / 43.0	59	Q259Q (CAA→CAG)	guaB	IMP dehydrogenase
Reads supporting (aligned to +/- strand): ref base T (7/27); new base C (8/17); total (15/44)
Fisher's exact test for biased strand distribution p-value = 3.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.19e-01

GCTGTTGCCACGTTGCCGCCGATAATTTGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGCAGTACACCTTCTGAGTGACCGTGGGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAC  >  minE/1665915‑1666045

gCTGTTGCCACGTTGCCGCCGATAATTTGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGc                                                                <  1:161073/69‑1 (MQ=255)

 cTGTTGCCACGTTGCCGCCGATAATTTGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGCa                                                               <  1:1222750/69‑1 (MQ=255)

 cTGTTGCCACGTTGCCGCCGATAATTTGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGCa                                                               <  1:2497408/69‑1 (MQ=255)

 cTGTTGCCACGTTGCCGCCGATAATTTGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGCa                                                               <  1:628635/69‑1 (MQ=255)

 cTGTTGCCACGTTGCCGCCGATAATTTGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGCa                                                               <  1:3305362/69‑1 (MQ=255)

 cTGTTGCCACGTTGCCGCCGATAATTTGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGCa                                                               <  1:2588796/69‑1 (MQ=255)

 cTGTTGCCACGTTGCCGCCGATAATTTGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGCa                                                               <  1:1404789/69‑1 (MQ=255)

 cTGTTGCCACGTTGCCGCCGATAATTTGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGCa                                                               <  1:2343865/69‑1 (MQ=255)

 cTGTTGCCACGTTGCCGCCGATAATTTGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGCa                                                               <  1:1970578/69‑1 (MQ=255)

 cTGTTGCCACGTTGCCGCCGATAATTTGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGCa                                                               <  1:2163986/69‑1 (MQ=255)

 cTGTTGCCACGTTGCCGCCGATAATTTGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGCa                                                               <  1:1762655/69‑1 (MQ=255)

 cTGTTGCCACGTTGCCGCCGATAATTTGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGCa                                                               <  1:3373365/69‑1 (MQ=255)

 cTGTTGCCACGGTGCCGCCGATAATTTGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGCa                                                               <  1:1888529/69‑1 (MQ=255)

  tGTTGCCACGTTGCCGCCGATAATTTGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGCAg                                                              <  1:19303/69‑1 (MQ=255)

  tGTTGCCACGTTGCCGCCGATAATTTGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGCAg                                                              <  1:3310222/69‑1 (MQ=255)

  tGTTGCCACGTTGCCGCCGATAATTTGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGCAg                                                              <  1:1354960/69‑1 (MQ=255)

  tGTTGCCACGTTGCCGCCGATAATTTGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGCAg                                                              <  1:2011805/69‑1 (MQ=255)

   gTTGCCACGTTGCCGCCGATAATTTGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGCAg                                                              <  1:355172/68‑1 (MQ=255)

      gCCACGTTGCCGCCGATAATTTGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGCAGTa                                                            <  1:2929705/67‑1 (MQ=255)

      gCCACGTTGCCGCCGATAATTTGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGCAGTa                                                            <  1:2545189/67‑1 (MQ=255)

      gCCACGTTGCCGCCGATAATTTGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGCAGTa                                                            <  1:2510776/67‑1 (MQ=255)

      gCCACGTTGCCGCCGATAATTTGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGCAGTa                                                            <  1:1732689/67‑1 (MQ=255)

      gCCACGTTGCCGCCGATAATTTGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGCAGTa                                                            <  1:620743/67‑1 (MQ=255)

      gCCACGTTGCCGCCGATAATTTGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGCAGTa                                                            <  1:840657/67‑1 (MQ=255)

                       gaTTTGCAGGTCAGGATATTTAGCGCGAGTTTCACGGATACGCTGCAGAACACCTTCAGAGTGAcc                                            >  1:3211438/2‑66 (MQ=255)

                       gaTTTGCAGGTCAGGATATTTAGCGCGAGTTTCACGGATACGCTGCAGAACACCTTCAGAGTGAcc                                            >  1:2878472/2‑66 (MQ=255)

                       gaTTTGCAGGTCAGGATATTTAGCGCGAGTTTCACGGATACGCTGCAGAACACCTTCAGAGTGAcc                                            >  1:485990/2‑66 (MQ=255)

                       gaTTTGCAGGTCAGGATATTTAGCGCGAGTTTCACGGATACGCTGCAGAACACCTTCAGAGTGAcc                                            >  1:858483/2‑66 (MQ=255)

                       gaTTTGCAGGTCAGGATATTTAGCGCGAGTTTCACGGATACGCTGCAGAACACCTTCAGAGTGAcc                                            >  1:997326/2‑66 (MQ=255)

                       gaTTTGCAGGTCAGGATATTTAGCGCGAGTTTCACGGATACGCTGCAGAACACCTTCAGAGTGAcc                                            >  1:1356041/2‑66 (MQ=255)

                       gaTTTGCAGGTCAGGATATTTAGCGCGAGTTTCACGGATACGCTGCAGAACACCTTCAGAGTGAcc                                            >  1:1704295/2‑66 (MQ=255)

                       gaTTTGCAGGTCAGGATATTTAGCGCGAGTTTCACGGATACGCTGCAGAACACCTTCAGAGTGAcc                                            >  1:1834286/2‑66 (MQ=255)

                          ttGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGCAGTACACCTTCTGAGTGACCGTggg                                       <  1:2114374/68‑1 (MQ=255)

                                  ccGGATATTTAGCACGGGTTTCACGGATACGTTGCAGTACACCTTCTGAGTGACCGTGGGAGGAGTCg                               >  1:2786754/1‑68 (MQ=255)

                                         tttAGCACGGGTTTCACGGATACGTTGCAGTACACCTTCTGAGTGACCGTGGGAGGAGTCGATcaggag                       >  1:2963779/1‑69 (MQ=255)

                                         tttAGCACGGGTTTCACGGATACGTTGCAGTACACCTTCTGAGTGACCGTGGGAGGAGTCGATcagcag                       >  1:1066890/1‑69 (MQ=255)

                                         tttAGCACGGGTTTCACGGATACGTTGCAGTACACCTTCTGAGTGACCGTGGGAGGAGTCGATcagcag                       >  1:2591836/1‑69 (MQ=255)

                                         tttAGCACGGGTTTCACGGATACGTTGCAGTACACCTTCTGAGTGACCGTGGGAGGAGTCGATcagcag                       >  1:612765/1‑69 (MQ=255)

                                         tttAGCACGGGTTTCACGGATACGTTGCAGTACACCTTCTGAGTGACCGTGGGAGGAGTCGATcagcag                       >  1:378671/1‑69 (MQ=255)

                                         tttAGCACGGGTTTCACGGATACGTTGCAGTACACCTTCTGAGTGACCGTGGGAGGAGTCGATcagcag                       >  1:367638/1‑69 (MQ=255)

                                               aCGGGTTTCACGGATACGTTGCAGTACACCTTCTGAGTGACCGTGGGAGGAGTCGATCAGCAGAACGTc                 <  1:3001855/69‑1 (MQ=255)

                                                cGGGTTTCACGGATACGTTGCAGTACACCTTCTGAGTGACCGTGGGAGGAGTCGATCAGCAGAACGTc                 <  1:1158920/68‑1 (MQ=255)

                                                              aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTTCGGCAAc  <  1:2369249/69‑1 (MQ=255)

                                                              aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc  <  1:3092124/69‑1 (MQ=255)

                                                              aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc  <  1:1897878/69‑1 (MQ=255)

                                                              aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc  <  1:3071324/69‑1 (MQ=255)

                                                              aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc  <  1:1816757/69‑1 (MQ=255)

                                                              aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc  <  1:1813324/69‑1 (MQ=255)

                                                              aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc  <  1:3068670/69‑1 (MQ=255)

                                                              aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc  <  1:3014893/69‑1 (MQ=255)

                                                              aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc  <  1:3001844/69‑1 (MQ=255)

                                                              aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc  <  1:1648556/69‑1 (MQ=255)

                                                              aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc  <  1:2058640/69‑1 (MQ=255)

                                                              aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc  <  1:2635992/69‑1 (MQ=255)

                                                              aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc  <  1:1350574/69‑1 (MQ=255)

                                                              aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc  <  1:1248622/69‑1 (MQ=255)

                                                              aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc  <  1:2606107/69‑1 (MQ=255)

                                                              aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc  <  1:941699/69‑1 (MQ=255)

                                                              aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc  <  1:2529081/69‑1 (MQ=255)

GCTGTTGCCACGTTGCCGCCGATAATTTGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGCAGTACACCTTCTGAGTGACCGTGGGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAC  >  minE/1665915‑1666045

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