BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,355,216	T→C	14.9%	Y328C (TAT→TGT)	yihE ←	predicted kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,355,216	0	T	C	14.9%	118.5 / 13.7	47	Y328C (TAT→TGT)	yihE	predicted kinase
Reads supporting (aligned to +/- strand): ref base T (19/21); new base C (5/2); total (24/23)
Fisher's exact test for biased strand distribution p-value = 4.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

TAAAACTAAACCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACTTTTGCCTGTTCTATAAAAGTCGCCGTCT  >  minE/2355157‑2355282

taaaACTAAACCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACa                                                                 >  1:281605/1‑63 (MQ=255)

  aaaCTAAACCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGtt                                                         <  1:3389168/69‑1 (MQ=255)

        aaCCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCaa                                                 <  1:761408/71‑1 (MQ=255)

        aaCCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCaa                                                 <  1:3595646/71‑1 (MQ=255)

        aaCCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCaa                                                 <  1:2958280/71‑1 (MQ=255)

         aCCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCaaa                                                <  1:3377020/71‑1 (MQ=255)

              ccagcgccagcCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAg                                                              >  1:4310701/1‑52 (MQ=255)

              ccagcgccagcCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAg                                                              >  1:3298566/1‑52 (MQ=255)

               cagcgccagcCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGt                                          >  1:1212992/1‑71 (MQ=255)

                     cagccaAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTg                                     >  1:2396406/1‑70 (MQ=255)

                       gccaAATCTTTTTCATGATCTACTCTCTCCGATTAACACCTa                                                               <  1:675536/42‑1 (MQ=38)

                       gacaAATCTTTTTCATGATCTACTCTCTCCGATTAACACCTa                                                               <  1:965228/40‑1 (MQ=38)

                        ccaAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTa                                                        <  1:3989157/48‑1 (MQ=255)

                        ccaAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTa                                                        <  1:3849032/48‑1 (MQ=255)

                        ccaAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTa                                                        <  1:2919393/48‑1 (MQ=255)

                        ccaAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCa                                                  <  1:3343559/54‑1 (MQ=255)

                         caAATCTTTTTCATGATCTACTCTCTTCGATTAACACCTatg                                                             >  1:3189067/1‑40 (MQ=37)

                         caAATCTTTTTCATGATCTACTCTCTCCGATTAATACa                                                                 >  1:3149841/1‑38 (MQ=255)

                         caAATCTTTTTCATGATCTACTCTCTCCGATTAATACa                                                                 >  1:1133073/1‑38 (MQ=255)

                         caAATCTTTTTCATGATCTACTCTCTCCGATTAATACa                                                                 >  1:1194804/1‑38 (MQ=255)

                         caAATCTTTTTCATGATCTACTCTCTCCGATTAACACCTatg                                                             >  1:2375038/1‑40 (MQ=38)

                         caAATCTTTTTCATGATCTACTCTCTCCGATTAACACCTatg                                                             >  1:774990/1‑40 (MQ=38)

                            aTCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACtt                             <  1:2633648/71‑1 (MQ=255)

                            aTCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACtt                             <  1:909682/71‑1 (MQ=255)

                            aTCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACtt                             <  1:630926/71‑1 (MQ=255)

                             tCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGtt                                                         <  1:2778344/42‑1 (MQ=255)

                             tCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCaaa                                                >  1:3124640/1‑51 (MQ=255)

                             tCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCaaa                                                >  1:416352/1‑51 (MQ=255)

                             tCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACttt                            >  1:2658563/1‑71 (MQ=255)

                             tCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACttt                            >  1:814380/1‑71 (MQ=255)

                             tCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACttt                            >  1:4343580/1‑71 (MQ=255)

                             tCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACttt                            >  1:1784110/1‑71 (MQ=255)

                             tCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACtt                             <  1:1791323/70‑1 (MQ=255)

                               tttttCATGATCTACTCTCTCCGATTAATACATAGGTGTTAAtt                                                     >  1:3828928/1‑44 (MQ=255)

                               tttttCATGATCTACTCTCTCCGATTAATACATAGGTGTTAAtt                                                     >  1:1063275/1‑44 (MQ=255)

                                 tttCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTg                                     <  1:119827/58‑1 (MQ=255)

                                  ttCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTg                                     <  1:4044288/57‑1 (MQ=255)

                                  ttCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTGGAACTTTTGCCt                       >  1:3605216/1‑71 (MQ=255)

                                     aTGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACTTTTGCCTGtt                    <  1:3476139/71‑1 (MQ=255)

                                             ctctctCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACTTTTGCCTGTTCtata               <  1:2268380/68‑1 (MQ=255)

                                             ctctctCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACTTTTGCCTGTTCta                 <  1:4363735/66‑1 (MQ=255)

                                             ctctctCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACTTTTGCCTGTTCta                 <  1:250299/66‑1 (MQ=255)

                                                  tcCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACTTTTGCCTGTTCTATAAAAGTCGc       >  1:3949315/1‑71 (MQ=255)

                                                    cGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACTTTTGCCTGTTCTATAAAAGTCGc       <  1:96321/69‑1 (MQ=255)

                                                    cGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACTTTTGCCTGTTCTATAAAAGTCGCCg     >  1:1336467/1‑71 (MQ=255)

                                                       ttAACACCTATGTATTAATTGCAAAGGGGGTTCTTGTAGAACTTTTGCCTGTTCTATAAAAGTCGCCGTCt  >  1:2014098/1‑71 (MQ=255)

                                                       ttAACACCTATGTATTAATTGCAAAGGGGGTTCTTGTAGAACTTTTGCCTGTTCTATAAAAGTAGCCGTCt  >  1:3388958/1‑71 (MQ=255)

TAAAACTAAACCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACTTTTGCCTGTTCTATAAAAGTCGCCGTCT  >  minE/2355157‑2355282

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