BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	61,959	A→G	12.5%	Y500Y (TAT→TAC)	polB ←	DNA polymerase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	61,959	0	A	G	12.5%	145.0 / 11.6	56	Y500Y (TAT→TAC)	polB	DNA polymerase II
Reads supporting (aligned to +/- strand): ref base A (31/18); new base G (2/5); total (33/23)
Fisher's exact test for biased strand distribution p-value = 1.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

TGGTGATCGACGATGCCAGCCGCGGATCGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGTTTGGCTTCA  >  minE/61898‑62025

tGGTGATCGACGATGCCAGCCGCGGATCGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCa                                                           <  1:4050852/71‑1 (MQ=255)

tGGTGATCGACGATGCCAGCCGCGGATCGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCa                                                           <  1:4895270/71‑1 (MQ=255)

     aTCGACGATGCCAGCCGCGGATCGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCat                                                      <  1:529948/71‑1 (MQ=255)

     aTCGACGATGCCAGCCGCGGATCGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCat                                                      <  1:4960578/71‑1 (MQ=255)

     aTCGACGATGCCAGCCGCGGATCGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCat                                                      <  1:254312/71‑1 (MQ=255)

     aTCGACGATGCCAGCCGCGGATCGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCat                                                      <  1:1848242/71‑1 (MQ=255)

       cgacgaTGCCAGCCGCGGATCGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCatga                                                    >  1:848268/1‑71 (MQ=255)

       cgacgaTGCCAGCCGCGGATCGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCatga                                                    >  1:2847933/1‑71 (MQ=255)

        gacgaTGCCAGCCGCGGATCGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCatga                                                    <  1:2108599/70‑1 (MQ=255)

              gCCAGCCGCGGATCGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGAtttt                                             >  1:4742047/1‑71 (MQ=255)

              gCCAGCCGCGGATCGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGAtttt                                             >  1:358330/1‑71 (MQ=255)

              gCCAGCCGCGGATCGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGAtttt                                             >  1:2344276/1‑71 (MQ=255)

               ccAGGCGCGGGTCGAAGAAACGGCAGCCGGTGGTGCCGAGCACGCCGTAGAAGGCGTTCATGATGATCTTc                                            <  1:2149633/71‑1 (MQ=255)

               ccAGGCGCGGGTCGAAGAAACGGCAGCCGGTGGTGCCGAGCACGCCGTAGAAGGCGTTCATGATGATCTTc                                            <  1:3155193/71‑1 (MQ=255)

               ccAGGCGCGGGTCGAAGAAACGGCAGCCGGTGGTGCCGAGCACGCCGTAGAAGGCGTTCATGATGATCTTc                                            <  1:3149026/71‑1 (MQ=255)

               ccAGGCGCGGGTCGAAGAAACGGCAGCCGGTGGTGCCGAGCACGCCGTAGAAGGCGTTCATGATGATCTTc                                            <  1:4216314/71‑1 (MQ=255)

               ccAGGCGCGGGTCGAAGAAACGGCAGCCGGTGGTGCCGAGCACGCCGTAGAAGGCGTTCATGATGATCTTc                                            <  1:4210221/71‑1 (MQ=255)

                  gCCGCGGATCGAAGAAGCGTCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAgc                                         >  1:2545034/1‑71 (MQ=255)

                  gCCGCGGATCGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAgc                                         >  1:3147449/1‑71 (MQ=255)

                  gCCGCGGATCGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAgc                                         >  1:2699352/1‑71 (MQ=255)

                      cgGATCGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAGCGCCt                                     <  1:4821529/71‑1 (MQ=255)

                      cgGATCGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAGCGCCt                                     <  1:2998960/71‑1 (MQ=255)

                      cgGATCGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAGCGCCt                                     <  1:2774116/71‑1 (MQ=255)

                        gATCGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAGCGCCTGc                                   <  1:2332710/71‑1 (MQ=255)

                         aTCGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCat                                                      <  1:2889953/51‑1 (MQ=255)

                           cGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAGCGCCTGCGa                                 >  1:1225392/1‑70 (MQ=255)

                           cGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAGCGCCTGCGa                                 >  1:1743065/1‑70 (MQ=255)

                           cGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAGCGCCTGCGa                                 >  1:4422257/1‑70 (MQ=255)

                           cGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAGCGCCTGCGa                                 >  1:3991775/1‑70 (MQ=255)

                           cGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTAATGATGATTTTCAGCGCCTGCGa                                 >  1:3824803/1‑70 (MQ=255)

                             aagaagCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCAtgatga                                                 >  1:4069394/1‑52 (MQ=255)

                             aagaagCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCAtgatga                                                 >  1:1185120/1‑52 (MQ=255)

                             aagaagCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCAtgatga                                                 >  1:4431695/1‑52 (MQ=255)

                             aagaagCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCAtgatga                                                 >  1:3820330/1‑52 (MQ=255)

                             aagaagCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCAtgatga                                                 >  1:1290761/1‑52 (MQ=255)

                             aagaagCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCAtgatga                                                 >  1:1553632/1‑52 (MQ=255)

                             aagaagCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCAtgatga                                                 >  1:2522947/1‑52 (MQ=255)

                             aagaagCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCAtgatga                                                 >  1:2297980/1‑52 (MQ=255)

                               gaagCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCatgat                                                   <  1:252617/48‑1 (MQ=255)

                                 agCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAGCGCCTGCGACAGCGGt                          >  1:381994/1‑71 (MQ=255)

                                  gCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATTATTTTCATCGCCt                                     <  1:4741629/59‑1 (MQ=255)

                                     gCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAGCGCCTGCGACAGCGGTTTGt                      <  1:2124368/71‑1 (MQ=255)

                                          cGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAg                                          >  1:3792336/1‑46 (MQ=255)

                                          cGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAg                                          >  1:101698/1‑46 (MQ=255)

                                          cGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAg                                          >  1:814777/1‑46 (MQ=255)

                                          cGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAg                                          >  1:2733028/1‑46 (MQ=255)

                                          cGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAg                                          >  1:3137110/1‑46 (MQ=255)

                                          cGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAg                                          >  1:3415326/1‑46 (MQ=255)

                                          cGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAg                                          >  1:3657449/1‑46 (MQ=255)

                                                tgCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAGCGCCTTCGACAGCGGTTTGTTACCCTGGCGt           <  1:615969/71‑1 (MQ=255)

                                                 gCCGAGCACGCCGTAGAAGGCGTTCATGATTATCTTCAGCGCCTGCGa                                 >  1:531144/1‑48 (MQ=255)

                                                 gCCGAGCACGCCGTAGAAGGCGTTCATGATTATCTTCAGCGCCTGCGa                                 >  1:922626/1‑48 (MQ=255)

                                                 gCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCa                                           >  1:2499409/1‑38 (MQ=255)

                                                   cGAGCACGCCATAAAAGGCATTCATGATGATTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGTTTg        >  1:1223149/1‑71 (MQ=255)

                                                     aGCACGCCATAAAAGGCATTCATGATGATTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGTTTGGc      <  1:4997668/71‑1 (MQ=255)

                                                         cGCCATAAAAGGCATTCATGATGATTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGTTTGGCTTCa  <  1:3058615/71‑1 (MQ=255)

TGGTGATCGACGATGCCAGCCGCGGATCGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGTTTGGCTTCA  >  minE/61898‑62025

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