BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	196,451	C→T	37.9%	V384I (GTT→ATT)	proS ←	prolyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	196,451	0	C	T	37.9%	63.8 / 68.6	108	V384I (GTT→ATT)	proS	prolyl‑tRNA synthetase
Reads supporting (aligned to +/- strand): ref base C (14/53); new base T (22/19); total (36/72)
Fisher's exact test for biased strand distribution p-value = 6.98e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 108 total aligned reads displayed.

TGATCAGCAGCCTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAGTGTTTACCATCGATGTTAGC  >  minE/196385‑196514

tGACCAGCAGCCTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACtt                                                               >  1:457008/1‑69 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                         cGGGCTTGGATCGCCAGCCACCACGTTGCGGATATCAGCAATTTCCggag                                                         >  1:1519186/1‑48 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                ggATCGCCAGCCACCACGTTACGGATATCTGCAACTTcc                                                             >  1:450035/1‑39 (MQ=255)

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

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

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

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

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

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

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

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

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

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

                                             ccacGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAGTGtt                  <  1:3300884/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:2967524/50‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

                                                       gATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAGTGTTTACCATCGAt        <  1:382114/69‑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:160516/47‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

TGATCAGCAGCCTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAGTGTTTACCATCGATGTTAGC  >  minE/196385‑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: `‑`