BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	196,458	G→A	38.0%	T381T (ACC→ACT)	proS ←	prolyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	196,458	0	G	A	38.0%	62.1 / 66.6	105	T381T (ACC→ACT)	proS	prolyl‑tRNA synthetase
Reads supporting (aligned to +/- strand): ref base G (12/53); new base A (21/19); total (33/72)
Fisher's exact test for biased strand distribution p-value = 4.47e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.95e-01

Only 100 of 106 total aligned reads displayed.

AGCCTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAGTGTTTACCATCGATGTTAGC  >  minE/196393‑196514

agcCTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGtc                                                       <  1:3258542/69‑2 (MQ=255)

  ccTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGc                                                     <  1:723463/69‑1 (MQ=255)

  ccTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGc                                                     <  1:1427781/69‑1 (MQ=255)

  ccTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGc                                                     <  1:1578947/69‑1 (MQ=255)

  ccTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGc                                                     <  1:1962766/69‑1 (MQ=255)

  ccTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGc                                                     <  1:2970881/69‑1 (MQ=255)

  ccTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGc                                                     <  1:2545680/69‑1 (MQ=255)

             cATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGa                                           <  1:323423/68‑1 (MQ=255)

             cATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGa                                           <  1:3332740/68‑1 (MQ=255)

             cATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGa                                           <  1:1634278/68‑1 (MQ=255)

             cATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGa                                           <  1:3029387/68‑1 (MQ=255)

             cATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGa                                           <  1:840895/68‑1 (MQ=255)

             cATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGa                                           <  1:1219812/68‑1 (MQ=255)

             cATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGa                                           <  1:2821701/68‑1 (MQ=255)

             cATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGa                                           <  1:704513/68‑1 (MQ=255)

             cATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGa                                           <  1:1357938/68‑1 (MQ=255)

              aTCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGa                                           <  1:720098/67‑1 (MQ=255)

                 cGGGCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATccc                                       >  1:2809779/1‑68 (MQ=255)

                 cGGGCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATccc                                       >  1:2599746/1‑68 (MQ=255)

                 cGGGCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCa                                      >  1:679728/1‑69 (MQ=255)

                 cGGGCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCa                                      >  1:814584/1‑69 (MQ=255)

                 cGGGCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCa                                      >  1:2182912/1‑69 (MQ=255)

                 cGGGCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCa                                      >  1:2658913/1‑69 (MQ=255)

                 cGGGCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCa                                      >  1:2334160/1‑69 (MQ=255)

                 cGGGCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCa                                      >  1:2536775/1‑69 (MQ=255)

                 cGGGCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCa                                      >  1:1818568/1‑69 (MQ=255)

                 cGGGCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCa                                      >  1:2492017/1‑69 (MQ=255)

                 cGGGCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCa                                      >  1:992969/1‑69 (MQ=255)

                 cGGGCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCa                                      >  1:965421/1‑69 (MQ=255)

                 cGGGCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCa                                      >  1:3408442/1‑69 (MQ=255)

                 cGGGCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCa                                      >  1:152378/1‑69 (MQ=255)

                 cGGGCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCa                                      >  1:3435222/1‑69 (MQ=255)

                 cGGGCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCa                                      >  1:1284012/1‑69 (MQ=255)

                 cGGGCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCa                                      >  1:1272162/1‑69 (MQ=255)

                 cGGGCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCa                                      >  1:173710/1‑69 (MQ=255)

                 cGGGCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCa                                      >  1:1214967/1‑69 (MQ=255)

                 cGGGCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCa                                      >  1:3228723/1‑69 (MQ=255)

                 cGGGCTTGGATCGCCAGCCACCACGTTGCGGATATAAGCAATTTCCGGAGTAGCGACATCGCGATCCCa                                      >  1:2655389/1‑69 (MQ=255)

                    gCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCAGtt                                   <  1:2710709/69‑1 (MQ=255)

                    gCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCAGtt                                   <  1:3434536/69‑1 (MQ=255)

                    gCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCAGtt                                   <  1:3131164/69‑1 (MQ=255)

                    gCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCAGtt                                   <  1:3376567/69‑1 (MQ=255)

                    gCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCAGtt                                   <  1:397888/69‑1 (MQ=255)

                    gCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCAGtt                                   <  1:2921566/69‑1 (MQ=255)

                    gCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCAGtt                                   <  1:964855/69‑1 (MQ=255)

                    gCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCAGtt                                   <  1:1124207/69‑1 (MQ=255)

                    gCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCAGtt                                   <  1:1162671/69‑1 (MQ=255)

                    gCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCAGtt                                   <  1:1170213/69‑1 (MQ=255)

                    gCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCAGtt                                   <  1:1179920/69‑1 (MQ=255)

                    gCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCAGtt                                   <  1:859910/69‑1 (MQ=255)

                    gCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCAGtt                                   <  1:732566/69‑1 (MQ=255)

                    gCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCAGtt                                   <  1:1693512/69‑1 (MQ=255)

                    gCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCAGtt                                   <  1:713857/69‑1 (MQ=255)

                    gCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCAGtt                                   <  1:632957/69‑1 (MQ=255)

                    gCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCAGtt                                   <  1:215814/69‑1 (MQ=255)

                    gCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCAGtt                                   <  1:610998/69‑1 (MQ=255)

                    gCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCGGAGTAGCGACATCGCGATCCCAGtt                                   <  1:499460/69‑1 (MQ=255)

                        ggATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGtt                                   >  1:2000880/1‑65 (MQ=255)

                        ggATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGtt                                   >  1:347297/1‑65 (MQ=255)

                        ggATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGtt                                   >  1:417483/1‑65 (MQ=255)

                        ggATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGtt                                   >  1:959093/1‑65 (MQ=255)

                        ggATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGtt                                   >  1:1095272/1‑65 (MQ=255)

                        ggATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATg                               >  1:3121327/1‑69 (MQ=255)

                        ggATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATg                               >  1:522291/1‑69 (MQ=255)

                        ggATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATg                               >  1:1039401/1‑69 (MQ=255)

                                     ccacGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAGTGtt                  <  1:3300884/69‑1 (MQ=255)

                                     ccacGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAGTGtt                  <  1:481848/69‑1 (MQ=255)

                                     ccacGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAGTGtt                  <  1:386499/69‑1 (MQ=255)

                                     ccacGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAGTGtt                  <  1:2364150/69‑1 (MQ=255)

                                       acGTTACGTATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGtt                                   <  1:1016259/50‑1 (MQ=255)

                                       acGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCGCAGtt                                   <  1:2020108/50‑1 (MQ=255)

                                       acGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGtt                                   <  1:1027269/50‑1 (MQ=255)

                                       acGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGtt                                   <  1:2967524/50‑1 (MQ=255)

                                       acGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGtt                                   <  1:1692720/50‑1 (MQ=255)

                                        cGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGtt                                   <  1:532426/49‑1 (MQ=255)

                                          ttACGGATATCTGCTACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAGTGtt                  <  1:1015335/64‑1 (MQ=255)

                                               gATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAGTGTTTACCATCGAt        <  1:2359629/69‑1 (MQ=255)

                                               gATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAGTGTTTACCATCGAt        <  1:1927264/69‑1 (MQ=255)

                                               gATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAGTGTTTACCATCGAt        <  1:2202242/69‑1 (MQ=255)

                                               gATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAGTGTTTACCATCGAt        <  1:1675957/69‑1 (MQ=255)

                                               gATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAGTGTTTACCATCGAt        <  1:2263108/69‑1 (MQ=255)

                                               gATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAGTGTTTACCATCGAt        <  1:2327740/69‑1 (MQ=255)

                                               gATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAGTGTTTACCATCGAt        <  1:382114/69‑1 (MQ=255)

                                               gATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAGTGTTTACCATCGAt        <  1:3400747/69‑1 (MQ=255)

                                                     tGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAagt                        <  1:3243472/47‑1 (MQ=255)

                                                     tGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAagt                        <  1:160516/47‑1 (MQ=255)

                                                     tGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAagt                        <  1:1107255/47‑1 (MQ=255)

                                                     tGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAagt                        <  1:1036330/47‑1 (MQ=255)

                                                     tGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAagt                        <  1:1909001/47‑1 (MQ=255)

                                                     tGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAagt                        <  1:637967/47‑1 (MQ=255)

                                                     tGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAagt                        <  1:2866204/47‑1 (MQ=255)

                                                     tGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAagt                        <  1:302501/47‑1 (MQ=255)

                                                     tGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAGTGTTTACCATCGATGTTAg   <  1:539879/68‑1 (MQ=255)

                                                     tGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAGTGTTTACCATCGATGTTAg   <  1:2792603/68‑1 (MQ=255)

                                                     tGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAGTGTTTACCATCGATGTTAg   <  1:2824866/68‑1 (MQ=255)

                                                     tGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAGTGTTTACCATCGATGTTAg   <  1:3226568/68‑1 (MQ=255)

                                                     tGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAGTGTTTACCATCGATGTTAGc  >  1:2086872/1‑69 (MQ=255)

                                                     tGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAGTGTTTACCATCGATGTTAGc  >  1:209050/1‑69 (MQ=255)

                                                     tGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAGTGTTTACCATCGATGTTAGc  >  1:1354208/1‑69 (MQ=255)

                                                               ggggTAGCGACATCGCGATCCCAGTTGATGCCGAagt                        <  1:1798680/37‑1 (MQ=255)

AGCCTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAGTGTTTACCATCGATGTTAGC  >  minE/196393‑196514

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