BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	61,980	T→C	12.3%	K493K (AAA→AAG)	polB ←	DNA polymerase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	61,980	0	T	C	12.3%	149.3 / 11.0	57	K493K (AAA→AAG)	polB	DNA polymerase II
Reads supporting (aligned to +/- strand): ref base T (35/15); new base C (2/5); total (37/20)
Fisher's exact test for biased strand distribution p-value = 8.37e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GCCAGCCGCGGATCGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGTTTGGCTTCATCGCGCCCGTGCCAGATGT  >  minE/61912‑62044

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

             cGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTAATGATGATTTTCAGCGCCTGCGa                                                    >  1:3824803/1‑70 (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:3137110/1‑46 (MQ=255)

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

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

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

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

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

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

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

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

                                   gCCGAGCACGCCGTAGAAGGCGTTCATGATTATCTTCAGCGCCTGCGa                                                    >  1:531144/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)

                                                tAAAAGGCATTCATGATGATTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGTCGTTTg                           <  1:670793/60‑1 (MQ=255)

                                                  aaaGGCATTCATGATGATTTTCAGCGCCTGCGACAGCGGTTTGTTa                                       >  1:3349075/1‑46 (MQ=255)

                                                  aaaGGCATTCATGATGATTTTCAGCGCCTGCGACAGCGGTTTGTTAc                                      >  1:3105215/1‑47 (MQ=255)

                                                  aaaGGCATTCATGATGATTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGtt                             >  1:3171794/1‑56 (MQ=255)

                                                  aaaGGCATTCATGATGATTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGtt                             >  1:4124873/1‑56 (MQ=255)

                                                  aaaGGCATTCATGATGATTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGtt                             >  1:319388/1‑56 (MQ=255)

                                                  aaaGGCATTCATGATGATTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGTTTGGc                         >  1:2630115/1‑60 (MQ=255)

                                                  aaaGGCATTCATGATGATTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGTTTGGc                         >  1:2624123/1‑60 (MQ=255)

                                                  aaaGGCATTCATGATGATTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGTTTGGc                         >  1:1547859/1‑60 (MQ=255)

                                                   aaGGCATTCATGATGATTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGTTTGGCTTCATcgcg                <  1:5143955/68‑1 (MQ=255)

                                                   aaGGCATTCATGATGATTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGTTTGGCTTCATCGCTccc             <  1:2903264/71‑1 (MQ=255)

                                                       cATTCATGATGATTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGTTTGGCTTCATCGCGCCCGTGc         <  1:4586408/71‑1 (MQ=255)

                                                         ttCATGATGATTTTCAGCGCCTGCGACAGCGGttt                                           <  1:3007911/35‑1 (MQ=255)

                                                           cATGATGATTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGTTTGGCTTCATCGCGCCCGTGc         >  1:754400/1‑67 (MQ=255)

                                                            atgatgATTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGTTTg                           >  1:1227919/1‑48 (MQ=255)

                                                            atgatgATTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGTTTg                           >  1:2781189/1‑48 (MQ=255)

                                                            atgatgATTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGTTTg                           >  1:4932427/1‑48 (MQ=255)

                                                            atgatgATTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGTTTGGCTTCATCGCGCCCGTGCCAGa     <  1:716386/70‑1 (MQ=255)

                                                               atgatTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGTTTgg                          >  1:1221174/1‑46 (MQ=255)

                                                               atgatTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGTTTGGCTTCa                     >  1:1405235/1‑51 (MQ=255)

                                                               atgatTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGTTTGGCTTCATCGCGCCCGTGCCAGATGt  >  1:486699/1‑70 (MQ=255)

                                                               atgatTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGTTTGGCTTCATCGCGCCCGTGCCAGATGt  >  1:2953691/1‑70 (MQ=255)

                                                               atgatTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGTTTGGCTTCATCGCGCCCGTGCCAGATGt  >  1:3711000/1‑70 (MQ=255)

GCCAGCCGCGGATCGAAGAAGCGGCAGGCGGTGGTGCCGAGCACGCCATAAAAGGCATTCATGATGATTTTCAGCGCCTGCGACAGCGGTTTGTTACCCTGGCGTTTGGCTTCATCGCGCCCGTGCCAGATGT  >  minE/61912‑62044

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