BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	4,501,102	T→C	69.1%	noncoding (16/85 nt)	leuX →	tRNA‑Leu

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	4,501,102	0	T	C	69.1%	77.1 / 44.5	55	noncoding (16/85 nt)	leuX	tRNA‑Leu
Reads supporting (aligned to +/- strand): ref base T (11/6); new base C (27/11); total (38/17)
Fisher's exact test for biased strand distribution p-value = 7.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

GGTGGCGTGCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCT  >  W3110S.gb/4501036‑4501162

ggtggCGTGCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAACCgg                                                           >  1:2176762/1‑70 (MQ=255)

ggtggCGTGCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAACCgg                                                           >  1:850732/1‑70 (MQ=255)

ggtggCGTGCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAACCgg                                                           >  1:3198224/1‑70 (MQ=255)

 gtggCGTGCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAACCGGTa                                                         <  1:1642951/71‑1 (MQ=255)

      gTGCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGc                                                    <  1:2540625/71‑1 (MQ=255)

      gTGCGACAGGTATAATCCACAACGTTTCCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGc                                                    <  1:726539/71‑1 (MQ=255)

       tGCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGc                                                    <  1:303504/70‑1 (MQ=255)

          gACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGtt                                                <  1:2217070/71‑1 (MQ=255)

          gACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGtt                                                <  1:3204591/71‑1 (MQ=255)

          gACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAACCGGTAGACGCAGtt                                                <  1:3170420/71‑1 (MQ=255)

                 ataATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAACCgg                                                           >  1:3083779/1‑53 (MQ=255)

                 ataATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAACCgg                                                           >  1:21489/1‑53 (MQ=255)

                 ataATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAACCg                                                            >  1:1113849/1‑52 (MQ=255)

                 ataATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAACCGGTAGACGCAGTTGATTCa                                          >  1:113269/1‑70 (MQ=255)

                  taATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAACCGGTAGACGCAGTTGATTCaa                                         >  1:494392/1‑70 (MQ=255)

                  taATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAACCGGTAGACGCAGTTGATTCaa                                         >  1:2616829/1‑70 (MQ=255)

                  taATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAACCGGTAGACGCAGTTGATTCaa                                         >  1:861267/1‑70 (MQ=255)

                    aTCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAACCGGTAGACGCAGTTGATTCAAAAt                                      >  1:226305/1‑71 (MQ=255)

                    aTCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAACCGGTAGACGCAGTTGATTCAAAAt                                      >  1:2744235/1‑71 (MQ=255)

                    aTCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAACCGGTAGACGCAGTTGATTCAAAAt                                      >  1:1787356/1‑71 (MQ=255)

                      ccACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATc                                     >  1:2342939/1‑70 (MQ=255)

                      ccACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATc                                     >  1:2766501/1‑70 (MQ=255)

                       cacaACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAACCGGTAGACGCAGTTGATTCAAAATCaa                                   >  1:4988/1‑71 (MQ=255)

                       cacaACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAACCGGTAGACGCAGTTGATTCAAAATCaa                                   >  1:3078785/1‑71 (MQ=255)

                       cacaACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAACCGGTAGACGCAGTTGATTCAAAATCaa                                   >  1:715506/1‑71 (MQ=255)

                           aCGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAACCGGTAGACGCAGTTGATTCaaaa                                       >  1:2415748/1‑63 (MQ=255)

                           aCGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAACCGGTAGACGCAGTTGATTCaaaa                                       >  1:643997/1‑63 (MQ=255)

                           aCGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAACCGGTAGACGCAGTTGATTCaaaa                                       >  1:1578039/1‑63 (MQ=255)

                           aCGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAACCGGTAGACGCAGTTGATTCaaa                                        >  1:3688987/1‑62 (MQ=255)

                              ttttCCGCATACCTCTTCAGTGCCGAAGTGGCGAAACCGGTAGACGCAGTTGATTCAAAATCAACCGTAGa                            <  1:3334219/71‑1 (MQ=255)

                                   cGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATa                        <  1:2780747/70‑1 (MQ=255)

                                             ttCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGa                                              >  1:1368168/1‑38 (MQ=255)

                                             ttCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGa                                              >  1:3655689/1‑38 (MQ=255)

                                             ttCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGa                                              >  1:3022729/1‑38 (MQ=255)

                                             ttCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGa                                              >  1:1941115/1‑38 (MQ=255)

                                             ttCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGa                                              >  1:737353/1‑38 (MQ=255)

                                             ttCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGa                                              >  1:1263650/1‑38 (MQ=255)

                                             ttCAGTGCCGAAGTGGCGAAACCGGTAGACGCAGtt                                                <  1:3432917/36‑1 (MQ=255)

                                               cAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTc             >  1:2923149/1‑69 (MQ=255)

                                               cAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTc             >  1:2473079/1‑69 (MQ=255)

                                               cAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTc             >  1:3693003/1‑69 (MQ=255)

                                               cAGTGCCGAAGTGGCGAAACCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTc             <  1:843680/69‑1 (MQ=255)

                                               cAGTGCCGAAGTGGCGAAACCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTc             <  1:2397308/69‑1 (MQ=255)

                                               cAGTGCCGAAGTGGCGAAACCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTc             <  1:1372668/69‑1 (MQ=255)

                                                aGTGCCGAAGTGGCGAAACCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTc             <  1:1353952/68‑1 (MQ=255)

                                                  tGCCGAAGTGGCGAAACCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTc        <  1:1740046/71‑1 (MQ=255)

                                                  tGCCGAAGTGGCGAAACCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTc        <  1:2436705/71‑1 (MQ=255)

                                                  tGCCGAAGTGGCGAAACCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTc        <  1:2306914/71‑1 (MQ=255)

                                                    ccGAAGTGGCGAAACCGGTAGACGCAGTTGATTCAAAATCAAc                                  >  1:540535/1‑43 (MQ=255)

                                                    ccGAAGTGGCGAAACCGGTAGACGCAGTTGATTCAAAATCAAc                                  >  1:707083/1‑43 (MQ=255)

                                                    ccGAAGTGGCGAAACCGGTAGACGCAGTTGATTCAAAATCAAc                                  >  1:1909318/1‑43 (MQ=255)

                                                    ccGAAGTGGCGAAACCGGTAGACGCAGTTGATTCAAAATCAAc                                  >  1:2981315/1‑43 (MQ=255)

                                                    ccGAAGTGGCGAAACCGGTAGACGCAGTTGATTCAAAATCAAc                                  >  1:2038410/1‑43 (MQ=255)

                                                    ccGAAGTGGCGAAACCGGTAGACGCAGTTGATTCAAAATCAAc                                  >  1:2394713/1‑43 (MQ=255)

                                                         gTGGCGAAACCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCt  >  1:1479978/1‑70 (MQ=255)

GGTGGCGTGCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCT  >  W3110S.gb/4501036‑4501162

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